libname my "C:\Users\Bartek\Desktop\AHZ";

data colon;
    set my.colon;
run;
/*
proc contents data=colon;
run;
*/
data colon;
  set colon;
  rx=tranwrd(rx,"Obs","0");
  rx=tranwrd(rx,"Lev+5FU","1");
  rx=tranwrd(rx,"Lev","2");
  nodes=tranwrd(nodes,"NA","0");
  differ=tranwrd(differ,"NA","0");
run;

data colon;
  set colon;
  rx1=input(rx, best12.);
  nodes1 = input(nodes, best12.);
  differ1 = input(differ, best12.);
  drop rx;
run;
DATA colon;
    SET colon;
    IF arrest=0 THEN week=53;
RUN;

data dead remission ;
    set  colon;
if etype = '2' then output dead;
else if etype = '1' then output remission;

run;


proc freq data=remission;
    table status;
run;
proc freq data=dead;
    table status;
run;

proc means data=remission
    n nmiss mean std var stderr q1 median q3 range qrange skewness kurtosis;
var age time;
run;

proc univariate data=remission;
    var time;
    histogram time;
run;

ODS GRAPHICS ON;
PROC LIFETEST DATA=dead PLOTS=S(TEST);
    TIME time*status(0);
    STRATA rx1/ ADJUST=TUKEY;
RUN;
ODS GRAPHICS OFF;

PROC LIFETEST DATA=colon PLOTS=(S, LS, LLS);
    TIME time*status(0);
    test rx1 ;
RUN;