Predicting the Future Course of Opioid Overdose Mortality: An Example From Two US States

The rapid growth of opioid abuse and the related mortality across the United States has spurred the development of predictive models for the allocation of public health resources. These models should characterize heterogeneous growth across states using a drug epidemic framework that enables assessments of epidemic onset, rates of growth, and limited capacities for epidemic growth.

We used opioid overdose mortality data for 146 North and South Carolina counties from 2001 through 2014 to compare the retrodictive and predictive performance of a logistic growth model that parameterizes onsets, growth, and carrying capacity within a traditional Bayesian Poisson space-time model.

In fitting the models to past data, the performance of the logistic growth model was superior to the standard Bayesian Poisson space-time model (deviance information criterion: 8,088 vs. 8,256), with reduced spatial and independent errors. Predictively, the logistic model more accurately estimated fatality rates 1, 2, and 3 years in the future (root mean squared error medians were lower for 95.7% of counties from 2012 to 2014). Capacity limits were higher in counties with greater population size, percent population age 45-64, and percent white population. Epidemic onset was associated with greater same-year and past-year incidence of overdose hospitalizations.

Growth in annual rates of opioid fatalities was capacity limited, heterogeneous across counties, and spatially correlated, requiring spatial epidemic models for the accurate and reliable prediction of future outcomes related to opioid abuse. Indicators of risk are identifiable and can be used to predict future mortality outcomes.