您好,我今天按照代码来跑,怎么横轴都是0.1了,您能帮我改下代码么,拜托了。还有加垂直线的代码是xline(0.2142)我直接加在代码里面就好了么,麻烦您给看下,救救小白吧。
search permute //先安装此命令
permute did beta = _b[did] se = _se[did] df = e(df_r), seed(10000) ///
reps(500) saving("安7检验.dta"): ///
reghdfe y did $c ,a(id year)
use "安7检验.dta" , clear
gen t_value = beta / se
gen p_value = 2 * ttail(df, abs(beta/se))
replace beta = beta
#delimit ;
twoway (kdensity beta, yaxis(1) lc(black) ) (scatter p_value beta, yaxis(2) m(Oh) mc(black)) ,
xline(0.2142)
xline(0, lc(black*0.5) lp(solid))
yline(0.05, lc(black*0.5) lp(dash) axis(2))
xtitle("Estimator", size(*0.8)) xlabel(, format(%4.1f) labsize(small))
ytitle("Density", size(*0.8)) ylabel(, nogrid format(%4.1f) labsize(small))
ytitle("P Value", size(*0.8) axis(2)) ///
ylabel(, nogrid format(%4.1f) labsize(small) axis(2))
legend(ring(2) order(1 "Estimator" 2 "P Value") pos(6) rows(1) )
graphregion(color(white)) scheme(s1mono) ;
#delimit cr
graph export "安检验.png", width(1000) replace