cuminc2: Cumulative incidence analysis

In raredd/cmprsk2: cmprsk2

Cumulative incidence analysis

Description

Estimate cumulative incidence functions from competing risks data and test equality across groups.

Usage

cuminc2(
  formula,
  data,
  rho = 0,
  cencode = NULL,
  subset = NULL,
  na.action = getOption("na.action")
)

Arguments

formula	a survival object formula, Surv(time, status(censor)) ~ response, where censor is a unique value of the status variable indicating the censoring code; note all other unique values of status will be treated as competing risks
data	a data frame in which to interpret the variables named in formula
rho, cencode, subset, na.action	passed to cuminc; the censoring indicator will be guessed from formula but may be overridden by cencode

See Also

summary.cuminc2; plot.cuminc2; cuminc;

Examples

tp <- within(transplant, {
  event_ind <- as.integer(event) - 1L
})

## no group, no strata
form <- Surv(futime, event_ind(0)) ~ 1
identical(
  cuminc2(form, tp)$cuminc,
  with(tp, cuminc(futime, event_ind))
)

## group-only
form <- Surv(futime, event_ind(0)) ~ sex
identical(
  cuminc2(form, tp)$cuminc,
  with(tp, cuminc(futime, event_ind, sex))
)

## strata-only
form <- Surv(futime, event_ind(0) == 1) ~ strata(abo)
identical(
  cuminc2(form, tp)$cuminc,
  with(tp, cuminc(futime, event_ind, strata = sex))
)

## group and strata
form <- Surv(futime, event_ind(0) == 1) ~ sex + strata(abo)
identical(
  cuminc2(form, tp)$cuminc,
  with(tp, cuminc(futime, event_ind, sex, abo))
)

raredd/cmprsk2 documentation built on March 29, 2024, 5:34 a.m.