Robust analyzes for longitudinal clinical trials with missing and non-normal continuous outcomes.

Missing data is unavoidable in longitudinal clinical trials, and outcomes are not always normally distributed. In the presence of outliers or heavy-tailed distributions, the conventional multiple imputation with the mixed model with repeated measures analysis of the average treatment effect (ATE) based on the multivariate normal assumption may produce bias and power loss. Control-based imputation (CBI) is an approach for evaluating the treatment effect under the assumption that participants in both the test and control groups with missing outcome data have a similar outcome profile as those with an identical history in the control group. We develop a robust framework to handle non-normal outcomes under CBI without imposing any parametric modeling assumptions. Under the proposed framework, sequential weighted robust regressions are applied to protect the constructed imputation model against non-normality in the covariates and the response variables. Accompanied by the subsequent mean imputation and robust model analysis, the resulting ATE estimator has good theoretical properties in terms of consistency and asymptotic normality. Moreover, our proposed method guarantees the analysis model robust-ness of the ATE estimation in the sense that its asymptotic results remain intact even when the analysis model is misspecified. The superiority of the proposed robust method is demonstrated by comprehensive simulation studies and an AIDS clinical trial data application.

Multiply robust estimators in longitudinal studies with missing data under control-based imputation.

Longitudinal studies are often subject to missing data. The recent guidance from regulatory agencies, such as the ICH E9(R1) addendum addresses the importance of defining a treatment effect estimand with the consideration of intercurrent events. Jump-to-reference (J2R) is one classical control-based scenario for the treatment effect evaluation, where the participants in the treatment group after intercurrent events are assumed to have the same disease progress as those with identical covariates in the control group. We establish new estimators to assess the average treatment effect based on a proposed potential outcomes framework under J2R. Various identification formulas are constructed, motivating estimators that rely on different parts of the observed data distribution. Moreover, we obtain a novel estimator inspired by the efficient influence function, with multiple robustness in the sense that it achieves n1/2-consistency if any pairs of multiple nuisance functions are correctly specified, or if the nuisance functions converge at a rate not slower than n-1/4 when using flexible modeling approaches. The finite-sample performance of the proposed estimators is validated in simulation studies and an antidepressant clinical trial.

A novel power prior approach for borrowing historical control data in clinical trials.

There has been an increased interest in borrowing information from historical control data to improve the statistical power for hypothesis testing, therefore reducing the required sample sizes in clinical trials. To account for the heterogeneity between the historical and current trials, power priors are often considered to discount the information borrowed from the historical data. However, it can be challenging to choose a fixed power prior parameter in the application. The modified power prior approach, which defines a random power parameter with initial prior to control the amount of historical information borrowed, may not directly account for heterogeneity between the trials. In this paper, we propose a novel approach to pick a power prior based on some direct measures of distributional differences between historical control data and current control data under normal assumptions. Simulations are conducted to investigate the performance of the proposed approach compared with current approaches (e.g. commensurate prior, meta-analytic-predictive, and modified power prior). The results show that the proposed power prior improves the study power while controlling the type I error within a tolerable limit when the distribution of the historical control data is similar to that of the current control data. The method is developed for both superiority and non-inferiority trials and is illustrated with an example from vaccine clinical trials.

SMIM: A unified framework of survival sensitivity analysis using multiple imputation and martingale.

Censored survival data are common in clinical trial studies. We propose a unified framework for sensitivity analysis to censoring at random in survival data using multiple imputation and martingale, called SMIM. The proposed framework adopts the δ-adjusted and control-based models, indexed by the sensitivity parameter, entailing censoring at random and a wide collection of censoring not at random assumptions. Also, it targets a broad class of treatment effect estimands defined as functionals of treatment-specific survival functions, taking into account missing data due to censoring. Multiple imputation facilitates the use of simple full-sample estimation; however, the standard Rubin's combining rule may overestimate the variance for inference in the sensitivity analysis framework. We decompose the multiple imputation estimator into a martingale series based on the sequential construction of the estimator and propose the wild bootstrap inference by resampling the martingale series. The new bootstrap inference has a theoretical guarantee for consistency and is computationally efficient compared to the nonparametric bootstrap counterpart. We evaluate the finite-sample performance of the proposed SMIM through simulation and an application on an HIV clinical trial.

Sensitivity analyses in longitudinal clinical trials via distributional imputation.

Missing data is inevitable in longitudinal clinical trials. Conventionally, the missing at random assumption is assumed to handle missingness, which however is unverifiable empirically. Thus, sensitivity analyses are critically important to assess the robustness of the study conclusions against untestable assumptions. Toward this end, regulatory agencies and the pharmaceutical industry use sensitivity models such as return-to-baseline, control-based, and washout imputation, following the ICH E9(R1) guidance. Multiple imputation is popular in sensitivity analyses; however, it may be inefficient and result in an unsatisfying interval estimation by Rubin's combining rule. We propose distributional imputation in sensitivity analysis, which imputes each missing value by samples from its target imputation model given the observed data. Drawn on the idea of Monte Carlo integration, the distributional imputation estimator solves the mean estimating equations of the imputed dataset. It is fully efficient with theoretical guarantees. Moreover, we propose weighted bootstrap to obtain a consistent variance estimator, taking into account the variabilities due to model parameter estimation and target parameter estimation. The superiority of the distributional imputation framework is validated in the simulation study and an antidepressant longitudinal clinical trial.

Analysis of time-to-event data using a flexible mixture model under a constraint of proportional hazards.

Cox proportional hazards (PH) model evaluates the effects of interested covariates under PH assumption without specified the baseline hazard. In clinical trial applications, however, the explicitly estimated hazard or cumulative survival function for each treatment group helps to assess and interpret the meaning of treatment difference. In this paper, we propose to use a flexible mixture model under the PH constraint to fit the underline survival functions. Simulations are conducted to evaluate its performance and show that the proposed mixture PH model is very similar to the Cox PH model in terms of estimating the hazard ratio, bias, confidence interval coverage, type-I error and testing power. Application to several real clinical trial examples demonstrates that the results from this approach are almost identical to the results from Cox PH model. The explicitly estimated hazard function for each treatment group provides additional useful information and helps the interpretation of hazard comparisons.

A flexible parametric survival model for fitting time to event data in clinical trials.

Time-to-event data are common in clinical trials to evaluate survival benefit of a new drug, biological product, or device. The commonly used parametric models including exponential, Weibull, Gompertz, log-logistic, log-normal, are simply not flexible enough to capture complex survival curves observed in clinical and medical research studies. On the other hand, the nonparametric Kaplan Meier (KM) method is very flexible and successful on catching the various shapes in the survival curves but lacks ability in predicting the future events such as the time for certain number of events and the number of events at certain time and predicting the risk of events (eg, death) over time beyond the span of the available data from clinical trials. It is obvious that neither the nonparametric KM method nor the current parametric distributions can fulfill the needs in fitting survival curves with the useful characteristics for predicting. In this paper, a full parametric distribution constructed as a mixture of three components of Weibull distribution is explored and recommended to fit the survival data, which is as flexible as KM for the observed data but have the nice features beyond the trial time, such as predicting future events, survival probability, and hazard function.