/*  SAS code for using lifereg to
    estimate parameters in a Weibull
    model for survial times of breast
    cancer patients  */

data set1;
  input d y;
  datalines;
  1   5
  1   8
  1  10
  1  13
  1  18
  1  24
  1  26
  1  26
  1  31
  1  35
  1  40
  1  41
  1  48
  1  50
  1  59
  1  61
  1  68
  1  71
  0  76
  0 105
  0 107
  0 109
  1 113
  0 116
  1 118
  1 143
  0 154
  0 162
  0 188
  0 212
  0 217
  0 225
  run;

  proc lifereg covout outest=set2;
    model y*d(0)= / dist=weibull covb itprint
                    maxiter=30 converge=.0001;
    output out=set3 censored=c cdf=cdf predicted=p
           quantiles = .1 to .9 by .1 std=std;
   run;

proc print data=set2; run;

proc print data=set3; run;