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

survival: Improve plots

Browse Source 
Display a risk table, median survival line and adjust theme.
This commit is contained in:
Jens Sauer 2020-11-17 13:28:06 +01:00
parent 12ffcd4e38
commit c0c890c78c
1 changed files with 14 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
survival.R
View File

@ -17,7 +17,10 @@ surv_asct <- Surv(time = secmal$event_time_asct, event = secmal$event_status)
sma_plot_surv_dx <- function() {
ggsurvplot(survfit(surv_dx ~ 1), data = secmal, xscale = "d_y",
title = "Survival after diagnosis",
break.time.by = 5 * 365.25)
break.time.by = 5 * 365.25,
surv.median.line = "hv",
risk.table = "abs_pct",
ggtheme = theme_bw())
}
# plot survival after diagnosis per sex
@ -26,6 +29,9 @@ sma_plot_surv_dx_sex <- function() {
xscale = "d_y", title = "Survival after diagnosis",
legend.labs = c("Female", "Male"),
break.time.by = 5 * 365.25,
surv.median.line = "hv",
risk.table = "abs_pct",
ggtheme = theme_bw(),
pval = TRUE)
}
@ -33,7 +39,10 @@ sma_plot_surv_dx_sex <- function() {
sma_plot_surv_asct <- function() {
ggsurvplot(survfit(surv_asct ~ 1), data = secmal, xscale = "d_y",
title = "Survival after transplantation",
break.time.by = 5 * 365.25)
break.time.by = 5 * 365.25,
surv.median.line = "hv",
risk.table = "abs_pct",
ggtheme = theme_bw())
}
# plot survival after transplantation per sex
@ -42,6 +51,9 @@ sma_plot_surv_asct_sex <- function() {
xscale = "d_y", title = "Survival after transplantation",
legend.labs = c("Female", "Male"),
break.time.by = 5 * 365.25,
surv.median.line = "hv",
risk.table = "abs_pct",
ggtheme = theme_bw(),
pval = TRUE)
}