Compliance subsampling designs for comparative research: estimation and optimal planning

For studies with treatment noncompliance, analyses have been developed recently to better estimate treatment efficacy. However, the advantage and cost of measuring compliance data have implications on the study design that have not been as systematically explored. In order to estimate better treatment efficacy with lower cost, we propose a new class of compliance subsampling (CSS) designs where, after subjects are assigned treatment, compliance behavior is measured for only subgroups of subjects. The sizes of the subsamples are allowed to relate to the treatment assignment, the assignment probability, the total sample size, the anticipated distributions of outcome and compliance, and the cost parameters of the study. The CSS design methods relate to prior work (i) on two-phase designs in which a covariate is subsampled and (ii) on causal inference because the subsampled postrandomization compliance behavior is not the true covariate of interest. For each CSS design, we develop efficient estimation of treatment efficacy under binary outcome and all-or-none observed compliance. Then we derive a minimal cost CSS design that achieves a required precision for estimating treatment efficacy. We compare the properties of the CSS design to those of conventional protocols in a study of patient choices for medical care at the end of life.     

The compliance score as a regressor in randomized trials

The compliance score in randomized trials is a measure of the effect of randomization on treatment received. It is in principle a group-level pretreatment variable and so can be used where individual-level measures of treatment received can produce misleading inferences. The interpretation of models with the compliance score as a regressor of interest depends on the link function. Using the identity link can lead to valid inference about the effects of treatment received even in the presence of nonrandom noncompliance; such inference is more problematic for nonlinear links. We illustrate these points with data from two randomized trials.     

Likelihood methods for treatment noncompliance and subsequent nonresponse in randomized trials

While several new methods that account for noncompliance or missing data in randomized trials have been proposed, the dual effects of noncompliance and nonresponse are rarely dealt with simultaneously. We construct a maximum likelihood estimator (MLE) of the causal effect of treatment assignment for a two-armed randomized trial assuming all-or-none treatment noncompliance and allowing for subsequent nonresponse. The EM algorithm is used for parameter estimation. Our likelihood procedure relies on a latent compliance state covariate that describes the behavior of a subject under all possible treatment assignments and characterizes the missing data mechanism as in Frangakis and Rubin (1999, Biometrika 86, 365-379). Using simulated data, we show that the MLE for normal outcomes compares favorably to the method-of-moments (MOM) and the standard intention-to-treat (ITT) estimators under (1) both normal and non-normal data, and (2) departures from the latent ignorability and compound exclusion restriction assumptions. We illustrate methods using data from a trial to compare the efficacy of two antipsychotics for adults with refractory schizophrenia.     

Multiple imputation methods for treatment noncompliance and nonresponse in randomized clinical trials

Summary .  Randomized clinical trials are a powerful tool for investigating causal treatment effects, but in human trials there are oftentimes problems of noncompliance which standard analyses, such as the intention-to-treat or as-treated analysis, either ignore or incorporate in such a way that the resulting estimand is no longer a causal effect. One alternative to these analyses is the complier average causal effect (CACE) which estimates the average causal treatment effect among a subpopulation that would comply under any treatment assigned. We focus on the setting of a randomized clinical trial with crossover treatment noncompliance (e.g., control subjects could receive the intervention and intervention subjects could receive the control) and outcome nonresponse. In this article, we develop estimators for the CACE using multiple imputation methods, which have been successfully applied to a wide variety of missing data problems, but have not yet been applied to the potential outcomes setting of causal inference. Using simulated data we investigate the finite sample properties of these estimators as well as of competing procedures in a simple setting. Finally we illustrate our methods using a real randomized encouragement design study on the effectiveness of the influenza vaccine.     

Analysis of times to repeated events in two-arm randomized trials with noncompliance and dependent censoring

This article develops randomization-based methods for times to repeated events in two-arm randomized trials with noncompliance and dependent censoring. Structural accelerated failure time models are assumed to capture causal effects on repeated event times and dependent censoring time, but the dependence structure among repeated event times and dependent censoring time is unspecified. Artificial censoring techniques to accommodate nonrandom noncompliance and dependent censoring are proposed. Estimation of the acceleration parameters are based on rank-based estimating functions. A simulation study is conducted to evaluate the performance of the developed methods. An illustration of the methods using data from an acute myeloid leukemia trial is provided.     

Estimating causal treatment effects from longitudinal HIV natural history studies using marginal structural models

Several recently completed and ongoing studies of the natural history of HIV infection have generated a wealth of information about its clinical progression and how this progression is altered by therepeutic interventions and environmental factors. Natural history studies typically follow prospective cohort designs, and enroll large numbers of participants for long-term prospective follow-up (up to several years). Using data from the HIV Epidemiology Research Study (HERS), a six-year natural history study that enrolled 871 HIV-infected women starting in 1993, we investigate the therapeutic effect of highly active antiretroviral therapy regimens (HAART) on CD4 cell count using the marginal structural modeling framework and associated estimation procedures based on inverse-probability weighting (developed by Robins and colleagues). To evaluate treatment effects from a natural history study, specialized methods are needed because treatments are not randomly prescribed and, in particular, the treatment-response relationship can be confounded by variables that are time-varying. Our analysis uses CD4 data on all follow-up visits over a two-year period, and includes sensitivity analyses to investigate potential biases attributable to unmeasured confounding. Strategies for selecting ranges of a sensitivity parameter are given, as are intervals for treatment effect that reflect uncertainty attributable both to sampling and to lack of knowledge about the nature and existence of unmeasured confounding. To our knowledge, this is the first use in "real data" of Robins's sensitivity analysis for unmeasured confounding (Robins, 1999a, Synthese 121, 151-179). The findings from our analysis are consistent with recent treatment guidelines set by the U.S. Panel of the International AIDS Society (Carpenter et al., 2000, Journal of the American Medical Association 280, 381-391).     

Assessing the effect of an influenza vaccine in an encouragement design

Many randomized experiments suffer from noncompliance. Some of these experiments, so-called encouragement designs, can be expected to have especially large amounts of noncompliance, because encouragement to take the treatment rather than the treatment itself is randomly assigned to individuals. We present an extended framework for the analysis of data from such experiments with a binary treatment, binary encouragement, and background covariates. There are two key features of this framework: we use an instrumental variables approach to link intention-to-treat effects to treatment effects and we adopt a Bayesian approach for inference and sensitivity analysis. This framework is illustrated in a medical example concerning the effects of inoculation for influenza. In this example, the analyses suggest that positive estimates of the intention-to-treat effect need not be due to the treatment itself, but rather to the encouragement to take the treatment: the intention-to-treat effect for the subpopulation who would be inoculated whether or not encouraged is estimated to be approximately as large as the intention-to-treat effect for the subpopulation whose inoculation status would agree with their (randomized) encouragement status whether or not encouraged. Thus, our methods suggest that global intention-to-treat estimates, although often regarded as conservative, can be too coarse and even misleading when taken as summarizing the evidence in the data for the effects of treatments.     

Application of the principal stratification approach to the Faenza randomized experiment on breast self-examination

In this article we present an extended framework based on the principal stratification approach (Frangakis and Rubin, 2002, Biometrics 58, 21-29), for the analysis of data from randomized experiments which suffer from treatment noncompliance, missing outcomes following treatment noncompliance, and "truncation by death." We are not aware of any previous work that addresses all these complications jointly. This framework is illustrated in the context of a randomized trial of breast self-examination.     

Stratified analysis in randomized trials with noncompliance

This article develops methods for stratified analyses of additive or multiplicative causal effect on binary outcomes in randomized trials with noncompliance. The methods are based on a weighted estimating function for an unbiased estimating function under randomization in each stratum. When known weights are used, the derived estimator is a natural extension of the instrumental variable estimator for stratified analyses, and test-based confidence limits are solutions of a quadratic equation in the causal parameter. Optimal weights that maximize asymptotic efficiency incorporate variability in compliance aspects across strata. An assessment based on asymptotic relative efficiency shows that a substantial enhancement in efficiency can be gained by using optimal weights instead of conventional ones, which do not incorporate the variability in compliance aspects across strata. Application to a field trial for coronary heart disease is provided.     

Semiparametric regression estimation in the presence of dependent censoring

We propose a semiparametric estimation procedure for estimatingthe regression of an outcome Y, measured at the end of a fixedfollow-up period, on baseline explanatory variables X, measuredprior to start of follow-up, in the presence of dependent censoringgiven X. The proposed estimators are consistent when the dataare ‘missing at random’ but not ‘missing completelyat random’ (Rubin, 1976), and do not require full specificationof the complete data likelihood. Specifically, we assume thatthe probability of censoring at time t is independent of theoutcome Y conditional on the recorded history up to t of a vectorof time-dependent covariates that are correlated with Y. Ourestimators can be used to adjust for dependent censoring andnonrandom noncompliance in randomised trials studying the effectof a treatment on the mean of a response variable of interest.Even with independent censoring, our methods allow the investigatorto increase efficiency by exploiting the correlation of theoutcome with a vector of time-dependent covariates.