[R] Tukey HSD

[Erin Hodgess]

Dear R People:

Here is a possible answer to Christophe's problem
>    summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks))
            Df Sum Sq Mean Sq F value   Pr(>F)   
wool         1  450.7   450.7  3.3393 0.073614 . 
tension      2 2034.3  1017.1  7.5367 0.001378 **
Residuals   50 6747.9   135.0                    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 
>      TukeyHSD(fm1, "tension", ordered = TRUE)
  Tukey multiple comparisons of means
    95% family-wise confidence level
    factor levels have been ordered

Fit: aov(formula = breaks ~ wool + tension, data = warpbreaks)

$tension
         diff        lwr      upr
M-H  4.722222 -4.6311985 14.07564
L-H 14.722222  5.3688015 24.07564
L-M 10.000000  0.6465793 19.35342

>
This is from the Example for TukeyHSD.

It gives you the confidence intervals for the level that you set.

Not p values.

Hope this helps!  Sorry to send to the entire list, but deleted original
message.

Sincerely,
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: hodgess at gator.uhd.edu