- Research article
- Open Access
- Open Peer Review
Drug-related deaths in Scotland 1979–2013: evidence of a vulnerable cohort of young men living in deprived areas
BMC Public Healthvolume 18, Article number: 357 (2018)
Even after accounting for deprivation, mortality rates are higher in Scotland relative to the rest of Western Europe. Higher mortality from alcohol- and drug-related deaths (DRDs), violence and suicide (particularly in young adults) contribute to this ‘excess’ mortality. Age-period and cohort effects help explain the trends in alcohol-related deaths and suicide, respectively. This study investigated whether age, period or cohort effects might explain recent trends in DRDs in Scotland and relate to exposure to the changing political context from the 1980s.
We analysed data on DRDs from 1979 to 2013 by sex and deprivation using shaded contour plots and intrinsic estimator regression modelling to identify and quantify relative age, period and cohort effects.
The peak age for DRDs fell around 1990, especially for males as rates increased for those aged 18 to 45 years. There was evidence of a cohort effect, especially among males living in the most deprived areas; those born between 1960 and 1980 had an increased risk of DRD, highest for those born 1970 to 1975. The cohort effect started around a decade earlier in the most deprived areas compared to the rest of the population.
Age-standardised rates for DRDs among young adults rose during the 1990s in Scotland due to an increased risk of DRD for the cohort born between 1960 and 1980, especially for males living in the most deprived areas. This cohort effect is consistent with the hypothesis that exposure to the changing social, economic and political contexts of the 1980s created a delayed negative health impact.
Since around 1950 Scotland has had higher, and more slowly improving, mortality rates than the rest of Western Europe [1,2,3]. From 1981 onwards less and less of the higher mortality in Scotland compared to England and Wales could be explained by deprivation [4, 5]. Higher mortality from alcohol- and drug-related deaths (DRDs), violence and suicide (particularly in young adults) and higher mortality from heart disease, stroke and cancer throughout adulthood contribute to this ‘excess’ mortality (i.e. higher mortality over and above that explained by differences in socioeconomic deprivation), the slower improvement in mortality compared with the rest of Europe and the greater health inequalities [6, 7]. It has been suggested that a third of Scotland’s ‘excess’ mortality over England among 15 to 54 year olds can be accounted for by a higher prevalence of problem drug use in Scotland, which rose rapidly in Scotland in the 1980s .
It has been proposed that these mortality phenomena are best explained in terms of the interaction of the neoliberal political and economic policies of the 1980s (which favoured monetarist economic policies and led to a rapid rise in income inequality and an increase in unemployment) with the underlying vulnerabilities in the Scottish population due to historical regional and urban policy (which meant that the population in Scotland had different housing, employment and social context, which was more badly affected) [6, 7, 9, 10]. These policy changes may have impacted on health in the same way as in other countries [11,12,13]. Part of the Scottish population, particularly the working class living in the deindustrialising regions, may have been particularly badly affected by these dramatic changes of the 1980s [6, 10, 11, 14,15,16,17].
If this impact of political change was partially responsible for the mortality phenomena in Scotland, then the worst affected cohort of individuals would be those most directly exposed (i.e. working class, young adults first experiencing working life during 1979–1990). While more males would have been impacted directly in terms of job losses, females would have been affected by the loss of household income and the associated effect of the changing political context. There is, however, ambiguity around this exposure cohort; the timing of the exposure is rather indistinct with deindustrialisation occurring earlier  and the change in political approach persisted after [11, 12].
Period or cohort effects (age effects would have been adjusted for in the age-standardised analyses) may be partially responsible for the trends in DRDs in Scotland, something already identified for alcohol-related mortality and suicide [19, 20]. DRDs have continued to increase in Scotland but with large differences between groups [21, 22]. For instance, males account for the majority of deaths and the average age of those dying has increased reflecting an ageing profile of problem drug users (problem drug use is defined by Information Services Division Scotland (ISD) Scotland as the problematic use of opiates (including illicit and prescribed methadone use) and/or the illicit use of benzodiazepines and implies routine and prolonged use as opposed to recreational and occasional drug use) . Consistent with this, the prevalence of drug use has fallen among adolescents, although the vast majority of the reported drug use at this age relates to drugs that do not carry a high risk of mortality [24, 25]. There are also clear differences in DRD rates by deprivation, with higher rates for those living in more deprived areas [22, 26].
In this paper we explore whether there are age, period or cohort (APC) effects in DRDs in Scotland corresponding to the changing political context from the 1980s, which might help explain these recent trends in DRDs in Scotland.
We obtained data on the number of DRDs by sex, single year of age at death and year of registration of death for Scotland from 1979 to 2013, and mid-year population estimates for Scotland by sex and single year of age, from age 0 to 89 years (with pooled data for those aged 90 years and over), from 1979 to 2013, from the National Records of Scotland (NRS).
DRDs from 1979 to 1999 were defined by the International Classification of Diseases, Ninth Revision 9 (ICD9) codes 292, 304, 305.2–305.9, E850-E858, E950.0-E950.5, E962.0, and E980.0-E980.5 and from 2000 onward by the ICD10 codes F11-F16, F18–19, X40-X44, X60-X64, X85, and Y10-Y14. This follows the Office for National Statistics’ ‘wide’ definition for deaths from drug poisoning, referred to in this paper as DRDs, which include both deaths from poisoning from any drug or medication as well as controlled drugs (see NRS  for a discussion of data availability including why only data on DRDs post-1979 are used and Additional file 1: Figure S1 to compare DRDs determined by other common definitions).
Analyses by deprivation were based on Carstairs deprivation categories and compared the most deprived quintile to the remaining quintiles grouped to create the less deprived four-fifths of the population (this allowed sufficient numbers for analysis) . For DRDs data, Carstairs look-up files were obtained from the ISD website; death data were matched to Carstairs quintiles using the appropriate output area (mean population c.115 in 2011). Data were then assigned to an appropriate Carstairs time period; cases with no associated Carstairs (2.4% of dataset) were removed from the analyses by deprivation. Population data by Carstairs, from age 0 to 84 years (with pooled data for those aged 85 years and over), were obtained from ISD with interpolation between the censuses.
Unless stated otherwise, analyses were restricted to those aged 15–89 years, or aged 15–84 years for deprivation analysis.
Descriptive analysis - rates
For presentation of trends, death rates for all ages from 1979 to 2013 by sex were age-standardised to the 2013 version of the European Standard Population. For age effects analysis, data from 1979 to 1996 and data from 1997 to 2013 were separately combined for presentation of crude death rates by five-year age groups (from 0 to 4 years to 90+ years (85+ years for deprivation analysis)) for age at death by sex.
To examine APC effects, DRD rates were plotted using data grouped by age at death (five-year age groups), year of death (five-year periods) and birth cohort (those born over five-year periods, estimated from age at death and year of registration of death. Additional file 1: Figure S2 shows the cohorts, and their corresponding age at death, covered by our dataset).
Descriptive analysis - age-year and age-cohort shaded contour plots
Shaded contour plots (SCPs) were produced in R version 3.3.0 (in these Lexis surface visualisations, lines and colour show population age groups at particular points in time with the same DRD rates) . Separate plots were produced by sex, and then separately by deprivation, created using rates for single age and year of death or for single age and birth year (for information on SCPs see Appendix 1 in Parkinson et al. (2016) ). In age-year SCPs, where age is on the vertical axis and year on the horizontal axis, period effects appear as vertical patterns, age effects as horizontal patterns and cohort effects as diagonal ‘disruptions’ running along diagonal lines from the bottom left to top right. In age-cohort SCPs the data are rearranged by birth cohort on the horizontal axis so cohorts instead run vertically from bottom to top of the figure and period effects appear as diagonal lines running from bottom right to top left.
Inferential analysis – Intrinsic estimator models
For APC regression modelling, DRD and population data were grouped into five-year age groups and five-year time periods stratified by sex and deprivation. The year for the cohorts was calculated as a function of the five-year age and period groups, as being equal to mean period minus mean age. Standard negative binomial regression models were initially run to ensure that age, period or cohort alone or any combination of two factors did not fit the data better than a full APC model. The intrinsic estimator (IE) command (apc_ie) in Stata version 13 (Stata Corp LP, College Station, Texas, USA) was used . As the count data were over-dispersed, negative binomial models were used. Models were fitted separately for each sex and deprivation specific strata. For the most deprived fifth of females, the 1899 cohort data were replaced by those of the 1904 cohort to obtain a suitable model. This equated to replacing data for age group 80–84 years and period 1979–1983 with that of age group 80–84 years and period 1984–1988. The replacement was made at sub-group level, which would not have significantly affected trends reported.
Age-sex-standardised trends over time
From 1979 until 1986, age-standardised DRD rates (Additional file 1: Figure S3) for both males and females were low and decreased following similar trends and levels. Subsequently, rates rose for both sexes but much faster for males, reaching a level in 2013 far greater than that of 1979, while those for females increased to a level similar to that of 1979.
These worsening trends in DRDs were sharpest for those living in the most deprived areas and for males (Additional file 1: Figure S4). Indeed, for those living in less deprived areas, rates for males increased slowly only after 1991 and those for females remained consistently low.
Age-year and age-cohort shaded contour plots
Rates for DRDs differed markedly over the life-course (Additional file 1: Figures S5 and S6) and the age profiles have changed over the period and differ by sex. The peak rate for females was initially older than that for males, but this fell as rates decreased around 1985 for females aged 50–60 years and then increased around the mid-1990s for young adult females aged 20–30 years resulting in the peak age becoming more comparable to that for males.
Around 1990, the peak age of DRDs fell for males as rates increased for those aged around 18 to 45 years, and especially for those < 30 years. Subsequently the peak age began to increase, suggesting a cohort effect with the most affected cohort born between 1970 and 1980, where DRDs peaked at about 45/100,000 for 25 year olds (around the 1975 cohort) and also 5 years later at 65/100,000 for 30 year olds (Fig. 1). This effect was obvious after calendar year 1990 (red line in Fig. 1). A similar but weaker cohort effect was seen in females born between 1965 and 1980.
The age profiles by deprivation (Additional file 1: Figures S7 and S8) showed clear differences by deprivation and over time for both males and females. The vertical pattern in Fig. 2 indicates a cohort effect of higher DRDs for those born between around 1960 and 1980, with men and those in deprived areas affected most. This effect was again obvious after calendar year 1990 (red line in Fig. 2). The cohort effect also began up to a decade earlier for males in the most deprived group compared to the less deprived group. For men in the most deprived group, the most affected cohort, born between around 1970 and 1975, experienced high DRD risk at two stages in their life course, both in their mid to late twenties and then in their mid to late thirties.
Analysis by period showed that the peak age for DRDs for males decreased from 1979-1983 to 1994–1998 before increasing again as the cohort born in the early 1970s aged, appearing to carry with them an increased DRD risk compared to previous cohorts (Additional file 1: Figure S9). The effect was less clear for females.
The cohort effect was clear for males living in both deprivation groups, with higher rates for those in the most deprived areas (Fig. 3). Rates also started to increase five years earlier for the most deprived group resulting in the peak age group for each period being five years older for this deprivation group. A weaker cohort effect was also evident for females living in the most deprived areas. The overall DRD trends are largely accounted for by the trends in the most deprived population fifth.
Regression analysis confirmed that no single or two factor combination of age, period or cohort provided a better fit than the full APC model (Additional file 2). Results from the APC models showed congruence between the plots of IE coefficients (Figs. 4 and 5) and the SCPs (Figs. 1 and 2).
The IE coefficients (see Additional file 2 for how to interpret IE coefficients) suggested that there is a strong age effect for both sexes, with a peak in middle age (younger and larger among males); a small period effect for males, increasing through the 1990s, but little effect for females; and a birth cohort effect for both sexes, with a low around the 1944 cohort (Fig. 4) (see Additional file 3 for the IE coefficient data with 95% confidence intervals). A larger cohort effect was observed for males than females in the successive birth cohorts 1959 to 1989.
The deprivation stratified data showed similar overall patterns to the whole population for males (Fig. 5) (see Additional file 3 for the IE coefficient data with 95% confidence intervals). The timing of the cohort effects differed between deprivation groups, with the effect since the 1960s increasing five years earlier in the most deprived population. For females, both deprivation populations showed age and cohort effects, similar to those for males, but with slightly different timings, while a period effect similar to that for males was only evident for the most deprived group. As for males, the cohort effect since the 1960s increased earlier in the most deprived population.
Between 1979 and 1986 the DRD rate in Scotland decreased and was similar for males and females. Subsequently, rates for men increased, due to a sharp rise among males living in the most deprived areas. In contrast, rates for females increased more slowly, also primarily due to a rise among females in the most deprived areas, leading to a growing gap between the sexes.
Around 1990, DRD rates increased noticeably (particularly for young males and in the most deprived areas), and the cohort born between around 1960 and 1980 experienced considerably higher DRD rates until the end of the time series and are primarily responsible for the increasing number of deaths over time. A similar although weaker cohort effect was observed for females. The cohort effect was more prominent for those living in the most deprived areas for both sexes, where it also started around a decade earlier than for the rest of the population. The IE analyses supported the descriptive findings of a cohort effect alongside the known age patterning and suspected period effect from the 1990s.
Strengths and weaknesses of the analysis
Defining a DRD is not straightforward and changes over the years make establishing trends problematic . To obtain the longest time trend, we used the ONS ‘wide’ definition rather than the ‘baseline’ definition for the UK Drugs Strategy used to report DRDs by NRS and many other statistical bodies in Scotland [21, 27]. The method for this ‘wide’ definition has not changed over our time series, but its use means that figures for DRDs are higher than those obtained using other definitions, as it includes both deaths from poisoning from any drug or medication as well as controlled drugs. However, long-term trends (if not the levels) are broadly the same using the different definitions (Additional file 1: Figure S1).
Although NRS has information on each death registered in Scotland since 1974, data on DRDs were restricted to post-1979 when NRS started to use ICD9 codes; ICD8 codes were used from 1974 to 1978 but there is no suitable cross-match between codes that may be used for DRDs . The numbers identified by ICD9 codes are considered to give a good indication of what happened between 1979 and 1999 . Thereafter, ICD10 codes have been used. There is some overlap in the ICD codes used to identify DRDs and suicides. These codes concern deaths due to intentional self-poisoning or undetermined intent poisoning by drugs, medicaments and biological substances. The number of DRDs that are a subset of suicides is greater for ICD9 than ICD10 coding and also for the ONS definition than for the NRS used definition, due to the way the data are obtained.
We used data with the year the death was registered rather than the year it occurred (chosen to match the reporting convention of NRS) but the impact on our conclusions is negligible . The disaggregation by age over time led to small numbers in each group and a substantial degree of random variation in the Lexis surface but the confidence intervals for the IE analyses and the rounding in the SCPs ensure that the patterns are very unlikely to be due to chance.
Modelling APC effects simultaneously to disentangle their effects is problematic due to their mathematical co-dependence, the ‘identification problem’. There is disagreement as to whether this is possible [31,32,33,34] and proposed methods to separate APC effects have limitations [35,36,37]. Imposing a constraint on the model, such as applying an assumption of equality between groups to break the co-dependence, is a common solution but the results are sensitive to the chosen assumptions and there is often little basis for making these . The IE approach used here assumes that the sum of the coefficients in the vector are zero thus avoiding the need to make such constraints; although it is not without limitations , it is used here to provide an inductive exploration of the hypothesis to complement the line plots and SCPs. This triangulation adds strength to the findings; there is congruence between the IE plots and SCPs.
Finally, we only have a partial view of each birth cohort as we do not have data for all years for all cohorts. This is always the case with APC analyses and it is mitigated through the comparison of age-specific rates for different cohorts and time periods. We can, however, only have confidence that our observations are valid for the period 1979–2013.
Comparison with other studies
Methodological differences in recording of DRDs make comparison of DRDs between countries problematic. However, a recent comparison indicates that the Scottish drug-death rate could be at least double that of the UK as a whole . The trends in Scottish DRDs over time and differences between males and females and between age groups have been previously described for Scotland . The substantially higher rates of DRDs among those living in deprived circumstances is also well-recognised as is the increase in the average age of individuals dying from DRDs from 1996 to 2014, seen as reflecting an ageing profile of problem drug users and evidenced here as a cohort effect [21, 22, 26].
The findings presented here echo those of recent APC research on suicide in Scotland where a cohort effect was identified for people born between 1960 and 1980 . Both cohort effects are largely driven by males and those living in the most deprived areas but the cohort most affected by suicide occurs slightly earlier than that for DRDs, 1965–1974 versus 1970–1980, respectively. Further, although for both health outcomes the timing of the cohort effect and cohort at highest risk differs by deprivation, being earlier for males in the most deprived group, this difference is not the same. For suicides there is around a five year difference while for DRDs it is around ten years. There are also similarities to the finding of an age-period interaction for alcohol-related deaths in Scotland which impacted on the same birth cohorts as those most affected by suicide and DRDs . These APC results are consistent with these health outcomes being important in explaining a growing relative mortality risk for young to middle-aged males since the mid-1990s .
Recently Walsh et al.  showed that ‘excess’ mortality in Scotland can be best explained by a greater vulnerability, due to historical processes of urban change and selective migration to new towns, to economic and social changes resulting from the political changes of the 1980s. The cohort effect identified here is consistent with exposure of young working-age adults to the changed socioeconomic context during the 1980s, an event which would have occurred primarily when young adults entered the job market and therefore earlier for those in the most deprived areas, and adds further evidence to this thesis. It seems likely that alcohol-related deaths, DRDs and suicide outcomes that contribute to Scotland’s ‘excess’ mortality could share a common causal pathway stemming from the changing social and economic policies of the 1980s. It has also been suggested recently that there is a similar patterning of external causes of death for particular cohorts in the USA . Another possibility is that there was a change in the type of drugs available for the cohorts with the highest DRDs as they moved into adulthood, as hard drugs became more available. However, we would have expected any change in DRDs due to the introduction of new drugs to society to have also been evident in subsequent cohorts.
DRDs negatively impact society [41, 42] and make an important contribution to the health inequalities and ‘excess’ mortality in Scotland . The rationale for our study was to investigate whether there is evidence of a cohort effect for those most likely to have been exposed to the negative consequences of the changing social, economic and political contexts of the 1980s. The cohort effect we identify here for DRDs, which appears for young adult males (and to a lesser extent for females) in the early 1990s, seems consistent with this hypothesis given that young adults would have entered the job market at this time and, particularly in the most deprived areas, been exposed to high unemployment levels and reduced support [43, 44].
The full impact of ‘excess’ mortality, and indeed morbidity, in cohorts with high DRDs is unlikely to be known for some time. As the cohort of people at greatest risk of DRDs continues to age, drugs services will need to adapt to their needs as co-morbidities from chronic conditions associated with ageing and drug use become more prevalent. There is also a risk that more recent exposures to a more ‘flexible’ labour market and greater conditionality and sanctions in the social security system, particularly for young working-age adults, may in time lead to another cohort at high risk [45, 46]. Continued surveillance of this population is therefore merited.
International Classification of Diseases
Information Services Division Scotland
National Records of Scotland
shaded contour plot.
McCartney G, Walsh D, Whyte B, Collins C. Has Scotland always been the 'sick man' of Europe? An observational study from 1855 to 2006. Eur J Pub Health. 2012;22(6):756–60.
Popham F, Boyle P. Assessing socio-economic inequalities in mortality and other health outcomes at the Scottish national level: Final report. Scottish Collaboration for Public Health Research and Policy; 2010. https://calls.ac.uk/wp-content/uploads/final_report_-_assessing_socio-economic_inequalities_in_mortality_and_other_health_outcomes_at_the_scottish_national_level..pdf.
Norman P, Boyle P, Exeter D, Feng Z, Popham F. Rising premature mortality in the UK’s persistently deprived areas: only a Scottish phenomenon? Soc Sci Med. 2011;73(11):1575–84.
Hanlon P, Lawder R, Buchanan D, Redpath A, Walsh D, Wood R. Why is mortality higher in Scotland than in England & Wales? Decreasing influence of socioeconomic deprivation between 1981 and 2001 supports the existence of a ‘Scottish effect’. J Public Health. 2005;27(2):199–204.
Schofield L, Walsh D, Munoz-Arroyo R, McCartney G, Buchanan D, Lawder R, et al. Dying younger in Scotland: trends in mortality and deprivation relative to England and Wales, 1981 –2011. Health Place. 2016;40:106–15.
McCartney G, Collins C, Walsh D, Batty GD. Accounting for Scotland’s excess mortality: towards a synthesis. Glasgow Centre for Population Health; 2011. http://www.gcph.co.uk/assets/0000/1080/GLA147851_Hypothesis_Report__2_.pdf.
Walsh D, McCartney G, Collins C, Taulbut M, Batty GD. History, politics and vulnerability: explaining excess mortality in Scotland and Glasgow. Glasgow Centre for Population Health; 2016. http://www.gcph.co.uk/assets/0000/5988/Excess_mortality_final_report_with_appendices.pdf.
Bloor M, Gannon M, Hay G, Jackson G, Leyland AH, McKeganey N. Contribution of problem drug users’ deaths to excess mortality in Scotland: secondary analysis of cohort study. BMJ. 2008;337:a478.
Galea S, Ahern J, Karpatic A. A model of underlying socioeconomic vulnerability in human populations: evidence from variability in population health and implications for public health. Soc Sci Med. 2005;60(11):2417–30.
McCartney G, Collins C, Walsh D, Batty GD. Why the scots die younger: synthesizing the evidence. Public Health. 2012;126(6):459–70.
Collins C, McCartney G. The impact of neoliberal "political attack" on health: the case of the "Scottish effect". Int J Health Ser. 2011;41(3):501–23.
Scott-Samuel A, Bambra C, Collins C, Hunter DJ, McCartney G, Smith K. The impact of Thatcherism on health and wellbeing in Britain. Int J Health Ser. 2014;44(1):53–71.
Beckfield J, Krieger N. Epi + demos + cracy. Linking political systems and priorities to the magnitude of health inequities—evidence, gaps, and a research agenda. Epidemiol Rev. 2009;31(1):152–77.
Mackenzie M, Collins C, Connolly J, Doyle D, McCartney G. Working-class discourses of politics, policy and health: ‘I don’t smoke; I don’t drink. The only thing wrong with me is my health’. Policy Politics. 2017;45(2):231–49.
Bell I. Still suffering Thatcher’s legacy 30 years on. Herald (Scotland); 2009.
MacWhirter I. That bloody woman. NewStatesman, 2009. https://www.newstatesman.com/uk-politics/2009/02/scottish-scotland-thatcher.
McCormack C. The wee yellow butterfly. Glasgow: Argyll Publishing; 2009.
Walsh D, Taulbut M, Hanlon P. The aftershock of deindustrialization—trends in mortality in Scotland and other parts of post-industrial Europe. Eur J Pub Health. 2010;20(1):58–64.
McCartney G, Bouttell J, Craig N, Craig P, Graham L, Lakha F, et al. Explaining trends in alcohol-related harms in Scotland, 1991-2011 (I): the role of incomes, effects of socio-economic and political adversity and demographic change. Public Health. 2016;132:13–23.
Parkinson J, Minton J, Lewsey J, Bouttell J, McCartney G. Recent cohort effects in suicide in Scotland: a legacy of the 1980s? J Epidemiol Community Health. 2017;71(2):194–200.
National Records of Scotland. Drug-related deaths in Scotland in 2016. Statistics of drug-related deaths in 2016 and earlier years, broken down by age, sex, selected drugs reported, underlying cause of death and NHS Board and Council areas. National Records of Scotland; 2017. https://www.nrscotland.gov.uk/files/statistics/drug-related-deaths/drd2016/16-drug-rel-deaths.pdf.
Barnsdale L, Gordon R, Graham L, Walker D, Elliott V, Graham B. The National Drug-Related Deaths Database (Scotland) Report: analysis of deaths occurring in 2014. Information Services Division Scotland; 2016. https://www.isdscotland.org/Health-Topics/Drugs-and-Alcohol-Misuse/Publications/2016-03-22/2016-03-22-NDRDD-Report.pdf.
ISD Scotland. Estimating the national and local prevalence of problem drug use in Scotland 2012/13. Information Services Division Scotland; Update 2016. http://www.isdscotland.org/Health-Topics/Drugs-and-Alcohol-Misuse/Publications/2014-10-28/2014-10-28-Drug-Prevalence-Report.pdf.
Scottish Government. Scottish Schools Adolescent Lifestyle and Substance Use Survey (SALSUS): Drug Use Report 2015. Scottish Government; 2016. http://www.gov.scot/Resource/0050/00508357.pdf.
Currie C, Van der Sluijs W, Whitehead R, Currie D, Rhodes G, Neville F, et al. HBSC 2014 survey in Scotland national report. Child and Adolescent Health Research Unit (CAHRU); 2015. http://www.cahru.org/content/03-publications/04-reports/hbsc_nr14_interactive_final.pdf.
Leyland AH, Dundas R, McLoone P, Boddy FA. Cause-specific inequalities in mortality in Scotland: two decades of change. A population-based study. BMC Public Health. 2007;7:172.
National Records of Scotland. Figures for drug-related deaths for Scotland for 1995 and earlier years. National Records of Scotland; 2016. https://www.nrscotland.gov.uk/files//statistics/drug-related-deaths/drd14/figures-for-drug-related-deaths-1995-and-earlier-years.pdf.
Carstairs V, Morris R. Deprivation. Explaining differences in mortality between Scotland and England and Wales. BMJ. 1989;299:886.
Minton J, Vanderbloemen L, Dorling D. Visualizing Europe’s demographic scars with coplots and contour plots. Int J Epidemiol. 2013;42(4):1164–76.
Yang Y, Schulhofer-Wohl S, Fu WJ, Land KC. The intrinsic estimator for age-period-cohort analysis: what it is and how to use it. Am J Sociol. 2008;113(6):1697–736.
Bell A, Jones K. The impossibility of separating age, period and cohort effects. Soc Sci Med. 2013;93:163–5.
Kramer MR, Casper M. Kramer and Casper respond to "A-P-C ... It's easy as 1-2-3!". Am J Epidemiol. 2015;182(4):318–9.
Kramer MR, Valderrama AL, Casper ML. Decomposing black-white disparities in heart disease mortality in the United States, 1973-2010: an age-period-cohort analysis. Am J Epidemiol. 2015;182(4):302–12.
Reither EN, Land KC, Jeon SY, Powers DA, Masters RK, Zheng H, et al. Clarifying hierachical age-period-cohort models: a rejoiner to bell and Jones. Soc Sci Med. 2015;145:125–8.
Yang Y, Land KC. Age-period-cohort analysis: new models, methods, and empirical applications: 1st ed. Boca Raton: CRC Press; 2013.
Reither EN, Masters RK, Yang YC, Powers DA, Zheng H, Land KC. Should age-period-cohort studies return to the methodologies of the 1970s? Soc Sci Med. 2015;128:356–65.
Land KC, Fu Q, Guo X, Jeon SY, Reither EN, Zang E. Playing with the rules and making misleading statements: a response to Luo, Hodges, Winship, and powers. Am J Sociol. 2016;122(3):962–73.
Keyes KM, Miech R. Age, period, and cohort effects in heavy episodic drinking in the US from 1985 to 2009. Drug Alcohol Depend. 2013;132:140–8.
Minton J, Shaw R, Green MA, Vanderbloemen L, Popham F, McCartney G. Visualising and quantifying 'excess deaths' in Scotland compared with the rest of the UK and the rest of Western Europe. J Epidemiol Community Health. 2017;71(5):461–7.
Case C, Deaton A. Mortality and morbidity in the 21st century. Brookings Papers on Economic Activity, Spring; 2017. p. 397–476. https://www.brookings.edu/wp-content/uploads/2017/08/casetextsp17bpea.pdf.
ScotPHO. Drug data: social harm and health harm. 2017; http://www.scotpho.org.uk/behaviour/drugs/data/health-harm (accessed March 2017).
The Adverse Childhood Experiences Study. http://acestudy.org/ (accessed March 2017).
Furlong A, Cartmel F, Biggart A. Choice biographies and transitional linearity: re-conceptualising modern youth transitions. Papers: Revista de Sociologia. 2006;79:225–39.
Furlong A, Cartmel F. Vulnerable Young Men in Fragile Labour Markets. Joseph Rowntree Foundation; 2004. https://www.jrf.org.uk/report/vulnerable-young-men-fragile-labour-markets.
Hills J, Cunliffe J, Gambaro L, Obolenskaya P. Who won and who lost as the economic crisis hit? Changing inequality in the UK, 2007–2010. Centre for Analysis of Social Exclusion; 2013. http://sticerd.lse.ac.uk/dps/case/spcc/SRR02.pdf.
Taulbut M, Hearty W, Myers F, Craig N, McCartney G. Pulling in different directions? The impact of economic recovery and continued changes to social security on health and health inequalities in Scotland. NHS Health Scotland; 2016. http://www.healthscotland.scot/media/1087/pulling-in-different-directions_march-2016.pdf.
We would like to thank Frank Dixon at NRS for the dataset and for advising on the best coding to use for obtaining long-term trends, and also David Readhead at NHS National Services Scotland for providing data by deprivation groupings.
JP and GM are employees of NHS Health Scotland which is a special NHS health board funded by the Scottish Government. JL and JB are employees of University of Glasgow. JM is funded by the ESRC (Code: ES/K006460/1).
Availability of data and materials
The authors have provided additional data in the Additional files 1, 2 and 3. The data used for this study are owned by NRS and were used under a confidentiality agreement to allow access to individual level data for the current study. Researchers can request the data from NRS.
Ethics approval and consent to participate
Not applicable. The data for the research was obtained from an existing database containing details of Scottish deaths from 1974 held by the National Records of Scotland. All requirements specified for the use of this dataset by the National Records of Scotland were adhered to.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Drug-related deaths in Scotland as determined by different definitions. Figure S2. Age and period ranges for birth cohorts available in the drug-related deaths dataset, 1979–2013. Figure S3. European age-standardised drug-related death rates per 100,000 population per year in Scotland, 1979–2013, by sex. Figure S4. European age-standardised drug-related death rates per 100,000 population per year in Scotland, 1979–2013, by sex and deprivation. Figure S5. Age distribution of crude drug-related death rates per 100,000 population per year in Scotland over periods 1979–1996 and 1997–2013 combined data by sex. Figure S6. Smoothed shaded contour plot of age-specific crude drug-related death rates per 100,000 population in Scotland for each single age from 18 to 86 years of age and each year from 1982 to 2010 stratified by sex. Figure S7. Age distribution of crude drug-related death rates per 100,000 population per year in Scotland over periods 1979–1996 and 1997–2013 combined data by sex and deprivation. Figure S8. Smoothed shaded contour plot of age-specific crude drug-related death rates per 100,000 population in Scotland for each single age from 18 to 81 years of age and each year from 1982 to 2010 stratified by sex and deprivation. Figure S9. Crude drug-related death rates per 100,000 population per year in Scotland by age (five-year age groups, for ages 15–89 years), period (five-year periods 1979–83 to 2009–13) and sex. (PDF 384 kb)
Regression models and interpretation of intrinsic estimator coefficients. Table S1. Goodness-of-fit statistics for single and two factor regression models and APC model. Table S2. Goodness-of-fit statistics for single and two factor regression models and APC model by Carstairs deprivation category. (PDF 305 kb)
Intrinsic estimator statistics. Table S1. Intrinsic estimator coefficient statistics for age, period and birth cohorts for drug-related deaths in Scotland for females. Table S2. Intrinsic estimator coefficient statistics for age, period and birth cohorts for drug-related deaths in Scotland for males. Figure S1. Intrinsic estimator coefficients with 95% confidence intervals for age, period and birth cohorts for drug-related deaths in Scotland stratified by sex. Table S3. Intrinsic estimator coefficient statistics for age, period and birth cohorts for drug-related deaths in Scotland for females living in the most deprived quintile. Table S4. Intrinsic estimator coefficient statistics for age, period and birth cohorts for drug-related deaths in Scotland for females living in the less deprived four quintiles. Table S5. Intrinsic estimator coefficient statistics for age, period and birth cohorts for drug-related deaths in Scotland for males living in the most deprived quintile. Table S6. Intrinsic estimator coefficient statistics for age, period and birth cohorts for drug-related deaths in Scotland for males living in the less deprived four quintiles. Figure S2. Intrinsic estimator coefficients with 95% confidence intervals for age, period and birth cohorts for drug-related deaths in Scotland for females stratified by deprivation. Figure S3. Intrinsic estimator coefficients with 95% confidence intervals for age, period and birth cohorts for drug-related deaths in Scotland for males stratified by deprivation. (PDF 1151 kb)