Tuckey校正与R中的emmeans和Pairs()进行计划的对比

时间:2019-01-15 15:54:22

标签: r anova posthoc emmeans

我有一个混合设计,其中因素(组:ASD,CTR)之间有1个,因素内有三个(时间:之前,之后;受托人:好,坏;步骤:1、2、3、4、5)。 / p>

我用aov_car进行了混合模型方差分析(III型):

anova <- aov_car(rating ~ time*trustee*step*group + Error(id/time*trustee*step), data = datatemp, anova_table = list(es = "ges"))

我进行了重要的三种交互:时间*受托人*步骤

我想通过进行计划的时间对比来分解交互(我只是对不同受托人和步骤的前后之间的对比感兴趣)。

这是我用过的姿势:

trustee_time_step <- emmeans(anova, ~ time | trustee | step) 
pairs(trustee_time_step)

与时间进行比较的结果

enter image description here

默认情况下,pairs()应该执行Tuckey校正。 但是,没有对对比度进行校正。

仅当我使用步骤来划分对比度时,Tuckey校正才会出现:

step_time_trustee <- emmeans(anova,  ~ step | trustee | time) 
pairs(step_time_trustee)

与step比较的结果

enter image description here

我想这是因为时间只有两个级别(每个对比只是后期对比之前),而步骤只有5个级别(请参见图片末尾的“用于比较5个估计量的tukey方法”)。 / p>

这是我的问题:为什么在最初的情况下,pairs()不应用校正?即使post vs. pre只是一个对比,它也进行了10次(对于受托人和步骤的每个组合一次),因此我认为应该予以纠正。 我该怎么办?我可以手动更正吗?如果是的话,怎么办?

我对此事非常困惑,如果这不是“统计上好的”,我也不想报告未更正的结果。

我希望我的问题足够清楚。 谢谢您的帮助!

0 个答案:

没有答案