协方差矩阵lm对象R

Question

这应该很容易，但是在R中是否有一个内置函数来得到lm对象的方差协方差模型的无偏估计？我的意思是$\epsilon \epsilon‘/(N 1)$？因为使用cov(lm$resid)不会调整变量$k$的数量，或者它是吗？

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

Answer 1

我相信，假设适合的模型lm1，你只想要summary(lm1)$sigma。

是的，你的直觉是正确的；cov(lm$resid)不会调整自由度。