在 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

使用摘要（rd）。这将 return 系数估计

Answer 2

不幸的是，我不相信 rdrobust() 允许您恢复通过 covs 选项引入的系数。

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

rd$coef

只会为您提供 rd 估计器的点估计。

Josh McCrain 有 written-up 一个很好的小插图 here to replicate rdrobust using lfe，它还允许您恢复协变量的系数。

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

这可能已经离题了，但 RD 设计中的期刊要求很奇怪。