Extract median survival from survfit objects

[kassambara]

No description provided.

[kassambara]

New function surv_median() added to easily extract median survivals from one or a list of survfit objects.

# Load required packages
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::
library("survival")
library("survminer")

# compute different survfits
#:::::::::::::::::::::::::::::::::::::::::::::::::::::::
fit.null <- surv_fit(Surv(time, status) ~ 1, data = colon)

fit1 <- surv_fit(Surv(time, status) ~ sex, data = colon)

fit2 <- surv_fit(Surv(time, status) ~ adhere, data = colon)

fit.list <- list(sex = fit1, adhere = fit2)

surv_median(fit.null)

strata median lower upper
1    All   2351  2018  2910

surv_median(fit2)

    strata median lower upper
1 adhere=0   2718  2213    NA
2 adhere=1   1272   997  1885

surv_median(fit.list)

$sex
  strata median lower upper
1  sex=0   2174  1752    NA
2  sex=1   2527  1976  2910

$adhere
    strata median lower upper
1 adhere=0   2718  2213    NA
2 adhere=1   1272   997  1885

surv_median(fit.list, combine = TRUE)

      id   strata median lower upper
1    sex    sex=0   2174  1752    NA
2    sex    sex=1   2527  1976  2910
3 adhere adhere=0   2718  2213    NA
4 adhere adhere=1   1272   997  1885