我创建了一个随机截距和斜率模型 ( lmer
),需要检查随机效应的协方差以查看哪个是显着的。
我已经看到 SAS 输出提供了一个“协方差参数估计”表作为模型摘要的一部分 - 请参阅 PROC MIXED 的“输出 56.2.6 重复测量分析”部分
在 R 中有没有办法做到这一点(并获得每个协方差的 p 值)?
最佳答案
是的,您可以在模型摘要的输出中找到:
x <- lm(formula)
y <- summary(x)
print(y$cov.unscaled)
从
?summary
的文档中:cov.unscaled a p x p matrix of (unscaled) covariances of the
coef[j], j=1, …, p
.
关于r - R中的协方差参数估计表,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/46172305/