Search

A difficult decision for patients in need of kidney–pancreas transplant is whether to seek a living kidney donor or wait to receive both organs from one deceased donor. The framework of dynamic treatment regimes (DTRs) can inform this choice, but a patient‐relevant strategy such as "wait for deceased‐donor transplant" is ill‐defined because there are multiple versions of treatment (i.e., wait times, organ qualities). Existing DTR methods average over the distribution of treatment versions in the data, estimating survival under a "representative intervention." This is undesirable if transporting inferences to a target population such as patients today, who experience shorter wait times thanks to evolutions in allocation policy. We, therefore, propose the concept of a generalized representative intervention (GRI): a random DTR that assigns treatment version by drawing from the distribution among strategy compliers in the target population (e.g., patients today). We describe an inverse‐probability‐weighted product‐limit estimator of survival under a GRI that performs well in simulations and can be implemented in standard statistical software. For continuous treatments (e.g., organ quality), weights are reformulated to depend on probabilities only, not densities. We apply our method to a national database of kidney–pancreas transplant candidates from 2001–2020 to illustrate that variability in transplant rate across years and centers results in qualitative differences in the optimal strategy for patient survival. [ABSTRACT FROM AUTHOR]

Individuals can vary drastically in their response to the same treatment, and this heterogeneity has driven the push for more personalized medicine. Accurate and interpretable methods to identify subgroups that respond to the treatment differently from the population average are necessary to achieving this goal. The Virtual Twins (VT) method is a highly cited and implemented method for subgroup identification because of its intuitive framework. However, since its initial publication, many researchers still rely heavily on the authors' initial modeling suggestions without examining newer and more powerful alternatives. This leaves much of the potential of the method untapped. We comprehensively evaluate the performance of VT with different combinations of methods in each of its component steps, under a collection of linear and nonlinear problem settings. Our simulations show that the method choice for Step 1 of VT, in which dense models with high predictive performance are fit for the potential outcomes, is highly influential in the overall accuracy of the method, and Superlearner is a promising choice. We illustrate our findings by using VT to identify subgroups with heterogeneous treatment effects in a randomized, double-blind trial of very low nicotine content cigarettes. [ABSTRACT FROM AUTHOR]

The increasing multiplicity of data sources offers exciting possibilities in estimating the effects of a treatment, intervention, or exposure, particularly if observational and experimental sources could be used simultaneously. Borrowing between sources can potentially result in more efficient estimators, but it must be done in a principled manner to mitigate increased bias and Type I error. Furthermore, when the effect of treatment is confounded, as in observational sources or in clinical trials with noncompliance, causal effect estimators are needed to simultaneously adjust for confounding and to estimate effects across data sources. We consider the problem of estimating causal effects from a primary source and borrowing from any number of supplemental sources. We propose using regression‐based estimators that borrow based on assuming exchangeability of the regression coefficients and parameters between data sources. Borrowing is accomplished with multisource exchangeability models and Bayesian model averaging. We show via simulation that a Bayesian linear model and Bayesian additive regression trees both have desirable properties and borrow under appropriate circumstances. We apply the estimators to recently completed trials of very low nicotine content cigarettes investigating their impact on smoking behavior. [ABSTRACT FROM AUTHOR]

Unbiased estimation of causal effects with observational data requires adjustment for confounding variables that are related to both the outcome and treatment assignment. Standard variable selection techniques aim to maximize predictive ability of the outcome model, but they ignore covariate associations with treatment and may not adjust for important confounders weakly associated to outcome. We propose a novel method for estimating causal effects that simultaneously considers models for both outcome and treatment, which we call the bilevel spike and slab causal estimator (BSSCE). By using a Bayesian formulation, BSSCE estimates the posterior distribution of all model parameters and provides straightforward and reliable inference. Spike and slab priors are used on each covariate coefficient which aim to minimize the mean squared error of the treatment effect estimator. Theoretical properties of the treatment effect estimator are derived justifying the prior used in BSSCE. Simulations show that BSSCE can substantially reduce mean squared error over numerous methods and performs especially well with large numbers of covariates, including situations where the number of covariates is greater than the sample size. We illustrate BSSCE by estimating the causal effect of vasoactive therapy vs. fluid resuscitation on hypotensive episode length for patients in the Multiparameter Intelligent Monitoring in Intensive Care III critical care database. [ABSTRACT FROM AUTHOR]

Estimating causal effects from randomized controlled trials is often complicated due to participant noncompliance to randomized treatment. Although there are a variety of methods to estimate causal effects in the presence of noncompliance, they generally make the assumption that noncompliance is measured without error. This is frequently an untenable assumption, particularly when noncompliance is based on participant self-report. To overcome this issue, we treat compliance as an unobserved variable and show how to estimate the probability of compliance given a biomarker of treatment and the other observed data. We present inverse probability weighted estimators, regression-based estimators, and a doubly-robust augmented estimator that rely on the estimated probability of compliance rather than an indicator of compliance. We investigate the finite-sample properties of the estimators and their efficiency and robustness under correctly specified or misspecified models, and we apply the estimators to a recently completed trial of very low nicotine content cigarettes. [ABSTRACT FROM AUTHOR]

Summary: Patients awaiting cadaveric organ transplantation face a difficult decision if offered a low‐quality organ: accept the organ or remain on the waiting list and hope a better organ is offered in the future. A dynamic treatment regime (DTR) for transplantation is a rule that determines whether a patient should decline an offered organ. Existing methods can estimate the effect of DTRs on survival outcomes, but these were developed for applications where treatment is abundantly available. For transplantation, organ availability is limited, and existing methods can only estimate the effect of a DTR assuming a single patient follows the DTR. We show for transplantation that the effect of a DTR depends on whether other patients follow the DTR. To estimate the anticipated survival if the entire population awaiting transplantation were to adopt a DTR, we develop a novel inverse probability weighted estimator (IPCW) which re‐weights patients based on the probability of following their transplant history in the counterfactual world in which all patients follow the DTR of interest. We estimate this counterfactual probability using hot deck imputation to fill in data that is not observed for patients who are artificially censored by IPCW once they no longer follow the DTR of interest. We show via simulation that our proposed method has good finite‐sample properties, and we apply our method to a lung transplantation observational registry. [ABSTRACT FROM AUTHOR]

Summary: The efficiency of doubly robust estimators of the average causal effect (ACE) of a treatment can be improved by including in the treatment and outcome models only those covariates which are related to both treatment and outcome (i.e., confounders) or related only to the outcome. However, it is often challenging to identify such covariates among the large number that may be measured in a given study. In this article, we propose GLiDeR (Group Lasso and Doubly Robust Estimation), a novel variable selection technique for identifying confounders and predictors of outcome using an adaptive group lasso approach that simultaneously performs coefficient selection, regularization, and estimation across the treatment and outcome models. The selected variables and corresponding coefficient estimates are used in a standard doubly robust ACE estimator. We provide asymptotic results showing that, for a broad class of data generating mechanisms, GLiDeR yields a consistent estimator of the ACE when either the outcome or treatment model is correctly specified. A comprehensive simulation study shows that GLiDeR is more efficient than doubly robust methods using standard variable selection techniques and has substantial computational advantages over a recently proposed doubly robust Bayesian model averaging method. We illustrate our method by estimating the causal treatment effect of bilateral versus single‐lung transplant on forced expiratory volume in one year after transplant using an observational registry. [ABSTRACT FROM AUTHOR]

Noncompliance or non-adherence to randomized treatment is a common challenge when interpreting data from randomized clinical trials. The effect of an intervention if all participants were forced to comply with the assigned treatment (i.e., the causal effect) is often of primary scientific interest. For example, in trials of very low nicotine content (VLNC) cigarettes, policymakers are interested in their effect on smoking behavior if their use were to be compelled by regulation. A variety of statistical methods to estimate the causal effect of an intervention have been proposed, but these methods, including inverse probability of compliance weighted (IPCW) estimators, assume that participants' compliance statuses are reported without error. This is an untenable assumption when compliance is based on self-report. Biomarkers (e.g., nicotine levels in the urine) may provide more reliable indicators of compliance but cannot perfectly discriminate between compliers and non-compliers. However, by modeling the distribution of the biomarker as a mixture distribution and writing the probability of compliance as a function of the mixture components, we show how the probability of compliance can be directly estimated from the data even when compliance status is unknown. To estimate the causal effect, we develop a new approach which re-weights participants by the product of their probability of compliance given the observed data and the inverse probability of compliance given confounders. We show that our proposed estimator is consistent and asymptotically normal and show that in some situations the proposed approach is more efficient than standard IPCW estimators. We demonstrate via simulation that the proposed estimator achieves smaller bias and greater efficiency than ad hoc approaches to estimating the causal effect when compliance is measured with error. We apply our method to data from a recently completed randomized trial of VLNC cigarettes. [ABSTRACT FROM AUTHOR]

Because the number of patients waiting for organ transplants exceeds the number of organs available, a better understanding of how transplantation affects the distribution of residual lifetime is needed to improve organ allocation. However, there has been little work to assess the survival benefit of transplantation from a causal perspective. Previous methods developed to estimate the causal effects of treatment in the presence of time-varying confounders have assumed that treatment assignment was independent across patients, which is not true for organ transplantation. We develop a version of G-estimation that accounts for the fact that treatment assignment is not independent across individuals to estimate the parameters of a structural nested failure time model. We derive the asymptotic properties of our estimator and confirm through simulation studies that our method leads to valid inference of the effect of transplantation on the distribution of residual lifetime. We demonstrate our method on the survival benefit of lung transplantation using data from the United Network for Organ Sharing. [ABSTRACT FROM AUTHOR]