我是单期did，研究时间框为2009-2019，研究2014年碳交易试点的影响，因为用随机抽取的方法安慰剂检验不太好看，所以改用permute，

这是我的代码：permute did beta = _b[did] se = _se[did] df = e(df_r),  ///          reps(1000) rseed(123) saving("simulations.dta"):  ///
          reg lnggdp did lnpgdp lnfdl lnop lntc lnis lnpop lne i.year i.province, robust
use "simulations.dta", clear
gen t_value = beta / se
gen p_value = 2 * ttail(df, abs(beta/se))
twoway (scatter p_value beta, msymbol(smcircle_hollow) mcolor(black))       (kdensity beta ,yaxis(2) lp(solid)) , ///
       title("") ///
       xlabel(-0.4(0.05)0.4 -0.05 "-0.05" -0.1 "-0.1" -0.15 "-0.15" -0.2"-0.2" -0.25"-0.25" -0.3"-0.3" -0.35"-0.35" -0.4"-0.4" 0.05 "0.05" 0.1 "0.1" 0.15 "0.15" 0.2"0.2" 0.25"0.25" 0.3"0.3" 0.35"0.35" 0.4"0.4" ,format(%7.1f) angle(0)) ///
       ylabel(0(0.2)1, format(%7.1f) angle(0) nogrid axis(1)) ///
       ylabel(0(2)10, format(%7.1f) angle(0)  nogrid axis(2)) ///
       xtitle("回归系数") ///
       ytitle("P" "值" ,orientation(horizontal) axis(1)) ///
       ytitle("核" "密" "度" ,orientation(horizontal)  axis(2)) ///
       xline(0, lwidth(0.2) lp(dash))  ///
       xline(0.3847, lwidth(0.3) lp(solid)) ///
       yline(0.05,lwidth(0.2) lp(dash)) ///
       legend(label(1 "P值") label( 2 "核密度")) ///
       plotregion(style(none)) ///
       graphregion(color(white)) ///
saving(安慰剂检验图.png,replace)

请教大佬我的代码有问题吗？这个permute命令得出的安慰剂检验结果图算通过吗？