A global logrank test for adaptive treatment strategies based on observational studies

In studying adaptive treatment strategies, a natural question that is of paramount interest is whether there is any significant difference among all possible treatment strategies. When the outcome variable of interest is time-to-event, we propose an inverse probability weighted logrank test for testing the equivalence of a fixed set of pre-specified adaptive treatment strategies based on data from an observational study. The weights take into account both the possible selection bias in an observational study and the fact that the same subject may be consistent with more than one treatment strategy. The asymptotic distribution of the weighted logrank statistic under the null hypothesis is obtained. We show that, in an observational study where the treatment selection probabilities need to be estimated, the estimation of these probabilities does not have an effect on the asymptotic distribution of the weighted logrank statistic, as long as the estimation of the parameters in the models for these probabilities is n-consistent. Finite sample performance of the test is assessed via a simulation study. We also show in the simulation that the test can be pretty robust to misspecification of the models for the probabilities of treatment selection. The method is applied to analyze data on antidepressant adherence time from an observational database maintained at the Department of Veterans Affairs' Serious Mental Illness Treatment Research and Evaluation Center.

[1]  J M Robins,et al.  Marginal Mean Models for Dynamic Regimes , 2001, Journal of the American Statistical Association.

[2]  J. Robins,et al.  Marginal Structural Models and Causal Inference in Epidemiology , 2000, Epidemiology.

[3]  Alena I. Oetting,et al.  Statistical Methodology for a SMART Design in the Development of Adaptive Treatment Strategies , 2011 .

[4]  Abdus S Wahed,et al.  Sample size for two‐stage studies with maintenance therapy , 2009, Statistics in medicine.

[5]  Marie Davidian,et al.  Estimation of Survival Distributions of Treatment Policies in Two‐Stage Randomization Designs in Clinical Trials , 2002, Biometrics.

[6]  J. Robins,et al.  Marginal structural models to estimate the causal effect of zidovudine on the survival of HIV-positive men. , 2000, Epidemiology.

[7]  Susan A Murphy,et al.  Sample size formulae for two-stage randomized trials with survival outcomes. , 2011, Biometrika.

[8]  A. Wahed,et al.  Supremum weighted log-rank test and sample size for comparing two-stage adaptive treatment strategies , 2008 .

[9]  C. Glymour,et al.  STATISTICS AND CAUSAL INFERENCE , 1985 .

[10]  Marina Schmid,et al.  Design And Analysis Of Experiments For Statistical Selection Screening And Multiple Comparisons , 2016 .

[11]  A. Wahed,et al.  Weighted log-rank statistic to compare shared-path adaptive treatment strategies. , 2013, Biostatistics.

[12]  Bill Ravens,et al.  An Introduction to Copulas , 2000, Technometrics.

[13]  J. Robins,et al.  Estimation and extrapolation of optimal treatment and testing strategies , 2008, Statistics in medicine.

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

[15]  D. Kupfer,et al.  Sequenced treatment alternatives to relieve depression (STAR*D): rationale and design. , 2004, Controlled clinical trials.

[16]  Norman E. Breslow,et al.  A Z‐theorem with Estimated Nuisance Parameters and Correction Note for ‘Weighted Likelihood for Semiparametric Models and Two‐phase Stratified Samples, with Application to Cox Regression’ , 2008 .

[17]  A. Tsiatis.  Semiparametric Theory and Missing Data , 2006 .

[18]  S. Murphy,et al.  Developing adaptive treatment strategies in substance abuse research. , 2007, Drug and alcohol dependence.

[19]  Bin Nan,et al.  Weighted Likelihood Method for Grouped Survival Data in Case–Cohort Studies with Application to HIV Vaccine Trials , 2008, Biometrics.

[20]  P. Gänssler.  Weak Convergence and Empirical Processes - A. W. van der Vaart; J. A. Wellner. , 1997 .

[21]  Yuliya Lokhnygina,et al.  Cox Regression Methods for Two‐Stage Randomization Designs , 2007, Biometrics.

[22]  A. Wahed,et al.  Weighted Kaplan–Meier estimators for two‐stage treatment regimes , 2010, Statistics in medicine.

[23]  J. Robins.  A new approach to causal inference in mortality studies with a sustained exposure period—application to control of the healthy worker survivor effect , 1986 .

[24]  Jon A. Wellner,et al.  Weak Convergence and Empirical Processes: With Applications to Statistics , 1996 .

[25]  Ree Dawson,et al.  Flexible treatment strategies in chronic disease: clinical and research implications , 2000, Biological Psychiatry.

[26]  James M. Robins,et al.  Causal Inference from Complex Longitudinal Data , 1997 .

[27]  Guo Xiang,et al.  A Weighted Risk Set Estimator for Survival Distributions in Two-Stage Randomization Designs with Censored Survival Data , 2005 .

[28]  Bin Nan,et al.  Relative risk regression for current status data in case‐cohort studies , 2011 .

[29]  S. Murphy,et al.  An experimental design for the development of adaptive treatment strategies , 2005, Statistics in medicine.

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

[31]  J. Robins,et al.  The International Journal of Biostatistics CAUSAL INFERENCE Dynamic Regime Marginal Structural Mean Models for Estimation of Optimal Dynamic Treatment Regimes , Part I : Main Content , 2011 .

[32]  J. Robins,et al.  Estimation of Regression Coefficients When Some Regressors are not Always Observed , 1994 .

[33]  Ian W. McKeague,et al.  Comparing Distribution Functions via Empirical Likelihood , 2006 .

[34]  Anastasios A. Tsiatis,et al.  Semiparametric efficient estimation of survival distributions in two-stage randomisation designs in clinical trials with censored data , 2006 .