困惑:为什么 lme4 将技术从 Wald F 测试更改为 Wald Chisquare?

机器算法验证 假设检验 混合模式 lme4-nlme 重复测量 平方和
2022-03-27 13:54:03

我已经构建了一些使用二分变量及其交互作为回归量的 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 *  

有人知道为什么吗?

1个回答

它不是

Anova(conflictcon1, type = 3, Test = 'F')

反而

Anova(conflictcon1, type = 3, test.statistic = 'F')