[R] Help with TukeyHSD

eugen pircalabelu eugen_pircalabelu at yahoo.com
Thu Apr 15 18:05:01 CEST 2010


Hi,
Is this what you want?
 

?TukeyHSD
summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks))
tk.fm1<-TukeyHSD(fm1, "tension", ordered = TRUE)
tk.fm1
(tk.fm1[[1]][,4])




 Eugen Pircalabelu
(0032)471 842 140
(0040)727 839 293



----- Original Message ----
From: Amit Patel <amitrhelp at yahoo.co.uk>
To: R-help at r-project.org
Sent: Thu, April 15, 2010 4:46:08 PM
Subject: [R] Help with TukeyHSD

Hi

I am conducting ANOVA using the aov function
I am also conducting TukeyHSD to obtain which of the groups show variance
How can I obtain the first three p values from the list below?

>zzz.aov <- aov(Intensity ~ Group, data = zzzanova)

> TukeyHSD(zzz.aov) 
  Tukey multiple comparisons of means
    95% family-wise confidence level

Fit: aov(formula = Intensity ~ Group, data = zzzanova)

$Group
                      diff       lwr      upr     p adj
Group2-Group1  0.778354812 -3.414233 4.970943 0.9607836
Group3-Group1 -0.734044848 -5.073786 3.605696 0.9698685
Group4-Group1 -0.000158625 -4.192747 4.192429 1.0000000
Group3-Group2 -1.512399661 -5.852140 2.827341 0.7933015
Group4-Group2 -0.778513438 -4.971101 3.414075 0.9607611
Group4-Group3  0.733886223 -3.605855 5.073627 0.9698870





______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list