* SAS code to fit tBHQ data from Example 6.5;
 options linesize=70;

 data tBHQ;
 input dose y @@;
      n = 2000;
 cards;
 0     8    120   9    240   9
 480  18    600  39    720  44

 proc logistic;
 model y/n = dose / link=logit  covb ;
 run;