Dear R-help,
I have been trying to use the TukeyHSDs function in the "pgirmess"
package to quickly extract all
significant pairwise comparisons in an aov object. However, it seems
that this function isn't working
as intended when only the two last populations means being tested are
significant.
An example of this can be seen below:
>numbers<-c(464,482,453,434,495,487)
>group<-gl(3,2,label=c("A","B","C"))
>testobject<-aov(numbers~group)
>result<-TukeyHSD(testobject,conf.level=0.95)
>error<-TukeyHSDs(result)
>error$group[,1]
Error in error$group[, 1] : incorrect number of dimensions
To illustrate the "normal" function, the data set below can be used.
>numbers<-c(845,829,682, 689,581,495)
Is there something wrong with this function? Or is there a better way
to extract significant comparisons
from a tukey test?
Cheers,
Rune
-------------------------------------------
Rune Vejen Petersen, M.Sc.Eng.
Statens Serum Institute
Dept. of Infectious Disease Immunology
Artillerivej 5, 81/344
DK-2300 Copenhagen
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
