R 中 rdrobust 包的回归输出表

Question

我对 R 比较陌生，目前正在尝试根据我的稳健结果制作回归表。这是我的目标的一个示例：

我知道 rdrobust modelsummary() 不支持，我必须使用 tidy.rdrobust 和 glance.rdrobust 才能使其适用。所以我的问题是如何在这两个命令（tidy，glance）中指定参数以使用 modelsummary 包获得所需的输出？ 如图所示，我希望将结果变量/模型作为行，将三种方法的系数（带*）、SE、z 值和置信区间作为列。

正如我所说，我对 R 比较陌生，而且对一般编程语言也完全陌生。因此，这个问题可能不够具体，并且没有进行任何编码试验。也许无论如何有人可以帮助我。

多谢！

Answer 1

即使它不能满足您的所有要求，您也可以利用 modelsummary 的功能进行并排模型：https://vincentarelbundock.github.io/modelsummary/articles/modelsummary.html

library(rdrobust)
library(modelsummary)

tidy.rdrobust <- function(model, ...) {
  ret <- data.frame(
    term = row.names(model$coef),
    std.error = model$se[, 1],
    p.value = model$pv[, 1]
  )
  row.names(ret) <- NULL
  ret
}

glance.rdrobust <- function(model, ...) {
  ret <- data.frame(
    Kernel = model$kernel,
    Bandwidth = model$bwselect
  )
  ret
}

如果我们用虚拟实现数据，它将显示以下内容：

x1 <- runif(1000, -1, 1)
x2 <- runif(1000, -1, 1)
y1 <- 5 + 3 * x1 + 2 * (x1 >= 0) + rnorm(1000)
y2 <- 5 + 3 * x2 + 2 * (x2 >= 0) + rnorm(1000)

fit1 <- rdrobust(y1, x1)
fit2 <- rdrobust(y2, x2)
models <- list(fit1, fit2)
modelsummary(models, statistic = "p.value")

	模型 1	模型 2
系数	X	X
标准误差	X	X