I'm performing a relatively simple one-way ANOVA as part of a class exercise. I'd like to get pairwise uncorrected p-values out of R, so I can do a sequential FDR test in another package (I realize there are FDR packages in R as well). I've set up my ANOVA as shown below, and it works fine, produces results, but I can't seem to figure out how to get the raw, uncorrected p-values back. Ultimately, I'd like to do the pairwise tests for both FDR and sequential Bonferroni in R, but this seems like a first step.
It looks like
multcomp might get me where I'm trying to go, but I'm having a difficult time figuring out which will do what I'm looking for.
R> head(d10) time breed 1 27.4 type.A 2 18.3 type.A 3 24.3 type.B 4 19.6 type.B 5 21.6 type.C 6 30.3 type.D a10 <- aov(time~breed,data=d10) # reports the overall significance, but nothing on the pairs summary(a10) # reports corrected values only TukeyHSD(a10)