twoway bar  proportion occupation , barwidth(.75)  ||	///
       rcap lb ub occupation ,				///
	xlabel(1(1)13 , valuelabel angle(45))		///
	graphregion(margin(l=12))