/*-----------------------------------------------------
| Filename : COXMOD.SAS                               |
| The Rockall Trial - Interpretation of Linear Model  |
|                     Results for the Cox Model       |
-----------------------------------------------------*/

options nodate nonumber ps = 64 ls = 78;

libname perm 'c:\user';
libname library 'c:\user';

title1 'PRACTICAL 4';

/* Calculation of Percentile Survival Times */

data explan;
input sii siii age treat;
cards;
1 0 30 0
1 0 30 1
;
proc phreg data = perm.lurgi;
  model days * cens(0) = sii siii age treat;
  baseline covariates = explan
           out = savedata
           survival = surv / nomean;
  title2 'Calculating Percentiles';
run;

proc print data = savedata;
  where surv <0.5 AND surv >0.49;
  title2 'data = savedata (output data set from PROC PHREG)';
  title3 'The survivor Function estimates used to Calculate Percentiles';
run;

proc print data=savedata;
  where surv <0.1 AND surv >0.095;
run;

  title2;
  title3;