如何在R中从modelsummary中的仪器变量中获取弱仪器测试？

mwg9r5ms 于 2023-04-03 发布在其他

关注(0)|答案(1)|浏览(133)

如何获得modelsummary中的弱仪器F检验统计量？

data(mtcars)
library(ivreg)

iv_model <- ivreg(mpg ~ qsec + cyl + drat | disp | wt, data = mtcars)
summary(iv_model, diagnostics = TRUE)

Call:
ivreg(formula = mpg ~ qsec + cyl + drat | disp | wt, data = mtcars)

Coefficients:
            Estimate Std. Error t value Pr(>|t|)   
(Intercept) 23.28560   20.84029   1.117  0.27370   
disp        -0.05730    0.02053  -2.791  0.00953 **
qsec         0.20443    0.59223   0.345  0.73263   
cyl          0.88477    1.52033   0.582  0.56542   
drat         0.25095    2.19015   0.115  0.90962   

Diagnostic tests:
                 df1 df2 statistic  p-value    
Weak instruments   1  27     19.96 0.000127 ***
Wu-Hausman         1  26     13.87 0.000956 ***
Sargan             0  NA        NA       NA    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

默认情况下，metrics = "all"选项显示Wu-Hausman和Sargan检验。

library(modelsummary)
modelsummary(iv_model, metrics = "all")

r

来源：https://stackoverflow.com/questions/75754280/how-to-get-weak-instrument-test-from-instrument-variables-in-modelsummary-in-r

1条答案

按热度按时间

kh212irz1#

modelsummary包含一个向表中添加任何拟合优度统计信息的简单机制。有关详细说明和示例，请在此页上搜索glance_custom：https://vincentarelbundock.github.io/modelsummary/articles/modelsummary.html
在ivreg的特定情况下，贡献者只是向performance和modelsummary包添加了功能，以使其更容易。安装这两个包的开发版本：

library(remotes)
install_github("easystats/performance")
install_github("vincentarelbundock/modelsummary")

完全重新启动R以使更改生效。然后，

library(ivreg)
library(modelsummary)

iv_model <- ivreg(mpg ~ qsec + cyl + drat | disp | wt, data = mtcars)

gm <- transform(modelsummary::gof_map, omit = FALSE)
modelsummary(iv_model, gof_map = gm, metrics = "all")

	(1)
（拦截）	二十三点二八六
	（20.840英镑）
分散	-0.057
	（0.021）
质量秒	零点二零四
	（0.592）
圆柱体	0.885
	（1.520）
德拉特	零点二五一
	（2.190美元）
Num.Obs.	三十二
R2	零点六五五
R2调整	零点六零四
AIC	一百八十二点七
BIC	一百九十一点五
均方根误差	三点四十八分
弱IV F状态	20.0年
弱IV p	千
西格玛	三千七百九十三
吴-豪斯曼	13.869
Wu-Hausman（p）	千分之一

赞(0）回复(0）举报 2023-04-03

我来回答

如何在R中从modelsummary中的仪器变量中获取弱仪器测试？

1条答案

相关问题

热门标签

最新问答