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

survival: Change scaling to days per year

Browse Source 
The base data was changed from month to years for more precision.
Therefor the scaling of the graphs must be adjusted to days per year
instead of month per year.
This commit is contained in:
Jens Sauer 2020-10-22 16:38:39 +02:00
parent e4cead5e99
commit 42aec82e36
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
survival.R
View File

@ -13,13 +13,13 @@ surv_asct <- Surv(time = secmal$event_time_asct, event = secmal$event_status)
# plot survival after diagnosis
# scaled to years
png(filename = "survival_dx.png", width = 3000, height = 3000, res = 300)
plot(survfit(surv_dx ~ 1), mark.time = TRUE, xscale = 12, xlab = "Years", ylab = "Survival")
plot(survfit(surv_dx ~ 1), mark.time = TRUE, xscale = 365.25, xlab = "Years", ylab = "Survival")
title("Kaplan-Meier estimate for\nsecMalASCT study", "Survival after diagnosis")
dev.off()
# One graph per sex
png(filename = "survival_dx_sex.png", width = 3000, height = 3000, res = 300)
plot(survfit(surv_dx ~ sex, data = secmal), mark.time = TRUE, xscale = 12, xlab = "Years", ylab = "Survival", lty = 2:3)
plot(survfit(surv_dx ~ sex, data = secmal), mark.time = TRUE, xscale = 365.25, xlab = "Years", ylab = "Survival", lty = 2:3)
title("Kaplan-Meier estimate for\nsecMalASCT study", "Survival after diagnosis")
legend(100, .9, c("Female", "Male"), lty = 2:3)
dev.off()
@ -27,13 +27,13 @@ dev.off()
# plot survival after diagnosis
# scaled to years
png(file = "survival_asct.png", width = 3000, height = 3000, res = 300)
plot(survfit(surv_asct ~ 1), mark.time = TRUE, xscale = 12, xlab = "Years", ylab = "Survival")
plot(survfit(surv_asct ~ 1), mark.time = TRUE, xscale = 365.25, xlab = "Years", ylab = "Survival")
title("Kaplan-Meier estimate for\nsecMalASCT study", "Survival after transplantation")
dev.off()
# One graph per sex
png(filename = "survival_asct_sex.png", width = 3000, height = 3000, res = 300)
plot(survfit(surv_asct ~ sex, data = secmal), mark.time = TRUE, xscale = 12, xlab = "Years", ylab = "Survival", lty = 2:3)
plot(survfit(surv_asct ~ sex, data = secmal), mark.time = TRUE, xscale = 365.25, xlab = "Years", ylab = "Survival", lty = 2:3)
title("Kaplan-Meier estimate for\nsecMalASCT study", "Survival after transplantation")
legend(100, .9, c("Female", "Male"), lty = 2:3)
dev.off()