Stratifying Subjects for Treatment Selection with Censored Event Time Data from a Comparative Study

The conventional approach to comparing a new treatment with a standard therapy is often based on a summary measure for the treatment difference over the entire study population. A positive trial with respect to such a global measure, however, does not mean that all individual future patients would benefit from the new treatment. On the other hand, a negative finding may not be sufficiently conclusive to claim that the new treatment is entirely futile. In this article, we propose a systematic approach to identify future patients who would benefit from the new treatment with respect to an event time outcome via a two-stage inference procedure. We first develop a scoring index to stratify study patients based on parametric or semiparametric survival models with the observed event times and covariates. We then use a nonparametric method to estimate the average treatment difference for each stratum defined by the score. Sampling variation of the resulting estimator is also provided across the entire spectrum of the score by controlling certain local and global error rates. With a numerical study, we show that the new proposal performs well under various practical settings. Our method is illustrated with the data from a recent clinical trial to evaluate whether a specific anti-hypertensive drug would prolong the lives for patients with stable coronary artery disease and normal or slightly reduced left ventricular function.

[1]  Z. Ying,et al.  Checking the Cox model with cumulative sums of martingale-based residuals , 1993 .

[2]  Jianqing Fan,et al.  Data‐Driven Bandwidth Selection in Local Polynomial Fitting: Variable Bandwidth and Spatial Adaptation , 1995 .

[3]  Stephen W Lagakos,et al.  Statistics in medicine--reporting of subgroup analyses in clinical trials. , 2007, The New England journal of medicine.

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

[5]  Nils Lid Hjort,et al.  On inference in parametric survival data models , 1992 .

[6]  M S Pepe,et al.  Weighted Kaplan-Meier statistics: a class of distance tests for censored survival data. , 1989, Biometrics.

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

[8]  Tianxi Cai,et al.  Semiparametric Box–Cox power transformation models for censored survival observations , 2005 .

[9]  Michael R Kosorok,et al.  Simultaneous Inferences on the Contrast of Two Hazard Functions with Censored Observations , 2002, Biometrics.

[10]  Enno Mammen,et al.  Bootstrap, wild bootstrap, and asymptotic normality , 1992 .

[11]  R. Gill Large Sample Behaviour of the Product-Limit Estimator on the Whole Line , 1983 .

[12]  Irène Gijbels,et al.  Local likelihood and local partial likelihood in hazard regression , 1997 .

[13]  John D. Kalbfleisch,et al.  Estimation of the average hazard ratio , 1981 .

[14]  M. Pfeffer,et al.  Angiotensin-converting-enzyme inhibition in stable coronary artery disease. , 2004, The New England journal of medicine.

[15]  M. Pepe The Statistical Evaluation of Medical Tests for Classification and Prediction , 2003 .

[16]  L. Tian,et al.  Analysis of randomized comparative clinical trial data for personalized treatment selections. , 2011, Biostatistics.

[17]  Changbao Wu,et al.  Jackknife, Bootstrap and Other Resampling Methods in Regression Analysis , 1986 .

[18]  Thomas R. Fleming,et al.  Weighted Kaplan‐Meier Statistics: Large Sample and Optimality Considerations , 1991 .

[19]  Robert Kohn,et al.  Simple Transformation Techniques for Improved Non‐parametric Regression , 1997 .

[20]  Margaret Sullivan Pepe,et al.  Evaluating Markers for Selecting a Patient's Treatment , 2004, Biometrics.

[21]  Yuhyun Park,et al.  Estimating subject-specific survival functions under the accelerated failure time model , 2003 .

[22]  M. Bonetti,et al.  A graphical method to assess treatment-covariate interactions using the Cox model on subsets of the data. , 2000, Statistics in medicine.

[23]  M. Pfeffer,et al.  Prevention of events with angiotensin-converting enzyme inhibition (the PEACE study design). Prevention of Events with Angiotensin-Converting Enzyme Inhibition. , 1998, The American journal of cardiology.

[24]  Hajime Uno,et al.  Calibrating parametric subject-specific risk estimation. , 2010, Biometrika.

[25]  James Stephen Marron,et al.  Transformations in Density Estimation , 1991 .

[26]  Marc A Pfeffer,et al.  The charisma of subgroups and the subgroups of CHARISMA. , 2006, The New England journal of medicine.

[27]  Hani Doss,et al.  An Approach to Nonparametric Regression for Life History Data Using Local Linear Fitting , 1995 .

[28]  A A Tsiatis,et al.  Sequential Methods for Comparing Years of Life Saved in the Two‐Sample Censored Data Problem , 1999, Biometrics.

[29]  L. J. Wei,et al.  On the Cox Model With Time-Varying Regression Coefficients , 2005 .

[30]  Lihui Zhao,et al.  Utilizing the integrated difference of two survival functions to quantify the treatment contrast for designing, monitoring, and analyzing a comparative clinical study , 2012, Clinical trials.

[31]  Yi Li,et al.  Estimating Subject‐Specific Dependent Competing Risk Profile with Censored Event Time Observations , 2011, Biometrics.

[32]  D. Pollard Empirical Processes: Theory and Applications , 1990 .

[33]  P. Rothwell,et al.  External validity of randomised controlled trials: “To whom do the results of this trial apply?” , 2005, The Lancet.

[34]  Els Goetghebeur,et al.  Model evaluation based on the sampling distribution of estimated absolute prediction error , 2007 .