【发布时间】:2021-03-22 21:26:48
【问题描述】:
我正在运行带有协变量选项的 rdrobust 包中的 rdrobust 函数。我正在使用作者的指南(“回归不连续设计的实用介绍”2019 Cattaneo,Idboro,Titiunik)。在 R sn-p 26 的代码是这样的:
# R Snippet 26
# Using rdrobust with covariates
Z = cbind(data$vshr_islam1994, data$partycount, data$lpop1994,data$merkezi, data$merkezp, data$subbuyuk, data$buyuk)
colnames(Z) = c("vshr_islam1994", "partycount", "lpop1994","merkezi", "merkezp", "subbuyuk", "buyuk")
out = rdrobust(Y, X, covs = Z, kernel = 'triangular', scaleregul = 1, p = 1, bwselect = 'mserd')
根据作者的说法,Z 上应该有一组系数 gamma(第 79 页指南中的方程 4.3)但是在输出返回中,每个变量有两个(所以它们是 p+1 的顺序)并且它们左右两侧也不是对称的,因为它们应该是:
> out$beta_p_r
buyuk merkezi merkezp lpop1994 subbuyuk partycount vshr_islam1994
[1,] 15.7202048 0.02413397 0.42115723 0.066234600 8.285038847 0.054458918 5.80939355 33.6301014
[2,] -0.1623748 -0.00151349 -0.00525247 -0.001610547 -0.005106353 0.003832462 0.01467446 0.7929476
> out$beta_p_l
buyuk merkezi merkezp lpop1994 subbuyuk partycount vshr_islam1994
[1,] 12.8148371 0.017154325 0.478279531 0.0300389537 8.24125708 0.049303686 5.854047932 33.1039297
[2,] -0.2058053 0.001048017 0.009266473 0.0006171706 0.02126247 0.001227102 0.003424324 0.7899774
这些系数是多少?我只知道第一个未命名的向量对应于回归不连续性的局部线性函数(例如 p=1)。但其余的都是一个谜,它们与上述指南不符。我担心协变量的估计不是它的描述方式,如果估计方程与指南中描述的不同,我想知道估计方程是什么。
【问题讨论】: