3

我正在使用混合模型进行重复测量方差分析。我想运行一个事后测试来查看交互 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")))
4

3 回答 3

2

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

于 2016-12-07T11:11:03.790 回答
0

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

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

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

于 2016-11-30T04:32:39.863 回答
0

对于混合模型，您可以使用aov_ez()afex 包中的函数代替lme()，然后使用 执行事后分析lsmeans()。

您将在此处找到详细教程：

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

于 2017-02-19T14:17:41.640 回答