1
这应该很容易,但是R中是否有一个内置函数来获得lm对象的方差协方差模型的无偏估计?然后我的意思是$ \ epsilon \ epsilon'/(n-k-1)$?因为使用cov(lm$resid)
不会针对变量$ k $的数量进行调整,还是会这样做?协方差矩阵lm对象R
我当然知道如何手动编程,但我只是好奇是否没有内置函数。提前致谢。
这应该很容易,但是R中是否有一个内置函数来获得lm对象的方差协方差模型的无偏估计?然后我的意思是$ \ epsilon \ epsilon'/(n-k-1)$?因为使用cov(lm$resid)
不会针对变量$ k $的数量进行调整,还是会这样做?协方差矩阵lm对象R
我当然知道如何手动编程,但我只是好奇是否没有内置函数。提前致谢。
我相信,假设拟合模型lm1
,你只需要summary(lm1)$sigma
。
是的,你的直觉是正确的; cov(lm$resid)
不会针对自由度进行调整。