Abstract:BACKGROUND: Overuse and misuse of prescription opioids have become significant public health burdens in the United States. About 11.5 million people are estimated to have misused prescription opioids for nonmedical purposes in 2016. This has led to a significant number of drug overdose deaths in the United States. Previous studies have examined spatiotemporal clusters of opioid misuse, but they have been restricted to circular shaped regions.OBJECTIVE: The goal of this study was to identify spatiotemporal hot spots of opioid users and opioid prescription claims using Medicare data.METHODS: We examined spatiotemporal clusters with significantly higher number of beneficiaries and rate of prescriptions for opioids using Medicare payment data from the Centers for Medicare & Medicaid Services. We used network scan statistics to detect significant clusters with arbitrary shapes, the Kulldorff scan statistic to examine the significant clusters for each year (2013, 2014, and 2015) and an expectation-based version to examine the significant clusters relative to past years. Regression analysis was used to characterize the demographics of the counties that are a part of any significant cluster, and data mining techniques were used to discover the specialties of the anomalous providers.RESULTS: We examined anomalous spatial clusters with respect to opioid prescription claims and beneficiary counts and found some common patterns across states: the counties in the most anomalous clusters were fairly stable in 2014 and 2015, but they have shrunk from 2013. In Virginia, a higher percentage of African Americans in a county lower the odds of the county being anomalous in terms of opioid beneficiary counts to about 0.96 in 2015. For opioid prescription claim counts, the odds were 0.92. This pattern was consistent across the 3 states and across the 3 years. A higher number of people in the county with access to Medicaid increased the odds of the county being in the anomalous cluster to 1.16 in both types of counts in Virginia. A higher number of people with access to direct purchase of insurance plans decreased the odds of a county being in an anomalous cluster to 0.85. The expectation-based scan statistic, which captures change over time, revealed different clusters than the Kulldorff statistic. Providers with an unusually high number of opioid beneficiaries and opioid claims include specialties such as physician's assistant, nurse practitioner, and family practice.CONCLUSIONS: Our analysis of the Medicare claims data provides characteristics of the counties and provider specialties that have higher odds of being anomalous. The empirical analysis identifies highly refined spatial hot spots that are likely to encounter prescription opioid misuse and overdose. The methodology is generic and can be applied to monitor providers and their prescription behaviors in regions that are at a high risk of abuse.

Estimation of opioid misuse prevalence in New York State counties, 2007-2018. A Bayesian spatio-temporal abundance model approach

An integrated abundance model for estimating county-level prevalence of opioid misuse in Ohio

Invited commentary: Motivating better methods—and better data collection—for measuring prevalence of drug misuse

Small Area Forecasting of Opioid-Related Mortality: Bayesian Spatiotemporal Dynamic Modeling Approach

A latent spatial factor approach for synthesizing opioid associated deaths and treatment admissions in Ohio counties

Large-scale Analysis of Opioid Poisoning Related Hospital Visits in New York State

Massachusetts Prevalence of Opioid Use Disorder Estimation Revisited: Comparing a Bayesian Approach to Standard Capture-Recapture Methods

Discovering heterogeneous subpopulations for fine-grained analysis of opioid use and opioid use disorders

Detection of Spatiotemporal Prescription Opioid Hot Spots With Network Scan Statistics: Multistate Analysis

Simulating the simultaneous impact of medication for opioid use disorder and naloxone on opioid overdose death in eight New York counties

Predicting U.S. county opioid poisoning mortality from multi-modal social media and psychological self-report data

Spatiotemporal Forecasting of Opioid-Related Fatal Overdoses: Towards Best Practices for Modeling and Evaluation

Investigating opioid vulnerability profiles through a spatial Kalman filter: insights from social and environmental factors

Intersectional inequities and longitudinal prevalence estimates of opioid use disorder in Massachusetts 2014-2020: a multi-sample capture-recapture analysis

Investigating the importance of social vulnerability in opioid-related mortality across the United States

Firsthand Opiates Abuse on Social Media: Monitoring Geospatial Patterns of Interest Through a Digital Cohort

Spatial heterogeneity and socioeconomic determinants of opioid prescribing in England between 2015 and 2018

Network-Based Discovery of Opioid Use Vulnerability in Rats Using the Bayesian Stochastic Block Model

Agent-based Dynamics of a SPAHR Opioid Model on Social Network Structures

Risk factors for persistent fatal opioid-involved overdose clusters in Massachusetts 2011-2021: a spatial statistical analysis with socio-economic, accessibility, and prescription factors