> # Print out the survival times, censored values have a "+" next to them (Surv must be spelled exactly this way) > > Surv(time, censor==0) [1] 7 7 7 7 7 7 7 7 7 7 7 7 7 13 13 13 13 13 [19] 13 13 13 13 13 13 13 13 13 13 19 19 19 19 19 19 19 19 [37] 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 25 25 25 [55] 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 [73] 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 [91] 25 25 25 25 25 31 31 31 31 31 31 31 31 31 31 31 31 31 [109] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [127] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [145] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [163] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [181] 31 31 31 31 31 37 37 37 37 37 37 37 37 37 37 37 37 37 [199] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [217] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [235] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [253] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [271] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 43 [289] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [307] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [325] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [343] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [361] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [379] 43 43 43 43 43 43 43 43 43 43 43 49 49 49 49 49 49 49 [397] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [415] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [433] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [451] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [469] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [487] 49 49 49 49 49 49 55 55 55 55 55 55 55 55 55 55 55 55 [505] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [523] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [541] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [559] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [577] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [595] 55 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [613] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [631] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [649] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [667] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 67 [685] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [703] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [721] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [739] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [757] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 73 73 73 [775] 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [793] 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [811] 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [829] 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [847] 73 73 73 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 [865] 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 [883] 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 [901] 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 [919] 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 85 [937] 85 85 85 85 85 85 85 85 85 85 7 7 7 7 7 7 7 7 [955] 7 7 7 7 7 7 7 13 13 13 13 13 13 13 13 13 13 13 [973] 13 13 13 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 [991] 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 19 [1009] 19 19 19 19 25 25 25 25 25 25 25 25 25 25 25 25 25 25 [1027] 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 [1045] 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 [1063] 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 25 [1081] 25 25 25 25 25 25 25 25 25 31 31 31 31 31 31 31 31 31 [1099] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [1117] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [1135] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [1153] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [1171] 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 31 [1189] 31 31 31 31 31 31 31 31 31 37 37 37 37 37 37 37 37 37 [1207] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [1225] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [1243] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [1261] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [1279] 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 37 [1297] 37 37 37 37 37 37 37 37 43 43 43 43 43 43 43 43 43 43 [1315] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [1333] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [1351] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [1369] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [1387] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [1405] 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 43 [1423] 43 43 43 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [1441] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [1459] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [1477] 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 49 [1495] 49 49 49 49 49 49 49 49 49 49 49 49 49 55 55 55 55 55 [1513] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [1531] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [1549] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [1567] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [1585] 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 55 [1603] 55 55 55 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [1621] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [1639] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [1657] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [1675] 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 61 [1693] 61 61 61 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [1711] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [1729] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [1747] 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 67 [1765] 67 67 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [1783] 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [1801] 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [1819] 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 73 [1837] 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 [1855] 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 79 [1873] 79 79 79 85 85 > > # Obtain survival tables and curves by regimen > > surv.byreg <- survfit(Surv(time, censor==0)~trt) > > summary(surv.byreg) Call: survfit(formula = Surv(time, censor == 0) ~ trt) trt=1 time n.risk n.event survival std.err lower 95% CI upper 95% CI 7 946 13 0.9863 0.00379 0.97887 0.9937 13 933 15 0.9704 0.00551 0.95966 0.9813 19 918 23 0.9461 0.00734 0.93181 0.9606 25 895 44 0.8996 0.00977 0.88063 0.9189 31 851 90 0.8044 0.01290 0.77956 0.8301 37 761 102 0.6966 0.01495 0.66793 0.7265 43 659 102 0.5888 0.01600 0.55826 0.6210 49 557 103 0.4799 0.01624 0.44911 0.5128 55 454 103 0.3710 0.01571 0.34149 0.4031 61 351 88 0.2780 0.01457 0.25088 0.3081 67 263 88 0.1850 0.01262 0.16183 0.2115 73 175 78 0.1025 0.00986 0.08492 0.1238 79 97 86 0.0116 0.00349 0.00646 0.0209 85 11 11 0.0000 NaN NA NA trt=2 time n.risk n.event survival std.err lower 95% CI upper 95% CI 7 931 15 0.98389 0.00413 0.975834 0.99201 13 916 14 0.96885 0.00569 0.957756 0.98007 19 902 37 0.92911 0.00841 0.912768 0.94574 25 865 77 0.84640 0.01182 0.823555 0.86988 31 788 108 0.73040 0.01454 0.702442 0.75947 37 680 107 0.61547 0.01594 0.584998 0.64752 43 573 121 0.48550 0.01638 0.454434 0.51869 49 452 82 0.39742 0.01604 0.367199 0.43013 55 370 98 0.29216 0.01490 0.264361 0.32288 61 272 90 0.19549 0.01300 0.171605 0.22270 67 182 71 0.11923 0.01062 0.100127 0.14197 73 111 70 0.04404 0.00672 0.032648 0.05940 79 41 39 0.00215 0.00152 0.000538 0.00858 85 2 2 0.00000 NaN NA NA > > plot(surv.byreg) > > # Conduct the Log-Rank test to compare survival curves > > survdiff(Surv(time,censor==0)~trt) Call: survdiff(formula = Surv(time, censor == 0) ~ trt) N Observed Expected (O-E)^2/E (O-E)^2/V trt=1 946 946 1050 10.2 31.5 trt=2 931 931 827 13.0 31.5 Chisq= 31.5 on 1 degrees of freedom, p= 2.04e-08 >