On the propensity score weighting analysis with survival outcome: Estimands, estimation, and inference

Propensity score analysis is widely used in observational studies to adjust for confounding and estimate the causal effect of a treatment on the outcome. When the outcome is survival time, there are special considerations on the definition of the causal estimand, point, and variance estimation that have not been thoroughly studied in the literature. We investigate propensity score analysis of survival data with a class of weighting methods. We consider the following estimands in the two-sample context: average survival time, restricted average survival time, survival probability, survival quantile, and the marginal hazard ratio. We propose a unified analytic framework to obtain the point and variance estimators. The proposed methodology properly adjusts for the sampling variability in the estimated propensity scores. Extensive simulations show that the point and variance estimators possess desired finite sample properties and demonstrate better numerical performance than some existing weighting and matching methods commonly used in the literature. The proposed methodology is illustrated with data from a breast cancer study.

[1]  Liang Li,et al.  A Weighting Analogue to Pair Matching in Propensity Score Analysis , 2013, The international journal of biostatistics.

[2]  E. Stuart,et al.  The performance of inverse probability of treatment weighting and full matching on the propensity score in the presence of model misspecification when estimating the effect of treatment on survival outcomes , 2015, Statistical methods in medical research.

[3]  P. Austin.  The use of propensity score methods with survival or time-to-event outcomes: reporting measures of effect similar to those used in randomized experiments , 2013, Statistics in medicine.

[4]  Elizabeth A. Stuart,et al.  Optimal full matching for survival outcomes: a method that merits more widespread use , 2015, Statistics in medicine.

[5]  Zhiliang Ying,et al.  Semiparametric transformation models for causal inference in time‐to‐event studies with all‐or‐nothing compliance , 2015, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[6]  Yoshiaki Uyama,et al.  Moving beyond the hazard ratio in quantifying the between-group difference in survival analysis. , 2014, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[7]  Gary King,et al.  MatchIt: Nonparametric Preprocessing for Parametric Causal Inference , 2011 .

[8]  Chaofeng Liu,et al.  Adjusted Kaplan–Meier estimator and log‐rank test with inverse probability of treatment weighting for survival data , 2005, Statistics in medicine.

[9]  P. Royston,et al.  The use of restricted mean survival time to estimate the treatment effect in randomized clinical trials when the proportional hazards assumption is in doubt , 2011, Statistics in medicine.

[10]  H. Uno,et al.  Alternatives to Hazard Ratios for Comparing the Efficacy or Safety of Therapies in Noninferiority Studies. , 2015, Annals of internal medicine.

[11]  J. Qin,et al.  Semiparametric estimation and inference for distributional and general treatment effects , 2009 .

[12]  Jing Cheng,et al.  Inference for the Effect of Treatment on Survival Probability in Randomized Trials with Noncompliance and Administrative Censoring , 2011, Biometrics.

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

[14]  Miguel A Hernán,et al.  The hazards of hazard ratios. , 2010, Epidemiology.

[15]  G. Imbens,et al.  Estimation of Causal Effects using Propensity Score Weighting: An Application to Data on Right Heart Catheterization , 2001, Health Services and Outcomes Research Methodology.

[16]  Patrick Royston,et al.  Restricted mean survival time: an alternative to the hazard ratio for the design and analysis of randomized trials with a time-to-event outcome , 2013, BMC Medical Research Methodology.

[17]  Stephen R Cole,et al.  Adjusted survival curves with inverse probability weights , 2004, Comput. Methods Programs Biomed..

[18]  Richard K. Crump,et al.  Moving the Goalposts: Addressing Limited Overlap in Estimation of Average Treatment Effects by Changing the Estimand , 2006, SSRN Electronic Journal.

[19]  Peter C Austin,et al.  Conditioning on the propensity score can result in biased estimation of common measures of treatment effect: a Monte Carlo study , 2007, Statistics in medicine.

[20]  D. Rubin.  Statistics and Causal Inference: Comment: Which Ifs Have Causal Answers , 1986 .

[21]  M. Sugihara.  Survival analysis using inverse probability of treatment weighted methods based on the generalized propensity score , 2010, Pharmaceutical statistics.

[22]  R. Porcher,et al.  Propensity score applied to survival data analysis through proportional hazards models: a Monte Carlo study , 2012, Pharmaceutical statistics.

[23]  D.,et al.  Regression Models and Life-Tables , 2022 .

[24]  G. Imbens,et al.  On the Failure of the Bootstrap for Matching Estimators , 2006 .

[25]  J. Pearl,et al.  Confounding and Collapsibility in Causal Inference , 1999 .

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

[27]  P. Austin,et al.  The performance of different propensity score methods for estimating absolute effects of treatments on survival outcomes: A simulation study , 2014, Statistical methods in medical research.

[28]  D. Rubin,et al.  Constructing a Control Group Using Multivariate Matched Sampling Methods That Incorporate the Propensity Score , 1985 .

[29]  Richard J. Cook,et al.  Does Cox analysis of a randomized survival study yield a causal treatment effect? , 2015, Lifetime data analysis.

[30]  D. Rubin,et al.  Reducing Bias in Observational Studies Using Subclassification on the Propensity Score , 1984 .

[31]  Erik T Parner,et al.  Causal inference in survival analysis using pseudo‐observations , 2017, Statistics in medicine.

[32]  G. Imbens,et al.  Matching on the Estimated Propensity Score , 2009 .

[33]  Elizabeth A Stuart,et al.  Matching methods for causal inference: A review and a look forward. , 2010, Statistical science : a review journal of the Institute of Mathematical Statistics.

[34]  Taisuke Otsu,et al.  Bootstrap Inference of Matching Estimators for Average Treatment Effects , 2017 .

[35]  G. Imbens,et al.  Efficient Estimation of Average Treatment Effects Using the Estimated Propensity Score , 2000 .

[36]  Peter C Austin,et al.  The performance of different propensity score methods for estimating marginal hazard ratios , 2007, Statistics in medicine.

[37]  Richard K. Crump,et al.  Dealing with limited overlap in estimation of average treatment effects , 2009 .

[38]  P. Ho,et al.  Extension of Kaplan-Meier methods in observational studies with time-varying treatment. , 2012, Value in health : the journal of the International Society for Pharmacoeconomics and Outcomes Research.

[39]  R. D'Agostino.  Adjustment Methods: Propensity Score Methods for Bias Reduction in the Comparison of a Treatment to a Non‐Randomized Control Group , 2005 .

[40]  Kari Lock Morgan,et al.  Balancing Covariates via Propensity Score Weighting , 2014, 1609.07494.

[41]  Peter C. Austin,et al.  Variance estimation when using inverse probability of treatment weighting (IPTW) with survival analysis , 2016, Statistics in medicine.

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

[43]  P. Holland.  Statistics and Causal Inference , 1985 .

[44]  Peter C Austin,et al.  A critical appraisal of propensity‐score matching in the medical literature between 1996 and 2003 , 2008, Statistics in medicine.

[45]  Rob J Hyndman,et al.  Mixed Model-Based Hazard Estimation , 2002 .