0

我有一个混合设计,其中因子(组: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 次(受托人和步骤的每种组合一次),因此我认为应该更正它。我应该怎么办?我要手动更正它们吗?如果是,如何?

我对这件事很困惑,如果这不是“统计良好”,我不想报告未更正的结果。

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

4

0 回答 0