http://michaelewens.com/2009/03/ ... -one-plot-in-stata/
sysuse auto,clear

qui sum price if foreign==0
global mean1=r(mean)
qui sum price if foreign==1
global mean2=r(mean)
#d;
kdensity price if foreign==0, xline($mean1)
addplot(kdensity price if foreign==1,  
xline($mean2,lpattern("-")) )
legend(ring(0) pos(2) label(1 "Domestic")
label(2 "Foreign")
region(lcolor(white)))
scheme(s2mono)
;
#d cr

谢谢！

不客气

学习了

学习