A multivariate spatiotemporal change-point model of opioid overdose deaths in Ohio

Ohio is one of the states most impacted by the opioid epidemic and experienced the second highest age-adjusted fatal drug overdose rate in 2017. Initially it was believed prescription opioids were driving the opioid crisis in Ohio. However, as the epidemic evolved, opioid overdose deaths due to fentanyl have drastically increased. In this work we develop a Bayesian multivariate spatiotemporal model for Ohio county overdose death rates from 2007 to 2018 due to different types of opioids. The log-odds are assumed to follow a spatially varying change point regression model. By assuming the regression coefficients are a multivariate conditional autoregressive process, we capture spatial dependence within each drug type and also dependence across drug types. The proposed model allows us to not only study spatiotemporal trends in overdose death rates but also to detect county-level shifts in these trends over time for various types of opioids.