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

R 从lme4输出中提取二级残差的协方差,r,lme4,R,Lme4,我需要知道如何从一个随机物体中提取二级剩余协方差 library(lme4) fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy) summary(fm1) 我知道Corr下的值0.07是协方差相关系数,我可以用它来计算实际的二级协方差。然而,我不是100%确定我提取的是正确的。我试过了 as.data.frame(VarCorr(model)) 得到 grp var1 var2

我需要知道如何从一个随机物体中提取二级剩余协方差

library(lme4)
fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)
summary(fm1)
我知道
Corr
下的值
0.07
是协方差相关系数,我可以用它来计算实际的二级协方差。然而,我不是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
grp var1 var2 vcov sdcor
1主题(截取)612.089751 24.74044768
2受试者日35.071662 5.92213326
3受试者(截距)天9.604335 0.06555134
4剩余654.941040 25.59181589
我要查找的是
sdcor
列第3行中的值,即协方差相关系数吗?这似乎很明显(0.065四舍五入到0.07),但我想确定一下


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

经过阅读,我发现vcov列第3行中的值实际上是2级残差的协方差。如果其他人需要的话,我会保留这个。读了一些资料后,我发现vcov列第3行的值实际上是2级残差的协方差。我会把这个留着,以防别人需要。
       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