stata实证代码答疑

Question

如何解决该有关RDD回归不连续设计stata代码中的报错，代码如下：

* set directories 
clear all
set maxvar  30000
set matsize 11000 
set more off
cap log close


use "$data/abatement_QJE_final.dta",clear
 use "/Users/a12420/Documents/sysu/卫生经济/Replication Materials/T4_Abatement/abatement_QJE_final.dta"
cap erase $result/T_abatement.txt 
cap erase $result/T_abatement.xml     


********Abatement in production process********
** reduced production hours 
rdrobust resid1_hours distance_new,masspoints(off) kernel(tri) all vce(cluster site_id)
esttab using "$result/T_abatement", type=excel cell(d) label replace varlist(resid prod hours Kernel Type) ///
    nonumber format(%9.0f %9.0g %9.0g %s) addnotes("Kernel Type: " e(kernel))
outreg2 using $result/T_abatement, excel addtext(Kernel Type, `e(kernel)') ctitle(resid prod hours) dec(0) replace
rdrobust resid1_hours distance_new,masspoints(off) kernel(epa) all vce(cluster site_id)
outreg2 using $result/T_abatement, excel addtext(Kernel Type, `e(kernel)') ctitle(resid prod hours) dec(0) append
rdrobust resid1_hours distance_new,masspoints(off) kernel(uni) all vce(cluster site_id)
outreg2 using $result/T_abatement, excel addtext(Kernel Type, `e(kernel)') ctitle(resid prod hours) dec(0) append

报错如下

. do "/var/folders/db/t4jw7f151271gdx8hy_sxj2m0000gn/T//SD31847.000000"

. rdrobust resid1_hours distance_new,masspoints(off) kernel(tri) all vce(cluster site_id)

Sharp RD estimates using local polynomial regression.

      Cutoff c = 0 | Left of c  Right of c            Number of obs =       7302
-------------------+----------------------            BW type       =      mserd
     Number of obs |      4357        2945            Kernel        = Triangular
Eff. Number of obs |      1223         965            VCE method    =    Cluster
    Order est. (p) |         1           1
    Order bias (q) |         2           2
       BW est. (h) |     5.222       5.222
       BW bias (b) |     8.999       8.999
         rho (h/b) |     0.580       0.580
Number of clusters |       141         113

Outcome: resid1_hours. Running variable: distance_new.
--------------------------------------------------------------------------------
            Method |   Coef.    Std. Err.    z     P>|z|    [95% Conf. Interval]
-------------------+------------------------------------------------------------
      Conventional |   288.3     101.38   2.8436   0.004    89.5892      487.003
    Bias-corrected |  331.83     101.38   3.2730   0.001    133.123      530.537
            Robust |  331.83     115.03   2.8848   0.004    106.383      557.277
--------------------------------------------------------------------------------
Std. Err. adjusted for clusters in site_id

. 
end of do-file

. do "/var/folders/db/t4jw7f151271gdx8hy_sxj2m0000gn/T//SD31847.000000"

. outreg2 using $result/T_abatement, excel addtext(Kernel Type, `e(kernel)') ctitle(resid prod hours) dec(
> 0) replace
file /T_abatement.txt could not be opened
r(603);

end of do-file

r(603);

Answer 1

阿里嘎多学长整理AIGC生成，因移动端显示问题导致当前答案未能完全显示，请使用PC端查看更加详细的解答过程

stata实证代码答疑

你遇到了Stata代码中的报错，可能是由于RDD回归不连续设计中的一些问题。下面是一些可能的解决方案：

1. 检查数据是否有缺失值：RDD回归需要数据的连续性，如果数据中存在缺失值，可能会导致报错。你可以使用summarize命令来检查数据的缺失值。
2. 检查变量的数据类型：RDD回归需要变量的连续性，如果变量的数据类型不是连续的，可能会导致报错。你可以使用describe命令来检查变量的数据类型。
3. 检查回归模型的 specification：RDD回归需要正确的回归模型 specification，如果模型 specification 不正确，可能会导致报错。你可以检查模型 specification 是否正确。

以下是一个简单的示例代码，用于解决RDD回归不连续设计中的报错：

* set directories 
clear all
set maxvar  30000
set matsize 11000

* 读取数据
use "your_data.dta", clear

* 检查数据是否有缺失值
summarize *

* 检查变量的数据类型
describe *

* 进行RDD回归
rdrobust y x1 x2, cluster(id)

请注意，这只是一个简单的示例代码，你需要根据你的实际情况进行修改。

如果你仍然遇到报错，可以提供更多的代码和数据信息，我将尽力帮助你解决问题。