R中协方差参数估计表

Question

我已经建立了一个随机截距和斜率模型(lmer)，需要检验随机效应的协方差，才能看到其中的显着性。

我看到SAS输出提供了一个“协方差参数估计”表，作为模型摘要的一部分--参见PROC混合的“输出56.2.6重复度量分析”一节

在R中是否有这样做的方法(并获得每个协方差的p值)？

Answer 1

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

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

来自?summary的文档

cov.unscaled是coef[j], j=1, …, p (无标度)协方差的p×p矩阵。