선형 혼합 효과 모델의 레버리지 추출 ​​및 쿡 거리 추출 방법

누구든지 mer클래스 객체 ( lme4패키지 를 통해 얻은 )를 활용하고 Cook의 거리를 계산 (또는 추출)하는 방법을 알고 있습니까? 잔차 분석을 위해 이것을 플롯하고 싶습니다.

R 패키지를 살펴 봐야합니다 influence.ME. 로 생성 된 혼합 효과 모델에 대한 영향력있는 데이터의 측정 값을 계산할 수 있습니다 lme4.

예제 모델 :

library(lme4)
model <- lmer(mpg ~ disp + (1 | cyl), mtcars)

이 기능 influence은 모든 추가 단계의 기초입니다.

library(influence.ME)
infl <- influence(model, obs = TRUE)

쿡의 거리를 계산하십시오 :

cooks.distance(infl)

쿡 거리 :

plot(infl, which = "cook")