options ls=80;
data belief;
input race $ belief $ count;
cards;
W Yes 1339
W No   300
B Yes  260
B No    55
O Yes   88
O No    22
;
proc genmod order=data;
class race belief;
model count = race belief / dist=poi link=log ;
run;




/*********************************************************************************
                     Criteria For Assessing Goodness Of Fit

          Criterion                 DF           Value        Value/DF

          Deviance                   2          0.3565          0.1782
          Scaled Deviance            2          0.3565          0.1782
          Pearson Chi-Square         2          0.3601          0.1800
          Scaled Pearson X2          2          0.3601          0.1800
          Log Likelihood                    11415.5149


  Algorithm converged.


                        Analysis Of Parameter Estimates

                              Standard       Wald 95%          Chi-
Parameter       DF  Estimate     Error   Confidence Limits   Square  Pr > ChiSq

Intercept        1    3.0003    0.1061    2.7924    3.2083   799.54      <.0001
race       W     1    2.7014    0.0985    2.5083    2.8944   752.22      <.0001
race       B     1    1.0521    0.1107    0.8350    1.2692    90.24      <.0001
race       O     0    0.0000    0.0000    0.0000    0.0000      .         .
belief     Yes   1    1.4985    0.0570    1.3868    1.6101   691.89      <.0001
belief     No    0    0.0000    0.0000    0.0000    0.0000      .         .
Scale            0    1.0000    0.0000    1.0000    1.0000

NOTE: The scale parameter was held fixed.
*********************************************************************************/





data belief2;
input R $ B $ count;
cards;
W Y  1339
W N   300
B Y   260
B N    55
O Y    88
O N    22
;

proc genmod order=data;
class R B;
model count = R B R*B/ dist=poi link=log ;
run;


/*********************************************************************************
                     Criteria For Assessing Goodness Of Fit

          Criterion                 DF           Value        Value/DF

          Deviance                   0          0.0000           .
          Scaled Deviance            0          0.0000           .
          Pearson Chi-Square         0          0.0000           .
          Scaled Pearson X2          0          0.0000           .
          Log Likelihood                    11415.6932


  Algorithm converged.


                        Analysis Of Parameter Estimates

                               Standard       Wald 95%          Chi-
Parameter        DF  Estimate     Error   Confidence Limits   Square  Pr > ChiSq

Intercept         1    3.0910    0.2132    2.6732    3.5089   210.20      <.0001
R          W      1    2.6127    0.2209    2.1798    3.0457   139.92      <.0001
R          B      1    0.9163    0.2523    0.4219    1.4107    13.19      0.0003
R          O      0    0.0000    0.0000    0.0000    0.0000      .         .
B          Y      1    1.3863    0.2384    0.9191    1.8535    33.82      <.0001
B          N      0    0.0000    0.0000    0.0000    0.0000      .         .
R*B        W  Y   1    0.1096    0.2468   -0.3741    0.5933     0.20      0.6569
R*B        W  N   0    0.0000    0.0000    0.0000    0.0000      .         .
R*B        B  Y   1    0.1671    0.2808   -0.3833    0.7174     0.35      0.5519
R*B        B  N   0    0.0000    0.0000    0.0000    0.0000      .         .
R*B        O  Y   0    0.0000    0.0000    0.0000    0.0000      .         .
R*B        O  N   0    0.0000    0.0000    0.0000    0.0000      .         .
Scale             0    1.0000    0.0000    1.0000    1.0000
*********************************************************************************/




proc freq data=belief order=data; 
weight count;
tables race*belief / chisq expected;
run;

/*********************************************************************************
                               The FREQ Procedure

                            Table of race by belief

                      race      belief

                      Frequency|
                      Expected |
                      Percent  |
                      Row Pct  |
                      Col Pct  |Yes     |No      |  Total
                      ---------+--------+--------+
                      W        |   1339 |    300 |   1639
                               | 1339.6 | 299.37 |
                               |  64.87 |  14.53 |  79.41
                               |  81.70 |  18.30 |
                               |  79.37 |  79.58 |
                      ---------+--------+--------+
                      B        |    260 |     55 |    315
                               | 257.46 | 57.536 |
                               |  12.60 |   2.66 |  15.26
                               |  82.54 |  17.46 |
                               |  15.41 |  14.59 |
                      ---------+--------+--------+
                      O        |     88 |     22 |    110
                               | 89.908 | 20.092 |
                               |   4.26 |   1.07 |   5.33
                               |  80.00 |  20.00 |
                               |   5.22 |   5.84 |
                      ---------+--------+--------+
                      Total        1687      377     2064
                                  81.73    18.27   100.00

                     Statistics for Table of race by belief

             Statistic                     DF       Value      Prob
             ------------------------------------------------------
             Chi-Square                     2      0.3601    0.8352
             Likelihood Ratio Chi-Square    2      0.3565    0.8367
             Mantel-Haenszel Chi-Square     1      0.0178    0.8939
             Phi Coefficient                       0.0132
             Contingency Coefficient               0.0132
             Cramer's V                            0.0132

                               Sample Size = 2064
*********************************************************************************/