如何提取r中稳健的标准误差？

Question

在使用lm()函数运行回归之后，我计算了健壮的标准错误。

# robust standard errors
cov2I         <- vcovHC(ols2I, type = "HC1")
robust_se2I    <- sqrt(diag(cov2I))
print(robust_se2I)

我想从结果矩阵中提取第二个值，并将其保存在一个新变量下。我尝试过下面的代码，但没有成功。

stderrorols2I <- (summary(robust_se2I))[2]

谢谢你的帮忙!

Answer 1

这就是如何手动将RSE添加到摘要输出中。另外，你可以看看coeftest()

library(sandwich)  
mod1 <- lm(mpg ~ cyl + disp, data = mtcars) 

# robust standard errors
cov2I         <- vcovHC(mod1, type = "HC1")
robust_se2I    <- sqrt(diag(cov2I)) 

mod1 %>% 
  broom::tidy() %>% 
  mutate(rse = robust_se2I)

Answer 2

我找到了我问题的答案。

为了保存特定的健壮标准错误，我应该编写以下代码：

stderrorols2I <- (robust_se2I)[2]

现在它运转得很好。无论如何，谢谢你的快速反馈！