options nodate pagesize=54 linesize=80;

/* Read in the input file   */

data ex0802;
* infile '../../../data/biostat/ex0802.dat';
infile 'C:\biostat\data\ex0802.dat';
input bldprs srmchl log10bp log10sc numdead numlive;
totcase=numdead+numlive;
run;

/* Run the logistic regression with dependent variable=death
  and independent variables = log10sc and log10bp  */
/* Proc Genmod is used for Generalized Linear Models  
  in this model, we are using a binomial distribution with a logit link function.
  Note the nature of the data, for each combination of SC and BP, we have
  numdead, numlive and totcase. we model the proportion numdead/totcase as dependent var */

proc genmod;
model numdead/totcase = log10sc log10bp / dist=binomial link=logit;
run;

quit;