/* chex97.sas                                            */
 /*   Example 9.7 -- Cox's logistic regression problem    */
 /*   y = number of ingots not ready for rolling          */
 /*   x = heating times                                   */
 /*   m = number of ingots with x heating time            */
 /*                                                       */
  options ls=80 ;
  data a ;
    input x m y ;
    cards ;
   7   55   0
  14  157   2
  27  159   7
  51   16   3
 ;
 run ;
 /*                                                       */
 /*                       write out the data              */
  proc print data=a ;
    title 'Example 9.7 -- Cox''s logistic regression problem' ;
 run ;
 /*                                                       */
 /*                       logistic regression             */
 /*                         y is Binomial(m,p)            */
 /*                         logit(p) is b0 + b1*x         */
  proc logistic data=a ;
    model y / m = x ;
 run ;