diff --git a/survival.R b/survival.R index 9b2b3e1..c81a560 100644 --- a/survival.R +++ b/survival.R @@ -73,6 +73,16 @@ sma_plot_surv_asct_dx <- function() { ggtheme = theme_bw()) } +# plot survival after SM diagnosis +sma_plot_surv_sm <- function() { + ggsurvplot(survfit(Surv(event_time_sm, event_status) ~ 1, secmal), + data = secmal, xscale = "d_y", + break.time.by = 365.25, + legend = "none", + xlab = "Years", + ggtheme = theme_bw()) +} + # # Write survival plots to files # @@ -83,4 +93,5 @@ sma_plot_file_surv <- function() { sma_plot_file("survival_asct_sex.png", png, sma_plot_surv_asct_sex) sma_plot_file("survival_dx_dx.png", png, sma_plot_surv_dx_dx) sma_plot_file("survival_asct_dx.png", png, sma_plot_surv_asct_dx) + sma_plot_file("survival_sm.png", png, sma_plot_surv_sm) }