Last active
January 15, 2017 17:48
-
-
Save MarcinKosinski/47b0a38a4c812b66b324255e4a9877ba to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
devtools::install_github("kassambara/survminer", build_vignettes = TRUE) | |
library("survminer") | |
library("survival") | |
fit<- survfit(Surv(time, status) ~ sex, data = lung) | |
ggsurvplot( | |
fit, # survfit object with calculated statistics. | |
risk.table = TRUE, # show risk table. | |
pval = TRUE, # show p-value of log-rank test. | |
conf.int = TRUE, # show confidence intervals for | |
# point estimaes of survival curves. | |
xlim = c(0,500), # present narrower X axis, but not affect | |
# survival estimates. | |
xlab = "Time in days", # customize X axis label. | |
break.time.by = 100, # break X axis in time intervals by 500. | |
ggtheme = theme_light(), # customize plot and risk table with a theme. | |
risk.table.y.text.col = T,# colour risk table text annotations. | |
risk.table.y.text = FALSE,# show bars instead of names in text annotations | |
# in legend of risk table. | |
conf.int.style = "step", # customize style of confidence intervals | |
surv.median.line = "hv", # add the median survival pointer. | |
legend.labs = | |
c("Male", "Female"), # change legend labels. | |
palette = | |
c("#E7B800", "#2E9FDF"),# custom color palettes. | |
main = "Survival curves", # specify the title of the plot | |
submain = "Based on Kaplan-Meier estimates", # the subtitle of the plot | |
caption = "created with survminer", # the caption of the plot | |
font.main = c(16, "bold", "darkblue"), # font for titles of the plot, the table and censor part | |
font.submain = c(15, "bold.italic", "purple"), # font for subtitles in the plot, the table and censor part | |
font.caption = c(14, "plain", "orange"), # font for captions in the plot, the table and censor part | |
font.x = c(14, "bold.italic", "red"), # font for x axises in the plot, the table and censor part | |
font.y = c(14, "bold.italic", "darkred"), # font for y axises in the plot, the table and censor part | |
font.tickslab = c(12, "plain", "darkgreen"), # font for ticklabs in the plot, the table and censor part | |
########## risk table #########, | |
risk.table.title = "Note the risk set sizes", # the title of the risk table | |
risk.table.subtitle = "and remember about censoring.", # the subtitle of the risk table | |
risk.table.caption = "source code: website.com", # the caption of the risk table | |
risk.table.height = 0.35, # the height of the risk table | |
########## ncensor plot ###### | |
ncensor.plot = TRUE, # plot the number of censored subjects at time t | |
ncensor.plot.title = "Number of censorings", # as above but for the censoring plot | |
ncensor.plot.subtitle = "over the time.", | |
ncensor.plot.caption = "data available at data.com", | |
ncensor.plot.height = 0.35 | |
) |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment