我有一个混合设计,其中因子(组:ASD、CTR)和三个因子(时间:前、后;受托人:好、坏;步骤:1、2、3、4、5)之间的混合设计。
我使用 aov_car 执行了混合模型 anova(III 型):
anova <- aov_car(rating ~ time*trustee*step*group + Error(id/time*trustee*step), data = datatemp, anova_table = list(es = "ges"))
我得到了显着的三向交互:时间 * 受托人 * 步骤
我想通过对时间进行有计划的对比来分解交互(我只对不同受托人和步骤的前后对比感兴趣)。
这是我用的cos:
trustee_time_step <- emmeans(anova, ~ time | trustee | step)
pairs(trustee_time_step)
与时间比较的结果
默认情况下,pairs() 应该执行 Tuckey 校正。然而,没有对对比度进行校正。
只有当我使用 step 来划分对比时,才会出现 Tuckey 校正:
step_time_trustee <- emmeans(anova, ~ step | trustee | time)
pairs(step_time_trustee)
与步骤比较的结果
我猜这是因为 Time 只有两个级别(每个对比只是 post vs. pre),而 step 有 5 个级别(请参阅图片末尾的“tukey method for comparison a family of 5estimates”)。
所以这是我的问题:为什么pairs() 在第一种情况下不应用更正?即使 post 与 pre 只是一种对比,它也会进行 10 次(受托人和步骤的每种组合一次),因此我认为应该更正它。我应该怎么办?我要手动更正它们吗?如果是,如何?
我对这件事很困惑,如果这不是“统计良好”,我不想报告未更正的结果。
我希望我的问题足够清楚。感谢您的帮助!