jsauer/secmalasct
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

survival: Add plots per diagnosis

Browse Source 
This adds survival curves per diagnosis. They still look a bit
overloaded.
This commit is contained in:
Jens Sauer 2020-11-17 14:33:56 +01:00
parent 3c160c37bb
commit 95be046985
1 changed files with 20 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
survival.R
View File

@ -57,6 +57,24 @@ sma_plot_surv_asct_sex <- function() {
pval = TRUE)
}
# plot survival after diagnosis per diagnosis
sma_plot_surv_dx_dx <- function() {
ggsurvplot(survfit(surv_dx ~ diagnosis, data = secmal), data = secmal,
xscale = "d_y", title = "Survival after diagnosis",
break.time.by = sma_break.time.by,
risk.table = "abs_pct",
ggtheme = theme_bw())
}
# plot survival after diagnosis per diagnosis
sma_plot_surv_asct_dx <- function() {
ggsurvplot(survfit(surv_asct ~ diagnosis, data = secmal), data = secmal,
xscale = "d_y", title = "Survival after transplantation",
break.time.by = sma_break.time.by,
risk.table = "abs_pct",
ggtheme = theme_bw())
}
#
# Write survival plots to files
#
@ -65,4 +83,6 @@ sma_plot_file_surv <- function() {
sma_plot_file("survival_dx_sex.png", png, sma_plot_surv_dx_sex)
sma_plot_file("survival_asct.png", png, sma_plot_surv_asct)
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)
}