/*------------ CHD data ------------*/ options ls=80 nocenter; data chd; input CHD $ serum $ count xCHD yserum @@; datalines; chd 0-199 12 1 1 chd 200-199 8 1 2 chd 220-259 31 1 3 chd 260+ 41 1 4 nochd 0-199 307 2 1 nochd 200-199 246 2 2 nochd 220-259 439 2 3 nochd 260+ 245 2 4 ; proc freq; weight count; tables CHD*serum/chm1 scores=ridits; title 'CHM statstic for linear association'; run; proc genmod data=chd order=data; class CHD serum; model count=CHD serum/link=log dist=poisson lrci type3 obstats; title 'Indepedence Model'; run; proc genmod data=chd order=data; class CHD serum; model count=CHD serum xCHD*yserum/link=log dist=poisson lrci type3 obstats; title 'Linear by Linear Association Model'; run;