plotting survival curve using survfit2() and ggsurvfit(), clearly picking up numbers from somewhere else but I’m not sure.

surv_plot_sd <- survfit2(Surv(time_yrs_sd, status_sd) ~ drug_grp, data = final_matched) %>%
ggsurvfit() +
labs(
    x = "years",
    y = "overall KR/HR probability"
) +
add_confidence_interval() +
add_risktable() +
scale_x_continuous(breaks = seq(0, 25, by = 5), limits = c(0, 25)) +
coord_cartesian(ylim = c(0, 1))

The overlap only shows up for my first time for the risk table.