在使用DID方法时，如何在Stata软件中用命令画出类似如下所示的图？

http://s5/bmiddle/003z4oOvty6HP0ggaTa44&690

以下是博主所编的一个Stata程序，希望对读者有用。程序的基本思路是：先生成一个10乘4的矩阵，使用循环语句将期数、回归系数、上下95%的置信区间输出到矩阵中，再将该矩阵输出到txt文档中，最后将txt文档导入Stata软件中进行画图操作。

************以下程序可以将回归系数及上下95%的置信区间输出并绘图***************

  matrix A=J(10,4,0)

  forvalues i=1/10{

    local j=`i'+1

    local k=`i'-5

    scalar a`i'=_b[_IntrXfkin_`j'] 

    scalar b`i'=(_b[_IntrXfkin_`j']+invttail(e(df_r),0.025)*_se[_IntrXfkin_`j'])

    scalar c`i'=(_b[_IntrXfkin_`j']-invttail(e(df_r),0.025)*_se[_IntrXfkin_`j'])

    mat A[`i',1]=`k'

    mat A[`i',2]=a`i'

    mat A[`i',3]=b`i'

    mat A[`i',4]=c`i'

}

  mat2txt,matrix(A) saving(DID_coe_ci) replace 

  insheet using DID_coe_ci.txt,clear

  keep c1 c2 c3 c4

  rename c1 year

  rename c2 coe

  rename c3 up

  rename c4 low

  save DID_coe_ci.dta,replace

  twoway (connected coe year,msymbol(T)) (line up year,lpattern(shortdash)) (line low /// year,lpattern(longdash)), xline(0,lc(black)lp(dash)) yline(0,lc(black)lp(solid)lw(medthick)) /// xscale(range(-4 5)) xlabel(-4(1)5) xtitle("距离政策实施的时间") ytitle("回归系数:XX变动的 ///

百分比") legend(label(1 "估计值") label(2 "95%的上界") label(3 "95%的下界")) scale(0.8) /// saving(DID_coe_ci,replace)

*******************************************************************************

说明：

    本程序中，有-5、-4......0......4、5共11期，以-5期为基准组，共生成从_IntrXfkin_2到_IntrXfkin_11的10个交互项，实际操作中只需将其修改成自己回归中使用的变量名即可。

    本程序需在跑完DID的回归后使用。

mat2txt 是外部命令。

图片出处：Wang, J., 2013, "The Economic Impact of Special Economic Zones: Evidence from Chinese Municipalities", Journal of Development Economics, 101, 133-147.