R工具变量回归（2SLS）诊断结果的模型报告工具问询

• texreg包
  支持直接提取ivreg模型的诊断结果，可通过screenreg()（纯文本输出）或texreg()（LaTeX输出）生成报告，配合自定义参数整合诊断统计量。示例代码：

  library(texreg)
  # 先获取带诊断结果的模型摘要
  sum_iv <- summary(model_iv, diagnostics = TRUE)
  # 提取关键诊断值
  weak_instr_f <- sum_iv$diagnostics["Weak instruments", "F statistic"]
  wu_hausman_p <- sum_iv$diagnostics["Wu-Hausman", "Pr(>F)"]
  sargan_p <- sum_iv$diagnostics["Sargan", "Pr(>chisq)"]
  
  # 生成带诊断统计量的回归报告
  texreg(model_iv, 
         custom.gof.rows = list(
           "弱工具变量F统计量" = weak_instr_f,
           "Wu-Hausman检验p值" = wu_hausman_p,
           "Sargan检验p值" = sargan_p
         ),
         caption = "2SLS回归结果及诊断统计量")
  

  该包适配学术场景，输出格式规范，支持批量模型对比。

• huxtable包
  主打高度自定义表格构建，可分别提取模型系数和诊断结果，再自由组合成统一表格。示例代码：

  library(huxtable)
  sum_iv <- summary(model_iv, diagnostics = TRUE)
  
  # 提取系数表和诊断表
  coef_df <- as.data.frame(coef(summary(model_iv)))
  diag_df <- as.data.frame(sum_iv$diagnostics)
  
  # 构建并合并表格
  coef_hux <- hux(coef_df, add_colnames = TRUE)
  diag_hux <- hux(diag_df, add_colnames = TRUE)
  final_table <- rbind(coef_hux, diag_hux)
  
  # 设置表格样式和标题
  final_table <- final_table %>%
    set_caption("2SLS回归结果及诊断") %>%
    set_all_borders(1)
  print(final_table)
  

  支持输出HTML、LaTeX、Word等多种格式，适合需要精细化调整表格样式的场景。

• gtsummary包
  虽默认不直接适配ivreg的诊断模块，但可通过add_stat()函数自定义添加诊断统计量，输出表格美观易读。示例代码：

  library(gtsummary)
  sum_iv <- summary(model_iv, diagnostics = TRUE)
  
  # 生成基础回归表格
  tbl <- tbl_regression(model_iv)
  # 添加诊断统计量到表格底部
  tbl <- tbl %>%
    add_stat(
      fns = list(
        "弱工具变量F统计量" = ~ sum_iv$diagnostics["Weak instruments", "F statistic"],
        "Wu-Hausman检验p值" = ~ sum_iv$diagnostics["Wu-Hausman", "Pr(>F)"],
        "Sargan检验p值" = ~ sum_iv$diagnostics["Sargan", "Pr(>chisq)"]
      ),
      location = "bottom"
    )
  # 输出表格
  tbl
  

  生成的表格支持一键导出为Word、HTML等格式，适合快速生成报告。

内容的提问来源于stack exchange，提问作者Marco