R : NLME

Question

library(nlme) 
fm1 <- nlme(height ~ SSasymp(age, Asym, R0, lrc), 
      data = Loblolly, 
      fixed = Asym + R0 + lrc ~ 1, 
      random = Asym ~ 1, 
      start = c(Asym = 103, R0 = -8.5, lrc = -3.3)) 

나는 R.에서 nlme 패키지를 사용하여 비선형 혼합 효과 모델을 피팅하고 그리고는 1) 오류가 정규 분포를 그 모델 진단을 수행하고 가정을 확인하고 싶은 그와 확인 모델 가정 2) 무작위 효과가 정상적으로 분배됩니다.

경우 1

), 난 그냥 잔차

qqnorm(fm1$residuals) 

의 간단한 산점도을 수행 할 수 있습니다하지만 난 임의 효과도 정규 분포되어 있는지 어떻게 확인합니까?

Answer 1

일반 함수 ranef (또는 random.effects)을 사용하여 임의 효과를 추출 할 수 있습니다.

qqnorm(oo) 
qqline(oo) 

: 특히, 당신이 원하는 효과를

oo <- ranef(fm1)$Asym 

지금 당신이 정상을 확인하는 QQ 플롯을 사용할 수있다