事后测试双向混合模型ANOV

我正在做一个混合模型重复措施ANOVA。我想运行事后测试来查看交互TREAT * TIME的p值，但是我只设法使用下面的ghlt Tukey测试，它不会给我所寻找的交互。事后测试双向混合模型ANOV

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")))

来源

2016-11-27 kumbu

使用snk.test(model, term="TREAT*TIME", among="TREAT", within="TIME")从包装GAD如果你有从LSMEANS均衡模型和summary(lsmeans(oi, pairwise ~ TIME*TREAT), infer=TRUE)如果你的模式是不平衡

来源

2016-12-07 11:11:03 kumbu

我也有这个问题。似乎不存在双向ANOVA的直接事后检验。但是，您可能想尝试引导，这是一种双向ANOVA的稳健估计形式。我发现以下链接非常有帮助。

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

它包含使用rcompanion，WRS2，psych，并multcompView封装来执行你的自举方差分析和跟进事后一步一步的教程。祝你好运。

来源

2016-11-30 04:32:39

嗨！感谢您的回答。不过，我更深入地了解了这一点，并问了一些人，实际上这些功能是存在的。如果你有一个平衡模型和'summary（lsmeans（oi，pairwise）'，从包GAD中使用'snk.test（model，term =“TREAT * TIME”，在=“TREAT”中，在=“TIME” 〜TIME * TREAT），推断= TRUE）'如果你的模型不平衡，请享用！ – kumbu

@ kumbu，这太好了！我没有问这个问题，但现在我有一个答案。 :) –

因为你可以从AFEX包，而不是lme()找到与aov_ez()功能的替代混合模式，然后使用lsmeans()进行事后分析。

您将在这里找到一个详细的教程：

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

来源

2017-02-19 14:17:41

事后测试双向混合模型ANOV

回答

相关问题