Use of power-law analysis to predict abuse or diversion of prescribed medications: proof-of-concept mathematical exploration
BMC Research Notes volume 11, Article number: 523 (2018)
To conduct a proof-of-concept study comparing Lorenz-curve analysis (LCA) with power-law (exponential function) analysis (PLA), by applying segmented regression modeling to 1-year prescription claims data for three medications—alprazolam, opioids, and gabapentin—to predict abuse and/or diversion using power-law zone (PLZ) classification.
In 1-year baseline observation, patients classified into the top PLZ groups (PLGs) were demographically and diagnostically similar to those in Lorenz-1 (top 1% of utilizers) and Lorenz-25 (top 25%). For prediction of follow-up (6-month post-baseline) Lorenz-1 use of alprazolam and opioids (i.e., potential abuse/diversion), PLA had somewhat lower sensitivity compared with LCA (83.5–95.4% vs. 99.5–99.9%, respectively) but better specificity (98.2–98.8% vs. 75.5%) and much better positive predictive value (PPV; 34.5–45.3% vs. 4.0–4.6%). Of top-PLG alprazolam- and opioid-treated patients, respectively, 20.7 and 9.9% developed incident (new) Lorenz-1 in followup, compared with < 3% of Lorenz-25 patients. For gabapentin, neither PLA nor LCA predicted incident Lorenz-1 (PPV = 0.0–1.4%). For all three medications, PLA sensitivity for follow-up hospitalization was < 5%, but specificity was better for PLA (97.3–99.2%) than for LCA (74.3–75.4%). PLA better identified patients at risk of future controlled substance abuse/diversion than did LCA, but the technique needs refinement before widespread use.
For the United States (US) health care system, identifying and intervening on patients who are abusing or diverting controlled substances is a top priority because of high cost, morbidity, and mortality associated with fraud, waste, or abuse (FWA) [1, 2]. Numerous methods for identifying at-risk patients have been proposed [3,4,5]. Limiting their utility, most methods require integrated medical/pharmacy claims datasets, laboratory data, or complicated algorithms that may exceed the programming resources available in many health care organizations [3,4,5].
Various simpler analyses based solely on pharmacy claims data have been suggested [3, 6]. One common technique is Lorenz-curve analysis (LCA), which assesses the percentage of total medication supply dispensed to top utilizers, with Lorenz-1 indicating supply dispensed to the top 1% [6–8]. A Lorenz-1 of 15% or more indicates an abusable medication .
Although the Lorenz-1 metric is intuitively appealing and easily calculated, it indicates FWA that is already ongoing, perhaps at a dangerous level. One study of Lorenz-1 utilizers found dispensed dosages averaging 11,274 mg/day for gabapentin and 180 morphine-milligram equivalents (MME)/day for opioids , with both > 3 times the labeled/recommended dosage for any indication . A follow-up study found that concomitant high-dosage consumption/diversion of gabapentin and opioids approximately doubled the risk of inpatient hospitalization (IPH) and quadrupled the risk of respiratory depression-related emergency or IPH care . Thus, although important, Lorenz-1 may indicate damage that has already taken place. Ideally, health care systems would be able to predict medication FWA before it begins, acting on warning signs of future events.
One technique proposed to make predictions of this type, adapted from geophysical sciences, is power-law analysis (PLA), also known as fractal-scaling analysis , which was originally developed to predict catastrophic hurricanes, earthquakes, and other geophysical hazards . In PLA, logarithmically transformed 2-dimensional plots of cumulative frequency against event magnitude are segmented into zones, each characterized by a linear equation (power law). “Transitions” (i.e., points at which the power law changes) indicate fundamental transformations in the nature of the events . For example, a transition in the power law for hurricane wind speed magnitudes, identified by an early pioneer of the technique, indicated eyewall formation .
The present study was a proof-of-concept comparison of LCA with PLA, applied to prescription claims data for 3 medications: alprazolam, the benzodiazepine with the greatest abuse liability; [14, 15] opioids, which are a top public health concern because of an epidemic of misuse and overdose deaths; [1, 16] and gabapentin, a noncontrolled substance recently identified as a medication of abuse [10, 11].
Data source and sample
Study data were derived from the Truven Health Market Scan® Commercial Claims and Encounters Database for calendar years 2013–2015. The database, which comprises pharmacy claims, medical claims, and eligibility (enrollment) data for up to 50 million commercially insured enrollees annually, is commonly used in studies of US health care .
The study sample included enrollees aged 16–64 years with ≥ 2 claims for alprazolam, gabapentin, and/or any opioid during a 365-day period that began with the first observed pharmacy claim (baseline treatment year). To ensure accurate dosage calculations, patients using patches or fentanyl, or with missing or invalid dosages on any claims, were excluded. Patients using multiple study medications were not excluded; instead, each medication had its own baseline year, and concomitant uses of other study medications were measured.
Cohort group definitions
Based on utilization in the baseline treatment year, patients were categorized into groups using both PLA and LCA, including Lorenz-1 (top 1% of users) and Lorenz-25 (top 25%). These classifications were performed using SPSS v24.0 (IBM SPSS, Armonk, NY).
Specifically, for each patient and medication, total supply dispensed during the baseline year was summed as total mg for alprazolam and gabapentin and total MME for opioids, across all claims and strengths. For the LCA, again for each medication, patients were ranked according to total dispensed supply; and Lorenz-1 and quartile groups, including Lorenz-25, were identified.
The PLA also used rankings based on summed supply data, but patients were grouped differently: into 5-‰ bands from the 5th to 95th percentile, and into more granular (smaller) bands at the lowest and highest ends of the percentile distribution. Specifically, at the lowest end, patients were grouped as < 1st percentile and 1st < 5th percentile. At the highest end, patients were grouped into single-percentage increments (e.g., 96th, 97th, etc.) to 99th, then into bands representing one-tenth of 1 percent increments (e.g., 99.1, 99.2%, etc.) up to the 99.9th percentile. This approach was taken because PLA predicts rare events, making a granular analysis of top utilizers necessary to determine the power law for each linear segment.
Summed data were then loaded into Excel (Microsoft, Redmond, WA). Power-law calculations, including identification of transition points and power-law zones (PLZs), were performed in Excel using a method similar to that of geophysical sciences for event magnitude and frequency [12, 13]. Specifically, for each utilization band, ranked from highest to lowest magnitude (i.e., dosage), medication supply/patient was calculated (summed supply ÷ summed patient count). Then, again for each band, patient frequencies (counts) were accumulated, whereas dosage magnitudes (supply/patient) were not accumulated. Both cumulative frequencies and magnitudes for each band were then log-transformed (log10; Additional file 1: Appendix S1). The resulting values were plotted, with the x-axis representing log10-magnitude, and the y-axis representing log10-cumulative frequency. Using a segmented-regression approach , PLZs were identified by visual inspection coupled with model fitting using linear regression.
Cohort group analyses
Two types of analyses were performed in SPSS, each including comparisons of groups based on LCA quartiles and on PLZs. The first assessed patient characteristics in the baseline treatment year. The second, a criterion validity assessment, was limited to subgroups of patients continuously enrolled through the 6 months after the baseline year (follow-up). Follow-up outcomes included Lorenz-1 utilization (i.e., potential FWA), IPH, and dosages/day standardized as Z-scores (distance from the mean measured in standard deviation units). Top-PLZ groups (PLGs) and top-LCA quartiles at baseline were defined as “at risk”; and rates of sensitivity, specificity, positive predictive value (PPV), and negative predictive value were calculated based on follow-up outcomes.
For alprazolam, 4 separate power-law zones (linear segments) were identified, with R2 ranging from 0.990 to 0.999. For gabapentin and opioids, 3 zones were identified, with R2 = 0.961–0.998 (Additional file 1: Appendix S2, Additional file 1: Appendix S3).
Patient characteristics, baseline treatment year
For alprazolam and opioids, top PLGs were larger and used less medication compared with Lorenz-1 groups (Table 1; Additional file 1: Appendix S4). Threshold dosages for top-PLG and Lorenz-1, respectively, were 5.04 versus 7.40 mg/day alprazolam (i.e., 126 and 185% of maximum labeled/recommended dosage); and 130.2 versus 271.2 MME/day opioids (i.e., 260% vs. 542% of maximum labeled/recommended dosage). Of those in the top PLG, only 46.2% of alprazolam- and 36.9% of opioid-treated patients were Lorenz-1 in the baseline year. In all other respects, however, top PLG and Lorenz-1 patients were similar, demographically, diagnostically, and in proportions of claims exceeding recommended dosages. For example, comparing alprazolam-treated patients in PLG-4 and Lorenz-1, respectively, 63 and 64% were female; mean claims/month exceeding labeled/recommended dosage were 0.84 and 0.91; 21 and 20% were diagnosed with substance use disorder (SUD); and 59% in each group were diagnosed with anxiety. Similar patterns were observed in patients treated with opioids.
In contrast, for those treated with gabapentin, the threshold dosage for the top PLG was higher than the Lorenz-1 threshold, at 12,510 and 10,356 mg/day, respectively; and 100% of PLG-3 patients were Lorenz-1 (Table 1). However, gabapentin PLG-3 and Lorenz-1 patients were similar in other respects.
For all 3 medications, Lorenz-25 patient groups were characterized by relatively low dosage thresholds (ranging from 12 to 29% of maximum; Table 1). Baseline SUD prevalence rates were generally lower for these groups than for the top PLGs and Lorenz-1 groups, but otherwise, the groups were demographically and diagnostically similar.
Criterion validity analyses
From the 1-year baseline to 6-month follow-up periods for all 3 medications, mean daily dosages changed only modestly for most treated patients in PLG-1 and PLG-2 (Fig. 1). In contrast, for those in each top PLG, a distinctive splitting pattern occurred, in which a proportion of patients experienced substantial increases in dispensed dosage from baseline to follow-up. The LCA was dominated by Lorenz-25, with no visible distinctions among quartiles to characterize the threshold dosage at which the splitting pattern began (Additional file 1: Appendix S5).
For prediction of follow-up Lorenz-1 use of alprazolam and opioids, PLA had somewhat lower sensitivity compared with LCA (83.5–95.4% vs. 99.5–99.9%, respectively) but better specificity (98.2–98.8% vs. 75.5%) and much better PPV (34.5–45.3% vs. 4.0–4.6%; Table 2). Of those in the top PLG for alprazolam and opioids, respectively, 20.7 and 9.9% went from utilization < Lorenz-1 at baseline to Lorenz-1 in follow-up (i.e., incident Lorenz-1). For the same respective medications, only 2.5 and 1.3% in Lorenz-25 went on to incident Lorenz-1. For patients treated with gabapentin, PPV for incident Lorenz-1 was < 2% for both techniques.
For follow-up IPH, sensitivity was much less for PLA (1.0–4.6%) than for LCA (30.5–38.0%). However, both specificity and PPV were improved using PLA.
This proof-of-concept study applied PLA, which was originally developed to predict catastrophic events in physical systems, for the new purpose of predicting FWA development based solely on pharmacy claims data. For controlled substances, PLA performed better than did LCA at identifying a cohort of patients not currently engaging in, but at risk of developing, future FWA, with somewhat reduced sensitivity but better specificity and PPV. Notably, patients in the top PLGs, Lorenz-25, and Lorenz-1 were demographically and clinically similar at baseline in all respects other than pharmacy utilization. For a health care organization attempting to target FWA mitigation efforts to enrollees most in need of them, the improved PPV and specificity outcomes achieved with PLA are potentially important. However, in this preliminary analytic stage, PLA’s case-finding ability, especially for IPH, was not sufficient for widespread application. Several areas for future development are indicated.
Foremost, results suggest that patients in top PLGs may be at increased risk of future abuse/diversion; yet, top-PLG status is clearly not a perfect predictor of FWA. Future analyses should assess FWA predictors in subsamples limited to top-PLG patients. Bollinger-band analysis, which has been used to predict surges in use of intensive-care services , might provide additional insights into medication demand changes in top PLGs. Perhaps the lower PPVs for IPH compared with Lorenz-1 indicate detection of diversion, rather than consumption, of medication.
Additionally, the most common use of PLA in the physical sciences is performed at the regional level, predicting federal expenditures on catastrophic events for specific geographic areas . The increasing availability of data from prescription drug monitoring programs (PDMP)  may represent an opportunity to apply PLA to regional utilization and mortality data, thereby providing health care payers with better information about which regions are most at risk, and which PDMP features may most effectively mitigate FWA-related harms.
Application of PLA to pharmacy claims data is a promising new method for identifying patients at risk of FWA but needs additional refinement prior to widespread use. Potential future applications may include proactive management of known or emerging medications of abuse .
If PLA predicts fundamental system change, it is not clear why only a portion of patients in each top PLG went on to experience large dosage increases, indicating potential FWA. Our application of PLA to all treated patients with ≥ 2 claims may have been overly broad—perhaps analogous to measuring a catastrophic event of any type, rather than the technique’s original purpose of measuring only one type of catastrophic event. Future applications might limit samples to more restricted subgroups and outcomes specific to key clinical situations; for example, development of FWA in patients with SUD, or of SUD in patients treated with opioids for chronic noncancer pain . Similarly, since recent work has suggested that gabapentin misuse is likely only in patients who either have a SUD diagnosis or are concomitantly using opioids , a future sample might be limited to patients with concomitant gabapentin/opioid use to improve PLA’s PPV for gabapentin misuse. Moreover, the 6-month follow-up period, the maximum possible for this dataset and design, was not long enough to measure long-term developments. Additional analyses following patients for longer time periods may yield more information about FWA development.
fraud, waste, or abuse
positive predictive value
substance use disorder
prescription drug monitoring program
current procedural terminology
Healthcare Common Procedure Coding System
International Classification of Diseases
Rudd RA, Seth P, David F, Scholl L. Increases in drug and opioid-involved overdose deaths—United States, 2010–2015. MMWR Morb Mortal Wkly Rep. 2016;65(5051):1445–52.
Oderda GM, Lake J, Rudell K, Roland CL, Masters ET. Economic burden of prescription opioid misuse and abuse: a systematic review. J Pain Palliat Care Pharmacother. 2015;29(4):388–400.
Canan C, Polinski JM, Alexander GC, Kowal MK, Brennan TA, Shrank WH. Automatable algorithms to identify nonmedical opioid use using electronic data: a systematic review. J Am Med Inform Assoc. 2017;24(6):1204–10.
Starner CI, Qiu Y, Karaca-Mandic P. Gleason PP. J Manag Care Spec Pharm. 2016;22(12):1403–10.
Cochran G, Woo B, Lo-Ciganic WH, Gordon AJ, Donohue JM, Gellad WF. Defining non-medical use of prescription opioids within health care claims: a systematic review. Subst Abus. 2015;36(2):192–202.
Birt J, Johnston J, Nelson D. Exploration of claims-based utilization measures for detecting potential nonmedical use of prescription drugs. J Manag Care Pharm. 2014;20(6):639–46.
Hallas J. Drug utilization statistics for individual-level pharmacy dispensing data. Pharmacoepidemiol Drug Saf. 2005;14:455–61.
Gjerden P, Bramness JG, Slordal L. The use and potential abuse of anticholinergic antiparkinson drugs in Norway: a pharmacoepidemiological study. Br J Clin Pharmacol. 2009;67(2):228–33.
US Centers for Disease Control and Prevention. Calculating total daily dose of opioids for safer dosage. https://www.cdc.gov/drugoverdose/pdf/calculating_total_daily_dose-a.pdf. Accessed 05 Jun 2018.
Peckham AM, Fairman KA, Sclar DA. Prevalence of gabapentin abuse: comparison with agents with known abuse potential in a commercially insured US population. Clin Drug Invest. 2017;37(8):763–73.
Peckham AM, Fairman KA, Sclar DA. All-cause and drug-related medical events associated with overuse of gabapentin and/or opioid medications: a retrospective cohort analysis of a commercially insured US population. Drug Saf. 2018;41(2):213–28.
United States Geological Survey. Natural disasters-forecasting economic and life losses. USGS fact sheet. https://pubs.usgs.gov/fs/natural-disasters/. Accessed 05 Jun 2018.
Fairman KA, Rucker ML. Fractal mathematics in managed care? How a simple and revealing analysis could improve the forecasting and management of medical costs and events. J Manag Care Pharm. 2009;15(4):351–8.
Shukla L, Bokka S, Shukla T, et al. Benzodiazepine and “Z-drug” dependence: data from a tertiary care center. Prim Care Companion CNS Disord. 2017. https://doi.org/10.4088/PCC.16br02025.
Substance Abuse and Mental Health Services Administration. Emergency department visits involving nonmedical use of the anti-anxiety medication alprazolam. 2014. https://archive.samhsa.gov/data/2k14/DAWN153/sr153-alprazolam-2014.pdf. Accessed 05 Jun 2018.
Kanouse AB, Compton P. The epidemic of prescription opioid abuse, the subsequent rising prevalence of heroin use, and the federal response. J Palliat Care Pharmacother. 2015;29(2):102–14.
Truven Health Analytics. MarketScan bibliography. http://truvenhealth.com/markets/life-sciences/products/data-tools/marketscan-bibliography. Accessed 05 Jun 2018.
Vieth E. Fitting piecewise linear regression functions to biological responses. J Appl Physiol. 1989;67(1):390–6.
Pagel C, Ramnarayan P, Ray S, Peters MJ. A novel method to identify the start and end of the winter surge in demand for pediatric intensive care in real time. Pediatr Crit Care Med. 2015;16(9):821–7.
Moyo P, Simoni-Wastila L, Griffin BA, Onukwugha E, Harrington D, Alexander GC, Palumbo F. Impact of prescription drug monitoring programs (PDMPs) on opioid utilization among Medicare beneficiaries in 10 US states. Addiction. 2017;112(10):1784–96.
Peckham AM, Fairman KA, Sclar DA. Policies to mitigate nonmedical use of prescription medications: how should emerging evidence of gabapentin misuse be addressed? Exp Opin Drug Saf. 2018;17(5):519–23.
Dowell D, Haegerich T, Chou R. CDC guideline for prescribing opioids for chronic pain—United States, 2016. JAMA. 2016;315(15):1624–45.
Peckham AM, Evoy KE, Covvey JR, Ochs L, Fairman KA, Sclar DA. Predictors of gabapentin overuse with or without concomitant opioids in a commercially insured US population. Pharmacotherapy. 2018;38(44):436–43.
KAF and MLR performed concept and design, assisted by AMP, JHR, and DAS. KAF, MLR, and JHR performed data analyses and constructed the data presentations. MLR was the primary methodological expert, and AMP was the primary clinical expert, in data interpretation. KAF drafted the manuscript, assisted by AMP. All authors reviewed the manuscript and contributed to intellectual content. All authors read and approved the final manuscript.
AMP, DAS, KAF and JHR declare that they have no competing interests. MLR is a consulting engineer whose services include the development and application of the techniques discussed in this article. KAF and JHR are MLR’s spouse and son, respectively.
Availability of data and materials
The data that support the findings of this study are available from Truven Health Market Scan® Commercial Claims and Encounters Database but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Aggregated data that were used to define the power-law zones are available from the authors upon reasonable request.
Consent to publish
Ethics approval and consent to participate
The Institutional Review Board (IRB) Committee of Midwestern University (Glendale, Arizona) found that this study fulfilled criteria for IRB Exemption per Code of Federal Regulations, Title 45, Part 46, section 101 (45 CFR 26.101).
This study was not funded.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Power-law file preparation, alprazolam, baseline treatment year. Example calculation to illustrate how aggregated data are prepared for power-law analysis. CF cumulative frequency, mg milligrams. 540,752 = total n of alprazolam-treated patients meeting sample criteria. Appendix S2. Power-Law Curves. Excel graphics showing (a) logarithmically transformed 2-dimensional plots of cumulative frequency against event magnitude and (b) transition points for each power-law zone. Appendix S3. Linear Equations by Power-Law Zones and Medication. Slopes and R2 for each power-law zone. Appendix S4. Diagnoses. International classification of diseases, diagnosis-related group, place of service, and current procedural terminology codes for all diagnoses measured in the study. a Measured using both ICD-9 and ICD-10 coding because follow-up period included dates of service on and after October 1, 2015. CPT current procedural terminology, HCPCS healthcare common procedure coding system, ICD international classification of diseases. Appendix S5. Criterion validity analyses: standardized mean dosage/day, baseline treatment and 6-month follow-up, by utilization quartiles.
About this article
Cite this article
Fairman, K.A., Peckham, A.M., Rucker, M.L. et al. Use of power-law analysis to predict abuse or diversion of prescribed medications: proof-of-concept mathematical exploration. BMC Res Notes 11, 523 (2018). https://doi.org/10.1186/s13104-018-3632-y