statsby "ci union , binomial"				///
	proportion=(r(mean)) lb=(r(lb)) ub=(r(ub)) , 	///
	by(occupation) clear