- Research article
- Open Access
- Open Peer Review
The total sale of prescription drugs with an abuse potential predicts the number of excessive users: a national prescription database study
BMC Public Healthvolume 15, Article number: 288 (2015)
Prescription drug sales may vary considerably across regions and over time. This study aimed to assess whether there is an association between mean drug sales and prevalence of excessive use in a range of psychotropic prescription drugs with an abuse potential, and if so, whether the variation in mean drug sales mostly reflects variation in the prevalence of excessive use or mostly reflects variation in non-excessive use.
Data on all filled prescriptions taken from the Norwegian prescription database for 10 drugs with an abuse potential (pain relievers, anxiolytics, and hypnotics) during one calendar year (2005) in Norway (n = 4 053 624) included number of defined daily doses (DDD). These were aggregated to individual level (n = 815 836) and county level (n = 19).
Analyses of individual level data showed that the distribution of drug use was skewed; those who used more than 365 DDD per year accounted for almost half of the sales of both anxiolytics and hypnotics. At the county level, the mean sales per inhabitant and the prevalence of excessive users were closely correlated, but both prevalence of non-excessive use and prevalence of excessive drug use were associated with the county-wise variation in mean drug sales.
Despite a strong individual control of access to psychotropic drugs through health personnel’ prescribing, a small proportion of users account for a large fraction of the sales of these drugs. The sales vary significantly between regions and this variation is closely associated with the prevalence of excessive users. This suggests that sales figures as such may be used as an indicator to monitor variations in excessive use between regions and over time, and to evaluate interventions targeting over-prescription and excessive use.
Opioids, hypnotics, and anxiolytics are drugs mainly prescribed for the treatment of pain, insomnia or anxiety. In high income countries these drugs are widely used, however, they have a potential for abuse and are also used in excessive and non-therapeutic amounts by some patients [1-3]. Such excessive long-term use carries an elevated risk of many and diverse negative consequences. Beyond the risk of dependence and withdrawal symptoms that apply to all these drugs [4-6], the risks of long-term use of hypnotics and anxiolytics include also impaired cognitive functions and memory loss, impaired motor function and risk of falls and injuries, paradoxical effects such as anxiety, restlessness and insomnia, and aggravated dementia . In most countries such drugs are therefore subject to formal control, most often in terms of a prescription system , a system of scheduling, and guidelines advocating restricted use of these drugs. Nevertheless, excessive use or abuse of pharmaceutical drugs also occurs under such regimes . To some extent this can be attributed to diverted drugs and access from illegal markets , but some patients are also prescribed excessive amounts of these drugs . It is the latter aspect of excessive drug use that will be addressed here.
The prescribing of opioids, hypnotics and anxiolytics varies enormously between countries, but a significant variation within countries is also observed, both between regions and over time . For instance, prescriptions of benzodiazepines in Spain and the USA were four to five times higher than that in Germany  and in the USA, prescriptions of opioid analgesics varied by a factor of 12 between states . Most often a prescription system will allow for generating wholesale statistics of drug sales at the aggregate level, but it would be important to know whether – or to what extent – the aggregate variation in these sales is indicative of variation in the prevalence of excessive use and abuse of these drugs.
As the prescription system is meant to prevent excessive use and thereby interfere with a ‘natural’ distribution of drug use, we may assume that the prescription system’s – and mainly the prescribers’ – ability to obtain this goal, will affect the distribution of drug use. Very different distribution patterns of use could therefore – in principle - underlie an overall variation in total sales. The overall variation could e.g. be due to variation in the number of users and the amount of use among non-excessive users with little or no variation in the prevalence of excessive use. Or vice versa; the overall variation in sales could be due to variation in the number of excessive users and their consumption. Or the variation in sales could be due to some combination of the two.
Obviously, the nature of this underlying distribution is not only of theoretical relevance, but also important from a public health perspective. Thus, if the variation could be ascribed mainly to the number of users, but not excessive use, high sales figures could suggest a higher morbidity rate in the studied population or a higher coverage of the target population. This is illustrated with the hypothetical drug use distributions in Figures 1a and b; the higher total drug use in 1b compared to 1a solely reflects an increase in number of non-excessive users, whereas the number of excessive users is the same. If, on the other hand, the variation in sales mainly reflected variation in excessive use, an increase in total drug use would reflect a more right-skewed distribution, as illustrated when comparing Figures 1a to c. In the latter scenario, this would mean that sales figures as such could be used to monitor variations in excessive use. If this were true, wholesale figures, more readily available, could be used to monitor differences in excessive use between regions and over time, and to evaluate interventions targeting over-prescription and excessive use.
So far, the distribution pattern of the consumption of prescription drugs with an abuse potential has scarcely been addressed in empirical studies. Indeed, surveys in general population samples and in institutionalized and other patient samples have demonstrated some aspects of the distribution of consumption, such as differences in prevalence of use and excessive use with respect to gender and age [15,16]. However, few studies have addressed a possible association between total sales of psychoactive drugs and prevalence of excessive use of these drugs, yet there is some indirect evidence suggesting this type of association. In the USA, both sales of opioids and treatment admissions for opioid abuse have increased significantly and in a parallel manner from 1999 to 2009 . In a previous study , we examined the distribution pattern of one prescription drug with abuse potential (carisoprodol) and found that there was a strong and positive association between sales of carisoprodol and prevalence of excessive use. In other words, high total sales of this drug were indicative of high prevalence of excessive drug use, and vice versa. The aim of the present study was to pursue this finding and, by looking at the 10 most accessible prescription drugs with an abuse potential, to assess whether a close association between total drug sales and prevalence of excessive use could be found for a wider range of prescription drugs with an abuse potential, and if so, whether the variation in total drug sales mostly reflects variation in prevalence of excessive use or also reflects variation in non-excessive use.
Drug use data were retrieved from the Norwegian Prescription Database (NorPD). Access to the data was applied for and granted by the Norwegian Public Health Institute. The NorPD covers the entire Norwegian population (5.0 million inhabitants). From 1 January 2004, all pharmacies in Norway were obliged by law to submit electronic data every month on all prescriptions to the Norwegian Institute of Public Health. The NorPD contains information on all prescription drugs, reimbursed or not, dispensed at Norwegian pharmacies to individual patients who live outside institutions [19,20]. Over-the-counter drugs not requiring prescriptions are not registered. For each prescription the following data were collected and used in the present analyses: patients’ unique identifiers (encrypted), county of residence, and drug information; i.e. which drug based on the anatomical–therapeutic–chemical code (ATC-code) and the amount of the drug prescribed converted into number of defined daily doses (DDD). The DDD is the assumed average maintenance dose per day for the drug’s main indication in adults, it is assigned by the WHO collaborating centre and has been used in similar types of studies previously . The data collected for the present study comprised 10 of the most widely used prescription drugs with an abuse potential in Norway , but not drugs on the narcotics list, like morphine and other opioids which are heavily restricted. The drugs were as follows (ACT-code, DDD and main use in parentheses): carisoprodol (M03BA02, 1400 mg, centrally acting muscle relaxant), codeine (N02AA59, 150 mg, opioid pain reliever), diazepam (N05BA01, 10 mg, anxiolytic), oxazepam (N05BA04, 50 mg, anxiolytic), nitrazepam (N05CD02, 5 mg, hypnotic), flunitrazepam (N05CD03, 1 mg, hypnotic), zopiclone (N05CF01, 7,5 mg, hypnotic), and zolpidem (N05CF02, 10 mg, hypnotic), that were filled in Norwegian pharmacies during the calendar year 2005. The year 2005 was chosen as this was the last year before carisoprodol came under scrutiny in the Norwegian market. The number of prescriptions for alprazolam (N05BA12, anxiolytic) and clomethiazole (N05CM02, hypnotic) were too few to analyse separately, but data on prescriptions for these drugs were included in summary measures of anxiolytic drugs and hypnotic drugs, respectively. The data set comprised a total of 4,053,624 prescriptions. Data on population figures were obtained from Statistics Norway.
Excessive use of addictive drugs
Different definitions of excessive use of prescription drugs in general have been applied previously; for instance: a set amount of prescribed drug within a certain time period (e.g. more than 2 DDDs/day over a 30-day period) ; large prescriptions (e.g. more than twice the recommended amount in any prescription) ; and frequent prescribing (e.g. only a few days between prescriptions) . In Norway, the guidelines for prescribing hypnotics and sedatives are that these drugs should be used only temporarily; daily use should be limited to the shortest possible time and no more than two to four weeks. Intermittent treatment is recommended for patients with chronic insomnia . Carisoprodol was not to be used for more than 5–10 days. The use of opioids, including codeine use should be restricted to a minimum, but no strict time is indicated. Thus, we considered drug consumption exceeding daily use of a defined dose for a day (DDD) over an entire year as excessive use and we applied three indicators of excessive drug use for each drug; those who filled prescriptions during a year which amounted to: a) more than 365 DDDs, b) more than 730 DDDs and c) more than 1095 DDDs.
Strategy of analyses
First, all prescriptions were aggregated to the individual level (n = 815 836), summarizing the amount of prescribed drugs in DDDs during the calendar year for each drug and drug group (Figure 2). By applying individual level data for the country as a whole, we described for each drug first the prevalence of past year users, the mean amount of drug use per inhabitant and year in the overall population and the mean and median amount per year among the users. Then, for each drug we described the prevalence of excessive users (applying various thresholds for excessive use) and the fraction of drug sales that these users accounted for.
Next, the individual level data were aggregated to the county level (n = 19) (Figure 2). We assessed whether there was any significant variation in mean drug sales and prevalence of excessive users between counties, applying maximum-minimum value ratios. We then assessed whether drug sales and prevalence of excessive users correlated, applying Spearman’s rank correlation coefficient. Due to low numbers, the third indicator of excessive drug use (>1095 DDDs/year) could not be included in the analyses at the county level.
Finally, we assessed to what extent the variation in drug sales between counties could be attributed to the prevalence of excessive users and the prevalence of non-excessive users. Linear ordinary least square (OLS) regression models were estimated, and input variables included by model fit criteria. In these analyses input variables comprised prevalence figures for three mutually exclusive categories of drug users; non-excessive users (<365 DDDs/year), moderately excessive users (365–730 DDDs/year), and highly excessive users (>730 DDDs/year). Beta coefficients were obtained to assess the relative importance of the input variables to predict variation in drug sales and model R square was obtained to assess the overall predictive value of the input variables included in each model.
All analyses were first carried out for each drug. Then all analyses were carried out for two main groups of drugs; all anxiolytic drugs (diazepam, oxazepam, and alprazolam) and all hypnotic drugs (nitrazepam, flunitrazepam, zopiclone, zolpidem, and clomethiazole), and finally all analyses were carried out for all 10 drugs together (carisoprodol, codeine, diazepam, oxazepam, alprazolam, nitrazepam, flunitrazepam, zopiclone, zolpidem, and clomethiazole).
Distributions at the individual level
The prevalence of users and the mean sales per inhabitant varied significantly between the eight drugs; the highest prevalence of users was found for codeine (8.8% of the population) and zopiclone (6.3%) and the lowest for flunitrazepam (0.3%). Within one year, altogether 17.7% of the Norwegian population had filled one or several prescriptions for one or several drugs with an abuse potential. Mean sales per inhabitant per year were highest for zopiclone (10 DDD/inhabitant/year) and about equally low for carisoprodol, flunitrazepam and zolpidem (about 1 DDD/inhabitant/year) (Table 1). Mean consumption per user and year was particularly high for flunitrazepam and nitrazepam. For all drugs the mean consumption among the users exceeded the median considerably, implying a distribution heavily skewed to the right (Table 1).
In line with the skew distributions we observed that a substantial proportion of users filled prescriptions exceeding 365 DDD per year, this proportion was particularly high for hypnotics (Table 2). The proportion of users filling prescriptions exceeding 730 DDD per year was not very high, but far from negligible, i.e. it was mostly around 1%, yet among flunitrazepam users, 6% filled prescriptions exceeding 730 DDD per year (Table 2). These excessive users accounted for a disproportionately high fraction of the total drug sales. Those who filled prescriptions exceeding 365 DDDs accounted for almost half of the total sales of anxiolytics and hypnotics, and those who filled prescriptions exceeding 730 DDDs accounted for almost one fourth of all sales of anxiolytics and for 13% of the total sales of hypnotics (Table 2). Even the very few who filled prescriptions exceeding 1095 DDDs accounted for a non-negligible fraction of the total drug sales; more specifically, they accounted for 15.8% of all anxiolytics sales and 6.4% of all hypnotics sales (Table 2).
Distributions and associations at the county level
The mean sales – and assumingly consumption – per inhabitant varied considerably across counties; the maximum-minimum ratio was of the magnitude of 2 to 3 for the various drug categories, and the variation in the prevalence of excessive users tended to be higher (Table 3). For all drug categories the mean sales per inhabitant was closely associated with the prevalence of users receiving more than 365 DDDs/year (Spearman’s rank correlation coefficients ranging between 0.87 and 0.99) (Table 3). There was also a positive and statistically significant correlation between mean sales and prevalence of users receiving more than 730 DDDs/year, yet of a somewhat more moderate magnitude (ranging between 0.52 and 0.89) (Table 3).
Finally, we addressed whether the variation in drug sales across counties was merely reflecting variation in excessive use or whether it also reflected the prevalence of users. For all drugs the prevalence of non-excessive users (<365 DDDs/year) and the prevalence of excessive users both contributed significantly to explain the variation between counties in drug sales per inhabitant. The beta coefficients – and thus the relative importance – of prevalence of moderately excessive users (365–730 DDDs/year) varied somewhat across drug categories and were higher than that for non-excessive users and highly excessive users (>730 DDDs/year) with respect to codeine, anxiolytics and hypnotics (Table 4). For five of the eight drugs (diazepam, nitrazepam, flunitrazepam, zopiclone and zolpidem) and for both drug groups (anxiolytics and hypnotics), even the prevalence of highly excessive users contributed significantly to the variation in drug sales between counties (Table 4).
This study has demonstrated that across a wide a range of prescription drugs with an abuse potential the distribution of consumption as measured by filled prescriptions was heavily skewed, and that the relatively few excessive users accounted for a disproportionately high fraction of the drug sales. This was true for drugs like benzodiazepine anxiolytics and hypnotics, z-hypnotics and weak opioid analgesics. Moreover, the sales and prevalence of excessive users varied significantly between counties and there was a strong positive correlation between the two. Thus, the prevalence of excessive users accounted for much of the variation in drug sales, yet the prevalence on non-excessive users also contributed to the variation in drug sales between counties.
Our observations of a very skewed distribution of use of prescription drugs with an abuse potential resemble those from some previous studies of such drugs . Moreover, these findings add to a meagre literature on the distribution of prescription drug use and the association between drug sales on the one hand and indicators of excessive use on the other. It demonstrates that what we found earlier for carisoprodol alone  is applicable to a series of drugs with abuse potential, indicating that whole sales of various drugs with an abuse potential is closely related to excessive use. This is in line with previous observations of a correlation between wholesale of prescription opioids and treatment admissions , and relevant to the observed correlations between wholesale and overdoses of opioid pain relievers [17,28] and of carisoprodol [29,30].
Study strengths and limitations
The major strength of this study is that we have a complete prescription register providing data on all psychotropic drug sales in the non-institutionalized population. The data are not subject to attrition bias and reporting bias as in population surveys. The data are accurate with respect to drug, drug category and dispensed amounts as opposed to self-reported data. In this study we have large numbers of patients who have filled prescriptions, which allows for comparisons of distribution measures across geographical units, i.e. counties.
The study limitations include not having a clinical diagnosis of drug abuse, but having to rely on figures of filled prescriptions. We have attempted to compensate for this by introducing a conservative operationalization of excessive use. This implies that among those categorized as ‘non-excessive users’, a fraction have filled prescriptions exceeding the Norwegian national guidelines limits for prescriptions of opioids, benzodiazepines and z-hypnotics. For this reason, the prevalence of excessive use of these drugs may have been underestimated. Another limitation is that filled prescriptions may not necessarily have been entirely consumed by the recipients. Some drugs may have been diverted (sold or given away to others) or been lost, wasted or returned unused to the pharmacy and this is probably more likely to have occurred among excessive users. However, drug users may have supplemented their drug access with illegal or diverted drugs. Thus, the prescription register base data do not provide a perfect and accurate picture of excessive use (whatever the definition is).
Finally, the data set dates 10 years back in time. While figures for total sales and prevalence of excessive use for each drug may change over time, for instance in response to change in scheduling of particular drugs, the mechanisms underlying the strong associations between total sales and excessive use are probably less subject to change over time. In fact, the 10 drugs we studied differed hugely with respect to total sales and prevalence of users and for all these we found strong associations between sales and excessive use. This suggests that the strength of these associations is not very sensitive to the level of sales. Yet, further studies are needed to assess the generalizability of the findings both over time and across jurisdictions.
These findings are important as they show that changes – or differences – in excessive use of prescription drugs can be predicted with some accuracy on the basis of sales data only. Consequently, this may also apply to the prevalence of problems related to such use. In the absence of a programme for individual monitoring of prescription drug use, drug sales at the aggregate level thus seems to provide a useful proxy tool for assessment of variation in excessive drug use and may be a useful indicator of prevalence of excessive use in studies evaluating possible impacts of policy changes or other interventions targeting excessive psychotropic drug use. For instance, the increase in problems connected with prescription opioids in the USA  could have been predicted by the sharp increase in sales of these drugs alone .
Further research in this area is warranted and may take several directions. For instance, replications of the present study on prescription drugs with an abuse potential could address external validity of the present findings. It would be important to address whether these distributional patterns apply to other - often widely used - addictive prescription drugs, such as oxycodone and methylphenidate, and whether they apply also to non-addictive drugs. Moreover, it will be highly important to assess the prescribers’ role in excessive drug use. For instance, more knowledge is needed regarding the extent of ‘doctor shopping’ among excessive drug users, as well as a better understanding of the process and interaction between patient and prescriber that underlie the development of excessive use of prescribed drugs.
Despite that access to psychotropic drugs is subject to a strong individual control through clinicians’ prescribing, a non-negligible proportion of psychotropic drug users have filled prescriptions within a year clearly indicating excessive use, and these excessive users account for a disproportionately large fraction of the sales of these drugs. The sales vary significantly between regions and this variation is closely associated with the prevalence of excessive users. This suggests that sales figures as such may be used as an indicator to monitor variations in excessive use between regions and over time, and to evaluate interventions targeting over-prescription and excessive use.
Anatomical Therapeutic Chemical (Classification System)
Defined Daily Doses
Norwegian Prescription Database
Compton WM, Denisco R. Prescription Drug Abuse. In: Galanter M, Kleber HD, editors. The American Psychiatric Publishing Textbook of Substance Abuse Treatment. 4th ed. Washington, DC: American Psychiatric Publishing, Inc.; 2008. p. 595–607.
Casati A, Sedefov R, Pfeiffer-Gerschel T. Misuse of medicines in the European Union: a systematic review of the literature. Eur Addict Res. 2012;18:228–45.
Ghodse H. Pain, anxiety and insomnia–a global perspective on the relief of suffering Comparative review. Br J Psychiatry. 2003;183:15–21.
Longo LP, Johnson B. Addiction: Part I. Benzodiazepines - side effects, abuse risk and alternatives. Am Fam Physician. 2000;61:2121–8.
Griffiths RR, Johnson MW. Relative abuse liability of hypnotic drugs: a conceptual framework and algorithm for differentiating among compounds. J Clin Psychiatry. 2005;66(supplement 9):31–41.
Fishbain DA, Cole B, Lewis J, Rosomoff HL, Rosomoff RS. What Percentage of Chronic Nonmalignant Pain Patients Exposed to Chronic Opioid Analgesic Therapy Develop Abuse/Addiction and/or Aberrant Drug‐Related Behaviors? A Structured Evidence‐Based Review. Pain Med. 2008;9:444–59.
Lader M. Benzodiazepines revisited—will we ever learn? Addiction. 2011;106(12):2086–109.
Babor T, Caulkins J, Edwards G, Fischer B, Foxcroft DR, Humphreys K, et al. Drug Policy and the Public Good. Oxford: Oxford University Press; 2010.
Neutel CI, Skurtveit S, Berg C. What is the point of guidelines? Benzodiazepine and z-hypnotic use by an elderly population. Sleep Med. 2012;13:893–7.
Inciardi JA, Surratt HL, Cicero TJ, Kurtz SP, Martin SS, Parrino MW. The “black box” of prescription drug diversion. J Addict Dis. 2009;28:332–47.
Ashton H. The diagnosis and management of benzodiazepine dependence. Curr Opin Psychiatry. 2005;18:249–55.
Manchikanti L, Singh A. Therapeutic opioids: a ten-year perspective on the complexities and complications of the escalating use, abuse, and nonmedical use of opioids. Pain Physician. 2008;11(2 Suppl):S63–88.
Khong T, de Vries F, Goldenberg JSB, Klungel O, Robinson N, Ibáñez L, et al. Potential impact of benzodiazepine use on the rate of hip fractures in five large European countries and the United States. Calcif Tissue Int. 2012;91(1):24–31.
Curtis LH, Stoddard J, Radeva JI, Hutchison S, Dans PE, Wright A, et al. Geographic variation in the prescription of schedule II opioid analgesics among outpatients in the United States. Health Serv Res. 2006;41(3p1):837–55.
Donoghue J, Lader M. Usage of benzodiazepines: A review. Int J Psychiatry Clin Pract. 2010;14(2):78–87.
Fischer B, Ialomiteanu A, Boak A, Adlaf E, Rehm J, Mann RE. Prevalence and key covariates of non-medical prescription opioid use among the general secondary student and adult populations in Ontario, Canada. Drug Alcohol Review. 2013;32(3):276–87.
Paulozzi L, Jones C, Mack K, Rudd R. Centers for Disease Control and Prevention (CDC). Vital Signs: overdoses of prescription opioid pain relievers–United States, 1999–2008. MMWR Morb Mortal Wkly Rep. 2011;60(43):1487–92.
Bramness JG, Rossow I. Can the total consumption of a medicinal drug be used as an indicator of excessive use? The case of carisoprodol. Drugs Educ Prev Policy. 2010;17(2):168–80.
Hausken AM, Furu K, Skurtveit S, Engeland A, Bramness JG. Starting insomnia treatment: the use of benzodiazepines versus z-hypnotics. A prescription database study of predictors. Eur J Clin Pharmacol. 2009;65(3):295–301.
Furu K. Drug utilisation in a public health perspective: Establishing a national prescription register in Norway. Nor J Epidemiol. 2001;11(1):55–60.
Mason WA, Hitchings JE, Spoth RL. Emergence of delinquency and depressed mood throughout adolescence as predictors of late adolescent problem substance use. Psychol Addict Behav. 2007;21(1):13.
Eitle D, Taylor J, Eitle TM. Heavy episodic alcohol use in emerging adulthood: the role of early risk factors and young adult social roles. J Drug Issues. 2010;40(2):295–320.
Kampmann J. Sovemedicine og nervemedicine - storforbrugere. København: Lægemiddelstyrelsen; 1997.
Isacson D, Bingefors K, Wennberg M, Dahlström M. Factors associated with high-quantity prescriptions of benzodiazepines in Sweden. Soc Sci Med. 1993;36(3):343–51.
Varley H, Jordan J, Bascom P: Quantities of carisoprodol (Soma(R)) to be restricted. DUR Newsletter 2000. http://www.orpdl.org/durm/newsletter/articles/volume2/durv2i4.pdf.
Helsetilsynet. Vanedannende legemidler. Forskrivning og forsvarlighet. Oslo: Helsetilsynet; 2001.
Bramness JG, Furu K, Engeland A, Skurtveit S. Carisoprodol use and abuse in Norway. A pharmacoepidemiological study. Br J Clin Pharmacol. 2007;64:210–8.
Modarai F, Mack K, Hicks P, Benoit S, Park S, Jones C, et al. Relationship of opioid prescription sales and overdoses, North Carolina. Drug Alcohol Depend. 2013;132(1):81–6.
Høiseth G, Bramness J, Christophersen A, Mørland J. Carisoprodol intoxications: a retrospective study of forensic autopsy material from 1992–2003. Int J Legal Med. 2007;121(5):403–9.
Høiseth G, Karinen R, Sørlid HK, Bramness JG. The Effect of Scheduling and Withdrawal of Carisoprodol on Prevalence of Intoxications with the Drug. Basic Clin Pharmacol Toxicol. 2009;105(5):345–9.
Okie S. A flood of opioids, a rising tide of deaths. N Engl J Med. 2010;363:1981–5.
Manchikanti L. National drug control policy and prescription drug abuse: facts and fallacies. Pain Physician. 2007;10(3):399–424.
This study was funded by the Norwegian Institute for Alcohol and Drug Research and the Norwegian Centre for Addiction Research, University of Oslo. The authors are grateful for the constructive comments by three reviewers; Christian J. Teter; Joel Lexchin and Sebastien Haneuse.
The authors declare that they have no competing interests.
IR designed the study, analysed the data and drafted the manuscript; JGB contributed to acquisition of data, interpretation of data and revised the manuscript critically for intellectual content. Both authors have read and approved the final manuscript.