R 中的协方差参数估计表

OCT*_*OCT 5 r covariance lme4

我创建了一个随机截距和斜率模型 ( lmer),需要检查随机效应的协方差以确定哪个是显着的。

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

R 中有没有办法做到这一点(并获得每个协方差的 p 值)?

JAD*_*JAD 2

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

\n\n
x <- lm(formula)\ny <- summary(x)\nprint(y$cov.unscaled)\n
Run Code Online (Sandbox Code Playgroud)\n\n

来自以下文档?summary

\n\n
\n

cov.unscaled apxp 的(未缩放的)协方差矩阵coef[j],\n j=1, \xe2\x80\xa6, p

\n
\n