R中的第一级IV回归输出

Question

我想报告我的两阶段最小二乘回归法两个阶段的结果，但占星师的输出仅给我第二阶段的结果。

我已经用R中的ivreg命令计算了两阶段最小二乘回归。这是我的代码的样子：

ivmodel1 <- ivreg(Y ~ X + W1  + W2 + W3 + W4 | W1  + W2 + W3 + W4 + Z, data = df)

其中

Y =因变量（续）；

X =内生自变量（虚拟）；

W1-W4 =控制变量；

Z =外来工具（虚拟）

现在，我很难报告2SLS回归的第一阶段。当我使用普通的stargazer命令时：

stargazer(ivmodel1)

我只收到第二阶段的结果，但我也需要第一阶段的估计。有人知道在R中使用什么命令以接收两个阶段的结果吗？

Answer 1

library(AER)
library(stargazer)

y <- rnorm(100, 5, 10)
x <- rnorm(100, 3, 15)
z <- rnorm(100, 3, 7)
a <- rnorm(100, 1, 7)
b <- rnorm(100, 3, 5)

# Fitting IV models
fit1 <- ivreg(y ~ x + a  |
                a + z,
              model = TRUE)
fit2 <- ivreg(y ~ x + a  |
                a + b + z,
              model = TRUE)

# Create latex table
stargazer(fit1, fit2, type = "text")

来自：R: Robust SE's and model diagnostics in stargazer table