Skip to content

Instantly share code, notes, and snippets.

@MarcinKosinski
Last active January 15, 2017 17:48
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save MarcinKosinski/47b0a38a4c812b66b324255e4a9877ba to your computer and use it in GitHub Desktop.
Save MarcinKosinski/47b0a38a4c812b66b324255e4a9877ba to your computer and use it in GitHub Desktop.
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