options nodate pagesize=54 linesize=80;

/* Read in the input file   */

data ex0801;
* infile '../../../data/biostat/ex0801.dat';
infile 'C:\biostat\data\ex0801.dat';
input dose death;

/* Run the logistic regression with dependent variable=death
  and independent variable = dose */
/* Proc Genmod is used for Generalized Linear Models  
  in this model, we are using a binomial distribution with a logit link function  
  the DESCENDING option means to fit the probability of being in the "high"
  death category (1) as opposed to the "low" (0).
  SAS' default is to fit probability in "low" category*/

proc genmod descending;
model death = dose / dist=binomial link=logit;
run;

quit;