从lme4输出中提取2级残差的协方差

时间:2017-01-05 02:24:42

标签: r lme4

我需要知道如何从随机斜率的lmer对象中提取二级残差协方差。

library(lme4)
fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
summary(fm1)

这是随机效应表

Random effects:
 Groups   Name        Variance Std.Dev. Corr
 Subject  (Intercept) 612.09   24.740       
          Days         35.07    5.922   0.07
 Residual             654.94   25.592       
Number of obs: 180, groups:  Subject, 18

我知道0.07下的Corr值是协方差相关系数,我可以用它来计算实际的2级协方差。但是,我并不是100%确定我提取的是正确的。我试过了

as.data.frame(VarCorr(model))

得到了

       grp        var1 var2       vcov       sdcor
1  Subject (Intercept) <NA> 612.089751 24.74044768
2  Subject        Days <NA>  35.071662  5.92213326
3  Subject (Intercept) Days   9.604335  0.06555134
4 Residual        <NA> <NA> 654.941040 25.59181589

sdcor列中第3行的值是我正在寻找的,协方差相关系数?这似乎是显而易见的(0.065四舍五入到0.07),但我想确定一下。

vcov下第3行的值是多少?这不会出现在任何地方的摘要输出中。

0 个答案:

没有答案