sysuse nlsw88 , clear statsby "ci union , binomial" /// proportion=(r(mean)) lb=(r(lb)) ub=(r(ub)) , /// by(occupation) clear label variable lb "" label variable ub "" label variable proportion "" twoway bar proportion occupation twoway bar proportion occupation , /// barwidth(.75) /// xlabel(1(1)13 , valuelabel angle(45)) /// graphregion(margin(l=12)) twoway bar proportion occupation , barwidth(.75) || /// rcap lb ub occupation , /// xlabel(1(1)13 , valuelabel angle(45)) /// graphregion(margin(l=12)) twoway bar proportion occupation , horiz barwidth(.75) || /// rcap lb ub occupation , horiz /// ylabel(1(1)13 , valuelabel angle(0) nogrid) /// ytitle("") xlabel(,grid) /// legend(label(1 "proportion")) twoway scatter occupation proportion || /// rcap lb ub occupation , horiz /// ylabel(1(1)13 , valuelabel angle(0) nogrid) /// ytitle("") xlabel(,grid) /// legend(label(1 "proportion"))