我想报告feols()使用 IV 回归的第一阶段和第二阶段结果modelsummary()。我找不到方法（除了将第一阶段作为单独的模型运行）。
我可以使用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 …