* Effects of Hurricane Katrina on BLS; * Data from http://www.bls.gov/katrina/cpsesquestions.htm; data employ; input month $ employment $ residence $ count; count=count*1000; datalines; november employed same 204 november employed different 185 november unemployed same 29 november unemployed different 71 november notinLabor same 209 november notinLabor different 188 october employed same 153 october employed different 179 october unemployed same 18 october unemployed different 90 october notinLabor same 134 october notinLabor different 217 proc freq data=employ order=data; weight count; tables month*residence*employment /chisq nopercent cmh measures cmh2; run; proc freq data=employ order=data; weight count; tables employment*month*residence /chisq nopercent cmh measures cmh2; run; proc logistic order=data; weight count; class residence month; model employment (ref=first) = residence month residence*month/ link=glogit ; run; proc genmod data=employ order=data; class month employment residence ; model count = employment month residence employment*month employment*residence residence*employment / link=log dist=poi ; run;