quantileCumInc: Compute arbitrary quantile of a Cumulative Incidence estimate

In biostatUZH: Misc Tools of the Department of Biostatistics, EBPI, University of Zurich

Description

In survival analysis, often some quantiles (mainly the median) of an estimated cumulative incidence function are of interest. This function computes quantiles of cumulative incidence functions received via cuminc in the package cmprsk.

Usage

quantileCumInc(time, event, group, quant = 0.5)

Arguments

time	Event times, censored or observed.
event	Event indicator.
group	Quantiles can be computed for several survival curves, defined by group.
quant	Quantile to be computed. Real number in [0, 1].

Value

Dataframe with quantiles estimates for each event-group combination.

Author(s)

Kaspar Rufibach
[e-mail]

See Also

cuminc and timepoint in the package cmprsk.

Examples

if (require("cmprsk")) {
  ## illustrate on simulated data
  set.seed(1977)
  n <- 180
  time <- rgamma(n, 2, 1)
  status <- factor(sample(rep(0:2, n / 3)), levels = 0:2, labels = 
      c("censored", "Time to 2nd tumor", "Death"))
  x <- factor(round(runif(n, 1, 2)), levels = 1:2, labels = 
      c("Tmt A", "Tmt B"))
  
  ## plot
  plot(cuminc(time, status, x))

  ## median time to event
  quantileCumInc(time, status, group = x, quant = 0.25)
}

biostatUZH documentation built on May 2, 2019, 6:06 p.m.