* LIFETST1.SAS
* xref:  Cox & Oakes table 1.1
* input:   6-MP   DATA A
* output:
* does: - demonstrate proc lifetest - KM estimates and plots
        - one sample
***************************************************;
filename INF '6-MP.DAT';
options linesize=176;
***************************************************;
data A;
  infile INF;
  input time delta group;
  label
    time  = 'Follow-up time (wks)         ..'
    delta = 'Event: 0=censored 1=observed ..'
    group = 'Group: 0=control  1=drug     ..'
  ;
  if group;
  title2 'The treatment (6-MP) group';
  run;
***************************************************;
proc lifetest data = A  method=KM  outs=B alpha=0.05 plots=(s, ls, lls);
  time time * delta(0);
  run;
***************************************************;
data C;
  set B;
  if not _censor_;
  * compute intervals on the log-scale;
  sdf_lcl = survival * exp(sdf_lcl/survival-1);
  sdf_ucl = survival * exp(sdf_ucl/survival-1);
  if (sdf_ucl > 1) then  sdf_ucl = 1;
  drop _censor_;
  run;
***************************************************;
proc print data=C; run;
***************************************************;
proc plot data=C;
  plot (sdf_lcl survival sdf_ucl)  * time / overlay;
  run;