[R] TukeyHSDs function (pgirmess package)
ehlers at math.ucalgary.ca
Mon Apr 10 21:36:29 CEST 2006
Rune Vejen Petersen wrote:
> 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
> An example of this can be seen below:
> Error in error$group[, 1] : incorrect number of dimensions
> To illustrate the "normal" function, the data set below can be used.
> Is there something wrong with this function? Or is there a better way
> to extract significant comparisons
> from a tukey test?
It looks like the line
res[] <- x[y, ]
in TukeyHSDs() should be replaced with
res[] <- x[y, , drop = FALSE]
> 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]]
> R-help at stat.math.ethz.ch mailing list
> PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
More information about the R-help