next up previous

STAT 330 Lecture 23

Reading for Today's Lecture: 10.2, 10.3.

Goals of Today's Lecture:

Today's notes

Confidence Intervals

One at a time confidence intervals for differences between group means, that is, for tex2html_wrap_inline184 :

displaymath186

where tex2html_wrap_inline188 is the degrees of freedom used in computing the MSE and

displaymath190

Problem: If, for example, I=4 then there are tex2html_wrap_inline194 pairs so there are 6 confidence intervals. If each interval has a 5% chance of missing its target then

displaymath196

and

displaymath198

For purposes of interpretation we often pick out the ``most interesting'' looking interval but we want

displaymath200

Tukey's Studentized Range Procedure

Suppose that tex2html_wrap_inline202 are iid N(0,1) and tex2html_wrap_inline206 where tex2html_wrap_inline208 is independent of the Z's. Then we can choose constants tex2html_wrap_inline212 so that

displaymath214

Tukey applied this to

displaymath216

in the case where tex2html_wrap_inline218 and with

displaymath220

In this case we have that all r and t

displaymath226

with probability tex2html_wrap_inline228 . But

eqnarray56

but with the 2 sample pooled estimate of the standard deviation replaced by the root mean squared error, tex2html_wrap_inline232 .

If we solve the inequality

displaymath234

we get a confidence interval for tex2html_wrap_inline236 of the form

displaymath238

These are called simultaneous (for all r and t) confidence intervals because

displaymath244

Example: With I=4, J=6, n=24 the t type multiplier is tex2html_wrap_inline254 so that the one at a time confidence intervals would be

displaymath256

while the Tukey intervals have tex2html_wrap_inline258 (see page 711 and 712, Table A.8) and the intervals are

displaymath260

or

displaymath262

The Tukey multiplier of tex2html_wrap_inline264 which is much wider than the one at a time confidence intervals. This extra width is required to control the overall error rate.

Unequal Sample Sizes

When the sample sizes are not all equal we make an approximation and replace tex2html_wrap_inline266 by

displaymath268

and get intervals

displaymath270

The multiplier tex2html_wrap_inline272 is larger than the usual tex2html_wrap_inline274 for one at a time intervals.

Numerical Example: Coagulation Data

We have the 4 mean values 61, 66, 68 and 61 for diets A, B, C and D. The corresponding sample sizes are tex2html_wrap_inline276 , tex2html_wrap_inline278 and tex2html_wrap_inline280 . We have tex2html_wrap_inline282 , n=24, I=4, and tex2html_wrap_inline258 . SO our 95% Tukey intervals are

Comparison Diff of Means tex2html_wrap_inline290 tex2html_wrap_inline290
tex2html_wrap_inline294 -5 tex2html_wrap_inline296 4.28
tex2html_wrap_inline298 -7 tex2html_wrap_inline300 4.28
tex2html_wrap_inline302 0 tex2html_wrap_inline300 4.06
tex2html_wrap_inline306 -2 tex2html_wrap_inline300 3.82
tex2html_wrap_inline310 5 tex2html_wrap_inline300 3.58
tex2html_wrap_inline314 7 tex2html_wrap_inline300 3.58

Interpretation of Tukey Intervals

Typically we make a small dot plot of the means labelled by the group names and then underline those pairs of means where the difference is not significant, that is, where the Tukey interval includes 0. (The following example uses the numbers which will come out of an example in the next lecture.)

tex2html_wrap_inline324
tex2html_wrap_inline324 tex2html_wrap_inline324 tex2html_wrap_inline324
tex2html_wrap_inline332 tex2html_wrap_inline334 tex2html_wrap_inline336

The conclusion is that tex2html_wrap_inline338 is possible and tex2html_wrap_inline340 is possible but that tex2html_wrap_inline342 , tex2html_wrap_inline344 , tex2html_wrap_inline346 and tex2html_wrap_inline348 .

Warning of potential problems

1: The following picture can arise:

A B C D

so that you know tex2html_wrap_inline358 but think maybe tex2html_wrap_inline360 and maybe tex2html_wrap_inline362 .

2: It can also happen that when you do the F test of tex2html_wrap_inline366 you reject tex2html_wrap_inline368 but then when you look to see why you get the picture:

A B C D

Thus in this situation you are sure there are two different but not sure which two.


next up previous



Richard Lockhart
Fri Feb 13 15:33:43 PST 1998