lmer中混合效果模型的anova（）函数给出了错误的输出

Question

我正在将两种混合效果模型相互比较。

mod_1 <- lmer(firstfix ~ ps + cd + ps:cd + (1+ cd||subj) + (1+ ps||object), data = dat2, REML = FALSE)

mod_2 <- lmer(firstfix ~ ps  + ps:cd + (1+ cd||subj) + (1+ ps||object), data = dat2, REML = FALSE)

anova(mod_1, mod_2)

但是anova（）函数的输出是错误的。它看起来像是：

各个型号给出了正确的输出，所以我不确定为什么anova（）函数很难？感谢有关可能导致此问题或如何解决问题的任何建议。

Answer 1

我强烈怀疑您输入的是summary(anova(mod_1,mod_2))而非anova(mod_1,mod_2)。我可以通过

获得非常相似的输出

library(lme4)
mod_1 <- lmer(Reaction~Days+(Days|Subject),sleepstudy)
mod_2 <- lmer(Reaction~1+(Days|Subject),sleepstudy)
summary(anova(mod_1,mod_2))

而anova(mod_1,mod_2)没问题。

Answer 2

您可以尝试在拟合模型时使用*而不是分号对交互进行编码：

  

lmer（firstfix~ps + ps * cd +（1+ cd || subj）+（1+ ps || object），data = dat2，REML = FALSE）