如何将三明治调整后的 p 值摘要转换为数据框？
How to convert sandwich adjusted p values summary to data frame?
让我们考虑以下模型，它是调整后的 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)) 但在这里它似乎不起作用。我该怎么做？
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
让我们考虑以下模型，它是调整后的 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)) 但在这里它似乎不起作用。我该怎么做？
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