我可以使用以下命令从 nlme 摘要中提取固定效果summary(fm1)。却苦于如何获得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.916033210
Number of Observations: 108
Number of Groups: 27
任何帮助将不胜感激。谢谢
Use 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"
本文内容由网友自发贡献，版权归原作者所有，本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容，请联系:hwhale#tublm.com(使用前将#替换为@)
