Robust Estimation of Mean Functions and Treatment Effects for Recurrent Events Under Event-Dependent Censoring and Termination: Application to Skeletal Complications in Cancer Metastatic to Bone

In clinical trials featuring recurrent clinical events, the definition and estimation of treatment effects involves a number of interesting issues, especially when loss to follow-up may be event-related and when terminal events such as death preclude the occurrence of further events. This paper discusses a clinical trial of breast cancer patients with bone metastases where the recurrent events are skeletal complications, and where patients may die during the trial. We argue that treatment effects should be based on marginal rate and mean functions. When recurrent event data are subject to event-dependent censoring, however, ordinary marginal methods may yield inconsistent estimates. Incorporating correctly specified inverse probability of censoring weights into analyses can protect against dependent censoring and yield consistent estimates of marginal features. An alternative approach is to obtain estimates of rate and mean functions from models that involve some conditioning to render censoring conditionally independent. We consider three methods of estimating mean functions of recurrent event processes and examine the bias and efficiency of unweighted and inverse probability weighted versions of the methods with and without a terminating event. We compare the methods via simulation and use them to analyse the data from the breast cancer trial.

[1]  J. Robins,et al.  Recovery of Information and Adjustment for Dependent Censoring Using Surrogate Markers , 1992 .

[2]  P. Thall,et al.  Some covariance models for longitudinal count data with overdispersion. , 1990, Biometrics.

[3]  Jianwen Cai,et al.  Regression methods for gap time hazard functions of sequentially ordered multivariate failure time data , 2004 .

[4]  Niels Keiding,et al.  Statistical Models Based on Counting Processes , 1993 .

[5]  Gerda Claeskens,et al.  Nonparametric Estimation , 2011, International Encyclopedia of Statistical Science.

[6]  G. Satten,et al.  Estimation of Integrated Transition Hazards and Stage Occupation Probabilities for Non‐Markov Systems Under Dependent Censoring , 2002, Biometrics.

[7]  James M. Robins,et al.  Estimating the marginal survival function in the presence of time dependent covariates , 2001 .

[8]  M. Pepe,et al.  A qualifier Q for the survival function to describe the prevalence of a transient condition. , 1991, Statistics in medicine.

[9]  S. Lewis,et al.  Regression analysis , 2007, Practical Neurology.

[10]  R. Wolfe,et al.  Shared Frailty Models for Recurrent Events and a Terminal Event , 2004, Biometrics.

[11]  S. Keleş,et al.  Recurrent events analysis in the presence of time‐dependent covariates and dependent censoring , 2004 .

[12]  Laurence L. George,et al.  The Statistical Analysis of Failure Time Data , 2003, Technometrics.

[13]  Zhiliang Ying,et al.  Semiparametric regression for the mean and rate functions of recurrent events , 2000 .

[14]  James M Robins,et al.  Locally Efficient Estimation of a Multivariate Survival Function in Longitudinal Studies , 2002 .

[15]  Somnath Datta,et al.  Validity of the Aalen–Johansen estimators of stage occupation probabilities and Nelson–Aalen estimators of integrated transition hazards for non-Markov models , 2001 .

[16]  G. Hortobagyi,et al.  Long-term prevention of skeletal complications of metastatic breast cancer with pamidronate. Protocol 19 Aredia Breast Cancer Study Group. , 1998, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[17]  Virginie Rondeau,et al.  Joint frailty models for recurring events and death using maximum penalized likelihood estimation: application on cancer events. , 2006 .

[18]  P. Hougaard Shared frailty models , 2000 .

[19]  G. Hortobagyi,et al.  Efficacy of pamidronate in reducing skeletal complications in patients with breast cancer and lytic bone metastases. Protocol 19 Aredia Breast Cancer Study Group. , 1996, The New England journal of medicine.

[20]  B. Lin Nonparametric estimation of the gap time distributions for serial events with censored data , 1999 .

[21]  Jerald F. Lawless,et al.  Some Simple Robust Methods for the Analysis of Recurrent Events , 1995 .

[22]  J. Lawless Regression Methods for Poisson Process Data , 1987 .

[23]  J. Lawless Statistical Models and Methods for Lifetime Data , 2002 .

[24]  Luc Duchateau,et al.  Evolution of recurrent asthma event rate over time in frailty models , 2003 .

[25]  J. Robins,et al.  Semiparametric regression estimation in the presence of dependent censoring , 1995 .

[26]  M S Pepe,et al.  Modelling prevalence of a condition: chronic graft-versus-host disease after bone marrow transplantation. , 1997, Statistics in medicine.

[27]  T. Powles,et al.  Clodronate decreases the frequency of skeletal metastases in women with breast cancer. , 1996, Bone.

[28]  David V Glidden,et al.  Robust Inference for Event Probabilities with Non‐Markov Event Data , 2002, Biometrics.

[29]  T M Therneau,et al.  rhDNase as an example of recurrent event analysis. , 1997, Statistics in medicine.

[30]  L. J. Wei,et al.  Regression analysis of multivariate incomplete failure time data by modeling marginal distributions , 1989 .

[31]  R. Prentice,et al.  Commentary on Andersen and Gill's "Cox's Regression Model for Counting Processes: A Large Sample Study" , 1982 .

[32]  Zhiliang Ying,et al.  Nonparametric estimation of the gap time distributions for serial events with censored data , 1999 .

[33]  Margaret S. Pepe,et al.  Inference for Events with Dependent Risks in Multiple Endpoint Studies , 1991 .

[34]  J P Klein,et al.  Summary curves for patients transplanted for chronic myeloid leukaemia salvaged by a donor lymphocyte infusion: the current leukaemia‐free survival curve , 2000, British journal of haematology.

[35]  Debashis Ghosh,et al.  MARGINAL REGRESSION MODELS FOR RECURRENT AND TERMINAL EVENTS , 2002 .

[36]  D. Lin,et al.  Nonparametric Analysis of Recurrent Events and Death , 2000, Biometrics.

[37]  L J Wei,et al.  An overview of statistical methods for multiple failure time data in clinical trials. , 1997, Statistics in medicine.

[38]  Mei-Cheng Wang,et al.  Joint Modeling and Estimation for Recurrent Event Processes and Failure Time Data , 2004, Journal of the American Statistical Association.

[39]  O Borgan,et al.  Covariate Adjustment of Event Histories Estimated from Markov Chains: The Additive Approach , 2001, Biometrics.

[40]  Odd O. Aalen,et al.  Estimating stage occupation probabilities in non-Markov models , 2007, Lifetime data analysis.

[41]  Richard J. Cook,et al.  The Statistical Analysis of Recurrent Events , 2007 .

[42]  Analysis of Mean and Rate Functions for Recurrent Events , 1997 .

[43]  R J Cook,et al.  Marginal analysis of recurrent events and a terminating event. , 1997, Statistics in medicine.

[44]  D. Cox,et al.  The statistical analysis of series of events , 1966 .

[45]  J. Lawless Negative binomial and mixed Poisson regression , 1987 .

[46]  M. Crowder Classical Competing Risks , 2001 .

[47]  Donald Fraser,et al.  Nonparametric Estimation IV , 1951 .

[48]  G. Koch,et al.  Measures of Morbidity in Clinical Studies with Recurrent Skeletal Complications , 2004, Journal of biopharmaceutical statistics.

[49]  Jianwen Cai,et al.  Some Graphical Displays and Marginal Regression Analyses for Recurrent Failure Times and Time Dependent Covariates , 1993 .