clear twoway (function y = sqrt(x * (1 - x)), range(0 1)), /// ytitle("Sample standard deviation") xtitle("Proportion") /// xlabel(0 (.25) 1, format(%6.2f)) /// name(figure1, replace)