Adjusted Kaplan–Meier estimator and log‐rank test with inverse probability of treatment weighting for survival data

Estimation and group comparison of survival curves are two very common issues in survival analysis. In practice, the Kaplan-Meier estimates of survival functions may be biased due to unbalanced distribution of confounders. Here we develop an adjusted Kaplan-Meier estimator (AKME) to reduce confounding effects using inverse probability of treatment weighting (IPTW). Each observation is weighted by its inverse probability of being in a certain group. The AKME is shown to be a consistent estimate of the survival function, and the variance of the AKME is derived. A weighted log-rank test is proposed for comparing group differences of survival functions. Simulation studies are used to illustrate the performance of AKME and the weighted log-rank test. The method proposed here outperforms the Kaplan-Meier estimate, and it does better than or as well as other estimators based on stratification. The AKME and the weighted log-rank test are applied to two real examples: one is the study of times to reinfection of sexually transmitted diseases, and the other is the primary biliary cirrhosis (PBC) study.

[1]  James M. Robins,et al.  Commentary on ‘Using inverse weighting and predictive inference to estimate the effects of time‐varying treatments on the discrete‐time hazard’ , 2002 .

[2]  J. Coresh,et al.  Adjusting survival curves for confounders: a review and a new method. , 1996, American journal of epidemiology.

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

[4]  J F Lawless,et al.  Likelihood analysis of multi-state models for disease incidence and mortality. , 1988, Statistics in medicine.

[5]  J. Lunceford,et al.  Stratification and weighting via the propensity score in estimation of causal treatment effects: a comparative study , 2004, Statistics in medicine.

[6]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

[7]  O. Aalen.  A linear regression model for the analysis of life times. , 1989, Statistics in medicine.

[8]  J. Peto,et al.  Asymptotically Efficient Rank Invariant Test Procedures , 1972 .

[9]  P. Grambsch,et al.  Prognosis in primary biliary cirrhosis: Model for decision making , 1989, Hepatology.

[10]  B. Hankey,et al.  Evaluating differences in survival between two groups of patients. , 1971, Journal of chronic diseases.

[11]  J. Klein,et al.  Statistical Models Based On Counting Process , 1994 .

[12]  D G Altman,et al.  A note on the calculation of expected survival, illustrated by the survival of liver transplant patients. , 1991, Statistics in medicine.

[13]  P. Lavori,et al.  Using inverse weighting and predictive inference to estimate the effects of time‐varying treatments on the discrete‐time hazard , 2002, Statistics in medicine.

[14]  S Greenland,et al.  Analytic methods for two-stage case-control studies and other stratified designs. , 1991, Statistics in medicine.

[15]  D. Harrington,et al.  Counting Processes and Survival Analysis , 1991 .

[16]  L P Zhao,et al.  Designs and analysis of two-stage studies. , 1992, Statistics in medicine.

[17]  L. Zhao,et al.  Weighted Semiparametric Estimation in Regression Analysis with Missing Covariate Data , 1997 .

[18]  R B D'Agostino,et al.  Age-adjusted survival curves with application in the Framingham Study. , 1995, Statistics in medicine.

[19]  P. Sasieni,et al.  Adjusted Nelson–Aalen Estimates With Retrospective Matching , 2002 .

[20]  P. Sasieni,et al.  A weighted Kaplan–Meier estimator for matched data with application to the comparison of chemotherapy and bone‐marrow transplant in leukaemia , 2002, Statistics in medicine.

[21]  D. Amato,et al.  A generalized kaplan-meier estimator for heterogenous populations , 1988 .

[22]  J. Klein,et al.  Survival Analysis: Techniques for Censored and Truncated Data , 1997 .

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

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

[25]  J. Lunceford,et al.  Strati cation and weighting via the propensity score in estimation of causal treatment e ects : a comparative study , 2004 .