options ls=90 nocenter nodate; data concerns; input gender $ health $ count; datalines; male sex 6 male mens - male healthy 49 male none 77 female sex 16 female mens 12 female healthy 29 female none 102 ; /*without an indicator variable*/ proc genmod data=concerns order=data; class gender health; model count = gender health gender*health /link=log dist=poisson lrci type3 obstats; title 'Saturated model without indicator'; /*without an indicator variable*/ proc genmod data=concerns order=data; class gender health; model count = gender health /link=log dist=poisson lrci type3 obstats; title 'Independence model'; run;