Hello,
> require(graphics)
> summary(fm1 <- aov(breaks ~ wool + tension, data = warpbreaks))
> THSD <- TukeyHSD(fm1, "tension", ordered = TRUE)
> p.adj <- THSD$tension[,4]
HTH
Pascal
Le 13/12/2012 16:08, raz a écrit :
Hey,
I have this TukeyHSD output from which I would like to extract only the
P-values (p adj, last number).
The problem is that the test output is a character list.
How can I "break" this sentence to separate the Pv?
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = Fe1$Fe ~ Fe1$genotype)
$`Fe1$genotype`
diff lwr upr p adj
A1/D1-A1 -1.48 -12.660238 9.700238 0.9808395
D1-A1 2.72 -8.460238 13.900238 0.8970997
WT-A1 16.23 5.049762 27.410238 0.0037499
D1-A1/D1 4.20 -6.980238 15.380238 0.7092182
WT-A1/D1 17.71 6.529762 28.890238 0.0017396
WT-D1 13.51 2.329762 24.690238 0.0154169
Thanks
______________________________________________
R-help@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.
