双向混合模型ANOVA的事后tukey检验

时间:2015-04-22 17:34:46

标签: anova posthoc

我和我的一些学生已经在很多地方寻找解决方案,没有运气,几个月也没有。我继续被提到我不想使用的lme命令。提供的输出不是我的同事或我自己使用超过15年的输出。此外,鉴于我使用R作为教学工具,它在t测试之后不会流动,并且对于介绍统计学学生而言也是单向的anovas。我正在进行双向RM ANOVA,重复一次因子。我成功地让R复制了Sigmaplot为主要效果提供的内容。然而,R给出的事后分析与Sigmaplot中的相同事后显着不同。这是我使用的代码 - 附注(因为我也用它来教学生)。

#IV between: IVB1 - Independent variable - between subject factor
#IV within: IVW1  - Independent variable - within subject factor
#DV: DV           - Dependent variable.

aov1= aov(DV ~ IVB1*IVW1 + Error(Subject/IVW1)+(IVB1), data=objectL)
summary(aov1)

# post hoc analysis
ph1=TukeyHSD(aov(DV ~ IVB1*IVW1, data=objectL))
ph1

我希望有人可以提供帮助。 谢谢!

1 个答案:

答案 0 :(得分:0)

我也有这个问题,我找到了afex包中的aov_ez()函数而不是aov()的方便替代方法,然后使用lsmeans()而不是TukeyHSD()执行事后分析:

model <- aov_ez(data,                           
            id="SubjID",                     
            dv="DV",                         
            within=c("IVW1", "IVW2"), 
            between = "IVB1")

# Post hoc
comp = lsmeans(model,specs = ~ IVB1: IVW1: IVW2, adjust = "tukey")
contrast(comp,method="pairwise")

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

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