我已经构建了一些使用二分变量及其交互作为回归量的 LMEM,我对输出感到困惑。当我只评估单个交互时,使用下面的代码(这是我的统计培训提供的方法),lme4 执行 III 型 Wald F 测试:
value3 <- lmer(neurosynth_value ~ contrast*Teen_vs_Adult + (1 | participant_id),dCon)
Anova(value3, type= 3, test = 'F')
输出看起来像这样,这正是我想要的:
Analysis of Deviance Table (Type III Wald F tests with Kenward-Roger df)
Response: neurosynth_value
F Df Df.res
(Intercept) 5.9460 1 328.06
contrast 5.9198 1 511.34
Teen_vs_Adult 3.4148 1 329.88
contrast:Teen_vs_Adult 1.5901 1 509.75
Pr(>F)
(Intercept) 0.01528 *
contrast 0.01531 *
Teen_vs_Adult 0.06551 .
但是当我添加一个额外的二分回归器及其交互时,Anova 函数会改变方法来执行 III 型 Wald Chi-Square?该代码和结果如下。
代码:
valuecon1 <- lmer(neurosynth_value ~ contrast*valence*Teen_vs_Adult + (1 | participant_id),dConValue)
Anova(conflictcon1, type = 3, Test = 'F')
输出:
Analysis of Deviance Table (Type III Wald chisquare tests)
Response: neurosynth_conflict
Chisq Df Pr(>Chisq)
(Intercept) 95.4341 1 < 2.2e-16 ***
contrast 16.9597 1 3.818e-05 ***
valence 4.0655 1 0.043768 *
Teen_vs_Adult 32.8155 1 1.013e-08 ***
contrast:valence 8.0806 1 0.004474 **
contrast:Teen_vs_Adult 10.6639 1 0.001092 **
valence:Teen_vs_Adult 2.4147 1 0.120200
contrast:valence:Teen_vs_Adult 6.2718 1 0.012268 *
有人知道为什么吗?