最近我们在 SO 上讨论了如何使用 NeweyWest 标准误差更新标准线性回归摘要。我coeftest从sandwich包装中使用。有人告诉它使用 unclass 来更新我已经存在的摘要,如下所示:
library(sandwich)
library(lmtest)
temp.lm <- lm(runif(100) ~ rnorm(100))
temp.summ <- summary(temp.lm)
temp.summ$coefficients <- unclass(coeftest(temp.lm, vcov. = NeweyWest)
现在我想知道在使用 NeweyWest VC 矩阵时,摘要中显示的关节参数是否完全不受影响?我的意思是使用这段代码,它们显然不会受到影响——但这是正确的吗?请注意,这不是语法,而是统计问题:) 类似的东西
Residual standard error: 1.177 on 83 degrees of freedom
Multiple R-squared: 0.7265, Adjusted R-squared: 0.71
F-statistic: 44.1 on 5 and 83 DF, p-value: < 2.2e-16
保持不变。是否也有需要调整的情况?