cuminc {cmprsk}R Documentation

Cumulative Incidence Analysis

Description

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

Usage

cuminc(ftime, fstatus, group, strata, rho=0, cencode=0,
subset, na.action=na.omit)

Arguments

ftime failure time variable
fstatus variable with distinct codes for different causes of failure and also a distinct code for censored observations
group estimates will calculated within groups given by distinct values of this variable. Tests will compare these groups. If missing then treated as all one group (no test statistics)
strata stratification variable. Has no effect on estimates. Tests will be stratified on this variable. (all data in 1 stratum, if missing)
rho Power of the weight function used in the tests.
cencode value of fstatus variable which indicates the failure time is censored.
subset a logical vector specifying a subset of cases to include in the analysis
na.action a function specifying the action to take for any cases missing any of ftime, fstatus, group, strata, or subset.

Value

A list with components giving the subdistribution estimates for each cause in each group, and a component giving the test results (if the number of groups is >1). The components giving the estimates have names that are a combination of the group name and the cause code. The estimate for each group and cause combination is also a list.

$<groupcause>
$<groupcause>$time
the times where the estimates are calculated
$<groupcause>$est
the estimated sub-distribution functions. These are step functions (all corners of the steps given), so they can be plotted using ordinary lines() commands. Estimates at particular times can be located using the timepoints() function.
$<groupcause>$var
the estimated variance of the estimates, which are estimates of the asymptotic variance of Aalen (1978).
$Tests If the number of groups is >1, there is also a component at the top level giving the test statistics and p-values for comparing the subdistribution for each cause across groups. The test statistics are described in Gray (1988).

Author(s)

Robert Gray

References

Gray RJ (1988) A class of K-sample tests for comparing the cumulative incidence of a competing risk, ANNALS OF STATISTICS, 16:1141-1154.

Kalbfleisch and Prentice (1980) THE ANALYSIS OF FAILURE TIME DATA, p 168-9.

Aalen, O. (1978) Nonparametric estimation of partial transition probabilities in multiple decrement models, ANNALS OF STATISTICS, 6:534-545.

See Also

plot.cuminc timepoints print.cuminc

Examples

set.seed(2)
ss <- rexp(100)
gg <- factor(sample(1:3,100,replace=TRUE),1:3,c('a','b','c'))
cc <- sample(0:2,100,replace=TRUE)
strt <- sample(1:2,100,replace=TRUE)
print(xx <- cuminc(ss,cc,gg,strt))
plot(xx,lty=1,color=1:6)
# see also test.R, test.out

[Package cmprsk version 2.2-0 Index]