我创建了一个随机截距和斜率模型 ( lmer),需要检查随机效应的协方差以确定哪个是显着的。
我已经看到 SAS 输出提供了一个“协方差参数估计”表作为模型摘要的一部分 - 请参阅PROC MIXED的“输出 56.2.6 重复测量分析”部分
R 中有没有办法做到这一点(并获得每个协方差的 p 值)?
是的,您可以在模型摘要的输出中找到这一点:
\n\nx <- lm(formula)\ny <- summary(x)\nprint(y$cov.unscaled)\nRun Code Online (Sandbox Code Playgroud)\n\n来自以下文档?summary:
\n\ncov.unscaled apxp 的(未缩放的)协方差矩阵
\ncoef[j],\n j=1, \xe2\x80\xa6, p。