Skip to main content

Use of power-law analysis to predict abuse or diversion of prescribed medications: proof-of-concept mathematical exploration



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% [68]. A Lorenz-1 of 15% or more indicates an abusable medication [8].

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 [9], with both > 3 times the labeled/recommended dosage for any indication [10]. 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 [11]. 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 [12], which was originally developed to predict catastrophic hurricanes, earthquakes, and other geophysical hazards [13]. 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 [12]. For example, a transition in the power law for hurricane wind speed magnitudes, identified by an early pioneer of the technique, indicated eyewall formation [12].

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].

Main text


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 [17].

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 [18], 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.

Table 1 Patient characteristics and utilization patterns by medication and group, baseline treatment year

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).

Fig. 1

Criterion validity analyses: standardized mean dosage/day, baseline treatment and 6-month follow-up, by PLGs. PLG power-law group

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.

Table 2 Criterion validity assessment: Lorenz-1 status and inpatient hospital use

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 [19], 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 [12]. The increasing availability of data from prescription drug monitoring programs (PDMP) [20] 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 [21].


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 [22]. 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 [23], 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.



United States


fraud, waste, or abuse


Lorenz-curve analysis




morphine-milligram equivalents


inpatient hospitalization


power-law analysis


power-law zones


power-law groups


positive predictive value


substance use disorder


prescription drug monitoring program


cumulative frequency




current procedural terminology


Healthcare Common Procedure Coding System


International Classification of Diseases


  1. 1.

    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.

    Article  PubMed  Google Scholar 

  2. 2.

    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.

    Article  PubMed  Google Scholar 

  3. 3.

    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.

    Article  PubMed  Google Scholar 

  4. 4.

    Starner CI, Qiu Y, Karaca-Mandic P. Gleason PP. J Manag Care Spec Pharm. 2016;22(12):1403–10.

    Article  PubMed  Google Scholar 

  5. 5.

    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.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    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.

    Google Scholar 

  7. 7.

    Hallas J. Drug utilization statistics for individual-level pharmacy dispensing data. Pharmacoepidemiol Drug Saf. 2005;14:455–61.

    Article  PubMed  Google Scholar 

  8. 8.

    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.

    Article  PubMed  CAS  Google Scholar 

  9. 9.

    US Centers for Disease Control and Prevention. Calculating total daily dose of opioids for safer dosage. Accessed 05 Jun 2018.

  10. 10.

    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.

    Article  CAS  Google Scholar 

  11. 11.

    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.

    Article  PubMed  CAS  Google Scholar 

  12. 12.

    United States Geological Survey. Natural disasters-forecasting economic and life losses. USGS fact sheet. Accessed 05 Jun 2018.

  13. 13.

    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.

    PubMed  Google Scholar 

  14. 14.

    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.

    PubMed  Article  Google Scholar 

  15. 15.

    Substance Abuse and Mental Health Services Administration. Emergency department visits involving nonmedical use of the anti-anxiety medication alprazolam. 2014. Accessed 05 Jun 2018.

  16. 16.

    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.

    Article  Google Scholar 

  17. 17.

    Truven Health Analytics. MarketScan bibliography. Accessed 05 Jun 2018.

  18. 18.

    Vieth E. Fitting piecewise linear regression functions to biological responses. J Appl Physiol. 1989;67(1):390–6.

    Article  PubMed  CAS  Google Scholar 

  19. 19.

    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.

    Article  PubMed  Google Scholar 

  20. 20.

    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.

    Article  PubMed  Google Scholar 

  21. 21.

    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.

    Article  Google Scholar 

  22. 22.

    Dowell D, Haegerich T, Chou R. CDC guideline for prescribing opioids for chronic pain—United States, 2016. JAMA. 2016;315(15):1624–45.

    Article  PubMed  CAS  Google Scholar 

  23. 23.

    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.

    Article  PubMed  CAS  Google Scholar 

Download references

Authors’ contributions

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.


Not applicable.

Competing interests

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

Not applicable.

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.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information



Corresponding author

Correspondence to Kathleen A. Fairman.

Additional file

Additional file 1: Appendix S1. 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.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

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).

Download citation


  • Gabapentin
  • Alprazolam
  • Opioids
  • Abuse
  • Diversion
  • Power-law analysis
  • Lorenz-curve analysis