DATA TABLE12_2; INPUT nodes time censored; /***************************/ /* Create Design Variables */ /***************************/ X1=(nodes=0); X2=(nodes=2); X3=(nodes=4); /*************************************/ /* Create Labels for Survival Curves */ /*************************************/ IF nodes=0 then LABEL="0 Nodes "; IF nodes=2 then LABEL="1-3 Nodes"; IF nodes=4 then LABEL="4 Nodes "; DATALINES; 0 17 0 0 18 0 0 28 0 0 42 0 0 51 0 0 100 1 0 101 0 0 106 0 0 116 1 0 117 0 0 121 0 0 125 1 0 126 1 0 140 1 0 153 0 0 160 1 0 164 1 0 185 1 0 186 1 0 194 1 0 195 1 0 197 1 0 205 1 0 212 1 0 219 1 0 221 0 0 221 1 0 252 1 0 252 1 0 252 1 0 252 1 0 252 1 0 252 1 0 252 1 2 1 0 2 4 0 2 6 0 2 6 0 2 9 0 2 12 0 2 13 0 2 16 0 2 50 0 2 53 0 2 58 0 2 63 0 2 65 0 2 80 0 2 88 0 2 97 0 2 123 1 2 128 0 2 147 1 2 149 1 2 154 1 2 216 1 2 217 1 2 223 1 2 229 1 2 231 1 2 239 1 2 242 1 4 1 0 4 3 0 4 3 0 4 5 0 4 9 0 4 10 0 4 11 0 4 11 0 4 23 0 4 23 0 4 36 0 4 38 0 4 49 0 4 65 0 4 121 0 4 128 0 4 139 1 4 152 0 4 152 1 4 157 0 4 229 1 4 252 1 4 252 1 RUN; /***************************/ /* Request File */ /***************************/ DATA requests; LENGTH label $ 9; /* Reserve 9 characters for the label */ INPUT nodes; X1=(nodes=0); X2=(nodes=2); X3=(nodes=4); IF nodes=0 then LABEL="0 Nodes "; IF nodes=2 then LABEL="1-3 Nodes"; IF nodes=4 then LABEL="4 Nodes "; DATALINES; 0 2 4 RUN; /***************************/ /* Compute hazard ratios */ /* and survival curves */ /***************************/ PROC PHREG DATA=table12_2; MODEL time*censored(1) = X1 X2 X3 / RISKLIMITS; ID label; BASELINE COVARIATES=requests OUT=survivalcurves SURVIVAL=S LOWER=LCL UPPER=UCL / NOMEAN; RUN; /******************************/ /* Graph the survival curves */ /******************************/ SYMBOL1 VALUE=none COLOR=RED I=STEPJ; SYMBOL2 VALUE=none COLOR=BLACK I=STEPJ; SYMBOL3 VALUE=none COLOR=BLUE I=STEPJ; PROC GPLOT DATA=survivalcurves; PLOT S*time=label; RUN;