A.10  SAS file for loglinear model (AM, CM) with Table 6.3

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
;
proc genmod; class a c m;
   model  count = a c m a*m  c*m / dist=poi link=log obstats;
proc catmod; weight count;
   model a*c*m = _response_; loglin a|m c|m;
run;
              The SAS System           15:16 Tuesday, May 18, 1999   1


                        The GENMOD Procedure

                         Model Information

             Description                     Value

             Data Set                        WORK.DRUGS
             Distribution                    POISSON
             Link Function                   LOG
             Dependent Variable              COUNT
             Observations Used               8


                      Class Level Information

                      Class     Levels  Values

                      A              2  1 2
                      C              2  1 2
                      M              2  1 2


               Criteria For Assessing Goodness Of Fit

        Criterion             DF         Value      Value/DF

        Deviance               2      187.7543       93.8772
        Scaled Deviance        2      187.7543       93.8772
        Pearson Chi-Square     2      177.6144       88.8072
        Scaled Pearson X2      2      177.6144       88.8072
        Log Likelihood         .    11916.9222             .


                   Analysis Of Parameter Estimates

  Parameter          DF    Estimate     Std Err   ChiSquare  Pr>Chi

  INTERCEPT           1      5.1921      0.0609   7273.5388  0.0001
  A          1        1      1.1272      0.0641    309.0133  0.0001
  A          2        0      0.0000      0.0000           .       .
  C          1        1     -0.2351      0.0555     17.9385  0.0001
  C          2        0      0.0000      0.0000           .       .
  M          1        1     -6.6209      0.4737    195.3511  0.0001
  M          2        0      0.0000      0.0000           .       .
  A*M        1  1     1      4.1251      0.4529     82.9437  0.0001
  A*M        1  2     0      0.0000      0.0000           .       .
  A*M        2  1     0      0.0000      0.0000           .       .
  A*M        2  2     0      0.0000      0.0000           .       .
  C*M        1  1     1      3.2243      0.1610    401.1651  0.0001
  C*M        1  2     0      0.0000      0.0000           .       .
  C*M        2  1     0      0.0000      0.0000           .       .
  C*M        2  2     0      0.0000      0.0000           .       .
  SCALE               0      1.0000      0.0000           .       .

NOTE:  The scale parameter was held fixed.


                                     Observation Statistics

    COUNT       Pred      Xbeta        Std    HessWgt      Lower      Upper     Resraw     Reschi

      911   909.2396     6.8126     0.0332   909.2396   852.0263   970.2947     1.7604     0.0584
      538   438.8404     6.0841     0.0444   438.8404   402.3018   478.6977    99.1596     4.7335
       44    45.7604     3.8234     0.1475    45.7604    34.2745    61.0954    -1.7604    -0.2602
      456   555.1596     6.3193     0.0401   555.1596   513.2140   600.5333   -99.1596    -4.2085
        3     4.7606     1.5604     0.4473     4.7606     1.9813    11.4386    -1.7606    -0.8069
       43   142.1596     4.9570     0.0638   142.1596   125.4565   161.0865   -99.1596    -8.3166
        The SAS System           15:16 Tuesday, May 18, 1999   2


                                     Observation Statistics

    COUNT       Pred      Xbeta        Std    HessWgt      Lower      Upper     Resraw     Reschi

        2     0.2396    -1.4288     0.4698     0.2396     0.0954     0.6017     1.7604     3.5965
      279   179.8405     5.1921     0.0609   179.8405   159.6126   202.6319    99.1595     7.3942