survexp_plot: Observed Kaplan-Meier, expected and relative survival curves

Description

Displays the observed Kaplan-Meier, expected and relative survival curves

Usage

survexp_plot(
  futime,
  status,
  age,
  sex,
  entry_date,
  ratetable = survexp.fr::survexp.fr,
  main = "Observed and expected survival",
  xlab = "Time (years)",
  ylab = "Survival",
  col.km = "black",
  lwd.km = 2,
  lty.km = 1,
  conf.int.km = TRUE,
  col.exp = "blue",
  lwd.exp = 2,
  lty.exp = 1,
  main.rel = "Relative survival",
  ylab.rel = "Relative survival",
  col.rel = "black",
  lwd.rel = 2,
  lty.rel = 1,
  times = seq(0, max(futime, na.rm = TRUE)/365.241, length = 6)[-1],
  alpha = 0.05,
  xscale = 365.241,
  ...
)

Value

A matrix containing the values of relative survivals and their confidence intervals for each time of times

Arguments

futime: follow-up time of the subjects in days
status: 0 if censored or 1 if dead at futime
age: age in days
sex: "male" or "female"
entry_date: entry date in the study
ratetable: a table of event rates, such as survexp.fr or survexp.us
main: main title of the Kaplan-Meier and expected survivals plot
xlab: x-label of the plot
ylab: y-label of the plot
col.km: color of the observed survival curve
lwd.km: line width of the observed survival curve
lty.km: line type of the observed survival curve
conf.int.km: TRUE to display the confidence interval of the observed survival
col.exp: color of the expected survival curve
lwd.exp: line width of the expected survival curve
lty.exp: line type of the expected survival curve
main.rel: main title of the relative survival plot
ylab.rel: y-label of the relative survival plot
col.rel: color of the relative survival curve
lwd.rel: line width of the relative survival curve
lty.rel: line type of the relative survival curve
times: times to draw the confidence intervals of the relative survival
alpha: determines the confidence level (1-alpha) of the confidence intervals for the relative survival
xscale: see the xscale argument in plot.survfit
...: other arguments to be passed in plot.survfit

Author

Hugo Varet

Details

This function displays the observed and expected survivals, and the relative survival which is defined as: $$r(t) = exp(-exp(\beta) \times t)$$ where $exp(\beta)$ is the excess risk by time unit estimated by an additive Poisson model.

References

M. Pohar and J. Stare, Making relative survival analysis relatively easy, Computers in Biology and Medicine, 2007

M. Pohar and J. Stare, Relative survival analysis in R, Computers Methods and Programs in Biomedicine, 2006

Examples

Run this code

attach(data.example)
survexp_plot(futime, status, age, sex, entry_date)

Run the code above in your browser using DataLab