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 -- \m/ [[alternative HTML version deleted]]
Pascal Oettli
2012-Dec-13 07:56 UTC
[R] Physically extracting P-value from TukeyHSD test output
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 >