Computationally Efficient Estimation of Localized Treatment Effects for Multi-Level, Multi-Component Interventions to Address the Opioid Crisis

2026/01/31 by Abdulrahman A. Ahmed, M. Amin Rahimian, Qiushi Chen +1
Mathematics · Computer Science · Physics and Astronomy · #stat.AP #cs.MA #cs.SI #physics.soc-ph

paper · pdf

repository link: https://github.com/abdulrahmanfci/gpr-metamodel/

arxiv created 2026/07/31 · arxiv updated 2026/08/03

Abstract

The opioid epidemic remains a major public health challenge in the United States, requiring a multi-pronged intervention approach to mitigate harms to communities. Given the heterogeneity of the epidemic, it is crucial for policymakers to understand localized treatment effects of different intervention components and utilize limited resources efficiently. While locally calibrated simulation models can project epidemic outcomes for any given intervention policy, collecting simulation results for all intervention combinations to estimate localized treatment effects for each community is impractical because the number of combinations grows exponentially with the number of interventions and the levels at which they are applied. To tackle this, we develop a two-stage metamodel framework with a two-step sequential design for efficient sampling. The metamodel consists of a response function linking health outcomes to each intervention component's treatment effect, and a Gaussian process regression (GPR) to learn spatial and socio-economic structures of the treatment effects based on locally-contextualized covariates. With two-step sequential sampling, we leverage spatial correlations and posterior uncertainty to sequentially sample the most informative counties and treatment conditions. We apply this framework to estimate the treatment effects of buprenorphine dispensing and naloxone distribution on overdose mortality rates using a calibrated agent-based opioid epidemic model in Pennsylvania counties. Our approach achieves less than 5% average relative error using fewer than 2% of the runs required for an exhaustive simulation. Our two-stage framework provides a computationally efficient approach to support policymakers, enabling an efficient evaluation of alternative resource-allocation strategies to mitigate the opioid epidemic in local communities.

Citations

• Optimized Model Selection for Estimating Treatment Effects from Costly Simulations of the US Opioid Epidemic
• Evaluating diversion and treatment policies for opioid use disorder
• Estimating Treatment Effects Using Costly Simulation Samples from a Population-Scale Model of Opioid Use Disorder
• Inferring epidemic dynamics using Gaussian process emulation of agent-based simulations
• Safe Active Learning for Multi-Output Gaussian Processes
• Maximizing acquisition functions for Bayesian optimization
• A General Framework for Vecchia Approximations of Gaussian Processes
• Improved prediction accuracy for disease risk mapping using Gaussian Process stacked generalisation
• A Tutorial on Bayesian Optimization of Expensive Cost Functions, with Application to Active User Modeling and Hierarchical Reinforcement Learning