/*---------------------------- | Example: Friench skiiers -----------------------------*/ data ski; input treatment $ response $ count; datalines; placebo cold 31 placebo nocold 109 absorbic cold 17 absorbic nocold 122 ; run; * independence model; proc catmod data=ski order=data; weight count; model treatment*response=_response_; loglin treatment response; run; proc genmod data=ski order=data; class treatment response; model count = treatment response /link=log dist=poisson lrci type3 obstats; run; * saturated model; proc catmod data=ski order=data; weight count; model treatment*response=_response_; loglin treatment|response; run; proc genmod data=ski order=data; class treatment response; model count = response*treatment /link=log dist=poisson lrci type3 obstats; run; proc freq; weight count; tables treatment*response/ chisq relrisk riskdiff expected; exact or; run; proc genmod data=ski order=data; class treatment response; model count = treatment response response*treatment/link=log dist=poisson lrci type3 obstats; run;