在 rdrobust 包中获取协变量的估计值

Question

我正在使用 rdrobust 来估计 RDD，并且对于在期刊上提交的文章，该期刊要求我报告包含协变量及其估计值的表格。我认为这些不应该在这样的设计中报告，并且真的不知道它们的信息量有多大，但无论如何：我在 rdrobust 调用的输出中找不到它们，所以我想知道是否有实际获得它们。

这是我的代码：

library(rdrobust)
rd <- rdrobust(y = full_data$share_female,
         x = full_data$running,
         c = 0,
         cluster = full_data$constituency.name,
         covs=cbind(full_data$income, full_data$year_fct,
                    full_data$population, as.factor(full_data$constituency.name)))

然后我调用该对象

rd

并得到：

Call: rdrobust

Number of Obs.                 1812
BW type                       mserd
Kernel                   Triangular
VCE method                       NN

Number of Obs.                 1452          360
Eff. Number of Obs.             566          170
Order est. (p)                    1            1
Order bias  (q)                   2            2
BW est. (h)                   0.145        0.145
BW bias (b)                   0.221        0.221
rho (h/b)                     0.655        0.655
Unique Obs.                    1452          360

因此，正如您所看到的，输出和函数调用的对象上似乎都没有关于此的信息。我真的不知道该怎么办。 谢谢！

Answer 1

不幸的是，我不相信

rdrobust()

允许您恢复通过

covs

选项引入的系数。

在您的情况下，运行您提供的代码，然后运行：

rd$coef

只会给你 rd 估计器的点估计。

Josh McCrain here 写了一个很好的小插图，使用

rdrobust

 复制 

lfe，这也允许您恢复协变量的系数。

它涉及您的一些修改，当然不是用户友好的，但确实允许恢复协变量。

现在这可能不是重点，但 RD 设计中的期刊要求很奇怪。

Answer 2

在当前版本的

rdrobust

(2.2) 中，您可以通过输出的

beta_covs

属性恢复协变量的系数：

rd$beta_covs

Answer 3

使用摘要(rd)。这将返回系数估计值