## dat$trt
## group1 group2 group3
## 17 44 67
all time
Survival by STS (yes, no)
## Call: survfit(formula = Surv(time2death, death) ~ trt, data = dat)
##
## n events median 0.95LCL 0.95UCL
## trt=group1 17 6 2.256 1.864 NA
## trt=group2 44 27 0.723 0.356 NA
## trt=group3 67 49 0.742 0.490 2.19
## Call: survfit(formula = Surv(time2death, death) ~ trt, data = dat)
##
## trt=group1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 17 0 1.000 0.0000 1.0000 1
## 0.5 12 2 0.857 0.0935 0.6921 1
## 1.0 9 1 0.786 0.1097 0.5977 1
## 2.0 4 1 0.655 0.1505 0.4173 1
## 3.0 2 1 0.491 0.1812 0.2383 1
## 4.0 1 1 0.246 0.1958 0.0514 1
##
## trt=group2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 44 0 1.000 0.0000 1.0000 1.000
## 0.5 23 15 0.632 0.0759 0.4999 0.800
## 1.0 12 7 0.432 0.0821 0.2972 0.626
## 2.0 4 3 0.320 0.0826 0.1927 0.530
## 3.0 2 1 0.213 0.1029 0.0827 0.549
## 4.0 1 0 0.213 0.1029 0.0827 0.549
##
## trt=group3
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 67 0 1.000 0.0000 1.000 1.000
## 0.5 35 27 0.578 0.0621 0.469 0.714
## 1.0 29 6 0.479 0.0633 0.370 0.621
## 2.0 21 6 0.379 0.0620 0.275 0.522
## 3.0 13 3 0.322 0.0607 0.223 0.466
## 4.0 10 2 0.269 0.0613 0.172 0.420
## 5.0 9 0 0.269 0.0613 0.172 0.420
## Call:
## coxph(formula = Surv(time2death, death) ~ trt, data = dat)
##
## n= 128, number of events= 82
##
## coef exp(coef) se(coef) z Pr(>|z|)
## trtgroup2 0.9138 2.4937 0.4526 2.019 0.0435 *
## trtgroup3 0.7474 2.1116 0.4370 1.711 0.0872 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## trtgroup2 2.494 0.4010 1.0270 6.055
## trtgroup3 2.112 0.4736 0.8967 4.972
##
## Concordance= 0.548 (se = 0.033 )
## Rsquare= 0.038 (max possible= 0.994 )
## Likelihood ratio test= 4.97 on 2 df, p=0.08334
## Wald test = 4.08 on 2 df, p=0.1303
## Score (logrank) test = 4.29 on 2 df, p=0.1169
Death or Hospitalization by STS (yes, no)
## Call: survfit(formula = Surv(time2death, death.or.hosp) ~ trt, data = dat)
##
## n events median 0.95LCL 0.95UCL
## trt=group1 17 7 2.256 1.292 NA
## trt=group2 44 32 0.723 0.356 1.16
## trt=group3 67 59 0.742 0.400 1.88
## Call: survfit(formula = Surv(time2death, death.or.hosp) ~ trt, data = dat)
##
## trt=group1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 17 0 1.000 0.0000 1.000 1.000
## 0.5 12 2 0.857 0.0935 0.692 1.000
## 1.0 9 1 0.786 0.1097 0.598 1.000
## 2.0 4 2 0.573 0.1523 0.340 0.965
## 3.0 2 1 0.430 0.1686 0.199 0.927
## 4.0 1 1 0.215 0.1737 0.044 1.000
##
## trt=group2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 44 0 1.0000 0.0000 1.0000 1.000
## 0.5 23 16 0.6061 0.0772 0.4722 0.778
## 1.0 12 8 0.3840 0.0801 0.2552 0.578
## 2.0 4 5 0.2032 0.0736 0.0999 0.413
## 3.0 2 1 0.1354 0.0739 0.0465 0.395
## 4.0 1 1 0.0677 0.0605 0.0118 0.390
##
## trt=group3
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 67 0 1.000 0.0000 1.0000 1.000
## 0.5 35 28 0.565 0.0620 0.4559 0.701
## 1.0 29 6 0.468 0.0628 0.3602 0.609
## 2.0 21 7 0.355 0.0605 0.2538 0.495
## 3.0 13 8 0.220 0.0531 0.1367 0.353
## 4.0 10 3 0.169 0.0482 0.0965 0.295
## 5.0 9 1 0.152 0.0462 0.0837 0.276
## Call:
## coxph(formula = Surv(time2death, death.or.hosp) ~ trt, data = dat)
##
## n= 128, number of events= 98
##
## coef exp(coef) se(coef) z Pr(>|z|)
## trtgroup2 0.9756 2.6529 0.4187 2.33 0.0198 *
## trtgroup3 0.7037 2.0213 0.4044 1.74 0.0818 .
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## trtgroup2 2.653 0.3770 1.1677 6.027
## trtgroup3 2.021 0.4947 0.9149 4.465
##
## Concordance= 0.554 (se = 0.031 )
## Rsquare= 0.05 (max possible= 0.997 )
## Likelihood ratio test= 6.54 on 2 df, p=0.03795
## Wald test = 5.6 on 2 df, p=0.06095
## Score (logrank) test = 5.88 on 2 df, p=0.05294
Two year survival
Survival by STS (yes, no)
## Call: survfit(formula = Surv(fu2y, death2y) ~ trt, data = dat)
##
## n events median 0.95LCL 0.95UCL
## trt=group1 17 4 NA 1.864 NA
## trt=group2 44 25 0.723 0.356 NA
## trt=group3 67 39 0.742 0.490 NA
## Call: survfit(formula = Surv(fu2y, death2y) ~ trt, data = dat)
##
## trt=group1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 17 0 1.000 0.0000 1.000 1
## 0.5 12 2 0.857 0.0935 0.692 1
## 1.0 9 1 0.786 0.1097 0.598 1
## 2.0 4 1 0.655 0.1505 0.417 1
##
## trt=group2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 44 0 1.000 0.0000 1.000 1.000
## 0.5 23 15 0.632 0.0759 0.500 0.800
## 1.0 12 7 0.432 0.0821 0.297 0.626
## 2.0 4 3 0.320 0.0826 0.193 0.530
##
## trt=group3
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 67 0 1.000 0.0000 1.000 1.000
## 0.5 35 27 0.578 0.0621 0.469 0.714
## 1.0 29 6 0.479 0.0633 0.370 0.621
## 2.0 21 6 0.379 0.0620 0.275 0.522
## Call:
## coxph(formula = Surv(fu2y, death2y) ~ trt, data = dat)
##
## n= 128, number of events= 68
##
## coef exp(coef) se(coef) z Pr(>|z|)
## trtgroup2 1.2073 3.3445 0.5394 2.238 0.0252 *
## trtgroup3 1.1405 3.1282 0.5252 2.171 0.0299 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## trtgroup2 3.345 0.2990 1.162 9.627
## trtgroup3 3.128 0.3197 1.117 8.757
##
## Concordance= 0.551 (se = 0.034 )
## Rsquare= 0.056 (max possible= 0.99 )
## Likelihood ratio test= 7.34 on 2 df, p=0.02545
## Wald test = 5.18 on 2 df, p=0.07485
## Score (logrank) test = 5.79 on 2 df, p=0.0552
Death or Hospitalization by STS (yes, no)
## Call: survfit(formula = Surv(fu2y, death.or.hosp.2y) ~ trt, data = dat)
##
## n events median 0.95LCL 0.95UCL
## trt=group1 17 5 NA 1.292 NA
## trt=group2 44 29 0.723 0.356 1.16
## trt=group3 67 41 0.742 0.400 1.88
## Call: survfit(formula = Surv(fu2y, death.or.hosp.2y) ~ trt, data = dat)
##
## trt=group1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 17 0 1.000 0.0000 1.000 1.000
## 0.5 12 2 0.857 0.0935 0.692 1.000
## 1.0 9 1 0.786 0.1097 0.598 1.000
## 2.0 4 2 0.573 0.1523 0.340 0.965
##
## trt=group2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 44 0 1.000 0.0000 1.0000 1.000
## 0.5 23 16 0.606 0.0772 0.4722 0.778
## 1.0 12 8 0.384 0.0801 0.2552 0.578
## 2.0 4 5 0.203 0.0736 0.0999 0.413
##
## trt=group3
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 67 0 1.000 0.0000 1.000 1.000
## 0.5 35 28 0.565 0.0620 0.456 0.701
## 1.0 29 6 0.468 0.0628 0.360 0.609
## 2.0 21 7 0.355 0.0605 0.254 0.495
## Call:
## coxph(formula = Surv(fu2y, death.or.hosp.2y) ~ trt, data = dat)
##
## n= 128, number of events= 75
##
## coef exp(coef) se(coef) z Pr(>|z|)
## trtgroup2 1.1670 3.2122 0.4854 2.404 0.0162 *
## trtgroup3 0.9581 2.6068 0.4739 2.022 0.0432 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## trtgroup2 3.212 0.3113 1.241 8.317
## trtgroup3 2.607 0.3836 1.030 6.599
##
## Concordance= 0.556 (se = 0.033 )
## Rsquare= 0.057 (max possible= 0.994 )
## Likelihood ratio test= 7.54 on 2 df, p=0.02303
## Wald test = 5.8 on 2 df, p=0.05507
## Score (logrank) test = 6.31 on 2 df, p=0.04273
One year survival
Survival by STS (yes, no)
## Call: survfit(formula = Surv(fu1y, death1y) ~ trt, data = dat)
##
## n events median 0.95LCL 0.95UCL
## trt=group1 17 3 NA NA NA
## trt=group2 44 22 0.723 0.356 NA
## trt=group3 67 33 0.742 0.490 NA
## Call: survfit(formula = Surv(fu1y, death1y) ~ trt, data = dat)
##
## trt=group1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 17 0 1.000 0.0000 1.000 1
## 0.5 12 2 0.857 0.0935 0.692 1
## 1.0 9 1 0.786 0.1097 0.598 1
##
## trt=group2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 44 0 1.000 0.0000 1.000 1.000
## 0.5 23 15 0.632 0.0759 0.500 0.800
## 1.0 12 7 0.432 0.0821 0.297 0.626
##
## trt=group3
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 67 0 1.000 0.0000 1.000 1.000
## 0.5 35 27 0.578 0.0621 0.469 0.714
## 1.0 29 6 0.479 0.0633 0.370 0.621
## Call:
## coxph(formula = Surv(fu1y, death1y) ~ trt, data = dat)
##
## n= 128, number of events= 58
##
## coef exp(coef) se(coef) z Pr(>|z|)
## trtgroup2 1.2811 3.6008 0.6158 2.081 0.0375 *
## trtgroup3 1.2895 3.6308 0.6032 2.138 0.0325 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## trtgroup2 3.601 0.2777 1.077 12.04
## trtgroup3 3.631 0.2754 1.113 11.84
##
## Concordance= 0.575 (se = 0.036 )
## Rsquare= 0.053 (max possible= 0.983 )
## Likelihood ratio test= 7.01 on 2 df, p=0.03004
## Wald test = 4.7 on 2 df, p=0.09518
## Score (logrank) test = 5.39 on 2 df, p=0.0676
Death or Hospitalization by STS (yes, no)
## Call: survfit(formula = Surv(fu1y, death.or.hosp.1y) ~ trt, data = dat)
##
## n events median 0.95LCL 0.95UCL
## trt=group1 17 3 NA NA NA
## trt=group2 44 24 0.723 0.356 NA
## trt=group3 67 34 0.742 0.400 NA
## Call: survfit(formula = Surv(fu1y, death.or.hosp.1y) ~ trt, data = dat)
##
## trt=group1
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 17 0 1.000 0.0000 1.000 1
## 0.5 12 2 0.857 0.0935 0.692 1
## 1.0 9 1 0.786 0.1097 0.598 1
##
## trt=group2
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 44 0 1.000 0.0000 1.000 1.000
## 0.5 23 16 0.606 0.0772 0.472 0.778
## 1.0 12 8 0.384 0.0801 0.255 0.578
##
## trt=group3
## time n.risk n.event survival std.err lower 95% CI upper 95% CI
## 0.0 67 0 1.000 0.0000 1.000 1.000
## 0.5 35 28 0.565 0.0620 0.456 0.701
## 1.0 29 6 0.468 0.0628 0.360 0.609
## Call:
## coxph(formula = Surv(fu1y, death.or.hosp.1y) ~ trt, data = dat)
##
## n= 128, number of events= 61
##
## coef exp(coef) se(coef) z Pr(>|z|)
## trtgroup2 1.3742 3.9520 0.6127 2.243 0.0249 *
## trtgroup3 1.3185 3.7378 0.6024 2.189 0.0286 *
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## exp(coef) exp(-coef) lower .95 upper .95
## trtgroup2 3.952 0.2530 1.189 13.13
## trtgroup3 3.738 0.2675 1.148 12.17
##
## Concordance= 0.551 (se = 0.035 )
## Rsquare= 0.059 (max possible= 0.986 )
## Likelihood ratio test= 7.85 on 2 df, p=0.01975
## Wald test = 5.17 on 2 df, p=0.07525
## Score (logrank) test = 5.99 on 2 df, p=0.04998