2017-09-12 3 views
1

무작위 절편 및 기울기 모델 (lmer)을 만들었으며 임의 효과의 공분산을 검사하여 어느 것이 중요한 지 확인할 필요가 있습니다.공변 매개 변수 추정 테이블 R

내가 SAS 출력은 모델 요약의 일환으로 '공분산 매개 변수 추정' 테이블을 제공하는 것을 본 적이 - 절 참조 '출력 56.2.6 반복 측정 분석'

PROC MIXED의인가 거기 R에서이 일을하는 방법 (그리고 각 공분산에 대한 p 값을 얻는 것)? ?summary의 문서에서

x <- lm(formula) 
y <- summary(x) 
print(y$cov.unscaled) 

:

+0

'lme4 :: lmer'을 사용한다고 가정 할 때, 당신은 어떤 공분산 행렬을 원하나요? 'summary (yourmodel)'은 고정 효과 ('cov (yourmodel) '와 corv2cor (vcov (yourmodel))를 사용하여 cov 행렬을 추출 할 수 있음을 보여줍니다 - 아마 lme4에서 후자의 단계를 수행하는 방법입니다. VARCorr (yourmodel)'을 사용하여 RE 분산을 추출합니다. – user20650

+0

예, lme4. 무작위 효과의 공분산을 찾고 있는데,> as.data.frame (Varcorr (mymodel))을 사용하여 데이터 프레임을 추출하고 만들었습니다. – TOC

+0

그래, 각 RE var-cov가 0이 아닌 의미를 테스트하고 싶습니까? lme4가 직접 지원한다고 생각하지 않습니다. 다른 방법으로 다른 모델을 비교할 수 있습니다. RE는 조건부 AIC를 사용합니다 (ps ive는 attentoin을 끌기 위해 lme4 태그를 추가했습니다) – user20650

답변

1

예, 모델의 요약의 출력이 찾을 수

(스케일 없음) 공분산의의 cov.unscaled apxp 매트릭스 coef[j], j=1, …, p.