Trends and Relationships in Opioid Prescribing Rates and Overdose Rates in the United States 2013-2019 ====================================================================================================== * Luke F. Enthoven ## Abstract The opioid epidemic began more than two decades ago and continues to increase in severity with the rate of overdose death increasing steadily. The epidemic was born in problematic claims about opioids and aggressive prescribing. After initial analysis of the first decade of the epidemic it became clear that actions would need to be taken to curtail overprescribing and many states passed legislation tracking and limiting opioid prescribing. Since then opioid prescribing has declined, but overdose has continued to increase at an alarming rate. This analysis is focused on evaluating more recent trends in opioid prescribing and overdose to see if the close positive correlation that existed between these values prior to 2010 holds true for the following decade. My findings indicate that between 2013 and 2019 that there is no relationship, or possibly a negative relationship, between opioid prescribing and overdose. Between 2013 and 2019 there was a 25% decrease in opioid prescribing rate yet a 96% increase in opioid overdose deaths. The close relationship that existed between prescribing and overdose appears to have diverged to the point where overdose increases as prescribing decreases. Linear regression on both national and state data revealed a significant negative relationship between overdose and prescribing for national data between 2013 and 2019 (R2 value of 0.84 and a p-value of 0.003) and no relationship between overdose and prescribing for state data between 2014 and 2019 (R2 value of 0.04 and a p-value of < 0.001). This questions the value of further efforts to reduce prescribing in preventing overdose and supports the hypothesis that restriction of prescription opioids can contribute to increased use of dangerous illicit opioids. ## Introduction The modern opioid crisis has its roots in problematic prescribing practices during a period of aggressive introduction of new opioid medications and a push for increased awareness and treatment of patient pain. Opioid prescribing ramped up heavily from 1995 to 2010 with billions of prescriptions dispensed all over the United States during this period [1]. With this increased medical opioid availability came an increase in diversion and illicit use which led to addiction and physical dependence among both medical and non-medical users [2,3]. As opioid drugs became more popular and the market more saturated, opioid overdoses began to rise in the United States. The first decade or so of the epidemic (1995-2009) was unique in that a majority of overdoses resulted from misuse of these new and widely available potent prescription opioids. These drugs were to blame for the almost four fold increase in overdose deaths in this time span, far outpacing any other group of illicitly used substances [4]. It is also relevant that during these years there was little to no change in rates of heroin use, and fentanyl was not yet widely available in illicit markets. Prescription opioids were clearly preferred by opioid users [5,6]. Due to the majority of opioid users using mainly prescription opioids prior to 2010, the availability of these drugs was highly linked to the rate of overdose. This relationship can most clearly be seen in the 2015 publication by Kolodny et. al. [7] which shows the close relationship between opioid sales, opioid overdose deaths, and opioid treatment admissions from 1999 to 2010. The authors suggest more conservative prescribing as one of many ways to reduce overdose mortality [7]. This same idea guided policy in response to the crisis. By the end of 2011, 49 states had passed legislation establishing some sort of prescription monitoring program. In 2015 the U.S Department of Health and Human Services called for “strong steps” to prevent addiction and overdose, including training, prescription monitoring programs, naloxone availability, and expanded use of medication-assisted treatment for addiction [8]. The goal of this investigation is to evaluate the more recent relationship between opioid prescribing and overdose using publicly available data from the National Institute on Drug Abuse (NIDA) [9], the Centers for Medicare and Medicaid Services (CMS) [10], and the Centers for Disease Control Wonder Databases (CDC) [11]. I assess if the correlation that was highlighted in Kolodny et. al. [7] is still what is observed when looking at the relationship between prescribing and overdose. The analysis also serves to evaluate the merit of continued efforts to curtail opioid prescribing in an attempt to reduce opioid mortality that are still popular among the medical community [12,13,14,15,16] ## Results This investigation utilized data that is publicly available through NIDA, CMS, and the CDC. Overdose data was acquired through the NIDA and the CDC while opioid prescribing data was obtained from CMS. While CMS only provides data on prescribing under Medicare Part D, almost 50 million Americans are currently enrolled and therefore the prescribing rates provided by the CMS for Part D will be used to approximate national prescribing rates. Prescribing rates and overdose rates were compared in multiple ways including nationally by year and by state by year. Data for years 2013 to 2019 was used, but not all analyses cover this whole span. CMS data showed the national opioid prescribing rate decreased 25% between 2013 and 2019 and an AMA release reported a 44.4% decrease in prescribing between 2011 and 2020 [17]. During this same period CDC data shows a 96% increase in the rate of overdose deaths. NIDA data showed the overdose rate to increase in every state from 2013-2019 except for Arkansas, Oklahoma, and Wyoming while CMS data showed a decrease in the prescription rate in every state. All three comparisons show either no relationship between opioid prescribing and overdose or a negative relationship between opioid prescribing and overdose. *Figure 1* and *Figure 2* use opioid overdose rate while *Figure 3* uses total overdose rate as state level opioid overdose data is not readily available. Due to such a large proportion of total overdoses involving opioids, total overdose rate is used as an approximation for opioid overdose rate in analysis using state level data in *Figure 3*. The percentage of overdoses involving opioids in the United States increased over time from 57% in 2013 to 70% in 2019 [9]. *Figure 1* shows a comparison between opioid prescribing and opioid overdose over time and highlights the increase in overdose rate and decrease in prescription rate that occurred between 2013 and 2019. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/06/2022.12.04.22283072/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2022/12/06/2022.12.04.22283072/F1) Figure 1. Opioid overdose rate is calculated as age-adjusted overdose deaths per 100,000 total population. Prescription rate is calculated as total Medicare Part D opioid claims divided by total Medicare Part D prescription claims. Data obtained from CMS and NIDA. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/06/2022.12.04.22283072/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2022/12/06/2022.12.04.22283072/F2) Figure 2. Opioid overdose rate is calculated as overdose deaths per 100,000 population. Prescription rate is calculated as total Medicare Part D opioid claims divided by total Medicare Part D prescription claims. Each data point represents the prescribing and overdose values for the United States at a given year between 2013 and 2019 for a total of 7 national snapshots. The red line is the result of a linear regression with an R2 value of 0.84 and a p-value of 0.003. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/12/06/2022.12.04.22283072/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2022/12/06/2022.12.04.22283072/F3) Figure 3. Overdose rate is calculated as overdose deaths per 100,000 population. Prescription rate is calculated as total Medicare Part D opioid claims divided by total Medicare Part D prescription claims. Each data point represents the prescribing and overdose values for a state at a given year between 2014 and 2019 for a total of 300 state snapshots. The red line is the result of a linear regression with an R2 value of 0.04 and a p-value of < 0.001. Overall overdose rate is used not opioid overdose rate specifically in this figure and analysis. In order to test the direct relationship between prescribing and overdose, they were compared directly using state overdose and state prescribing data. The first of these is a reapplication of data from *Figure 1* that compares prescribing rate and overdose rate on the x and y axes respectively with each data point representing national overdose and prescribing in a given year between 2013 and 2019, *Figure 2*. This scatter was then subject to linear regression resulting in an R2 value of 0.84 and a p-value of 0.003. This indicates a significant negative relationship between prescribing rates and overdose rate in these limited annual national snapshots. The trend shown in the regression for national annual data was strong, but limited by the small number of data points. This led to a comparison of prescribing rate and overdose rate with data points representing a single state in a single year. This takes away the temporal dimension of the data and simply investigates if a state’s prescription rate is related to its overdose rate based on state annual snapshots. The resulting scatter is upgraded from the 7 data points available in *Figure 2* to the 300 shown in *Figure 3*. This scatter was then subjected to a linear regression that gave an R2 value of 0.04 and a p-value of < 0.001. This would indicate that there is no relationship between overdose rate and prescribing rate in these annual state snapshots. ## Discussion It is clear from both our data and similar wide and narrow scope analysis that efforts to reduce opioid prescribing in the United States have been extremely successful [18,19]. Prescribing rates are now as low as they were 20 years ago and continuing to trend downwards [20]. The goal of these changes, as previously stated, was to reduce diversion and abuse of prescription opioids which has also been very successful [21,22,23]. This success ultimately led to a decrease in deaths from prescription opioids in the United States, but appears to have had no effect, if not the opposite effect, on overall overdose mortality [23]. This supports other reports and analyses performed on this subject [17,24]. When the results of this analysis for 2013-2019 are compared to those found by Kolodny et. al., [7] it is clear that the close positive relationship that was once observed between opioid prescribing and overdose has diverged with more recent years showing no relationship or even a negative relationship. Despite the intended decrease in prescription misuse and prescription opioid overdose that inspired regulations, overdose mortality appears to be unaffected by efforts to limit prescription opioid availability [23]. Overall since implicating measures to restrict opioid prescriptions the rate of overdose has only increased in the United States with little indication of slowing [9]. These results continue to support the idea that those who abuse prescription opioids will not cease use when their drug of choice is not available, but instead will move to obtaining non pharmaceutical opioids from illicit sources. Such behavior increases risk of overdose and is reflected in the increased overdose risk for those whose opioid doses are reduced while on long term opioid therapy [25,26]. In the face of this evidence it does not seem efficacious to continue aiming efforts to solve the opioid epidemic at reducing prescribing. Prescribing has already been substantially reduced in the last decade and the changes have done little to impact trends in overdose mortality. Strategies such as patient education, harm reduction, accessible treatment, and physician training have been suggested as alternatives and are substantially supported in helping reduce overdose in the era of fentanyl [27,28,29]. ## Limitations This investigation also had its limitations. This was a very broad look at prescribing and overdose that did not dive deeper into understanding the many factors that influence both overdose and opioid prescribing. Data for study was not available for both metrics prior to 2013 and data was only available up until 2019. There has also been a rise in poly substance abuse. Our data for opioid overdose does not specify if multiple drugs were involved in mortality, making it more difficult to determine the exact role of opioids alone in contributing to overdose deaths. ## Methods As previously stated all data used in this analysis is publicly available through NIDA, CMS, and the CDC. Data was obtained as CSV files that were used in R Studio in order to create figures and perform regression analysis [30]. In R Studio CDC overdose data and CMS prescribing data were combined in a data table for use in creating *Figure 1* and *Figure 2* using national data. NIDA overdose data and CMS prescribing data were combined for use in creating *Figure 3* using state data. *Figure 1* used the ggplot function in the ggplot2 package in order to create the plot [31]. *Figure 2* and *Figure 3* used the native R plot function to create the plots and display regression lines. Regression lines were calculated using the native R abline function. Lines were summarized in R to obtain R2 and p-values for the regression lines. ## Data Availability All data produced are available online at [https://nida.nih.gov/research-topics/trends-statistics/overdose-death-rates](https://nida.nih.gov/research-topics/trends-statistics/overdose-death-rates) [https://data.cms.gov/summary-statistics-on-use-and-payments/medicare-medicaid-opioid-prescribing-rates/medicare-part-d-opioid-prescribing-rates-by-geography](https://data.cms.gov/summary-statistics-on-use-and-payments/medicare-medicaid-opioid-prescribing-rates/medicare-part-d-opioid-prescribing-rates-by-geography) [https://www.cdc.gov/nchs/pressroom/sosmap/drug\_poisoning\_mortality/drug\_poisoning.htm](https://www.cdc.gov/nchs/pressroom/sosmap/drug_poisoning_mortality/drug_poisoning.htm) * Received December 4, 2022. * Revision received December 4, 2022. * Accepted December 6, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Lyden J, Binswanger IA. The United States opioid epidemic. Seminars in perinatology. 2019;43(3):123–131. doi:10.1053/j.semperi.2019.01.001 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.semperi.2019.01.001&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 2. 2.Cicero TJ, Surratt H, Inciardi JA, Munoz A. Relationship between therapeutic use and abuse of opioid analgesics in rural, suburban, and urban locations in the United States. Pharmacoepidemiology and drug safety. 2007;16(8):827–840. doi:10.1002/pds.1452 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pds.1452&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17636553&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000248776700001&link_type=ISI) 3. 3.Compton WM, Volkow ND. Abuse of prescription drugs and the risk of addiction. Drug and alcohol dependence. 2006;83(Supplement 1):S4–S7. doi:10.1016/j.drugalcdep.2005.10.020 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.drugalcdep.2005.10.020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16563663&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000237598400002&link_type=ISI) 4. 4.Calcaterra S, Glanz J, Binswanger IA. National trends in pharmaceutical opioid related overdose deaths compared to other substance related overdose deaths: 1999–2009. Drug and alcohol dependence. 2012;131(3):263–270. doi:10.1016/j.drugalcdep.2012.11.018 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.drugalcdep.2012.11.018&link_type=DOI) 5. 5.Lipari RN, Hughes A. Trends in Heroin Use in the United States: 2002 to 2013. U.S. Dept. of Health & Human Services, Substance Abuse & Mental Health Services Administration, Center for Behavioral Health Statistics and Quality; 2015. 6. 6.Tuazon E, Kunins HV, Allen B, Paone D. Examining opioid-involved overdose mortality trends prior to fentanyl: New York City, 2000–2015. Drug and alcohol dependence. 2019;205:107614–107614. doi:10.1016/j.drugalcdep.2019.107614 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.drugalcdep.2019.107614&link_type=DOI) 7. 7.Kolodny A, Courtwright DT, Hwang CS, et al. The Prescription Opioid and Heroin Crisis: A Public Health Approach to an Epidemic of Addiction. Annual review of public health. 2015;36(1):559–574. doi:10.1146/annurev-publhealth-031914-122957 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-publhealth-031914-122957&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25581144&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 8. 8.Webster LR, Grabois M. Current Regulations Related to Opioid Prescribing. PM & R. 2015;7(11):S236–S247. doi:10.1016/j.pmrj.2015.08.011 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.pmrj.2015.08.011&link_type=DOI) 9. 9.National Center for Biotechnology Information (NCBI)[Internet]. Bethesda (MD): National Library of Medicine (US), National Institute on Drug Abuse; [2022] – [cited 2022 September 26]. Available from: [https://nida.nih.gov/research-topics/trends-statistics/overdose-death-rates](https://nida.nih.gov/research-topics/trends-statistics/overdose-death-rates) 10. 10.Centers for Medicare and Medicaid Services (CMS)[Internet]. Medicare Part D Opioid Prescribing Rates-by Geography; [2020] – [cited 2022 September 26]. Available from: [https://data.cms.gov/summary-statistics-on-use-and-payments/medicare-medicaid-opioid-prescribing-rates/medicare-part-d-opioid-prescribing-rates-by-geography](https://data.cms.gov/summary-statistics-on-use-and-payments/medicare-medicaid-opioid-prescribing-rates/medicare-part-d-opioid-prescribing-rates-by-geography) 11. 11.Centers for Disease Control and Prevention (CDC)[Internet]. National Center for Health Statistics (US), Drug Overdose Mortality by State; [2022] – [cited 2022 September 26]. Available from: [https://www.cdc.gov/nchs/pressroom/sosmap/drug\_poisoning\_mortality/drug\_poisoning.htm](https://www.cdc.gov/nchs/pressroom/sosmap/drug_poisoning_mortality/drug_poisoning.htm) 12. 12.Daoust R, Paquet J, Marquis M, et al. Evaluation of Interventions to Reduce Opioid Prescribing for Patients Discharged From the Emergency Department: A Systematic Review and Meta-analysis. JAMA network open. 2022;5(1):e2143425–e2143425. doi:10.1001/jamanetworkopen.2021.43425 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2021.43425&link_type=DOI) 13. 13.Thiels CA, Ubl DS, Habermann EB. Finding the Balance Between Reduced Opioid Prescribing and Patient-Reported Pain Management among General Surgery Patients. Journal of the American College of Surgeons. 2020;231(4):S97-. doi:10.1016/j.jamcollsurg.2020.07.183 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jamcollsurg.2020.07.183&link_type=DOI) 14. 14.Stanley B, Collins L, Norman A, Bonomo A, Bonomo Y. Reducing opioid prescribing on discharge after orthopaedic surgery: does a guideline and education improve prescribing practice 1 year later? ANZ journal of surgery. 2022;92(5):1171–1177. doi:10.1111/ans.17557 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ans.17557&link_type=DOI) 15. 15.Carnes KM, Singh Z, Ata A, Mian BM. Interventions to Reduce Opioid Prescriptions following Urological Surgery: A Systematic Review and Meta-Analysis. The Journal of urology. 2022;207(5):969–981. doi:10.1097/JU.0000000000002447 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/JU.0000000000002447&link_type=DOI) 16. 16.Muench J, Hoopes M, Mayhew M, et al. Reduction of Long-Term Opioid Prescribing for Back Pain in Community Health Centers After a Medicaid Policy Change. Journal of the American Board of Family Medicine. 2022;35(2):352–369. doi:10.3122/jabfm.2022.02.210306 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiamFiZnAiO3M6NToicmVzaWQiO3M6ODoiMzUvMi8zNTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8xMi8wNi8yMDIyLjEyLjA0LjIyMjgzMDcyLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 17. 17.REPORT SHOWS DECREASES IN OPIOID PRESCRIBING, INCREASE IN OVERDOSES. States News Service. 2021. 18. 18.Jairam V, Yang DX, Pasha S, et al. Temporal Trends in Opioid Prescribing Patterns Among Oncologists in the Medicare Population. JNCI: Journal of the National Cancer Institute. 2021;113(3):274–281. doi:10.1093/jnci/djaa110 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jnci/djaa110&link_type=DOI) 19. 19.Enzinger AC, Ghosh K, Keating NL, Cutler DM, Landrum MB, Wright AA. US Trends in Opioid Access Among Patients With Poor Prognosis Cancer Near the End-of-Life. Journal of clinical oncology. 2021;39(26):2948–2958. doi:10.1200/JCO.21.00476 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1200/JCO.21.00476&link_type=DOI) 20. 20.IQVIA Institute for Human Data Science Issues Report, “Prescription Opioid Trends in the United States.” Professional Services Close-Up. Published online 2020. 21. 21.Reid DBC, Shah KN, Ruddell JH, et al. Effect of narcotic prescription limiting legislation on opioid utilization following lumbar spine surgery. The spine journal. 2019;19(4):717–725. doi:10.1016/j.spinee.2018.09.007 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.spinee.2018.09.007&link_type=DOI) 22. 22.Surratt HL, O’Grady C, Kurtz SP, et al. Reductions in prescription opioid diversion following recent legislative interventions in Florida. Pharmacoepidemiology and drug safety. 2014;23(3):314–320. doi:10.1002/pds.3553 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pds.3553&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24677496&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 23. 23.Dart RC, Surratt HL, Cicero TJ, et al. Trends in Opioid Analgesic Abuse and Mortality in the United States. The New England journal of medicine. 2015;372(3):241–248. doi:10.1056/NEJMsa1406143 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMsa1406143&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25587948&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 24. 24.Champagne-Langabeer T, Madu R, Giri S, Stotts AL, Langabeer JR. Opioid prescribing patterns and overdose deaths in Texas. Substance abuse. 2021;42(2):161–167. doi:10.1080/08897077.2019.1675114 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/08897077.2019.1675114&link_type=DOI) 25. 25.Cicero TJ, Ellis MS, Harney J. Shifting Patterns of Prescription Opioid and Heroin Abuse in the United States. The New England journal of medicine. 2015;373(18):1789–1790. doi:10.1056/NEJMc1505541 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc1505541&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26510045&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 26. 26.Glanz JM, Binswanger IA, Shetterly SM, Narwaney KJ, Xu S. Association Between Opioid Dose Variability and Opioid Overdose Among Adults Prescribed Long-term Opioid Therapy. JAMA network open. 2019;2(4):e192613–e192613. doi:10.1001/jamanetworkopen.2019.2613 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2019.2613&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31002325&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 27. 27.Barglow P. Commentary: The opioid overdose epidemic: Evidence-based interventions. The American journal on addictions. 2018;27(8):605–607. doi:10.1111/ajad.12823 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ajad.12823&link_type=DOI) 28. 28.Volkow ND, Jones EB, Einstein EB, Wargo EM. Prevention and Treatment of Opioid Misuse and Addiction: A Review. Archives of general psychiatry. 2019;76(2):208–216. doi:10.1001/jamapsychiatry.2018.3126 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapsychiatry.2018.3126&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30516809&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 29. 29.Barry CL. Fentanyl and the Evolving Opioid Epidemic: What Strategies Should Policy Makers Consider? Psychiatric services (Washington, DC). 2018;69(1):100–103. doi:10.1176/appi.ps.201700235 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ps.201700235&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28967324&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F12%2F06%2F2022.12.04.22283072.atom) 30. 30.RStudio Team. RStudio: Integrated Development Environment for R. Published online 2020. [http://www.rstudio.com/](http://www.rstudio.com/) 31. 31.Wickham H. ggplot2: Elegant Graphics for Data Analysis. Published online 2016. [https://ggplot2.tidyverse.org](https://ggplot2.tidyverse.org)