/* GLM handout */
/* Partial program for getting Table 7.4, Table 7.6 */

data drugs;
input a c m count @@; 
cards;
1 1 1 911   1 1 2 538   1 2 1 44   1 2 2 456
2 1 1   3   2 1 2 43   2 2 1 2   2 2 2 279
;
/* model is (AM, AC, CM)  */
proc genmod; class a c m;
   model  count = a c m a*m a*c c*m / dist=poi link=log;
output out=expected  p=pred;

proc print  data=expected;  /* c.f. Table 7.6, p. 211 */


/* the estimated partial odds ratios, Table 7.5, p. 210 */   
proc freq data=expected; weight predicted;
   tables m*a*c / chisq expected measures;
   tables c*a*m / chisq expected measures;
   tables a*c*m / chisq expected measures;
/* the estimated marginal odds ratios, Table 7.5, p. 210 */
   tables a*c / chisq expected measures;
   tables a*m / chisq expected measures;
   tables c*m / chisq expected measures;

run;