协方差矩阵lm对象R

Question

这应该很容易，但是R中是否有一个内置函数来获得lm对象的方差协方差模型的无偏估计？然后我的意思是$ \ epsilon \ epsilon'/（n-k-1）$？因为使用cov(lm$resid)不会针对变量$ k $的数量进行调整，还是会这样做？

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

Answer 1

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

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