anova.merMod使用REML = FALSE重新计算模型，在哪里找到参考？

Question

我一直在使用anova.merMod包中的lme4函数，通过科学出版物的似然比检验获得固定效应的p值（大多数评论者仍然要求我的领域中的p值）。我注意到anova.merMod函数使用lmer重新计算REML = FALSE函数（参见下面的示例），这是一个非常好的功能，迫使不太熟悉的用户正确地进行测试。但是，我一直在尝试阅读lme4包的大部分文档，但找不到此功能的表示法（例如，请参阅?anova.merMod，将用户定向到?vcov.merMod） 。这让我感到困惑。

问题：为什么文档中没有明确提及此功能？或许我理解错了吗？

聚苯乙烯。 R-mailing lists似乎对此有疑问，但答案让我更加困惑。

library(lme4)
data(sleepstudy)

reml <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
noreml <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy, REML = F)

reml0 <- lmer(Reaction ~ (Days | Subject), sleepstudy)
noreml0 <- lmer(Reaction ~ (Days | Subject), sleepstudy, REML = F)

## Returns similar likelihood ratio test statistics:
(a <- anova(reml, reml0))
(b <- anova(noreml, noreml0))

## Not identical though
identical(a, b)
[1] FALSE

编辑： sessionInfo：R version 3.0.2 (2013-09-25)，lme4_1.0-5

Answer 1

感谢Roland做腿部工作，我发表评论作为答案。

我不相信答案是不一样的：如果任何浮点数不完全相同或任何变量的任何名称不同，则相同将返回FALSE。您能看看返回的感兴趣元素的实际值，看看它们之间的差异是否超过机器精度？ -

Roland进行了测试，发现唯一的区别是name属性。

Answer 2

至少从lme4的最新版本开始，{em} 在?anova.merMod中记录了（虽然详细信息）（强调添加）：

  

'anova'：返回贡献的顺序分解             固定效应术语，或多个参数，模型             比较统计。 对于'lmerMod'类的对象             默认行为是使用ML重新安装模型（如果已安装）             'REML = TRUE'，这可以通过'改装'来控制             论点。另见'anova'。