2017-12-01 13 views
0

나는 https://www.cmm.bris.ac.uk/lemma/login/index.php에서 혼합 모델을 배우고,이 문제는 모듈 5.은 [R]을

데이터는 508 개 학교에서 33,988 관측치가의 운동에 관한됩니다. 학생 데이터는 첫 번째 수준이고 학교는 두 번째 수준입니다. 내가 널 모델을 실행

nullmodel <- lmer(score ~ (1|schoolid), data = mydata, REML = FALSE) 

(I 질문이 더 개념 때문이다. 나는 그것이 문제가되지해야 추측 데이터 세트에 대한 오픈 액세스를 찾을 수 없습니다)

해야하는 우리가 실행하여 얻을 절편 추정 null 모델의 요약은 '점수'의 평균과 같지 않습니까? 절편 추정값을 30.60으로, 평균값 (평균값 (mydata $ score))을 31.09로 계산합니다. 값은 비슷하지만 같지 않습니다.

또한 학교 수준의 잔차 인 u0을 직접 얻을 수있는 방법이 있습니까?

답변

0

임의 요격 효과 u0는 평균 0과 분산 σ_u0^2로 정규 분포됩니다.

σ_u0^2는 학교 간 차이입니다.

직접적인 방법으로 σ_u0^2를 계산할 수 있습니다.

분산은이 2^σ_u0에 가까운 값이 각 학교

  • 계산 분산의 평균 점수를 계산한다.

  • +0

    REML 추정치가 무작위 효과에 대해 편협하기 때문에 최종 모델 적합성을 위해 REML 방법을 권장합니다. –