Multiple survival curves

This is an instruction of how to plot multiple survival curves

library(survival)
# Create a demo data set
#::::::::::::::::::::::::::::::::::::::::::::::::::::::::
set.seed(123)
demo.data <- data.frame(
  os.time = colon$time,
  os.status = colon$status,
  pfs.time = sample(colon$time),
  pfs.status = colon$status,
  sex = colon$sex, rx = colon$rx, adhere = colon$adhere
)
demo.data

# Ex1: Combine null models
#::::::::::::::::::::::::::::::::::::::::::::::::::::::
# Fit
library(survminer)
training <- survfit( Surv(os.time, os.status) ~ adhere, data = demo.data)
testing <- survfit( Surv(os.time, os.status) ~ sex, data = demo.data)

# Combine on the same plot
fit <- list(training,testing)
ggsurvplot_combine(fit, demo.data,pval = TRUE,palette = "jco",risk.table = FALSE,
                   legend.labs = c("Group 1","Group 2","Group 3","Group 4"))