协方差矩阵lm对象R

Question

这应该很简单，但R中是否有内置函数来获得lm对象的方差协方差模型的无偏估计？然后我的意思是$ \ epsilon \ epsilon＆＃39; /（n-k-1）$？因为使用cov(lm$resid)不会调整变量数量$ k $，还是它？

我当然知道如何手动编程，但我只是好奇是否没有内置函数。提前谢谢。

Answer 1

我相信，假设拟合模型lm1，您只需要summary(lm1)$sigma。

是的，你的直觉是正确的; cov(lm$resid)不会根据自由度进行调整。