나는 현재 일부 작업을 검토하고 있으며 다음과 같은 문제를 겪었습니다. 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
명령문은 일반적인 자유도를 가진 카이 제곱을 사용하여 모델 간의 차이를 테스트합니다. 이것은 나에게 맞지 않는 것 같습니다. 맞다면 누구든지 이것을 정당화하는 참조를 알고 있습니까? 시뮬레이션에 의존하는 방법 (Paper by Lewis et al., 2011)과 Vuong (1989)이 개발 한 접근 방식을 알고 있지만 이것이 여기서 생성 된 것이라고 생각하지 않습니다. 나는 그 anova
진술 의 사용 이 정확 하다고 생각하지 않습니다 .
anova()
하여 R 의 함수는 REML 하에서 적합 화 된 두 모델을 비교 하지 않습니다 . ML을 사용하여 다시 장착 한 다음 테스트를 수행합니다.lme4:::anova.merMod
행을 포함 하는를 참조하십시오mods <- lapply(mods, refitML)
. (하지만anova()
중첩되지 않았기 때문에 두 모델을 비교하는 데 사용할 수없는 것은 여전히 옳습니다 .)