/* SAS code to enter the recurrence of infection
   data for patients using a dialysis machine */

data set1;
 infile "c:\stat565\data\kidney_infection.txt";
 input ID age sex gn an pkd time status recur;
 age10=age-10;
run;

proc print data=set1; run;


%include 'c:\documents and settings\kkoehler.IASTATE\my documents\courses\st565\sas\therneau\phlev.sas';

%phlev(data=set1, time=time, event=status,
       xvars = age10 sex gn an pkd, id=id, outlev=phlev,
       outvar= phvar );


proc phreg data=set1 cov;
 model time*status(0)= age10 sex gn an pkd /ties=efron;
run;


proc phreg data=set1 covm covs;
 model time*status(0)= age10 sex gn an pkd /ties=efron;
run;