Analyses of cumulative incidence functions via non-parametric multiple imputation

We describe a non-parametric multiple imputation method that recovers the missing potential censoring information from competing risks failure times for the analysis of cumulative incidence functions. The method can be applied in the settings of stratified analyses, time-varying covariates, weighted analysis of case-cohort samples and clustered survival data analysis, where no current available methods can be readily implemented. The method uses a Kaplan-Meier imputation method for the censoring times to form an imputed data set, so cumulative incidence can be analyzed using techniques and software developed for ordinary right censored survival data. We discuss the methodology and show from both simulations and real data examples that the method yields valid estimates and performs well. The method can be easily implemented via available software with a minor programming requirement (for the imputation step). It provides a practical, alternative analysis tool for otherwise complicated analyses of cumulative incidence of competing risks data.