Displaying 1 result from an estimated 1 matches for "mcalinfct".
2011 Dec 11
2
multiple comparison of interaction of ANCOVA
...ion effect is significant. I hope to use
TukeyHSD() or glht() to do multiple comparison on Trt:year. However, for
variable year is not a factor, they all give error messages.
I try to follow the demo("MMC.WoodEnergy") in HH package, as follwoing:
> library(HH)
> mca.1993 <- mcalinfct(mydata.aov, "Trt")
> non.zero <- mca.1993[,5:6] != 0
> mca.1993[,5:6][non.zero] <- 1993 * sign(mca.1993[,5:6][non.zero])
> summary(glht(mydata.aov, linfct=mca.1993))
Simultaneous Tests for General Linear Hypotheses
Fit: aov(formula = y ~ Trt * year, data = my...