From 388e20673cbd73166fbf694acabe9e2319b00320 Mon Sep 17 00:00:00 2001
From: Jens Sauer <jsauer65@gmail.com>
Date: Sun, 17 Jan 2021 15:27:15 +0100
Subject: [PATCH] survival: Add SM survival plot

Plot survival after SM diagnosis.
---
 survival.R | 11 +++++++++++
 1 file changed, 11 insertions(+)

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)
 }