我目前正在审查一些工作，并且遇到了以下内容，这对我来说似乎是错误的。使用 lmer 拟合（在 R 中）两个混合模型。这些模型是非嵌套的，并通过似然比检验进行比较。简而言之，这是我所拥有的一个可重现的示例：

set.seed(105)
Resp = rnorm(100)
A = factor(rep(1:5,each=20))
B = factor(rep(1:2,times=50))
C = rep(1:4, times=25)
m1 = lmer(Resp ~ A + (1|C), REML = TRUE)
m2 = lmer(Resp ~ B + (1|C), REML = TRUE)
anova(m1,m2)

据我所知，lmer用于计算对数似然，该anova语句使用具有通常自由度的卡方检验模型之间的差异。这对我来说似乎不正确。如果它是正确的，有没有人知道任何证明这一点的参考资料？我知道依赖于模拟的方法（Lewis 等人的论文，2011 年）和 Vuong（1989 年）开发的方法，但我不认为这就是这里产生的。我不认为该anova语句的使用是正确的。