如何解决如何将夹心调整后的 p 值摘要转换为数据框?
让我们考虑以下模型,它是调整后的 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
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。