R

Question

我创建了一个随机拦截和斜率模型（lmer），需要检查随机效应的协方差，看看哪个是显着的。

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

在R？

中是否有办法（并获得每个协方差的p值）

Answer 1

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

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

来自?summary的文档：

  

cov.nnscaled coef[j], j=1, …, p的（未缩放的）协方差的p x p矩阵。