options formdlim = '-';

data donner;
  infile '/home/pdixon/WWW/stat401b/data/donner.txt';
  input age male survive;
  if male = 1 then sex = 'M';
    else sex = 'F';

  ageclass = floor((age-14)/10)+1;
  
proc freq;
  table sex*survive /norow nocol nopercent chisq measures;
  table ageclass*sex /norow nocol nopercent;
  table ageclass*survive /norow nocol nopercent;
  
proc sort;
  by sex;
         
proc freq;
  by sex;
  table ageclass*survive /norow nocol nopercent;
  
proc logistic descending;
  where sex = 'M';
  model survive = age;
  output out = pred p = phat;

options noovp;
      
proc plot;
  plot phat*age = '+' survive*age /overlay;
              
proc logistic descending data = donner;
  class sex;
  
  model survive = age sex;
  title 'Estimates for no interaction model';