lme中系数的协方差（或相关）矩阵

Question

我使用nlme包的lme函数拟合了线性混合模型。 ＆＃39;总结＆＃39;命令显示系数，它们的标准误差和系数的相关矩阵。 如何在R数据集形式中获得此相关矩阵（协方差矩阵也很好）？

lm命令（对于线性模型）确实在对象＆＃39;摘要（lm5）$ cov.nnscaled＆＃39;中提供协方差矩阵。我想在lme函数中找到这个对应物。

Answer 1

如果需要固定效果的方差 - 协方差矩阵，请使用

vcov(fitted.model)

如果您需要相关矩阵，请使用

cov2cor(vcov(fitted.model))

您也可以使用

summary(fitted.model)$corFixed

（使用str(summary(fitted.model))查找您需要的位数），但上面的访问器更好，因为它们不会使用结果的（不一定是稳定的）内部结构。