代码如下:
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