首页 文章

lme中系数的协方差(或相关)矩阵[重复]

提问于
浏览
0

这个问题在这里已有答案:

我使用nlme包的lme函数拟合了线性混合模型 . 'Summary'命令显示系数,它们的标准误差和系数的相关矩阵 . 如何在R数据集形式中获得此相关矩阵(协方差矩阵也很好)?

lm命令(对于线性模型)确实在对象'summary(lm5)$ cov.unscaled'中提供协方差矩阵 . 我想在lme函数中找到这个对应物 .

1 回答

  • 3

    如果您想要固定效果的方差 - 协方差矩阵,请使用

    vcov(fitted.model)
    

    如果需要相关矩阵,请使用

    cov2cor(vcov(fitted.model))
    

    你也可以使用

    summary(fitted.model)$corFixed
    

    (使用 str(summary(fitted.model)) 来查找所需的位),但上面的访问器更好,因为它们没有使用结果的(不一定是稳定的)内部结构 .

相关问题