rename current d0 //先把current改个名称
(1)回归方程重新设定为:
xtreg y d_4 d_3 d_2 d0 d1 d2 d3 d4 d5 i.year 控制变量, fe r
gen t = invttail ( 自由度 ,0.05)
(2)手工绘图法 部分代码修改地方
*生成b_j的系数与置信区间(假设有4期,即distance最小值是-4;你可以手动调整下面命令的参数)
forvalues i=2/4{//跳过d_1
gen b_`i' = _b[d_`i']
gen se_b_`i' = _se[d_`i']
gen b_`i'LB = b_`i'- t * se_b_`i'
gen b_`i'UB = b_`i'+ t * se_b_`i'
}*生成bj的系数与置信区间(假设有5期,即distance最大值是5)
forvalues i=0/5{//循环起点改成0
gen b`i' = _b[d`i']
gen se_b`i' = _se[d`i']
gen b`i'LB = b`i'- t * se_b`i'
gen b`i'UB = b`i'+ t * se_b`i'
}
gen b =.
gen LB =.
gen UB =.*生成系数(政策前)
forvalues i=2/4{
replace b = b_`i' if distance == -`i'}*生成系数(政策后)
forvalues i=0/5{
replace b = b`i' if distance == `i'}*生成系数置信区间下限(政策前)
forvalues i=2/4{
replace LB = b_`i'LB if distance == -`i'}*生成系数置信区间下限(政策后)
forvalues i=0/5{
replace LB = b`i'LB if distance == `i'}*生成系数置信区间上限(政策前)
forvalues i=2/4{
replace UB = b_`i'UB if distance == -`i'}*生成系数置信区间上限(政策后)
forvalues i=0/5{
replace UB = b`i'UB if distance == `i'}