协方差参数估计表中的表

Question

我创建了一个随机截距和斜率模型（lmer），需要检查随机效应的协方差以查看哪个很重要。

我已经看到，SAS输出提供了'协方差参数估计'表作为模型摘要的一部分 - 请参见'输出56.2.6重复测量分析'Proc混合

是否有一种方法（并在R？

Answer 1

是的，您可以在模型摘要的输出中找到它：

x <- lm(formula)
y <- summary(x)
print(y$cov.unscaled)

来自?summary的文档：

cov.unscaled coef[j], j=1, …, p的（未量化）协方差的P X P矩阵。