我想知道是否可以使用函数“TukeyHSD”来执行“aov()”模型与一个因子(例如,GROUP)和一个连续协变量(例如,AGE)的所有成对比较。我做了例如:
library(multcomp)
data('litter', package = 'multcomp')
litter.aov <- aov(weight ~ gesttime + dose, data = litter)
TukeyHSD(litter.aov, which = 'dose')
我收到这样的警告消息:警告消息:在复制中(粘贴(“〜”,xx),数据= mf):非因素被忽略:gesttime
上面这个过程正确吗?警告信息的含义是什么?“TukeyHSD”是否适用于严重不平衡的设计?
另外,上面和下面的流程有区别吗?
litter.mc <- glht(litter.aov, linfct = mcp(dose = 'Tukey'))
summary(litter.mc)
最好的,苏