在lme4包中,使用lmer命令:
fm2 <- lmer(Reaction ~ Days + (1|Subject) + (0+Days|Subject), data = sleepstudy)
得到如下结果:
Linear mixed model fit by REML
Formula: Reaction ~ Days + (1 | Subject) + (0 + Days | Subject)
Data: sleepstudy
AIC BIC logLik deviance REMLdev
1754 1770 -871.8 1752 1744
Random effects:
Groups Name Variance Std.Dev.
Subject (Intercept) 627.568 25.0513
Subject Days 35.858 5.9882
Residual 653.584 25.5653
Number of obs: 180, groups: Subject, 18
Fixed effects:
Estimate Std. Error t value
(Intercept) 251.405 6.885 36.51
Days 10.467 1.559 6.71
请教各位专家:
如果我只想抽取随机效应的残差值653.584 ,应该如何操作呢?
如果我只想抽取随机效应的截距值627.568 , 应该如何操作呢?
如果我只想抽取随机效应的Days的值35.858 ,应该如何操作呢?