協方差矩陣lm對象R

Question

這應該很容易，但是R中是否有一個內置函數來獲得lm對象的方差協方差模型的無偏估計？然後我的意思是$ \ epsilon \ epsilon'/（n-k-1）$？因爲使用cov(lm$resid)不會針對變量$ k $的數量進行調整，還是會這樣做？

我當然知道如何手動編程，但我只是好奇是否沒有內置函數。提前致謝。

Answer 1

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

是的，你的直覺是正確的; cov(lm$resid)不會針對自由度進行調整。