如何在R中将rddensity函数的summary输出导出为LaTeX表格

由于modelsummary暂不支持rddensity对象的直接转换，我们可以手动提取summary(rdd)中的核心统计量，再用xtable或kableExtra生成规范的LaTeX表格。以下是具体实现步骤：

1. 从rddensity对象中提取数据 #

rddensity的返回对象包含了所有需要的统计信息，无需捕获打印输出，直接提取即可：

library(rdrobust)
library(rddensity)
library(xtable) # 用于生成LaTeX表格
set.seed(42)
x <- rnorm(2000, mean = -0.5)
rdd <- rddensity(X = x, vce = "jackknife")

# 提取基本信息
basic_info <- data.frame(
  "指标" = c("观测数", "模型类型", "核函数", "带宽方法", "方差估计方法"),
  "值" = c(rdd$N, rdd$model, rdd$kernel, rdd$bwselect, rdd$vce)
)

# 提取截断点两侧的估计统计量
side_info <- data.frame(
  "指标" = c("观测数", "有效观测数", "估计阶数(p)", "偏差阶数(q)", "带宽估计(h)"),
  "左侧" = c(rdd$N_left, rdd$N_left_eff, rdd$p, rdd$q, rdd$h_left),
  "右侧" = c(rdd$N_right, rdd$N_right_eff, rdd$p, rdd$q, rdd$h_right)
)

# 提取稳健检验结果
robust_test <- data.frame(
  "方法" = "Robust",
  "T统计量" = rdd$test$T,
  "P值" = rdd$test$p
)

# 提取二项式检验结果
binom_test <- as.data.frame(rdd$binomial)
colnames(binom_test) <- c("半窗口长度", "<c", ">=c", "P>|T|")

2. 生成并保存LaTeX表格 #

将各部分统计量生成独立表格，合并为完整的LaTeX文档并保存：

# 生成各模块的LaTeX表格代码
tab_basic <- xtable(basic_info, caption = "操纵性检验基本信息", align = "ll")
tab_side <- xtable(side_info, caption = "截断点两侧估计统计量", align = "lll")
tab_robust <- xtable(robust_test, caption = "稳健检验结果", align = "lll")
tab_binom <- xtable(binom_test, caption = "二项式检验结果(H0: p=0.5)", align = "lccc")

# 写入LaTeX文件
cat("\\documentclass{article}",
    "\\usepackage{booktabs}",
    "\\begin{document}",
    print(tab_basic, booktabs = TRUE, include.rownames = FALSE),
    print(tab_side, booktabs = TRUE, include.rownames = FALSE),
    print(tab_robust, booktabs = TRUE, include.rownames = FALSE),
    print(tab_binom, booktabs = TRUE, include.rownames = FALSE),
    "\\end{document}",
    file = "rddensity_summary.tex")

补充说明 #

• 运行代码后会生成rddensity_summary.tex文件，编译后即可得到格式化的表格。
• 若需要更灵活的样式调整，可替换xtable为kableExtra，例如用kable()生成表格后添加kable_styling(latex_options = "booktabs")优化格式。

内容的提问来源于stack exchange，提问作者John M. Riveros