这个问题在这里已有答案:
我使用nlme包的lme函数拟合了线性混合模型 . 'Summary'命令显示系数,它们的标准误差和系数的相关矩阵 . 如何在R数据集形式中获得此相关矩阵(协方差矩阵也很好)?
lm命令(对于线性模型)确实在对象'summary(lm5)$ cov.unscaled'中提供协方差矩阵 . 我想在lme函数中找到这个对应物 .
如果您想要固定效果的方差 - 协方差矩阵,请使用
vcov(fitted.model)
如果需要相关矩阵,请使用
cov2cor(vcov(fitted.model))
你也可以使用
summary(fitted.model)$corFixed
(使用 str(summary(fitted.model)) 来查找所需的位),但上面的访问器更好,因为它们没有使用结果的(不一定是稳定的)内部结构 .
str(summary(fitted.model))
1 回答
如果您想要固定效果的方差 - 协方差矩阵,请使用
如果需要相关矩阵,请使用
你也可以使用
(使用
str(summary(fitted.model))
来查找所需的位),但上面的访问器更好,因为它们没有使用结果的(不一定是稳定的)内部结构 .