Question

我正在使用混合模型重复测量anova。我想运行一个事后测试来查看交互TREAT * TIME的p值，但我只是设法使用以下ghlt Tukey测试，它没有给我正在寻找的交互。

library(multcomp)
library(nlme)
oi<-lme(total ~ TREAT * TIME, data=TURN, random = ~1|NO_UNIT)
anova(oi)
summary(glht(oi, linfct=mcp(TIME="Tukey", TREAT="Tukey")))

我想要的是：

summary(glht(oi, linfct=mcp(TIME="Tukey",TREAT="Tukey",TREAT*TIME="Tukey")))

Answer 1

如果你有一个平衡模型，请使用GAD包中的snk.test(model, term="TREAT*TIME", among="TREAT", within="TIME")，如果模型不平衡，请使用lsmeans中的summary( lsmeans( oi, pairwise ~ TIME*TREAT), infer=TRUE)

Answer 2

我也遇到过这个问题。似乎不存在针对双向ANOVA的直接事后检验。但是，您可能希望尝试自举，这是对双向ANOVA的稳健估计的一种形式。我发现以下链接非常有用。

http://rcompanion.org/rcompanion/d_08a.html

它包含使用rcompanion，WRS2，psych和multcompView程序包的分步教程，以执行您的自举ANOVA并跟进帖子特别。祝你好运。

Answer 3

对于混合模型，您可以使用afex包中的aov_ez()函数而不是lme()找到替代方法，然后使用lsmeans()进行事后分析。

您可以在此处找到详细的教程：

https://www.psychologie.uni-heidelberg.de/ae/meth/team/mertens/blog/anova_in_r_made_easy.nb.html

双向混合模型anova的事后检验

3 个答案: