似然比测试-LMER R-非嵌套模型


14

我目前正在审查一些工作,遇到了以下问题,这对我来说似乎是错误的。使用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陈述的使用不正确。

Answers:


8

在两个方面是不正确的:

  1. (普通)似然比检验只能用于比较嵌套模型;
  2. 我们无法在REML下比较均值模型。(这里不是这种情况,请参阅下面的@KarlOveHufthammer的评论。)

在使用ML的情况下,我知道使用AIC或BIC来比较非嵌套模型。


9
关于第2点,anova()R中的功能将REML下安装的两个模型进行比较;它使用ML对其进行改装,然后执行测试。请参阅lme4:::anova.merMod,其中包含该行mods <- lapply(mods, refitML)。(但是您仍然正确,anova()不能
将这

2
还请注意,在嵌套方面存在一些分歧:Brian Ripley说嵌套对于AIC比较至关重要(请参见链接文档的第20页以进行讨论),而安德森和伯纳姆(请参见第2页)不同意..
Ben Bolker

2
@BenBolker 另一个参考请参阅对于非嵌套模型使用AIC的 this),只要您考虑所有标准化常数和非病理模型即可。但是,在LMM中,您必须使用AIC的一些修改
LessFaceMoreBook 2014年

2
链接混乱:我认为stats.ox.ac.uk/~ripley/ModelChoice.pdf应该有效。
Ben Bolker 2014年

2
@BenBolker好吧,Brian Ripley很有主见。但是,他并未提出反对在非嵌套模型中使用AIC的毁灭性论据:)。很抱歉重复您的链接。
LessFaceMoreBook 2014年
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.