> leukemia time status group 1 9 1 Maintained 2 13 1 Maintained 3 13 0 Maintained 4 18 1 Maintained 5 23 1 Maintained 6 28 0 Maintained 7 31 1 Maintained 8 34 1 Maintained 9 45 0 Maintained 10 48 1 Maintained 11 161 0 Maintained 12 5 1 Nonmaintained 13 5 1 Nonmaintained 14 8 1 Nonmaintained 15 8 1 Nonmaintained 16 12 1 Nonmaintained 17 16 0 Nonmaintained 18 23 1 Nonmaintained 19 27 1 Nonmaintained 20 30 1 Nonmaintained 21 33 1 Nonmaintained 22 43 1 Nonmaintained 23 45 1 Nonmaintained > Surv(leukemia$time,leukemia$status) [1] 9 13 13+ 18 23 28+ 31 34 45+ 48 161+ 5 5 [14] 8 8 12 16+ 23 27 30 33 43 45 > summary(Surv(leukemia$time,leukemia$status)) Min. 1st Qu. Median Mean 3rd Qu. Max. 0 1 3 15.13 23 161 > survfit(Surv(leukemia$time,leukemia$status)) Call: survfit(formula = Surv(leukemia$time, leukemia$status)) n events mean se(mean) median 0.95LCL 0.95UCL 23 18 36.4 9.85 27 18 45 > survfit(Surv(time,status), data=leukemia) Call: survfit(formula = Surv(time, status), data = leukemia) n events mean se(mean) median 0.95LCL 0.95UCL 23 18 36.4 9.85 27 18 45 > summary(survfit(Surv(time,status), data=leukemia)) Call: survfit(formula = Surv(time, status), data = leukemia) time n.risk n.event survival std.err lower 95% CI upper 95% CI 5 23 2 0.9130 0.0588 0.8049 1.000 8 21 2 0.8261 0.0790 0.6848 0.996 9 19 1 0.7826 0.0860 0.6310 0.971 12 18 1 0.7391 0.0916 0.5798 0.942 13 17 1 0.6957 0.0959 0.5309 0.912 18 14 1 0.6460 0.1011 0.4753 0.878 23 13 2 0.5466 0.1073 0.3721 0.803 27 11 1 0.4969 0.1084 0.3240 0.762 30 9 1 0.4417 0.1095 0.2717 0.718 31 8 1 0.3865 0.1089 0.2225 0.671 33 7 1 0.3313 0.1064 0.1765 0.622 34 6 1 0.2761 0.1020 0.1338 0.569 43 5 1 0.2208 0.0954 0.0947 0.515 45 4 1 0.1656 0.0860 0.0598 0.458 48 2 1 0.0828 0.0727 0.0148 0.462 > survfit(Surv(time,status)~group, data=leukemia) Call: survfit(formula = Surv(time, status) ~ group, data = leukemia) n events mean se(mean) median 0.95LCL 0.95UCL group=Maintained 11 7 52.6 19.83 31 18 NA group=Nonmaintained 12 11 22.7 4.18 23 8 NA > summary(survfit(Surv(time,status)~group, data=leukemia)) Call: survfit(formula = Surv(time, status) ~ group, data = leukemia) group=Maintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 9 11 1 0.909 0.0867 0.7541 1.000 13 10 1 0.818 0.1163 0.6192 1.000 18 8 1 0.716 0.1397 0.4884 1.000 23 7 1 0.614 0.1526 0.3769 0.999 31 5 1 0.491 0.1642 0.2549 0.946 34 4 1 0.368 0.1627 0.1549 0.875 48 2 1 0.184 0.1535 0.0359 0.944 group=Nonmaintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 5 12 2 0.8333 0.1076 0.6470 1.000 8 10 2 0.6667 0.1361 0.4468 0.995 12 8 1 0.5833 0.1423 0.3616 0.941 23 6 1 0.4861 0.1481 0.2675 0.883 27 5 1 0.3889 0.1470 0.1854 0.816 30 4 1 0.2917 0.1387 0.1148 0.741 33 3 1 0.1944 0.1219 0.0569 0.664 43 2 1 0.0972 0.0919 0.0153 0.620 45 1 1 0.0000 NA NA NA > summary(survfit(Surv(time,status), data=leukemia,type="fleming-harrington")) Call: survfit(formula = Surv(time, status), data = leukemia, type = "fleming-harrington") time n.risk n.event survival std.err lower 95% CI upper 95% CI 5 23 2 0.917 0.0590 0.8081 1.000 8 21 2 0.833 0.0797 0.6909 1.000 9 19 1 0.791 0.0869 0.6375 0.981 12 18 1 0.748 0.0927 0.5867 0.954 13 17 1 0.705 0.0973 0.5382 0.924 18 14 1 0.657 0.1028 0.4831 0.892 23 13 2 0.563 0.1105 0.3833 0.827 27 11 1 0.514 0.1121 0.3352 0.788 30 9 1 0.460 0.1141 0.2830 0.748 31 8 1 0.406 0.1143 0.2337 0.705 33 7 1 0.352 0.1130 0.1875 0.660 34 6 1 0.298 0.1101 0.1444 0.615 43 5 1 0.244 0.1053 0.1046 0.569 45 4 1 0.190 0.0987 0.0686 0.526 48 2 1 0.115 0.1011 0.0206 0.643 > summary(survfit(Surv(time,status)~group, data=leukemia,type="fleming-harrington")) Call: survfit(formula = Surv(time, status) ~ group, data = leukemia, type = "fleming-harrington") group=Maintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 9 11 1 0.913 0.0871 0.7575 1.000 13 10 1 0.826 0.1174 0.6253 1.000 18 8 1 0.729 0.1422 0.4974 1.000 23 7 1 0.632 0.1572 0.3882 1.000 31 5 1 0.517 0.1731 0.2687 0.997 34 4 1 0.403 0.1781 0.1695 0.958 48 2 1 0.244 0.2038 0.0477 1.000 group=Nonmaintained time n.risk n.event survival std.err lower 95% CI upper 95% CI 5 12 2 0.8465 0.109 0.6572 1.000 8 10 2 0.6930 0.141 0.4645 1.000 12 8 1 0.6116 0.149 0.3791 0.987 23 6 1 0.5177 0.158 0.2849 0.941 27 5 1 0.4239 0.160 0.2021 0.889 30 4 1 0.3301 0.157 0.1300 0.838 33 3 1 0.2365 0.148 0.0692 0.808 43 2 1 0.1435 0.136 0.0225 0.914 45 1 1 0.0528 Inf 0.0000 1.000 > lm_survfit(Surv(time,status)~group, data=leukemia) > plot(lm[1,]) > ovarian futime fustat age residual.dz rx ecog.ps 1 59 1 72.3315 2 1 1 2 115 1 74.4932 2 1 1 3 156 1 66.4658 2 1 2 4 421 0 53.3644 2 2 1 5 431 1 50.3397 2 1 1 6 448 0 56.4301 1 1 2 7 464 1 56.9370 2 2 2 8 475 1 59.8548 2 2 2 9 477 0 64.1753 2 1 1 10 563 1 55.1781 1 2 2 11 638 1 56.7562 1 1 2 12 744 0 50.1096 1 2 1 13 769 0 59.6301 2 2 2 14 770 0 57.0521 2 2 1 15 803 0 39.2712 1 1 1 16 855 0 43.1233 1 1 2 17 1040 0 38.8932 2 1 2 18 1106 0 44.6000 1 1 1 19 1129 0 53.9068 1 2 1 20 1206 0 44.2055 2 2 1 21 1227 0 59.5890 1 2 2 22 268 1 74.5041 2 1 2 23 329 1 43.1370 2 1 1 24 353 1 63.2192 1 2 2 25 365 1 64.4247 2 2 1 26 377 0 58.3096 1 2 1