Assessing the effect of a partly unobserved, exogenous, binary time-dependent covariate on survival probabilities using generalised pseudo-values

BackgroundInvestigating the impact of a time-dependent intervention on the probability of long-term survival is statistically challenging. A typical example is stem-cell transplantation performed after successful donor identification from registered donors. Here, a suggested simple analysis based on the exogenous donor availability status according to registered donors would allow the estimation and comparison of survival probabilities. As donor search is usually ceased after a patient’s event, donor availability status is incompletely observed, so that this simple comparison is not possible and the waiting time to donor identification needs to be addressed in the analysis to avoid bias. It is methodologically unclear, how to directly address cumulative long-term treatment effects without relying on proportional hazards while avoiding waiting time bias.MethodsThe pseudo-value regression technique is able to handle the first two issues; a novel generalisation of this technique also avoids waiting time bias. Inverse-probability-of-censoring weighting is used to account for the partly unobserved exogenous covariate donor availability.ResultsSimulation studies demonstrate unbiasedness and satisfying coverage probabilities of the new method. A real data example demonstrates that study results based on generalised pseudo-values have a clear medical interpretation which supports the clinical decision making process.ConclusionsThe proposed generalisation of the pseudo-value regression technique enables to compare survival probabilities between two independent groups where group membership becomes known over time and remains partly unknown. Hence, cumulative long-term treatment effects are directly addressed without relying on proportional hazards while avoiding waiting time bias.

[1]  Ralf Bender,et al.  Generating survival times to simulate Cox proportional hazards models , 2005, Statistics in medicine.

[2]  M. Tallman,et al.  In adults with standard-risk acute lymphoblastic leukemia, the greatest benefit is achieved from a matched sibling allogeneic transplantation in first complete remission, and an autologous transplantation is less effective than conventional consolidation/maintenance chemotherapy in all patients: fin , 2006, Blood.

[3]  H. Heimpel,et al.  Survival with bone marrow transplantation versus hydroxyurea or interferon for chronic myelogenous leukemia. The German CML Study Group. , 1998, Blood.

[4]  R. Gray,et al.  How to avoid bias when comparing bone marrow transplantation with chemotherapy. , 1991, Bone marrow transplantation.

[5]  T. Martinussen,et al.  A Note on the Large Sample Properties of Estimators Based on Generalized Linear Models for Correlated Pseudo‐observations , 2016 .

[6]  H. Heinzl,et al.  Gaining more flexibility in Cox proportional hazards regression models with cubic spline functions. , 1997, Computer methods and programs in biomedicine.

[7]  John P. Klein,et al.  SAS and R functions to compute pseudo-values for censored data regression , 2008, Comput. Methods Programs Biomed..

[8]  B. Dong,et al.  Empirical Likelihood for Cumulative Hazard Ratio Estimation with Covariate Adjustment , 2012, Biometrics.

[9]  M. Schumacher,et al.  On pseudo-values for regression analysis in competing risks models , 2009, Lifetime data analysis.

[10]  P. Sasieni,et al.  Some new estimators for Cox regression , 1993 .

[11]  Maja Pohar Perme,et al.  Checking hazard regression models using pseudo‐observations , 2008, Statistics in medicine.

[12]  J R Anderson,et al.  Analysis of survival by tumor response. , 1983, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[13]  E Marubini,et al.  Analysing Survival Data from Clinical Trials and Observational Studies; E. Marubini & M. G. Valsecchi Published by John Wiley & Sons 414 pages ISBN 0-971-93987-0. , 1995, British journal of clinical pharmacology.

[14]  A. Perperoglou Reduced rank hazard regression with fixed and time‐varying effects of the covariates , 2013, Biometrical journal. Biometrische Zeitschrift.

[15]  Inverse Probability Weighting in Survival Analysis , 2014 .

[16]  James M. Robins,et al.  Inverse Probability Weighting in Survival Analysis , 2005 .

[17]  Richard Gray,et al.  Commentary: Mendelian randomization--an update on its use to evaluate allogeneic stem cell transplantation in leukaemia. , 2004, International journal of epidemiology.

[18]  Richard Sposto,et al.  Cure model analysis in cancer: an application to data from the Children's Cancer Group , 2002, Statistics in medicine.

[19]  J. Passweg,et al.  Outcome of patients with acquired aplastic anemia given first line bone marrow transplantation or immunosuppressive treatment in the last decade: a report from the European Group for Blood and Marrow Transplantation (EBMT). , 2007, Haematologica.

[20]  R Simon,et al.  A non-parametric graphical representation of the relationship between survival and the occurrence of an event: application to responder versus non-responder bias. , 1984, Statistics in medicine.

[21]  P. Royston,et al.  A New Proposal for Multivariable Modelling of Time‐Varying Effects in Survival Data Based on Fractional Polynomial Time‐Transformation , 2007, Biometrical journal. Biometrische Zeitschrift.

[22]  D. Schaubel,et al.  Estimating Cumulative Treatment Effects in the Presence of Nonproportional Hazards , 2008, Biometrics.

[23]  J O'Quigley,et al.  Estimating average regression effect under non-proportional hazards. , 2000, Biostatistics.

[24]  John P. Klein,et al.  Regression Analysis of Restricted Mean Survival Time Based on Pseudo-Observations , 2004, Lifetime data analysis.

[25]  Maja Pohar Perme,et al.  Pseudo-observations in survival analysis , 2010, Statistical methods in medical research.

[26]  Qi Jiang,et al.  Illustrating the Impact of a Time-Varying Covariate With an Extended Kaplan-Meier Estimator , 2005 .

[27]  Martin Schumacher,et al.  Multi-state model for studying an intermediate event using time-dependent covariates: application to breast cancer , 2013, BMC Medical Research Methodology.

[28]  Patrick Royston,et al.  Comparison between splines and fractional polynomials for multivariable model building with continuous covariates: a simulation study with continuous response , 2013, Statistics in medicine.

[29]  M. Schrappe,et al.  Chemotherapy versus allogeneic transplantation for very-high-risk childhood acute lymphoblastic leukaemia in first complete remission: comparison by genetic randomisation in an international prospective study , 2005, The Lancet.

[30]  Per Kragh Andersen,et al.  Regression Analysis for Multistate Models Based on a Pseudo‐value Approach, with Applications to Bone Marrow Transplantation Studies , 2007 .

[31]  J. Klein,et al.  Comparing Treatments in the Presence of Crossing Survival Curves: An Application to Bone Marrow Transplantation , 2008, Biometrics.

[32]  C. Pui,et al.  Clinical outcome of children with newly diagnosed Philadelphia chromosome-positive acute lymphoblastic leukemia treated between 1995 and 2005. , 2008, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[33]  James F. Watkins,et al.  Analysing Survival Data from Clinical Trials and Observational Studies. , 1995 .

[34]  M. Overgaard,et al.  Asymptotic theory of generalized estimating equations based on jack-knife pseudo-observations , 2017 .

[35]  Mei-Jie Zhang,et al.  Statistical Challenges in Comparing Chemotherapy and Bone Marrow Transplantation as a Treatment for Leukemia , 1996 .

[36]  J. Klein,et al.  Generalised linear models for correlated pseudo‐observations, with applications to multi‐state models , 2003 .

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

[38]  J. Kalbfleisch,et al.  The Statistical Analysis of Failure Time Data , 1980 .

[39]  Gilbert MacKenzie,et al.  The Statistical Analysis of Failure Time Data , 1982 .

[40]  D. Schaubel,et al.  Double Inverse‐Weighted Estimation of Cumulative Treatment Effects Under Nonproportional Hazards and Dependent Censoring , 2011, Biometrics.

[41]  M. Schemper,et al.  The estimation of average hazard ratios by weighted Cox regression , 2009, Statistics in medicine.

[42]  Martin Schumacher,et al.  Competing Risks and Multistate Models with R , 2011 .

[43]  H. Putter,et al.  Dynamic Prediction in Clinical Survival Analysis , 2011 .

[44]  John P Klein,et al.  Analyzing survival curves at a fixed point in time , 2007, Statistics in medicine.