Faculty Bibliography

Studies often report estimates of the average treatment effect (ATE). While the ATE summarizes the effect of a treatment on average, it does not provide any information about the effect of treatment within any individual. A treatment strategy that uses an individual's information to tailor treatment to maximize benefit is known as an optimal dynamic treatment rule (ODTR). Treatment, however, is typically not limited to a single point in time; consequently, learning an optimal rule for a time-varying treatment may involve not just learning the extent to which the comparative treatments' benefits vary across the characteristics of individuals, but also learning the extent to which the comparative treatments' benefits vary as relevant circumstances evolve within an individual. The goal of this paper is to provide a tutorial for estimating ODTR from longitudinal observational and clinical trial data for applied researchers. We describe an approach that uses a doubly-robust unbiased transformation of the conditional average treatment effect. We then learn a time-varying ODTR for when to increase buprenorphine-naloxone (BUP-NX) dose to minimize return-to-regular-opioid-use among patients with opioid use disorder. Our analysis highlights the utility of ODTRs in the context of sequential decision making: the learned ODTR outperforms a clinically defined strategy.

INTRODUCTION/BACKGROUND:People with chronic pain are at increased risk of opioid misuse. Less is known about the unique risk conferred by each pain management treatment, as treatments are typically implemented together, confounding their independent effects. This study estimated the extent to which pain management treatments were associated with risk of opioid use disorder (OUD) for those with chronic pain, controlling for baseline demographic and clinical confounding variables and holding other pain management treatments at their observed levels. METHODS:Data were analyzed in 2024 from 2 chronic pain subgroups within a cohort of non-pregnant Medicaid patients aged 35-64 years, 2016-2019, from 25 states: those with (1) chronic pain and physical disability (CPPD) (N=6,133) or (2) chronic pain without disability (CP) (N=67,438). Nine pain management treatments were considered: prescription opioid (1) dose and (2) duration; (3) number of opioid prescribers; opioid co-prescription with (4) benzo- diazepines, (5) muscle relaxants, and (6) gabapentinoids; (7) nonopioid pain prescription, (8) physical therapy, and (9) other pain treatment modality. The outcome was OUD risk. RESULTS:Having opioids co-prescribed with gabapentin or benzodiazepine was statistically significantly associated with a 37-45% increased OUD risk for the CP subgroup. Opioid dose and duration also were significantly associated with increased OUD risk in this subgroup. Physical therapy was significantly associated with an 18% decreased risk of OUD in the CP subgroup. DISCUSSION/CONCLUSIONS:Coprescription of opioids with either gabapentin or benzodiazepines may substantially increase OUD risk. More positively, physical therapy may be a relatively accessible and safe pain management strategy.

BACKGROUND:Determining a therapeutic window for maintaining antiretroviral drug concentrations within an appropriate range is required for identifying effective dosing regimens. The limits of this window are typically calculated using predictive models. We propose that target concentrations should instead be calculated based on counterfactual probabilities of relevant outcomes and describe a counterfactual framework for this. METHODS:The proposed framework is applied in an analysis including longitudinal observational data from 125 HIV-positive children treated with efavirenz-based regimens within the CHAPAS-3 trial, which enrolled children < 13 years in Zambia/Uganda. A directed acyclic graph was developed to visualize the mechanisms affecting antiretroviral concentrations. Causal concentration-response curves, adjusted for measured time-varying confounding of weight and adherence, are calculated using g-computation. RESULTS:The estimated curves show that higher concentrations during follow-up, 12/24 h after dose, lead to lower probabilities of viral failure (> 100 c/mL) at 96 weeks of follow-up. Estimated counterfactual failure probabilities under the current target range of 1-4 mg/L range from 24% to about 2%. The curves are almost identical for slow, intermediate and extensive metabolizers and show that a mid-dose concentration level of ≥ 3.5 mg/L would be required to achieve a failure probability of < 5%. CONCLUSIONS:Our analyses demonstrate that a causal approach may lead to different minimum concentration limits than analyses that are based on purely predictive models. Moreover, the approach highlights that indirect causes of failure, such as patients' metabolizing status, may predict patients' failure risk, but do not alter the threshold at which antiviral activity of efavirenz is severely reduced.

Mediation analysis is appealing for its ability to improve understanding of the mechanistic drivers of causal effects, but real-world data complexities challenge its successful implementation, including (i) the existence of post-exposure variables that also affect mediators and outcomes (thus, confounding the mediator-outcome relationship), that may also be (ii) multivariate, and (iii) the existence of multivariate mediators. All three challenges are present in the mediation analysis we consider here, where our goal is to estimate the indirect effects of receiving a Section 8 housing voucher as a young child on the risk of developing a psychiatric mood disorder in adolescence that operate through mediators related to neighborhood poverty, the school environment, and instability of the neighborhood and school environments, considered together and separately. Interventional direct and indirect effects (IDE/IIE) accommodate post-exposure variables that confound the mediator-outcome relationship, but currently, no readily implementable nonparametric estimator for IDE/IIE exists that allows for both multivariate mediators and multivariate post-exposure intermediate confounders. The absence of such an IDE/IIE estimator that can easily accommodate both multivariate mediators and post-exposure confounders represents a significant limitation for real-world analyses, because when considering each mediator subgroup separately, the remaining mediator subgroups (or a subset of them) become post-exposure intermediate confounders. We address this gap by extending a recently developed nonparametric estimator for the IDE/IIE to allow for easy incorporation of multivariate mediators and multivariate post-exposure confounders simultaneously. We apply the proposed estimation approach to our analysis, including walking through a strategy to account for other, possibly co-occurring intermediate variables when considering each mediator subgroup separately.

This tutorial discusses a methodology for causal inference using longitudinal modified treatment policies. This method facilitates the mathematical formalization, identification, and estimation of many novel parameters and mathematically generalizes many commonly used parameters, such as the average treatment effect. Longitudinal modified treatment policies apply to a wide variety of exposures, including binary, multivariate, and continuous, and can accommodate time-varying treatments and confounders, competing risks, loss to follow-up, as well as survival, binary, or continuous outcomes. Longitudinal modified treatment policies can be seen as an extension of static and dynamic interventions to involve the natural value of treatment and, like dynamic interventions, can be used to define alternative estimands with a positivity assumption that is more likely to be satisfied than estimands corresponding to static interventions. This tutorial aims to illustrate several practical uses of the longitudinal modified treatment policy methodology, including describing different estimation strategies and their corresponding advantages and disadvantages. We provide numerous examples of types of research questions that can be answered using longitudinal modified treatment policies. We go into more depth with one of these examples, specifically, estimating the effect of delaying intubation on critically ill COVID-19 patients' mortality. We demonstrate the use of the open-source R package lmtp to estimate the effects, and we provide code on https://github.com/kathoffman/lmtp-tutorial.

Time-varying confounding is a common challenge for causal inference in observational studies with time-varying treatments, long follow-up periods, and participant dropout. Confounder adjustment using traditional approaches can be limited by data sparsity, weight instability, and computational issues. The Nicotine Dependence in Teens Study is a prospective cohort study, and we used data from 21 data collection cycles carried out from 1999 to 2008 among 1294 students recruited from 10 high schools in Montreal, Quebec, Canada, including follow-up into adulthood. Our aim in this study was to estimate associations of timing of alcohol initiation and cumulative duration of alcohol use with depression symptoms in adulthood. Based on the target trials framework, we defined intention-to-treat and as-treated parameters in a marginal structural model with sex as a potential effect-modifier. We then used the observational data to emulate the trials. For estimation, we used pooled longitudinal target maximum likelihood estimation, a plug-in estimator with double-robust and local efficiency properties. We describe strategies for dealing with high-dimensional potential drinking patterns and practical positivity violations due to a long follow-up time, including modifying the effect of interest by removing sparsely observed drinking patterns from the loss function and applying longitudinal modified treatment policies to represent the effect of discouraging drinking.

Mediation analysis is a strategy for understanding the mechanisms by which interventions affect later outcomes. However, unobserved confounding concerns may be compounded in mediation analyses, as there may be unobserved exposure-outcome, exposure-mediator, and mediator-outcome confounders. Instrumental variables (IVs) are a popular identification strategy in the presence of unobserved confounding. However, in contrast to the rich literature on the use of IV methods to identify and estimate a total effect of a non-randomized exposure, there has been almost no research into using IV as an identification strategy to identify mediational indirect effects. In response, we define and nonparametrically identify novel estimands-double complier interventional direct and indirect effects-when 2, possibly related, IVs are available, one for the exposure and another for the mediator. We propose nonparametric, robust, efficient estimators for these effects and apply them to a housing voucher experiment.

Longitudinal modified treatment policies (LMTP) have been recently developed as a novel method to define and estimate causal parameters that depend on the natural value of treatment. LMTPs represent an important advancement in causal inference for longitudinal studies as they allow the non-parametric definition and estimation of the joint effect of multiple categorical, ordinal, or continuous treatments measured at several time points. We extend the LMTP methodology to problems in which the outcome is a time-to-event variable subject to a competing event that precludes observation of the event of interest. We present identification results and non-parametric locally efficient estimators that use flexible data-adaptive regression techniques to alleviate model misspecification bias, while retaining important asymptotic properties such as [Formula: see text]-consistency. We present an application to the estimation of the effect of the time-to-intubation on acute kidney injury amongst COVID-19 hospitalized patients, where death by other causes is taken to be the competing event.

We describe semiparametric estimation and inference for causal effects using observational data from a single social network. Our asymptotic results are the first to allow for dependence of each observation on a growing number of other units as sample size increases. In addition, while previous methods have implicitly permitted only one of two possible sources of dependence among social network observations, we allow for both dependence due to transmission of information across network ties and for dependence due to latent similarities among nodes sharing ties. We propose new causal effects that are specifically of interest in social network settings, such as interventions on network ties and network structure. We use our methods to reanalyze an influential and controversial study that estimated causal peer effects of obesity using social network data from the Framingham Heart Study; after accounting for network structure we find no evidence for causal peer effects.

Natural direct and indirect effects are mediational estimands that decompose the average treatment effect and describe how outcomes would be affected by contrasting levels of a treatment through changes induced in mediator values (in the case of the indirect effect) or not through induced changes in the mediator values (in the case of the direct effect). Natural direct and indirect effects are not generally point-identified in the presence of a treatment-induced confounder; however, they may be identified if one is willing to assume monotonicity between the treatment and the treatment-induced confounder. We argue that this assumption may be reasonable in the relatively common encouragement-design trial setting, where the intervention is randomized treatment assignment and the treatment-induced confounder is whether or not treatment was actually taken/adhered to. We develop efficiency theory for the natural direct and indirect effects under this monotonicity assumption, and use it to propose a nonparametric, multiply robust estimator. We demonstrate the finite sample properties of this estimator using a simulation study, and apply it to data from the Moving to Opportunity Study to estimate the natural direct and indirect effects of being randomly assigned to receive a Section 8 housing voucher-the most common form of federal housing assistance-on risk developing any mood or externalizing disorder among adolescent boys, possibly operating through various school and community characteristics.