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"))