r：如何读取lmer中的相关矩阵？

Question

我正在lme4中使用r包拟合线性混合效果模型，结果显示为：

m4 <- lmer(y ~ 0 + X + (0+ X|subject))

我想知道如何在绿色框中读取相关矩阵并将其用于以后的计算？

我尝试了

VarCorr(m4)

它给出了

我也尝试过

有没有办法获得Corr部分？

谢谢！

Answer 1

尝试类似这样的方法。保存data.frame对象，您将拥有一个易于操作的对象。 这样做：

m<-lmer(Sepal.Length~0+Species|Sepal.Length,iris)
res<-summary(m)[13][[1]]

res1<-res$Sepal.Length #replace with res$Subject as per OP's question
attr(res1,"correlation")


                      Speciessetosa     Speciesversicolor Speciesvirginica
Speciessetosa                 1                 1                1
Speciesversicolor             1                 1                1
Speciesvirginica              1                 1                1