nlme 요약에서 임의 효과 추출 요약

summary(fm1)을 사용하여 nlme 요약에서 고정 효과를 추출 할 수 있습니다. 그러나 분투하는 방법을 얻으려면 Random effects: 부분을 얻으십시오.nlme 요약에서 임의 효과 추출 요약

fm1 <- lme(distance ~ age, Orthodont, random = ~ age | Subject) 
summary(fm1) 
Linear mixed-effects model fit by REML 
Data: Orthodont 
     AIC  BIC logLik 
    454.6367 470.6173 -221.3183 

Random effects: 
Formula: ~age | Subject 
Structure: General positive-definite, Log-Cholesky parametrization 
      StdDev Corr 
(Intercept) 2.3270340 (Intr) 
age   0.2264278 -0.609 
Residual 1.3100397  

Fixed effects: distance ~ age 
       Value Std.Error DF t-value p-value 
(Intercept) 16.761111 0.7752460 80 21.620377  0 
age   0.660185 0.0712533 80 9.265333  0 
Correlation: 
    (Intr) 
age -0.848 

Standardized Within-Group Residuals: 
     Min   Q1   Med   Q3   Max 
-3.223106086 -0.493761144 0.007316631 0.472151121 3.

Number of Observations: 108 
Number of Groups: 27

도움이 될 것입니다. 감사합니다

출처

2012-01-28 MYaseen208

각 주제에 대해 추출하려면 ranef(fm1)을 사용하십시오.

요약 테이블에서 추출 코드를 제공하도록 업데이트 : 답변에 대한

>VarCorr(fm1) 
Subject = pdLogChol(age) 
      Variance StdDev Corr 
(Intercept) 5.41508758 2.3270341 (Intr) 
age   0.05126955 0.2264278 -0.609 
Residual 1.71620400 1.3100397 

> temp <- VarCorr(fm1) 
> temp[,2] 
(Intercept)   age Residual 
"2.3270341" "0.2264278" "1.3100397" 

> temp[1,2] 
[1] "2.3270341"

출처

2012-01-28 08:28:18 Michelle

감사합니다. 'ranef'는 랜덤 한 효과를 주지만'Random Effects ::'의'StdDev' 부분은'summary (fm1)'에서 가져옵니다. – MYaseen208

니스. 나는 lme 객체에 대한 print/summary 경로를 따르고 있었지만'print.modelStruct'에서 길을 잃었다. –

감사합니다. 내 답변에 대한 의견을 보내 주시면 감사하겠습니다. :) – Michelle

nlme 요약에서 임의 효과 추출 요약

답변

관련 문제