代码如下:
gen treat = 0
replace treat = 1 if 地区=="上海市" | 地区=="浙江省" | 地区=="广东省" | 地区=="深圳市" | 地区=="江苏省" | 地区=="山东省" | 地区=="北京市" | 地区=="江西省" | 地区=="宁夏回族自治区" | 地区=="青岛市"
gen pd=时间-2014
replace pd = -3 if pd <=-3
forvalues i = 3(-1)1{
gen pre_`i'=(pd== -`i' & treat==1)
}
gen current =(pd==0&treat==1)
forvalues j = 1(1)5{
gen las_`j'=(pd==`j'& treat==1)
}
xtset id 时间
xtreg industry pre_3 pre_2 current las_1 las_2 las_3 las_4 las_5 pre_1 $CV i.时间, fe r
est sto reg
coefplot,baselevels omitted keep(pre* current las*) vertical recast(connect) order( pre_2 pre_1 current las_1 las_2 las_3 las_4 ) yline(0,lp(solid) lc(black)) xline(3,lp(solid)) ytitle("政策效应") xtitle("政策时点") xlabel(1 "2012" 2 "2013" 3 "2014" 4 "2015" 5 "2016" 6 "2017" 7 "2018" )ciopts(recast(rcap) lc(black) lp(dash) lw(thin)) scale(1.0)


雷达卡







京公网安备 11010802022788号







