我使用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
因此,正如您所看到的,在输出和函数调用的对象上似乎都没有关于这方面的信息。我真的不知道该怎么办。谢谢
不幸的是,我不认为rdrobust()
允许您恢复通过covs
选项引入的系数。
在您的情况下,按照您提供的方式运行代码,然后运行:
rd$coef
只会给你rd估计器的点估计。
Josh McCrain在这里写了一个很好的小插曲,用lfe
复制rdrobust
,也可以恢复协变量上的系数。
它涉及到您的一些修改,当然不是用户友好的,但确实允许恢复协变量。
到目前为止,这可能已经不重要了,但研发设计中的日志要求很奇怪。
使用摘要(rd(。这将返回系数估计