如何将三明治调整后的 p 值摘要转换为数据框？

Question

让我们考虑以下模型，它是调整后 pvalues 的 HAC 版本：

library(lmtest)
library(sandwich)
set.seed(42)
y <- rnorm(100)
x <- data.frame(runif(100))
model <- lm(y~., data = x)
coeftest(model, vcov = vcovHAC)

t test of coefficients:

            Estimate Std. Error t value Pr(>|t|)
(Intercept) -0.09510    0.16125 -0.5898   0.5567
runif.100.   0.28952    0.30452  0.9508   0.3441

我想知道-如何将其转换为数据框？在标准 lm 摘要中，我们只使用命令as.data.frame(coef(summary))，但在这里它似乎不起作用。我该怎么做？

Answer 1

broom::tidy 可以在这里提供帮助。

data <- coeftest(model, vcov = vcovHAC)
broom::tidy(data)

#  term        estimate std.error statistic p.value
#  <chr>          <dbl>     <dbl>     <dbl>   <dbl>
#1 (Intercept)  -0.0951     0.161    -0.590   0.557
#2 runif.100.    0.290      0.305     0.951   0.344