* SAS code to fit two-group minnow mortality data;
 options linesize=70;

  data minnow2;
     input pop $ conc mort n @@;
     ic = pop="C";
     icconc = ic*conc;
     cards;
  C   31.0 60   60 C   14.5 60   60 C   11.8 46   60
  C   11.2 47   60 C    7.5 10   60 C    3.9  1   60
  C    3.3  0   60 C    1.4  2   60 C    0.8  1   60
  E   30.6 60   60 E   14.6 60   60 E   12.2 59   60
  E   10.9 51   60 E    7.4 24   60 E    3.9 15   60
  E    3.3 10   60 E    1.4 11   60 E    0.8  8   60
  ;

  proc sort; by pop;

  proc probit; by pop;         * See Fig. 7.25;
     model mort/n=conc / d=logistic inversecl;

  proc genmod;                 * See Fig. 7.26;
     model mort/n=conc ic icconc / dist=bin link=logit type1;

 run;