Hi,
For illustration purposes I would like to plot a skew normal distribution to add to the following plot:

Code:
twoway function normalden(x, 80, 10), ra(-10 200) || function normalden(x, 80, 20), ra(-10 200)  ///
, legend(order(1 "{&sigma} = 10" 2 "{&sigma} = 20") row(1)) ///
xline(80, lcolor(grey) lpattern(dot))
Is there a way to add a skewness parameter to the density function in order to allow for non-zero skewness? Or a handy workaround?

Any suggestions are much appreciated!