Compute arbitrary quantile of a Cumulative Incidence estimate

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

1
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
kaspar.rufibach@gmail.com

See Also

cuminc and timepoint in the package cmprsk.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
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)
}

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.