从 feols() IV 回归中提取 modelsummary() 的第一阶段结果
Extract first stage results from feols() IV regression for modelsummary()
我想使用 modelsummary()
报告 feols()
IV 回归的第一阶段和第二阶段结果。我找不到办法(除了 运行 第一阶段作为一个单独的模型)。
我可以像这样使用 etable()
显示第一阶段和第二阶段的结果:
library(fixest)
library(tidyverse)
library(modelsummary)
# create a toy dataset
base <- iris
names(base) <- c("y", "x1", "x_endo_1", "x_inst_1", "fe")
base$x_inst_2 <- 0.2 * base$y + 0.2 * base$x_endo_1 + rnorm(150, sd = 0.5)
base$x_endo_2 <- 0.2 * base$y - 0.2 * base$x_inst_1 + rnorm(150, sd = 0.5)
# estimate an instrumental variable model
mod <- feols(y ~ x1 | fe | x_endo_1 + x_endo_2 ~ x_inst_1 + x_inst_2, base)
# First and second stage results
etable(mod, stage = 1:2)
如有指点,我将不胜感激。
谢谢,
乌穆特
可以通过以下方式生成第一阶段 IV 回归模型摘要:
modelsummary::modelsummary(mod$iv_first_stage)
来自 fixest documentation:
要获得第一阶段和第二阶段,请使用 summary
和 stage = 1:2
。然后你可以把它喂给modelsummary
:
modelsummary(summary(mod, stage = 1:2))
我想使用 modelsummary()
报告 feols()
IV 回归的第一阶段和第二阶段结果。我找不到办法(除了 运行 第一阶段作为一个单独的模型)。
我可以像这样使用 etable()
显示第一阶段和第二阶段的结果:
library(fixest)
library(tidyverse)
library(modelsummary)
# create a toy dataset
base <- iris
names(base) <- c("y", "x1", "x_endo_1", "x_inst_1", "fe")
base$x_inst_2 <- 0.2 * base$y + 0.2 * base$x_endo_1 + rnorm(150, sd = 0.5)
base$x_endo_2 <- 0.2 * base$y - 0.2 * base$x_inst_1 + rnorm(150, sd = 0.5)
# estimate an instrumental variable model
mod <- feols(y ~ x1 | fe | x_endo_1 + x_endo_2 ~ x_inst_1 + x_inst_2, base)
# First and second stage results
etable(mod, stage = 1:2)
如有指点,我将不胜感激。
谢谢, 乌穆特
可以通过以下方式生成第一阶段 IV 回归模型摘要:
modelsummary::modelsummary(mod$iv_first_stage)
来自 fixest documentation:
要获得第一阶段和第二阶段,请使用 summary
和 stage = 1:2
。然后你可以把它喂给modelsummary
:
modelsummary(summary(mod, stage = 1:2))