Estimation of the marginal survival time in the presence of dependent competing risks using inverse probability of censoring weighted (IPCW) methods

In medical studies, there is interest in inferring the marginal distribution of a survival time subject to competing risks. The Kyushu Lipid Intervention Study (KLIS) was a clinical study for hypercholesterolemia, where pravastatin treatment was compared with conventional treatment. The primary endpoint was time to events of coronary heart disease (CHD). In this study, however, some subjects died from causes other than CHD or were censored due to loss to follow-up. Because the treatments were targeted to reduce CHD events, the investigators were interested in the effect of the treatment on CHD events in the absence of causes of death or events other than CHD. In this paper, we present a method for estimating treatment group-specific marginal survival curves of time-to-event data in the presence of dependent competing risks. The proposed method is a straightforward extension of the Inverse Probability of Censoring Weighted (IPCW) method to settings with more than one reason for censoring. The results of our analysis showed that the IPCW marginal incidence for CHD was almost the same as the lower bound for which subjects with competing events were assumed to be censored at the end of all follow-up. This result provided reassurance that the results in KLIS were robust to competing risks.

[1]  A. V. Peterson,et al.  Bounds for a joint distribution function with fixed sub-distribution functions: Application to competing risks. , 1976, Proceedings of the National Academy of Sciences of the United States of America.

[2]  D. Rubin,et al.  The central role of the propensity score in observational studies for causal effects , 1983 .

[3]  J. Sasaki,et al.  Pravastatin use and risk of coronary events and cerebral infarction in Japanese men with hypercholesterolemia: The Kyushu Lipid Intervention Study (KLIS) , 2000 .

[4]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data: Kalbfleisch/The Statistical , 2002 .

[5]  James M. Robins,et al.  Estimation of the failure time distribution in the presence of informative censoring , 2002 .

[6]  J P Klein,et al.  Bounds on net survival probabilities for dependent competing risks. , 1988, Biometrics.

[7]  E. Kaplan,et al.  Nonparametric Estimation from Incomplete Observations , 1958 .

[8]  M Gail,et al.  A review and critique of some models used in competing risk analysis. , 1975, Biometrics.

[9]  Susanne Rosthøj,et al.  Competing risks as a multi-state model , 2002, Statistical methods in medical research.

[10]  Katsuhiko Yano,et al.  Elevated Serum Cholesterol Is a Risk Factor for Both Coronary Heart Disease and Thromboembolic Stroke in Hawaiian Japanese Men: Implications of Shared Risk , 1994, Stroke.

[11]  J. Robins,et al.  Marginal Structural Models and Causal Inference in Epidemiology , 2000, Epidemiology.

[12]  D. Horvitz,et al.  A Generalization of Sampling Without Replacement from a Finite Universe , 1952 .

[13]  R. Gray.  A Class of $K$-Sample Tests for Comparing the Cumulative Incidence of a Competing Risk , 1988 .

[14]  D Scharfstein,et al.  Inference in Randomized Studies with Informative Censoring and Discrete Time‐to‐Event Endpoints , 2001, Biometrics.

[15]  Robert Gray,et al.  A Proportional Hazards Model for the Subdistribution of a Competing Risk , 1999 .

[16]  John P. Klein,et al.  Estimates of marginal survival for dependent competing risks based on an assumed copula , 1995 .

[17]  L. J. Wei,et al.  The Robust Inference for the Cox Proportional Hazards Model , 1989 .

[18]  J. Robins,et al.  Correcting for Noncompliance and Dependent Censoring in an AIDS Clinical Trial with Inverse Probability of Censoring Weighted (IPCW) Log‐Rank Tests , 2000, Biometrics.

[19]  Y. Matsuyama,et al.  Marginal Structural Models as a Tool for Standardization , 2003, Epidemiology.

[20]  James M. Robins,et al.  Causal Inference from Complex Longitudinal Data , 1997 .

[21]  Eric V. Slud,et al.  Dependent competing risks and summary survival curves , 1983 .

[22]  J. Robins,et al.  Analysis of semiparametric regression models for repeated outcomes in the presence of missing data , 1995 .

[23]  J. Robins,et al.  Estimating the effect of zidovudine on Kaposi's sarcoma from observational data using a rank preserving structural failure-time model. , 1998, Statistics in medicine.

[24]  Somnath Datta,et al.  The Kaplan–Meier Estimator as an Inverse-Probability-of-Censoring Weighted Average , 2001, The American statistician.

[25]  N. Keiding,et al.  Multi-state models for event history analysis , 2002, Statistical methods in medical research.