我需要知道如何从随机斜率的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行的值是多少?这不会出现在任何地方的摘要输出中。