* SAS code to fit log-q GLiM to data in Table 9.8;
 options linesize=70;

 data ecoli;
 input naaz chrom y n  a0 a1  b0 b1  w00 w01 w10 w11;
 cards;
 0     0    10  50  1 0  1 0  1 0 0 0
 0     5.2  10  50  1 0  0 1  0 1 0 0
 0.04  0    19  50  0 1  1 0  0 0 1 0
 0.04  5.2  37  50  0 1  0 1  0 0 0 1

 proc genmod;
       if _xbeta_ < 1.e-6  then  prob = 1.e-6;
       else prob = 1 - exp(-_xbeta_) ;
 fwdlink link = -log(1 - _mean_) ;
 invlink  ilink = prob ;
 model y/n = a1 b1 w11/ d=binomial type1;

 run;