Trích xuất hiệu ứng Ngẫu nhiên từ tóm tắt nlme

Question

Tôi có thể trích xuất Hiệu ứng cố định từ bản tóm tắt nlme bằng cách sử dụng summary(fm1). Nhưng phải vật lộn cách lấy phần Random effects:.

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 

Mọi trợ giúp sẽ được đánh giá cao. Cảm ơn

Answer 1

Sử dụng ranef(fm1) để trích xuất cho từng chủ đề.

cập nhật để cung cấp cho mã cho khai thác từ bảng tóm tắt:

>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"