cuminc_pairs: Pairwise 'cuminc' comparisons

View source: R/cuminc2.R

cuminc_pairsR Documentation

Pairwise cuminc comparisons

Description

Evaluate pairwise group differences in cumulative incidence with cuminc.

Usage

cuminc_pairs(
  object,
  data = NULL,
  rho = 0,
  cencode = NULL,
  method = p.adjust.methods,
  digits = 3L
)

Arguments

object

a cuminc2 formula or object of class cuminc2

data

a data frame to use (ignored if formula) is a cuminc2 object

rho, cencode

passed to cuminc

method

p-value correction method (default is 'holm'; see p.adjust

digits

number of digits past the decimal point to keep

Examples

## these two are equivalent since there are only two levels
form <- Surv(futime, event(censored)) ~ sex
cuminc_pairs(form, transplant)
cuminc2(form, transplant)$cuminc$Tests


## more useful in the case of three or more levels
form <- Surv(futime, event(censored)) ~ abo
cuminc_pairs(form, transplant)$p.value
cuminc2(form, transplant)$cuminc$Tests


tp <- transplant[transplant$abo %in% c('A', 'B', 'O'), ]
tp <- tp[tp$event %in% c('censored', 'ltx', 'death'), ]
ci <- cuminc2(form, tp)
ci[[1L]]$Tests
ciplot(ci)

cuminc_pairs(ci)$p.value$ltx


## unexported functions useful for extracting/formatting
cmprsk2:::pw_pval(ci, which = 2)
cmprsk2:::pw_text(ci, which = 'ltx')


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