Article Text


Prescribed opioids in primary care: cross-sectional and longitudinal analyses of influence of patient and practice characteristics
  1. Robbie Foy1,
  2. Ben Leaman2,
  3. Carolyn McCrorie1,
  4. Duncan Petty1,
  5. Allan House1,
  6. Michael Bennett1,
  7. Paul Carder3,
  8. Simon Faulkner4,
  9. Liz Glidewell1,
  10. Robert West1
  1. 1Leeds Institute of Health Sciences, University of Leeds, Leeds, UK
  2. 2Calderdale Metropolitan Borough Council, Halifax, UK
  3. 3Yorkshire & Humber Commissioning Support Unit, Bradford, UK
  4. 4Health and Social Care Information Centre, Leeds, UK
  1. Correspondence to Dr Robbie Foy; r.foy{at}


Objectives To examine trends in opioid prescribing in primary care, identify patient and general practice characteristics associated with long-term and stronger opioid prescribing, and identify associations with changes in opioid prescribing.

Design Trend, cross-sectional and longitudinal analyses of routinely recorded patient data.

Setting 111 primary care practices in Leeds and Bradford, UK.

Participants We observed 471 828 patient-years in which all patients represented had at least 1 opioid prescription between April 2005 and March 2012. A cross-sectional analysis included 99 847 patients prescribed opioids between April 2011 and March 2012. A longitudinal analysis included 49 065 patient-years between April 2008 and March 2012. We excluded patients with cancer or treated for substance misuse.

Main outcome measures Long-term opioid prescribing (4 or more prescriptions within 12 months), stronger opioid prescribing and stepping up to or down from stronger opioids.

Results Opioid prescribing in the adult population almost doubled for weaker opioids over 2005–2012 and rose over sixfold for stronger opioids. There was marked variation among general practices in the odds of patients stepping up to stronger opioids compared with those not stepping up (range 0.31–3.36), unexplained by practice-level variables. Stepping up to stronger opioids was most strongly associated with being underweight (adjusted OR 3.26, 1.49 to 7.17), increasing polypharmacy (4.15, 3.26 to 5.29 for 10 or more repeat prescriptions), increasing numbers of primary care appointments (3.04, 2.48 to 3.73 for over 12 appointments in the year) and referrals to specialist pain services (5.17, 4.37 to 6.12). Compared with women under 50 years, men under 50 were less likely to step down once prescribed stronger opioids (0.53, 0.37 to 0.75).

Conclusions While clinicians should be alert to patients at risk of escalated opioid prescribing, much prescribing variation may be attributable to clinical behaviour. Effective strategies targeting clinicians and patients are needed to curb rising prescribing, especially of stronger opioids.

Statistics from

Strengths and limitations of this study

  • Novel use of routine data to identify patient-level and practice-level characteristics associated with individual opioid prescribing trajectories.

  • Findings based on a large sample of primary care practices and patients from diverse metropolitan populations.

  • Identified associations cannot be assumed as casual.

  • Exploratory work involving multiple statistical testing will have increased risk of type 1 errors and hence falsely positive associations.

  • Detection and magnitude of associations will have been affected by imprecision and biases inherent in routinely recorded data.


There is international concern over rising trends in opioid prescribing,1–9 largely attributed to prescribing for chronic non-cancer pain.10 Despite known short-term efficacy of opioids in neuropathic and musculoskeletal pain, evidence is limited on longer term effects on pain, quality of life, functioning, tolerance and addiction.11–14 The degree of pain relief obtained for chronic pain is often of marginal clinical importance.15 Efficacy decreases with long-term opioid use, although dependence can make treatment withdrawal challenging.16 Prescribed opioids are associated with psychosocial problems, hospitalisation and, even for weaker opioids such as codeine, increased mortality.17–20 The ‘opioid epidemic’ in the USA has coincided with sharp increases in prescription opioid-related deaths and overdoses.21

There is concern that patients with chronic pain are inappropriately being moved up the WHO ‘analgesic ladder’, originally developed for cancer pain, without considering alternatives to medication.22–24 Over a third of one patient cohort eventually more than doubled their original dose of opioids while around 10% were prescribed potentially hazardous high doses.23

Previous research has identified patient characteristics associated with long-term opioid prescribing, which include female gender; increasing age; comorbidities; mental health disorders; and prescribed benzodiazepines, anticonvulsants or antidepressants.2 ,23 ,25–29 Strong or higher dose long-term opioid prescribing is associated with female gender, middle age, socioeconomic deprivation, care home residence, mood disorders, neuropathy, low back pain, polypharmacy, prior prescribing of weak opioids and nicotine dependence.9 ,23 ,30 ,31 Initiation of opioid prescribing is associated with higher baseline levels of pain, psychological distress, substance misuse and unhealthy lifestyles.2 ,32–37 However, most of these studies took place in selected patient groups, for example, trial participants or military veterans;34 ,37 less is known about factors associated with such trajectories in typical primary care populations.

Long-term prescribing may follow hospitalisation but usually starts in primary care.29 There are well-documented variations in wider prescribing practice and safety in UK general practice.38 Structural factors such as practice list size and postgraduate general practitioner (GP) training status only marginally explain variations in high-risk prescribing (potentially inappropriate prescribing of drugs to patients vulnerable to adverse drug events)39 and none in strong opioid prescribing.9

We aimed to improve understanding of patient trajectories towards long-term and stronger opioid prescribing in a UK primary care population. We examined trends and analysed routinely recorded primary care population data to identify patient and general practice characteristics associated with long-term prescribing, stronger opioid prescribing, and stepping up or down between weaker and stronger opioids.


Study design and setting

First, we examined trends in opioid prescribing between 1 April 2005 and 31 March 2012 using routinely collected data from general practices in Leeds and Bradford, UK. Second, we conducted a cross-sectional analysis to identify patient and practice characteristics associated with long-term and stronger opioid prescribing (April 2011 to March 2012). Third, we undertook a longitudinal analysis to examine associations with stepping up to or down from stronger opioids (April 2008 to March 2012). Leeds and Bradford have a combined population of almost 1.3 million. Leeds is typical of UK cities in terms of social deprivation indices and demographics. Bradford is more socially deprived and ethnically diverse; one-fifth of the population is of Pakistani ethnic origin.40

Study participants

There were 192 general practices in Leeds and Bradford in May 2012. We invited all 157 practices then using the SystmOne computerised record system (, which permitted centralised data collection, to participate. We used opt-out recruitment to reduce selection bias by facilitating general practices' agreement to share anonymised patient data. Practices not wishing to share data could opt out. Where practices had migrated from other record systems, we included only full year data from SystmOne to ensure data completeness for the periods covered.

We collected patient data in 12-month segments spanning April to March to ensure complete data for each Quality and Outcomes Framework (QOF) reporting period.41 The QOF is a performance management system whereby general practices are remunerated according to achievement of targets reflecting quality of care across four domains of clinical, organisational, patient experience and additional services. We extracted data on all patients registered with practices on 1st April who were prescribed an opioid at any point in that year. Within each 12-month segment, we excluded patients under age 18 at the start of each year, patients with a coded cancer diagnosis and patients with coded substance misuse or prescribed opioids normally used exclusively to treat drug addiction (methadone and sublingual buprenorphine).


Our dependent variables comprised long-term prescribing of any opioid and at least one prescription of a strong opioid at any time point. We defined long-term use as a record of four or more opioid prescriptions within a given study year, assuming an average duration of 4 weeks supply for each prescription.42–44 We categorised opioid strength according to potency.45 The ‘weaker’ group included weak-to-moderate strength opioids and comprised codeine (with or without paracetamol or ibuprofen), dihydrocodeine (with or without paracetamol), tramadol, pethidine, meptazinol and tapentadol. The ‘stronger’ opioids comprised diamorphine, morphine, oxycodone, fentanyl, hydromorphone, buprenorphine (excluding preparations used for substance misuse), pentazocine, dipipanone and papaveretum. We defined stronger opioid use as at least one prescription of the latter group within the year.

Independent variables included patient-level and practice-level factors we hypothesised or recognised from earlier work to be associated with opioid prescribing.2 ,18 ,25 ,30 ,32–37 At the patient level, we examined demographic, lifestyle, illness, prescribing and health service utilisation variables. Demographic variables comprised age, gender and ethnicity. Lifestyle variables comprised alcohol use, smoking and body mass index (BMI). We assigned patients with clinical codes associated with alcohol use to one of two groups—problem drinker and not a problem drinker. The problem drinker group consisted of males recorded as drinking 50 units or more and females drinking 35 units or more per week. We also included those with clinical codes indicating heavy use or referral to alcohol services. We used the most recently recorded BMI in each year and assigned patients accordingly as underweight (BMI under 18.5), normal (18.5–25), overweight (25–30) and obese (over 30).

Illness variables included pain-associated conditions and clinical presentations, mental illness and comorbidity. We reviewed diagnostic codes and used an iterative process (involving BL, AH, RF and DP) to categorise codes into three groups which we hypothesised might affect GPs' prescribing behaviour: definitively diagnosed disease (eg, rheumatoid arthritis), clinical presentations without a definite diagnosis (eg, ‘knee pain’) and pain-associated syndromes with stronger psychosocial features (eg, fibromyalgia).

As a reflection of patient comorbidity and complexity, we assessed polypharmacy using the number of different repeat prescription medicines within a study year (1–4, 5–9 and 10 or more). We excluded prescriptions for dressings, stoma equipment and medically indicated food.

We assessed health service utilisation using general practice consultations, missed appointments and specialist referrals. We could not distinguish which type of healthcare provider patients consulted (eg, GP, practice nurse). We summed general practice consultations and missed appointments separately for each patient per annum and categorised them as low (<5), medium (between 5 and 10) and high (>10).

We categorised referrals by coded groups of destination specialty: diagnostic test (eg, X-ray), medicine, musculoskeletal, neurology, specialist pain management and surgery. We excluded codes used <100 times by all participating general practices over the study period.

General practice-level variables comprised registered patient list size, numbers of partners and salaried GPs, and training status—all recorded as of 1 April 2011. We used practice-level Index of Multiple Deprivation (IMD) scores. The IMD measures area deprivation and is determined for each patient on the list, where available, and then averaged over the practice. We used overall achievement in the QOF clinical domain as a proxy measure for overall quality of care.

One consequence of collecting data coded within each 12-month period is that we only captured clinical codes entered within that year. For example, we could capture whether a patient was coded as having an alcohol or mental health problem within the year but not if either was coded in earlier years. Our analysis was therefore restricted to incident or rerecorded diagnoses but did not include prevalent or previous diagnoses.

Data sources

The West and South Yorkshire and Bassetlaw Commissioning Support Unit remotely extracted data from participating SystmOne practices. The Unit also supplied information on general practice characteristics.


For our cross-sectional study, we analysed only 2011–2012 data to examine short-term versus long-term and stronger versus weaker opioid prescribing. This was to avoid bias that would arise from studying later samples: when practices changed from one patient record system to another (SystmOne), they only migrated live patient records, so patients that were alive early in our study period but died before the practice changed records system would have been omitted from the study.

For our longitudinal analysis of stepping up or down between weaker and stronger opioids, we assembled a cohort of patients prescribed an opioid during each of 4 years (2008–2012) and prescribed at least one strong opioid during 1 or more but not all 4 years. We excluded patients prescribed stronger opioids throughout all 4 years from the step-up analysis because they had no opportunity to step up. We excluded patients prescribed only weaker opioids from the step-down analysis since they have no opportunity to step down. We did not examine stepping up from none to weaker opioids as we lacked complete records from prior years of patients not prescribed opioids. We used the final 4 years of data for the analysis on trajectories for pragmatic reasons as it delivered a sufficiently large data set of patients with complete data.

We initially calculated unadjusted ORs and then adjusted for confounders and clustering of patients within practices. We modelled associations with long-term and stronger opioid prescribing using multilevel logistic regression with a random effect for practice to account for clustering of patients within practices. We found that practices coded ethnicity into 79 different categories, making analysis potentially problematic. We therefore classified ethnicity into four crudely aggregated groups for cross-tabulation of ethnicity with opioid strength and duration: British, Pakistani, other and unknown. We included a random term for ethnicity for the logistic regression modelling, permitting use of all of the recorded ethnicity categories. We considered this to be an efficient and appropriate way to handle this variable. Thus, we considered patients clustered within practice and within ethnicity.

We sought a parsimonious multivariate model to ensure greater robustness of our findings. We fitted models to show effects adjusted for other covariates and excluded covariates if their effect was negligible and did not reach statistical significance at the 5% level. For adjusted ORs, if a covariate or factor was dropped to achieve a more parsimonious model, then no OR is reported.

Practice characteristics were strongly correlated with one other as expected, with the number of GPs increasing with list size. We also included practice IMD score in the model. Partly for simplicity of interpretation, and concerns over non-linearity, we treated all covariates as categorical variables (eg, practice IMD). We judged that any subsequent measurement error would be compensated for by greater ease of interpretation. We generally used an ‘unknown’ category for missing variables (eg, BMI). For variables with very few missing values (eg, age), we excluded patient records as increasing the number of parameters would not have been sufficiently compensated for by including larger numbers of patient records.

We anticipated an interaction between sex and age and accounted for this in the analysis. We handled age in three categories: 18–49, 50–64 and over 65 years.

All analyses conducted used R statistical software, specifically Revolution R Open V.3.2.0, with the lme4 library V.1.1.7 (R Core Team. R: A language and environment for statistical computing. Secondary R: A language and environment for statistical computing 2015.

Study size

The sample size was originally based on the assumption that we would recruit at least 70% of the then estimated 185 practices using SystmOne (later found to be only 157 practices), that is, 148 practices and around 800 000 registered patients in total. We assumed the prevalence of repeat prescribing of (weaker and stronger) opioids was around 5%, and that a sample size of ∼600 000 patients over 18 years of age, even after allowing for clustering within practices, would provide an accuracy within 0.5%. This was based on an estimated design effect of 50 and calculating two SDs. We anticipated that an effective patient sample size of 12 000 would provide scope to robustly estimate patient-level parameters for covariates and factors. With 148 practices, we also anticipated good scope to examine practice-level covariates, so that our models could robustly estimate up to 14 practice-level covariates using the ratio of 10 observations per covariate.



One hundred and eleven (70.7%) practices shared patient data. Participating practices had a higher median number of all registered patients compared with non-participants (7091 vs 5857; p=0.04), a higher mean number of GPs (5 vs 4; p=0.04) and similar median practice IMD (33.7 vs 33.1; p=0.77).

The number of practices with available data rose from 43 in 2005–2006 to 111 in 2011–2012 (table 1). The number of patients prescribed an opioid at least once rose from 26 249 (6.6% of total adult patients) to 99 847 (12.8%) over the same period, and covered a total of 471 828 patient-years. We observed 471 828 patient-years in which all patients included had received at least one opioid prescription over the 7-year period.

Table 1

Practices, patients and opioid prescribing over 2005–2012

For the 2011–2012 cross-sectional analysis, data for almost all patients were complete. Data on gender were unrecorded and treated as missing for three patients. Some practice-level variables were unavailable for one practice with 469 patients registered.

Of the sample, 60 776 (60.9%) were female, 41 016 (41.1%) aged 18–49 years at the start of the year, 26 460 (26.5%) 50–64 years and 32 371 (32.4%) 65 years or over. Most were categorised as British (55 150; 55.2%), while 11 464 (11.5%) were Pakistani, 10 802 (10.8%) other and 22 431 (22.5%) of unknown ethnicity. Table 2 summarises other patient and practice characteristics. Having earlier excluded patients with a code for cancer or treatment for substance misuse, we categorised 37 508 (37.6%) as prescribed long-term opioids, equivalent to 4.8% of the adult practice population. A total of 6605 (6.6%) of the sample were prescribed stronger opioids, or 0.85% of adults.

Table 2

Patient and practice characteristics for cross-sectional analysis of associations with long-term and stronger opioid prescribing

The longitudinal analysis included 17 165 patients prescribed any opioid in each of the 4 years considered. For the analysis of stepping up to stronger opioids, we excluded 810 patients prescribed stronger opioids for each of the first 3 years (2008–2011). Thus, there were 16 355 patients with three opportunities each to step up, or 49 065 patient-years. Of these, 32 223 (65.7%) were female, 12 354 (25.2%) aged 18–49 years, 16 343 (33.3%) 50–65 years and 20 866 (42.5%) 65 years or over. Most patient-years were British (25 821; 52.6%), while 8004 (16.3%) were Pakistani, 4844 (9.9%) other and 10 515 (21.4%) of unknown ethnicity.

For the analysis of stepping down to weaker opioids, we excluded those 14 918 patients prescribed weaker opioids in the first 3 years. After excluding 3 patients with incomplete data, there were 2244 patients, or 6732 patient-years, available for the step-down analysis. Of these, 4662 (69.3%) were female, 1583 (23.5%) were aged 18–49 years, 2142 (31.8%) were 50–64 years and 3007 (44.7%) were 65 years or over. Most patient-years were British (4160; 61.8%), while 480 (7.1%) were Pakistani, 4844 (9.9%) other and 1476 (21.9%) were of unknown ethnicity. Table 3 summarises sample characteristics for the longitudinal analyses.

Table 3

Patient and practice characteristics for longitudinal step-up and step-down analyses

Prescribing trends

The proportion of all patients prescribed a weaker opioid at least once increased over 7 years from 6.5% to 12.4%, while the proportion prescribed a stronger opioid at least once increased from 0.13% to 0.85% (table 1). For patients prescribed any opioids, the proportion prescribed stronger opioids increased from 4.0% to 13.7%.

Cross-sectional analyses

The odds of long-term against shorter term opioid prescribing varied from 0.59 to 1.59 for practices, while those for stronger opioid prescribing varied from 0.32 to 2.32, suggesting major influences of practice on prescribing. None of the practice-level characteristics examined was associated with prescribing, except for practice IMD for long-term prescribing (1.21; 1.09–1.35; table 4).

Table 4

Cross-sectional analysis of associations with long-term and stronger opioid prescribing

For ethnicity, the odds of long-term against shorter term prescribing varied from 0.65 to 1.64, based on 85 223 patients with complete data on covariates. Plots of practice residuals and ethnicity residuals were consistent with the assumption that the random effects followed Gaussian distributions for the long-term prescribing model. The ORs for stronger against weaker-only opioid prescribing ranged from 0.50 to 2.29, further indicating a strong association between ethnicity and prescribing.

Compared with women aged 18–49 years, long-term opioid prescribing was more likely in women aged 50–64 years (adjusted OR 1.93, 95% CI 1.83 to 2.02) as was stronger opioid prescribing (1.17, 1.06 to 1.29). These associations were respectively stronger in women 65 and over (2.39, 2.28 to 2.51; 1.46, 1.33 to 1.59).

Men below 50 years were more likely to be prescribed long-term opioids than their corresponding female age group (1.25, 1.18 to 1.32) with greater ORs for 50–64 years (1.75, 1.65 to 1.84) and 65 and over (1.76, 1.67 to 1.86). Males aged 18–49 were more likely (1.82, 1.64 to 2.02) and those 65 and over less likely to be prescribed stronger opioids (0.86, 0.78 to 0.96).

We found weak associations with coded diagnoses and problems. Long-term and stronger prescribing were associated with mental health problems (1.12, 1.03 to 1.23; 1.23, 1.06 to 1.23, respectively). Long-term and stronger prescribing were also associated with definitively diagnosed disease (1.09, 1.05 to 1.14; 1.04, 1.04 to 1.20, respectively), while long-term prescribing was inversely associated with clinical presentations lacking a definitive diagnosis (0.83, 0.81 to 0.86). Being overweight or obese was inversely associated (0.87, 0.79 to 0.95) whereas being underweight was associated with stronger opioid prescribing (1.89, 1.34 to 2.67).

Increasing polypharmacy was incrementally associated with long-term and stronger opioid prescribing, with adjusted ORs of 5.68 (5.35 to 6.03) and 6.55 (5.81 to 7.38), respectively, for those having 10 or more repeat prescriptions.

Relative to patients attending up to 6 primary care appointments over the year, those with 7–12 appointments were less likely to be prescribed long-term opioids (0.90, 0.86 to 0.93), while those attending 12 or more appointments were more likely to be prescribed long-term and stronger opioids (1.07, 1.03 to 1.12; 2.03, 1.86 to 2.21, respectively). Long-term and stronger prescribing were also both more likely if patients had missed booked primary care appointments, with the respective ORs rising from 1.08 (1.04 to 1.12) and 1.12 (1.06 to 1.20) for one to two missed appointments to 1.31 (1.24 to 1.39) and 1.68 (1.53 to 1.84) for three or more.

Long-term prescribing was less likely for patients referred for diagnostic tests (0.95, 0.90 to 0.99) and to general medicine (0.87, 0.83 to 0.92) and more likely for those referred to musculoskeletal services (1.19, 1.14 to 1.24). Long-term and stronger opioid prescribing were associated with referrals to neurology (1.27, 1.15 to 1.40; 1.52, 1.32 to 1.75, respectively) and specialist pain management services (3.70, 3.28 to 4.18; 5.74, 5.09 to 6.47, respectively).

Longitudinal analyses

The odds of stepping up against not stepping up to a stronger opioid varied markedly from 0.31 to 3.36 by practice but less so, from 0.91 to 1.06, for stepping down. Practice-level variables explained none of the variation (table 5). There were variations by ethnic category, with ORs ranging from 0.55 to 1.40 for stepping up and 0.87 to 1.24 for stepping down. Compared with women aged 18–49 years, men aged 65 and over were less likely to step up to a stronger opioid (0.77, 0.63 to 0.94), while those aged 18–49 years were less likely to step down (0.53, 0.37 to 0.75).

Table 5

Longitudinal step-up and step-down analyses

Stepping up was associated with coded definitive diagnoses (1.64, 1.46 to 1.85), clinical presentations lacking a definitive diagnosis (1.60, 1.43 to 1.79) and being underweight (3.26, 1.49 to 7.17). Increasing polypharmacy (4.15, 3.26 to 5.29 for 10 or more repeat prescriptions) was strongly associated with stepping up and inversely associated with stepping down (0.60, 0.46 to 0.77 for those with 10 or more prescriptions). Increasing numbers of primary care appointments were also associated with stepping up, with the OR reaching 3.04 (2.48 to 3.73) for over 12 appointments in the year.

Referrals to neurology and specialist pain management services were associated with stepping up (1.77, 1.42 to 2.21; 5.17, 4.37 to 6.12, respectively) and inversely associated with stepping down (0.59, 0.38 to 0.91; 0.61, 0.45 to 0.83, respectively). Referrals to general medicine (1.17, 1.02 to 1.34), musculoskeletal services (1.48, 1.31 to 1.66) and surgery (1.29, 1.09 to 1.52) were also associated with stepping up.

Table 6 provides an overall summary of our findings.

Table 6

Overall summary of associations with long-term and stronger opioid prescribing (+=adjusted OR of 1–2; ++=OR of 2–3; +++=OR over 3; italics=OR <1)


Opioid prescribing in primary care has risen markedly, even after excluding patients with a code for cancer or drug dependence. The proportion of all patients prescribed a weaker opioid at least once almost doubled over 7 years, while the proportion prescribed a stronger opioid increased over sixfold. Much of this prescribing was to patients with unspecified non-malignant pain rather than to those with recorded diagnoses of a specific disease. We identified substantial, poorly explained variation among general practices in long-term and stronger opioid prescribing, particularly in the likelihood of stepping up to stronger opioids. Prescribing patterns varied by age, gender and ethnicity and suggested different risk groups; almost a third of all long-term opioid prescribing and of strong opioid prescribing occurred in women 65 and over. Younger men were particularly more likely to experience long-term and stronger opioid prescribing and less likely to step down once prescribed stronger opioids.

Our findings are based on a large sample of relatively unselected primary care patients from diverse UK metropolitan populations. While there is now a sizeable literature on risk factors for long-term or stronger opioid prescribing,2 ,9 ,23 ,25–27 ,2931 less is known about patient trajectories. We found that stepping up from weaker to stronger opioids was most strongly associated with polypharmacy, more frequent primary care attendances, being underweight and pain service referrals. Patients were likely to step up regardless of whether or not painful conditions were definitively diagnosed, but those with no definitive physical diagnosis were less likely to be prescribed longer term opioids.

Our study has several limitations. First, we cannot imply causation from association. For example, patients with a coded mental health problem could be more likely to be prescribed an opioid or they could have developed a mental health problem as a result of chronic pain for which they were prescribed an opioid. Second, our exploratory work involved multiple statistical testing, increasing the risk of type 1 errors and hence falsely positive associations. Third, we recognised the potential for misclassification bias from our use of routinely recorded and coded data. For example, although musculoskeletal presentations to GPs are recorded in patient notes in 85% of such encounters, they are only clinically coded 32% of the time.47 Under-recording may have been more likely for diagnoses and clinical findings associated with stigmatisation (eg, mental illness, alcohol misuse). We anticipated that such imperfectly and incompletely coded data would reduce the identification and strength of any true associations. However, we sought typically recorded data in a typical sample of practices to enhance generalisability to other UK general practices and the utility of our findings. Any subsequent intervention to change prescribing practice based on searches of patient records would need to work in ‘real-world’ practices rather than those selected for better quality coding. Further imprecision was caused by our use of practice-level rather than patient-level deprivation. Fourth, we used data from only one practice computerised record system. Practice use of SystmOne is associated with higher QOF achievement; it is not known whether particular system characteristics facilitate higher quality of care, better data recording or both.48 This may be of limited relevance to our study as our dependent variables were not directly related to QOF targets. Fifth, the detection and magnitude of associations concerning diagnostic codes was reduced by ascertainment bias because we could only capture diagnoses and events recorded within each 12-month period. Subsequently, the strengths of associations related to mental health problems, alcohol problems and pain-related diagnoses were less than those found in other studies.2 ,27 ,30 ,33 ,35 ,37 Conversely, including ‘ever-coded’ problems would have overestimated current diagnoses. Sixth, our models excluded patients with missing covariates; patients with missing data may systematically differ from those with complete data. Seventh, model selection was achieved through backwards stepwise selection, with careful attention paid to collinearity. This may have led to a model with optimistic goodness-of-fit characteristics (the R2 value will be too high) and for exaggerated coefficients that would be improved by shrinkage. These issues are unlikely to substantially alter the conclusions of our exploratory work. Eighth, participating practices were likely to be larger and have more primary care physicians than non-participating practices. Although smaller practice list size is associated with higher levels of high-risk prescribing,39 our adjusted analyses did not find a relationship between practice size and opioid prescribing. Ninth, the smaller sample size limited the precision of the step-down analysis and ability to detect associations of modest magnitude.

The prescribing patterns we found reflect varying interactions between population need and one or more of patient, clinician and health system behaviour. In women, long-term and stronger opioid prescribing increased with age, probably reflecting increasing population burden and impact of musculoskeletal disorders. The relative risk of stronger opioid prescribing was highest in younger men whereas older men were less likely to step up to and be prescribed stronger opioids. Younger men are prone to riskier and more severe accidental injuries.49 However, our finding that they were less likely to step down from stronger to weaker opioids is also compatible with an explanation of under-recognised or under-recorded opioid dependence.

There were consistent and large variations in prescribing by ethnicity. Patients coded as Pakistani or other ethnicities were particularly less likely than those coded as British to be prescribed long-term and stronger opioids, less likely to step up and more likely to step down. This is harder to explain according to population need and suggests a range of patient and clinician influences, including patient–clinician interactions.

We examined different ways of coding painful conditions to explore the influence of clinician certainty about diagnosis. Many musculoskeletal presentations are self-limiting but those of longer durations eventually acquire a diagnostic label. This is reflected in our finding that long-term and stronger opioid prescribing and stepping up were all associated with definitively diagnosed conditions (eg, osteoarthritis). Long-term prescribing was less likely for clinical presentations without a definitive diagnosis (eg, knee pain); the association with stepping up in this group may partly reflect additional symptomatic coding for exacerbations of definitively diagnosed conditions. However, we found no such associations for the smaller group of patients categorised as having pain associated syndromes with stronger psychosocial features (eg, fibromyalgia). The wider CIs suggest type 2 error but this finding may also reflect clinician reluctance to prescribe opioids without an obvious physical explanation. We could not determine how much opioid prescribing occurred before or as a consequence of specialist referrals; however, the association with referrals to neurology suggest that patients and GPs may struggle to find explanations for pain symptoms.

The association with coded mental health problems is consistent with an explanation that opioid prescribing occurs in response to psychological distress while also causing further dysphoria.

Some findings can only invite more speculative explanations. Being underweight was associated with stepping up to stronger opioids. This may reflect deterioration of long-term conditions, anorexia associated with taking stronger opioids or omitted codes for end-of-life care or drug dependence.

We found strong associations between increasing polypharmacy, which reflects multimorbidity, and all of long-term, stronger and stepping up opioid prescribing. Multimorbidity affects over a fifth of patients in primary care;50 osteoarthritis and cardiometabolic conditions (eg, hypertension, diabetes) are the mostly commonly paired.51 This complicates clinical management, such as prescribing decisions,52 and presents challenges for planning care, especially given that an average consultation covers 2.5 problems in 11.9 min.47 While patterns of opioid prescribing may partly reflect perceived clinical need, decisions to prescribe stronger opioids are often made reactively without explicitly assessing addiction risk or clear action planning.53

Only practice-level deprivation partly accounted for marked variation in prescribing patterns between general practices. Differences among practices in approaches to prescribing cannot be attributed to routinely recorded organisational characteristics such as the number of GPs in each practice or QOF achievement. Practice team and individual clinician behaviours, which appear poorly represented by routinely available indicators, are likely to account for substantial levels of variations in opioid prescribing. These may include established norms and procedures for reviewing prescribing, policies promoting continuity of care and adherence to best practice guidance on opioid prescribing.54 We further recognise the likely influence of wider system and social factors on prescribing, such as local prescribing formularies, legislation on controlled substances and media coverage.

Current prescribing trends signal a need for a concerted response to what is likely to be ineffective and harmful treatment for a large proportion of primary care patients. Once patients with coded cancer or substance dependence are accounted for, around 1 in 20 of the adult practice population is prescribed long-term opioids, a mean of 338 patients per practice. Nearly 1 in 100, or 59 adults per practice, are prescribed stronger opioids. Healthcare systems need to respond to two related problems: how to improve care for patients with chronic non-cancer pain; and how best to curb and reverse current trends in opioid prescribing, especially of stronger opioids.

One response is to recognise and manage chronic pain as a long-term condition which requires coordinated care including the formulation of care plans agreed between professionals and patients, effective non-pharmacological and self-management strategies, and ongoing active monitoring. We have demonstrated that it is feasible to identify patients within primary care based on prescribing records, potentially stratifying those at highest risk of stepping up to stronger opioid treatment as a priority group. These patient groups may extend beyond those with recognised and more severe chronic non-cancer pain. Patients with multimorbidity and polypharmacy represent the most obvious and pressing target for ‘deprescribing’.55 While we would not advocate non-negotiated reductions of regular modest prescriptions of weaker opioids, patients should be alerted to potential risks of continuing treatment at medication reviews.

Another response is to reduce initiation of opioid prescribing. Given the accumulating evidence about their harms, reversing the current trend in opioid prescribing is likely to benefit a substantial at-risk population. There are at least two evidence-based approaches to address this at practice level. First, audit and feedback is effective at changing professional practices.56 Its effects may be relatively modest but audit and feedback can be conducted efficiently using routinely recollected data and its effects may be enhanced by targeting higher baseline non-compliance with standards (eg, higher prescribing practices) and incorporating action plans for practices. Second, computerised prompts triggered by combinations of patient characteristics associated with stepping up to stronger opioids could alert clinicians to reconsider opioid prescribing at the point of care. Computerised clinical decision support is more likely to be effective if clinicians are required to supply a reason for over-riding advice.57

There is a growing consensus on best practice when considering or initiating opioids.58–60 This includes recognising and dealing with psychosocial aspects of pain, managing patient expectations about the degree of pain relief likely to be achievable, starting with a therapeutic trial and an agreement to stop or reduce opioids when they do not work, and recording care plans agreed with patients to guide all subsequent prescribers in maintaining the plan. GPs should avoid escalation above 120 mg morphine-equivalent daily equivalent and refer to pain management services if this seems likely.60

Strategies to bring about large-scale change in healthcare are unlikely to succeed if they fail to address multiple barriers and enablers across different levels of healthcare systems.61 Therefore, these types of above responses need to be coupled with educational campaigns to challenge or change professional and patient beliefs and expectations about opioid prescribing, drawing on lessons learned so far from strategies to reduce inappropriate prescribing of medicines such as antibiotics and minor hypnotics.62

Future research efforts should address strategies to reduce opioid prescribing which target patients and clinicians. There is little evidence on the effectiveness of interventions to reduce opioid prescribing in patients with chronic non-cancer pain.63 Non-pharmacological alternatives to pain management may be modestly beneficial but require further development and evaluation.64 There are parallel needs to further understand practice-level variation and evaluate interventions primarily aimed at clinicians to discourage initiation or continuation of opioid prescribing.


The authors thank the practices which participated in this study. They are also grateful for support provided to the study by Emma Yeomans and Sandra Holliday, Leeds Institute of Health Sciences, and to Paul Lord for the analysis of participating versus non-participating practices.


View Abstract


  • Contributors AH, RF, LG, DP and CM conceived the original idea for the study. All authors contributed to further design and conduct. SF and PC collected data and RW conducted the analyses. All authors contributed to the interpretation of the analyses. RF, BL and RW drafted the manuscript, and all authors read and approved the final manuscript.

  • Funding This paper summarises independent research funded by the National Institute for Health Research (NIHR) under its Research for Patient Benefit Programme (Grant Reference Number PB-PG- 1010-23041).

  • Disclaimer The views expressed are those of the authors and not necessarily those of the NHS, the NIHR or the Department of Health.

  • Competing interests None declared.

  • Ethics approval The study was approved by Yorkshire and the Humber (Humber Bridge) NRES Committee 12/YH/0109.

  • Provenance and peer review Not commissioned; externally peer reviewed.

  • Data sharing statement No additional data are available.

Request permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.