2

我一直在使用 package 中的anova.merMod函数lme4通过科学出版物的似然比检验获得固定效应的 p 值（大多数审稿人仍然要求我所在领域的 p 值）。我注意到该函数使用（参见下面的示例）anova.merMod重新计算函数，这是一个非常好的功能，它迫使不太熟悉的用户正确地进行测试。但是，我一直在尝试阅读package 的大部分文档，但找不到此功能的注释（例如，请参阅，它会将用户引导至）。这让我很困惑。lmerREML = FALSElme4?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

4

2 回答 2

1
于 2017-11-07T23:30:31.400 回答
0

感谢 Roland 的努力，我将我的评论作为答案发表。

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

罗兰跑了测试，发现唯一的区别就是name属性。

于 2014-02-03T12:46:53.870 回答