Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Testing Persistence of Cohort Effects in the Epidemiology of Suicide: an Age-Period-Cohort Hysteresis Model

  • Louis Chauvel ,

    louis.chauvel@uni.lu

    Affiliation Institute for Research on Socio-Economic Inequality (IRSEI), University of Luxembourg, Esch-sur-Alzette, Luxembourg

  • Anja K. Leist,

    Affiliation Institute for Research on Socio-Economic Inequality (IRSEI), University of Luxembourg, Esch-sur-Alzette, Luxembourg

  • Valentina Ponomarenko

    Affiliation Institute for Research on Socio-Economic Inequality (IRSEI), University of Luxembourg, Esch-sur-Alzette, Luxembourg

Abstract

Birth cohort effects in suicide rates are well established, but to date there is no methodological approach or framework to test the temporal stability of these effects. We use the APC-Detrended (APCD) model to robustly estimate intensity of cohort effects identifying non-linear trends (or ‘detrended’ fluctuations) in suicide rates. The new APC-Hysteresis (APCH) model tests temporal stability of cohort effects. Analysing suicide rates in 25 WHO countries (periods 1970–74 to 2005–09; ages 20–24 to 70–79) with the APCD method, we find that country-specific birth cohort membership plays an important role in suicide rates. Among 25 countries, we detect 12 nations that show deep contrasts among cohort-specific suicide rates including Italy, Australia and the United States. The APCH method shows that cohort fluctuations are not stable across the life course but decline in Spain, France and Australia, whereas they remain stable in Italy, the United Kingdom and the Netherlands. We discuss the Spanish case with elevated suicide mortality of cohorts born 1965–1975 which declines with age, and the opposite case of the United States, where the identified cohort effects of those born around 1960 increase smoothly, but statistically significant across the life course.

Introduction

A growing body of research investigates birth cohort effects in suicide rates [16]. While there is large international evidence that certain cohorts experience a higher risk of suicide compared to other birth cohorts, there is no knowledge if this higher risk is stable across the cohorts’ life course or if the risk that elevated suicide rates are capturing is only temporary. Further, no method is able to resolve the linear age-period-cohort relation. Even modern statistical models [7, 8] suffer from intrinsic problems to disentangle linear age, period, cohort (APC) effects [912]. In order to estimate the contribution of birth cohort effects apart from period and age effects in suicide rates, it is thus necessary to investigate the non-linearities of cohort effects (colloquially speaking the ‘cohort bumps’) as stable long–term deviations from a linear trend of social change, similar to the ‘zero linear trend’ (ZLT) [9]. We will investigate these non-linearities of cohort effects with a recently introduced APC-Detrended (APCD) approach [13]. APCD is useful to detect, test and internationally compare cohort contrasts where more gifted birth cohorts are systematically above the inter age-period-cohort trend compared with birth cohorts below the trend (see further elaboration in the methodological section). Further, with our advancement of APCD, APC-Hysteresis (APCH) that will be introduced here, we will estimate if these non-linear cohort effects are temporally stable or follow differential trajectories.

The need to disentangle age, period and cohort effects on social, economic and health outcomes has led to greater focus on age-period-cohort models. Cohort membership matters indeed in explaining different social phenomena, as shown for mortality rates [14], political behaviour [15], income position [13], educational attainment [16], suicide [17, 18] and other health outcomes [19]. Increases in cohort mortality by suicide have been shown for cohorts born after World War II in Britain [1], the United States [2] and Spain [3]. Similar results were found for pre-war cohorts in Australia [4], Brazil [5] and Quebec [6]. However, these models lack ability to control for confounders or draw causal conclusions.

Methodological developments to overcome these issues included, for instance, multi-country comparisons [20] and mixed models with random cohort effects that confirmed cohort effects in suicide and homicide rates [21], and more specifically, that youth suicide rates of larger cohorts and those with higher non-marital birth rates are larger [22, 23]. However, those approaches were not able to distinguish whether observed cohort characteristics indeed produced the cohort effect or simply mirror an unobserved cohort characteristic. By applying age-period-cohort models to overcome limitations of age-period models, cohort effects for suicide have been investigated in several studies [24, 25], but could not test the stability of these effects.

Well-debated in the scientific community, for instance with regard to the APC-Intrinsic estimator (APC-IE) [8], is the inherent difficulty to extract linear patterns of age, period and cohort effects [911]. Other problems involve choice of model constraints, hierarchies or substitutions to solve the identification problem that may seem arbitrary or convenience-driven [10, 11, 13]. If the identification of age, period, cohort non-linear fluctuations are methodologically feasible [26], linear trends generally do not contain meaningful information. Consequently, the recently established APCD method [13] focuses on non-linear fluctuations, i.e., ‘detrended’ cohort effects above and below the linear trend. It enables not only the inclusion of independent variables, but also post-estimation statistics. To sum up, the first aim of the present study is to robustly estimate non-linear cohort effects determining intensity of cohort contrasts, as APCD provides stable validated results independent of subjective choice of constraints. Here, after providing large-scale comparative multi-country evidence on variations in suicide rates, the second aim is to analyse the temporal stability of the non-linear effects. Therefore, we differentiate four potential trajectories.

In the sociological tradition of investigating suicide not as individual behaviour, but rather via collective experience and behaviour, we study the development of suicide risk by cohorts over their life course. Durkheim hypothesised anomic suicide as result of lack of social integration and social regulation [27]. Merton referred to antisocial behaviour due to relative frustration of the perceived divergence of aspirations and reality [28]. From this perspective, a collective social trauma due to exposure to unemployment, economic downturn, or poverty may lead to long-term increases in suicide risk. However, those assumptions in APC analyses are hardly testable. What we can test though is the stability of these effects.

Most prominent among possible trajectories is stability or ‘hysteresis’, well-illustrated by the scarring effect of early unemployment generating negative long-term outcomes to its victims [29, 30]. Indeed, Clark and colleagues used the term hysteresis to describe the process of long-term habituation to an initially negative condition, even if that condition may change or even vanish with time. Hysteresis can in the same sense be used to describe cohort processes of habituating to initially negative conditions, i.e. holding behaviour constant over time, even under changing period influences.

Persistence of cohort-specific risks, such as unemployment and lower chances of future employment of affected cohorts, known as scarring effect [29, 31], has up to now rather been hypothesised than clearly demonstrated and we test stability of scarring effects against other plausible trajectories of cohort effects. The second suggested trajectory is that of increasing cohort contrasts across the life course drawing from the theory of ‘cumulative disadvantage’ [32], assuming that earlier disadvantage prevents future access to certain resources or positions leading to an even greater deprivation in future life, the so-called ‘Matthew effect’ [33]. This increase in contrasting cohort trends results in ever increasing inequalities across the life course. The third trajectory is ‘compensation’, where initially high levels of suicide are compensated and diminish over the life course [34]. Fourth, early disadvantage could theoretically even lead to later advantage over other birth cohorts in the sense of a complete reversal or ‘inversion’ [35]. Bell [36] presented evidence of this pattern of inversion for subjective well-being in a British sample. In this paper, we intend to investigate these four differential trajectories by analysing (extreme) health risks of birth cohorts.

Materials and Methods

Data

We use data on suicide mortality from the WHO mortality database (http://www.who.int/healthinfo/mortality_data/en) for years 1970 to 2009 and for 25 countries. No ethical consent for the study was necessary since the data were analysed anonymously (mortality by cause data). From the WHO data sets, cause of death, country, gender, and population size of the country were retained for the male population. All considered data of the mortality database code causes of mortality with ICD 8 to ICD 10 (although classifications of causes of death were restructured through the ICD versions, classification of suicide was consistent through the versions). With age at death i and year of death j, we construct birth cohort k. We build eight five-year periods and five-year age groups from 20–24 to 75–79. We provide the Stata do file (S1 Dofile), dataset (S1 Dataset), and codebook (S1 Codebook) to replicate the analyses presented here.

Strategy of data analysis

Relying on the tradition of age-period-cohort models, a dependent variable is explained by combination of age a (variable αa), cohort membership c (variable γc) and period of measurement p (variable πp). This leads to the equation:

Before elaborating on APC models in more detail, three hypothetical configurations of suicide rates across time and cohorts illustrate the decomposition of age, period, cohort effects. Table 1 shows a data structure with linear trends which are impossible to attribute to period or cohort. To this data structure in Table 2 a non-linear cohort effect is added, which can be detected by APCD. Finally, data structure of Table 3 reflects temporal change in the cohort effect, showing a cohort effect at the beginning of the observed time period which vanishes across the life course, detectable by APCH. We first elaborate on linear trends shown in Table 1. While many papers proposed empirical estimations of the three effects [26, 37], an “identification problem” besets all APC models [38, 39]: an infinite number of solutions can fit the same outcome with equal efficiency, resulting from the relation a = p—c. That is, each variable is a combination of the other two. While in Table 1 the increase in suicide rates is obvious, the different possible effects cannot be disentangled, as the increases can either be understood as combination of age×period (age effect: suicide rate increases by 5 units per 5-year age group; period effect: suicide increases by 1 unit per 5-year period) or combination of age×cohort (age effect: again an increase of 5 units per 5-year age group; cohort effect: 1 unit more for each cohort born five years later), or one of the infinite more complex APC combinations. Despite some attempts to identify these linear trends [19, 40], no statistical model can overcome this intrinsic indetermination [10, 41, 42].

thumbnail
Table 1. Hypothetical data structure of suicide rates across time and cohorts: APC linear cohort effect.

Note. a/y—age groups versus period.

https://doi.org/10.1371/journal.pone.0158538.t001

thumbnail
Table 2. Hypothetical data structure of suicide rates across time and cohorts: Non-linear cohort effect with persistent intensity.

Note. a/y—age groups versus period.

https://doi.org/10.1371/journal.pone.0158538.t002

thumbnail
Table 3. Hypothetical data structure of suicide rates across time and cohorts: Temporal change in intensity of cohort effect.

Note. a/y—age groups versus period.

https://doi.org/10.1371/journal.pone.0158538.t003

The second table (Table 2) superimposes a non-linear sinusoidal birth-cohort transformation. Those non-linear fluctuations had been constructed as waves that strictly follow a birth-cohort structure: some cohorts face higher risks at any age (e.g., cohorts born in 1945 and in 1965). APCD is able to extract cohort effect coefficients representing fluctuations above and below linear combination of age and period effects.

Table 3 illustrates hypothetical temporal change in the cohort effect by combining non-linear deviations in young cohorts (Table 2) and a linear trend in more senior age (Table 1), meaning that an initially observed cohort effect vanishes across the life course. In many empirical cases, e.g., for unemployment, unstable non-linearities can be observed when fluctuations are strong for young and weak for seniors. Initial cohort effects are thus completely absorbed over the analyzed lifespan. APCH is able to extract cohort fluctuations like APCD but extends the analysis by assessing the temporal stability or instability of the cohort effect with, but it also measures, with a hysteresis coefficient h (see section APCH and ‘hysteresis’ coefficient).

The methodology improvements here include three steps. The first is a Poisson specification of APCD originally built for ordinary least squares (OLS) [13] or logit specifications [43]. Second, a ‘cohortality coefficient’ is developed to express the degree to which birth cohort is a relevant expression of age×period. Third, APCH is built to assess cohort effect persistence.

APCD in Poisson specification

First, a former OLS-type APCD model [13] is transformed into Poisson specification. For each country, a dependent variable rapc denotes suicide rate (ratio between counts of suicides dapc in reference population of size napc) in age group a of period p belonging to cohort c = pa. Categorical time effect variables pertaining to age effects αa, period effects πp and cohort effects γc, are then indexed by age a, period p and cohort c. In more detail, a, p and c are entangled in a two-dimensional space so two trend variables (rescale a and c) are added to absorb these linear trends. Three vectors representing age, period, cohort are constrained to sum = zero and slope = zero to reflect non-linearities. The first (k = 1) and the last (k = a + p—1) cohorts are omitted as they appear just once and their introduction substantially increases standard errors. The resulting model with only one solution provides fits relatively close to those of the APC-Intrinsic Estimators of Yang et al. [8] in many different contexts, except that the cohort effect in APCD shows a slope zero by construction. APCD is expressed in the following:

(APCD) (1) where αa is the age effect vector indexed by age group a, πp is the period vector and γc is the cohort vector, β0 the constant. To provide a correct identification [14, 44], simple constraints such as Σaαa = Σpπp = Σcγc = 0 imply centred coefficients. We add three constraints, slopea(αa) = slopep(πp) = slopec(γc) = 0, where Slope is the linear function that gives the linear slope of the coefficients, so that αa, πp and γc are detrended. The terms α0 Rescale(a) and γ0 Rescale(c) absorb the two linear trends of APC. Rescale is a transformation that standardizes the coefficients α0 and γ0: it transforms age from initial amin to amax to the interval −1 to +1. The constraint Slopeaa) = 0 means the age coefficients αa, are trend-zero and is true only if Σa [(2aaminamaxa] = 0. This constraint is easily expressed as linear equation of coefficients. Coefficients πp and γc are built analogously. This model is thus a traditional Mason and Smith [45] APC model, following Holford’s notion [41] that cohort is detrended as a constraint imposes zero slope on A, P and C effects.

If cohort effects do not exist, i.e. cohorts do not deviate from age and period characteristics, detrended cohort effect coefficients γc are zero: APCD provides no improvement compared with a simple age×period model (AP) with first and last cohorts omitted:

(AP) (2)

However, if at least one γc coefficient is significantly different from zero (i.e., a significant non-linear effect deviating from the linear trend absorbed by AP), then AP is insufficient. Thus, to retain appropriate parsimonious models, comparing the Bayesian Information Criterion of (AP) and (APCD) is a diagnosis of relevance of non-linear cohort effects.

Under these specifications, APCD is a non-linear cohort effect detector of suicide, and γc is the detrended cohort effect (DCE) which assesses cohort non-linearities.

‘Cohortality’ coefficient

The cohortality coefficient (CC) assesses the degree to which the DCEs absorb (or do not) the full age×period interaction. CC is then a measure of relevance of cohort analysis. In an OLS specification, a natural index of appropriateness of cohort analysis consists of comparing (APC) and (AP), where age and period effects are summed. The full interaction (A×P) is also called “saturated model”, for example in sociology of mobility [46]. The gap between (AP) and (APC) is the addition of the diagonal interaction following birth cohorts; the one between (APC) and (A×P) is that the latter fits all cells of the interaction.

Then: (3.1) where R2 is the coefficient of the determination of the model. Otherwise, CC is the percentage of residuals of the simple model (AP) explained by cohort effects. In a Poisson specification, the equivalent of R2 is the residual deviance statistics [47], so CC can be expressed as (3.2) where: (3.3)

Again, in this Poisson compatible specification, CC measures how a cohort accurately summarizes the interaction between A and P. When CC is close to one, the cohort term perfectly fits the interaction. If CC is closer to zero, birth cohort is not relevant in explaining the phenomenon.

APCH and ‘hysteresis’ coefficient

Next, the APCH model is built, which assesses persistence of cohort effects over the life span. APCH is designed to assess the degree of stability of cohort effects, and thus the ‘scarring’ of cohorts, by assessing the interaction between cohort fluctuations and age. Interested readers are referred to the original ado commands, as both APCD and APCH have been developed as STATA ado commands that can be downloaded with the usual ‘ssc install apcd’ and ‘ssc install apch’ orders. Because they are based on STATA glm, they can handle OLS, logit and Poisson log-linear specifications.

The principle of APCH is to test the degree of stability across the (observed) life course of detected significantly nonzero.γc The solution is to estimate hysteresis coefficient h to differentiate trajectories via APCH:

(APCH) (4)

APCH is an advancement of APCD where the simple γc term is replaced by a variable term (1+h.rescale(a)) γc. Negative values of h represent diminishing cohort fluctuations detected by γc: the ‘cohort effect’ is changing across the life course. There, γc captures the middle-age value of the cohort effect, somehow an average cohort effect that culminates when rescale(a) is minimum (rescale(a) = −1 at the beginning of the observed life span), and the cohort effect diminishes (when h < 0) at the end of the observed life span (when rescale(a) = +1). More precisely, APCH offers calibration of h where, if h = −1, cohort effect γc is an average between its maximum at entry and a null effect at the end of the observed period. The detected cohort non-linearity can thus either relate to stable or instable trajectories, which correspond to the spectrum of possible h values as follows:

  • h > 0 represents increasing cohort effects across the life course (cumulative disadvantage);
  • h = 0 is cohort effect without specific trend and APCH is APCD (hysteresis/stability);
  • −1 < h < 0 describes compensation in that cohort fluctuations diminish over the life span; h can be understood as percentage measurement of the decline from rescale(a) = −1 at entry to rescale(a) = +1 at the end of observed lifespan (compensation);
  • h ≤ −1 is cohort inversion.

Because of the non-linear type of interaction involved here, estimation of γc (the DCE coefficients) and of the h coefficient cannot be simultaneous. The strategy here is to build an iterative process of estimation where at step (n), we alternate a first substep (n,1) with estimation of the h of step (n) as an interaction between rescaled age and DCE of step (n−1) that had been previously estimated. At second substep (n,2) where the DCE of step (n) is now estimated based on the estimate of h of substep (n,1) and the [h.rescale(a) γc] term estimated at step (n−1). At step (1), h(1) is estimated based on the results DCE(0) of the initial APCD.

APCH (n,1) will estimate the hysteresis coefficient .

(5)

APCH (n,2) will estimate the DCE(n) (6)

This iterative process of estimation converges rapidly towards a unique result, i.e., a DCE and a hysteresis coefficient h. Testing the models with hypothetical data structures of Tables 1 to 3 provide respective coefficients in Table 4.

thumbnail
Table 4. Goodness-of-fit and relevance of cohort analysis for hypothetical data structures in Tables 1 to 3.

https://doi.org/10.1371/journal.pone.0158538.t004

Resuming the data structures of Tables 1 to 3, Table 1 showed linear trends in suicide rates impossible to decompose, Table 2 an additional non-linear cohort effect detectable by APCD, and Table 3 temporal change in the cohort effect, i.e., a cohort effect at the beginning of the observed time period which vanishes across time, detectable by APCH. As depicted in Table 4, for linear trends (Table 1), CC = zero reflects irrelevance of cohort analysis; comparison of BICs between AP and APC show APC fits better. APCH does not improve APC. As DCE are not relevant in this example, no h coefficient is provided. For persistent intensity of cohort effect (Table 2), CC = 1 and the cohort fluctuation perfectly represents the a×p interaction. Comparing BICAP to BICAPC shows that cohort is relevant but APCH is not further improving fit of the model, i.e. temporal change of the trend is no issue here: The value of h is not significantly different to zero: the DCE is stable across the life span. For temporal change in intensity of cohort effect (Table 3), CC = 0.9 reflects temporal change in cohort structure, confirmed by BIC comparison with best fit of APCH. The h coefficient is almost -1, expressing the non-linear cohort effect vanishing across the life course, reflecting a pattern close to the earlier, theoretically assumed inversion across life span. On other sets of simulations (http://www.louischauvel.org/apchex), these methodological choices return appropriate values of h and DCE. Standard error of h and BIC differences of the models can assess if APCD is preferred over AP, and if APCH is preferred over APCD.

Results

Descriptive step

A descriptive analysis of the interaction between age and period reveals relevant changes in the age profile of suicide rates. Whereas suicide for men in old age was typically very high, the age ratio is changing towards younger males. Fig 1 displays the changing age distribution of suicide in 1985 and 2005 for each country across periods. In 2005, compared with rates in 1985, suicide increased for younger ages in many countries. In Australia, Belgium, Canada, France, Italy, the Netherlands, Spain, the United States and United Kingdom, suicide in old age is less common in 2005 than in 1985. The opposite holds for youth and middle-aged populations where an age transformation in suicide is observed, i.e., suicide rates vary for different ages.

thumbnail
Fig 1. Transformation of the age profile of logged suicide rate in 1985 (dashed line) and 2005 (solid line).

https://doi.org/10.1371/journal.pone.0158538.g001

To better understand interactions age×period in each country, residuals of the (AP) model show higher or lower suicide risks for each cohort (Fig 2), as opposed to the hypothesised null interaction of country-specific age×period (i.e., the age profile does not change over periods). Indeed, lines of the different age groups would be completely flat if no cohort effects existed, but we detect strong fluctuations which indicate cohort effects (Fig 2). For example, in Spain, old-age suicide decreased for cohorts born after 1920, but suicide at age 30 and 40 increased considerably for cohorts born after 1940 (see Fig 3 for an illustration with a hypothetical and the empirical dataset). The same development can be observed in Australia, Belgium, France, Ireland, New Zealand, the United Kingdom, and the United States (see Fig 4). These Figures simply reflect age and period descriptives, not a result based on age-period-cohort models. Fig 2 shows that specific cohorts are systematically above or below zero, showing specific cohort-based risks below and above the trend across the life course.

thumbnail
Fig 2. Residuals of logged suicide rates after control of (AP) by country.

https://doi.org/10.1371/journal.pone.0158538.g002

thumbnail
Fig 3. Hypothetical versus empirical suicide fluctuations: The Spanish case.

Note. We compare three datasets: The actual empirical dataset a1 (suicide rates per 100,000 population), a simulated dataset without changes in intensity of cohort effects (i.e. hysteresis/stability over the life course) a2, and the difference between cells a3. Now consider b1 and b2, the residuals of the logged values of datasets a1 and a2 in the model (AP) (we suppress the linear combination of age and period dummies). b1 shows some cohort effects, that are ‘imperfect’ compared to b2 the expected values under the hypothesis of a pure APC without h effect. In these tables, the red cells are the high values of suicide and the blue cells the relatively protected groups. b3 presents the differences between empirical b1 and theoretical b2, and following the diagonally changing colours, we can detect temporal stability/changes of cohort fluctuations.

https://doi.org/10.1371/journal.pone.0158538.g003

thumbnail
Fig 4. Hypothetical versus empirical suicide fluctuations: The U.S. case.

Note. See legend of Fig 3 for explanations, for reasons of space only b1-b3 are presented.

https://doi.org/10.1371/journal.pone.0158538.g004

APCD: Examining intensity of cohort effects

After detecting cohort differences, their shape and intensity are displayed with APCD as detrended APC coefficients in Fig 5 with 95% confidence intervals and zero mean. Horizontal lines indicate very small detrended cohort effects. Germany, Austria and countries in Northern and Eastern Europe are characterized by less intensive non-linear cohort trends. For instance, for Austria, the age profile does not change and thus the DCEs do not show a cohort bump. Conversely, several countries show similar configurations with strong downward oscillations for cohorts born in 1940 to 1950: strong cohort fluctuations are found for Australia, Belgium, France, Italy, the Netherlands, New Zealand, Spain and the United Kingdom. In these countries cohorts born around 1940 are ‘luckier’ cohorts that are less likely to commit suicide when age, period and overall cohort linear trend effects are controlled for. By directly comparing age- and period-dependent patterns in suicide rates of different countries, particularly important fluctuations are detected in Belgium, New Zealand and the United Kingdom. They show massive variations with +0.4 to −0.4 standard deviations in Belgium for cohorts born between 1900 and 1950. This means that relative suicide rates vary by a magnitude of 1 to 2 among these cohorts at the same age. For countries with flat horizontal lines cohort analysis is not relevant as most of variation in suicide rates is due to age and period effects.

thumbnail
Fig 5. Cohort APCD coefficients: structure of cohort ‘bumps’ in male suicide.

Note. 95% confidence intervals in grey.

https://doi.org/10.1371/journal.pone.0158538.g005

To conclude, strong cohort contrasts exist in Australia, Belgium, France, Italy, the Netherlands, New Zealand, Spain and the United Kingdom. In order to examine stability of these cohort effects, cohort difference intensity is measured by the standard deviation of the DCE coefficients, thus assessing ‘cohort contrasts’ in suicide rates and thus intensity of disadvantage or advantage of one birth cohort compared to others (that we will cross with the CC).

Cohortality coefficient (CC): Identifying countries with consistent cohort effects

After identifying shape and significance of cohort effects (Fig 5), the CC tool detects cases where cohort is a consistent explanation of suicide rate, i.e., the extent to which the age-period interaction is due to cohort effects (Fig 6). Although no explicit rule has been introduced for a cut-off of the CC coefficient, based on this sample of countries we choose 0.8 (see also the gap in Fig 6 between CCs per country). This cut-off helps to run analyses parsimoniously by pre-selecting only those constellations where cohort fluctuations are present and can be tested for temporal stability (compare with Fig 2). CC is above 0.8 for the 12 upper cases, with Poland at the lower end, scoring just above 0.8, and the United Kingdom has the strongest CC at 0.93. Countries with CC below 0.8 show either no cohort effect (Austria, Germany), more complicated interactions of age, period, cohort in the sense of cohort fluctuations that are influenced by simultaneously occurring period or age effects, with Japanese cohorts born in the 1940s starting with high suicide rates becoming better protected and converse patterns for those born in the 1920s (see Fig 2), or shaky variations due to smaller national populations (Israel, Luxembourg).

thumbnail
Fig 6. Cohortality coefficients of the 25 countries. The 12 countries with higher CC are retained for further analysis of hysteresis.

https://doi.org/10.1371/journal.pone.0158538.g006

APCH: Detecting temporal stability of cohort effects

The CC detects relevant cohort effects in 12 countries: Poland, New Zealand, Hungary, the United States, Canada, Australia, the Netherlands, Spain, Italy, Belgium, France, and the United Kingdom. However, high CC does not inform about the stability of the cohort effect. APCH detects degree of stability of cohort disadvantage and as such the development of cohort contrasts across time. Fig 7 show standard deviations of the DCEs plotted against hysteresis parameters of those countries with relevant cohort effects.

thumbnail
Fig 7. Standard deviation of the detrended cohort effects (SDDCE) and h coefficients in the 12 countries with higher Cohortality Coefficient.

https://doi.org/10.1371/journal.pone.0158538.g007

In Fig 7, standard deviation of the DCE coefficients and h plot shows the difference between countries with weak cohort effects (left) and stronger cohort effects (right). In case of the −1 < h < 0, the cohort effects decrease across the life course (lower side). In the case of h > 0, cohort effects increase, so the effect is cumulative (upper side). The spikes indicate the 95% confidence intervals on both variables. Below the baseline are countries where cohort effects decline across the life course questioning the implicit hypothesis of stability of cohort effects. Considering standard errors of estimates of h, the hysteresis/stability hypothesis is acceptable in only six of the 12 cases: Poland, Canada, Hungary, Italy, the Netherlands, and the United Kingdom.

Our analysis provides evidence for cumulative disadvantage in the United States as cohort fluctuations are increasing. On the contrary, the cohort effect of Spain can be interpreted as strong cohort compensation, and initially unlucky cohorts born 1965–1975 could partially recover with time. If h is significantly different from zero disadvantaged cohorts indeed remain disadvantaged however with only partial stability. For instance, in Belgium only 15% (h = -.15+/-.11) of disadvantage is absorbed across the life course. Additionally, one-quarter to one-third of cohort effects weaken in Australia, New Zealand and France. Therefore, rather than merely assuming persistence of cohort effects, the observed cases show the importance of considering the value of h to find evidence regarding actual stability of cohort contrasts.

Testing this APCH model against a random effects model as suggested by one reviewer, we arrive at similar results for a simulated data set (see syntax at http://www.louischauvel.org/simu_hapc_apch.do). For the empirical dataset used here, APCH provides the further advantage of providing robust results without requiring manual fine-tuning, arriving at stable and ready-to-use estimations of hysteresis (stability) coefficients of cohort effects via an automatized iterative process. Random effects models may be useful when year and cohort measured variables are of interest [49]. Bayesian estimation may solve convergence issues for random effects models.

Discussion

The present study introduced a technique to disentangle non-linear cohort effects from age and period effects and to test their persistence across the life course. The APCD and the newly developed APCH method focus on non-linear effects net of any trends and provide evidence that cohort effects in suicide rates are relevant for some of the observed countries of which several are stable across the life course. Stable non-linear cohort effects are found for the United Kingdom, Italy and the Netherlands, whereas France and Spain show cohort fluctuations with decreasing intensity across the life course. Cohort disadvantages are significantly stable in six of 12 countries, whereas in the United States, cohort contrasts are smoothly but significantly increasing across life span.

Explanation of findings

In the first part of this paper, intensity of cohort disadvantages changing across the life course was assessed by the h coefficient to distinguish trajectories of initial cohort disadvantages across the life course. Despite commonly hypothesised stability of cohort disadvantage, our analyses showed that among suicide rates of 25 countries, only six cases showed stable cohort fluctuations (Poland, Italy, the Netherlands, Canada, United Kingdom and Hungary). We can identify other configurations which prove the usefulness of our method of testing stability of cohort effects: cumulative disadvantage (the United States) and compensation (Spain).

In order to explain the strong suicide fluctuations for cohorts born 1965–1975 in Spain, with elevated levels in suicide mortality in the 1980s but those effects diminishing over time, we need to consider the specific economic and social context of Spain in the 1980s. Increases in suicide mortality during this time have been long documented [3]. The overall loss of life expectancy in young adults during the 1980s in Spain, rather than increase as would be expected, has been mainly attributed as consequence of HIV/AIDS, drugs and violent deaths, particularly traffic accidents [50]. This has been explained with very particular social and economic developments at that time in Spain: Whereas political and structural crises led to breaking with old norms and values, at the same time unemployment and labour market reforms prohibited the young from realizing their aspirations, which presumably led them to adopting unhealthy lifestyle behaviours such as alcohol and drug use, and driving under the influence of alcohol [50]. Unemployment rate has indeed been linked to elevated suicide mortality for Galicia, Spain [51]. This socioeconomic context coupled with wide availability of drugs and repeated use of injection needles led to the HIV outbreak, resulting in strong increase in AIDS mortality. Deaths classified as suicide at that time could also be attributable to the diagnosis and bad prognosis of AIDS [52]. The total loss of life expectancy for men and women due to HIV/AIDS, traffic accidents and suicide between 1980 and 1990 has been estimated to be 0.8 and 0.2 years, respectively, in Catalonia [53]. The negative h coefficient found for Spain mirrors the efforts of public health management to absorb those shocks specific to these generations after 1990, and afterwards became a public health issue that was not specific to those generations anymore.

The case of the United States [2, 18] is a specific example but it echoes a cohort structure relatively common in other Western countries; U.S. cohort fluctuations are significant but relatively modest compared to eight other nations. When the linear trend is taken into account, the luckier cohorts were less likely to commit suicide, with the luckiest cohort born 1920–1950. Hence, the ‘luckiest’ cohort members being those born around World War II or immediately after. As one reviewer points out, the positive h coefficient for the U.S. could have also been a chance finding due to multiple testing. However, we have replicated the U.S. findings with another data source: Merging Centers for Disease Control and Prevention (CDC) Multiple Cause of Death files matched with Current Population Survey (CPS) data of 1990–2010, we find the exactly same increase in mid-age mortality by suicide of the 1960-born cohort, more specifically those low educated, white, and not married [54]. We have extended those calculations by running the APC-H method on those new data as well and, mirroring the presented findings, found again a positive h coefficient in the CDC-CPS merged file. Although we cannot make use of further explanatory variables, it is likely that cohorts born 1920 to 1950 experienced post-war economic boom and context of full employment specific of affluent societies [55] in their youth and middle age. The strong positive h coefficient of the U.S. can be explained with rising socioeconomic stressors across the life course of the affected cohorts: The problems the baby boomers have been facing for some time were strongly increasing in the first decade of the 21st century. In comparison, cohorts born around 1960 could have been affected by economic slowdown of the 1980s. Those cohorts born around 1960, more specifically the non-Hispanic white male members of these cohorts, have recently been identified as especially vulnerable to morbidity and mortality [56], a finding which has been validated with age-period-cohort methodology for suicide, especially for low-educated and un-married members of the cohort [54]. Referring to theory of relative frustration, cohorts born around 1960 experienced a different economic situation in middle age than they expected. Where cohorts born 1900 to 1920 might have experienced material deprivation and scarcity of resources, cohorts born 1950 to 1975 could have perceived relative deprivation in comparison with their parents’ situation at their age. Higher suicide rates could be an extreme outcome of this divergence.

Limitations and suggestions for further research

The analyses here use the WHO mortality database with a sample of 25 countries for a period of almost 40 years (1970–2009). Data on further countries and longer time frames (for some countries since 1950, newer data up to 2013 at the moment of writing) are available in the database. Hence it would be desirable if further research extends the presented analyses and tests presence and temporal stability of cohort effects for a larger sample of countries and even longer time spans. One should also note that with longer time frames under investigations, cohort dynamics may appear different overall because extreme conditions in the newly considered time interval (before or after the initially considered interval) may change the long-term cohort trajectory in some cases. Furthermore, additional methodological elaborations with regard to the different interactions of age, period and cohort may be considered relevant depending on the conceptual framework.

In some of the investigated countries, birth cohorts showed initial disadvantage but partial or even full recovery of disadvantage across the life course (compensation), whereas others remained disadvantaged (stability). However, searching for reasons for differential disadvantage and differential trajectories was beyond the scope of the paper. In further studies, it would be desirable to investigate causes for temporary or stable disadvantage. Future research should include both micro-level information on education and other socioeconomic indicators to investigate how birth cohort characteristics may ‘shape their disadvantage’, see a review of papers related to social stratification and suicide [57]. As the present study only included men it is necessary to investigate possible gender effects and possible interactions of age, period, cohort effects with gender. One advancement would be to provide hysteresis coefficients for each cohort separately, when this paper only considered the hysteresis of the overall society under investigation. Another further development of our method could be to disentangle slowly evolving cohort fluctuations (e.g. New Zealand, Fig 2) versus abruptly occurring cohort fluctuations (e.g. French baby boomers). Rapidly significantly fluctuating cohort contrasts means a strong differentiation of close cohorts. Conversely, slow U-shaped cohort fluctuations might hide more complicated evolutions, and be less diagnostic in identifying disadvantaged cohorts. The test of temporal stability of those cohort fluctuations however stays the same.

Conclusions

Investigating temporal stability of suicide fluctuations in a set of 25 countries and between 1970 and 2009, we detect 12 nations where cohort fluctuations are unstable, in particular in Spain elevated suicide mortality of cohorts born 1965–75 declines with age, and in the U.S. cohort effects of those born around 1960 and protected at earlier ages increase across the life course.

Further research should include more information on historical economic situation to investigate contextual causes for temporary or stable disadvantage: Now that we can differentiate persistent and temporary cohort specificities in suicide rates, we will be able to study better the links between socio-economic context of socialization and temporal change of cohort fluctuations in suicidal behaviour.

Supporting Information

S1 Codebook. Codebook for WHO suicide mortality data.

https://doi.org/10.1371/journal.pone.0158538.s001

(DOCX)

S1 Dataset. Stata dataset with suicide rates of WHO mortality database 1970–2009 for 25 countries.

https://doi.org/10.1371/journal.pone.0158538.s002

(DTA)

S1 Dofile. Stata do file to replicate APCD and APCH models and figures.

https://doi.org/10.1371/journal.pone.0158538.s003

(DO)

Acknowledgments

We would like to thank Andrew Bell, Herbert L. Smith, and Vladeta Ajdacic-Gross for helpful suggestions on an earlier version of the manuscript.

Author Contributions

Conceived and designed the experiments: LC. Performed the experiments: LC VP. Analyzed the data: LC VP AL. Wrote the paper: AL LC VP. Data interpretation: LC VP AL.

References

  1. 1. Gunnell D, Middleton N, Whitley E, Dorling D, Frankel S. Influence of cohort effects on patterns of suicide in England and Wales, 1950–1999. The British Journal of Psychiatry. 2003;182(2):164–70.
  2. 2. Phillips JA, Robin AV, Nugent CN, Idler EL. Understanding recent changes in suicide rates among the middle-aged: period or cohort effects? Public Health Reports. 2010;125(5):680–88. pmid:20873284
  3. 3. Granizo JJ, Guallar E, Rodriguez-Artalejo F. Age-period-cohort analysis of suicide mortality rates in Spain, 1959–1991. International Journal of Epidemiology. 1996;25(4):814–20. pmid:8921461
  4. 4. Snowdon J, Hunt GE. Age, period and cohort effects on suicide rates in Australia, 1919–1999. Acta Psychiatrica Scandinavica. 2002;105(4):265–70. pmid:11942930
  5. 5. Rodrigues NC, Werneck GL. Age-period-cohort analysis of suicide rates in Rio de Janeiro, Brazil, 1979–1998. Social Psychiatry and Psychiatric Epidemiology. 2005;40(3):192–6. pmid:15742223
  6. 6. Légaré G, Hamel D. An age-period-cohort approach to analyzing trends in suicide in Quebec between 1950 and 2009. Can J Public Health. 2013;104(2):e118–e23. pmid:23618203
  7. 7. Yang Y, Fu WJ, Land KC. A methodological comparison of age‐period‐cohort models: The Intrinsic Estimator and conventional generalized linear models. Sociological methodology. 2004;34(1):75–110.
  8. 8. 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. American Journal of Sociology. 2008;113(6):1697–736.
  9. 9. O’Brien RM. Constrained estimators and age-period-cohort models. Sociological Methods & Research. 2011;40(3):419–52.
  10. 10. Luo L. Assessing validity and application scope of the intrinsic estimator approach to the age-period-cohort problem. Demography. 2013;50(6):1945–67. pmid:24072610
  11. 11. Bell A, Jones K. The impossibility of separating age, period and cohort effects. Social Science & Medicine. 2013;93:163–5.
  12. 12. Fienberg SE. Cohort analysis’ unholy quest: A discussion. Demography. 2013;50(6):1981–4. pmid:24132742
  13. 13. Chauvel L, Schröder M. Generational inequalities and welfare regimes. Social Forces. 2014:sot156.
  14. 14. Wilmoth JR. Les modèles âge-période-cohorte en démographie [Age-period-cohort models in demography] In: Caselli G, Vallin J, Wunsch G, editors. Démographie: Analyse et Synthèse I La Dynamique des Populations. Paris: INED; 2001. p. 379–97.
  15. 15. Chauvel L, Smits F. The endless Baby Boomer Generation: Cohort differences in participation in political discussions in nine European countries in the period 1976–2008. European Societies. 2015;17(2):242–78.
  16. 16. Chauvel L. Overeducation and social generations in France. In: Attewell P, Newman KS, editors. Growing gaps Educational inequality around the world. Oxford: Oxford University Press; 2010. p. 210–38.
  17. 17. Chauvel L. L'uniformisation du taux de suicide masculin selon l'âge: effet de génération ou recomposition du cycle de vie? [The standardization of the male suicide rate by age: effect of generation or reconstruction of the life cycle?] Revue française de sociologie. 1997:681–734.
  18. 18. Phillips JA. A changing epidemiology of suicide? The influence of birth cohorts on suicide rates in the United States. Social Science & Medicine. 2014;114:151–60.
  19. 19. Reither EN, Hauser RM, Yang Y. Do birth cohorts matter? Age-period-cohort analyses of the obesity epidemic in the United States. Social Science & Medicine. 2009;69(10):1439–48.
  20. 20. Stockard J, O'Brien RM. Cohort variations in suicide rates among families of nations. An analysis of cohorts born from 1875 through 1985. International Journal of Comparative Sociology. 2006;47(1):5–33.
  21. 21. O'Brien RM, Stockard J. A common explanation for the changing age distributions of suicide and homicide in the United States, 1930 to 2000. Social Forces. 2006;84(3):1539–57.
  22. 22. Stockard J, O'Brien RM. Cohort effects on suicide rates: International variations. American Sociological Review. 2002:854–72.
  23. 23. Stockard J, O'Brien RM. Cohort variations and changes in age-specific suicide rates over time: Explaining variations in youth suicide. Social Forces. 2002;81(2):605–42.
  24. 24. Allebeck P, Brandt L, Nordstrom P, Åsgård U. Are suicide trends among the young reversing? Age, period and cohort analyses of suicide rates in Sweden. Acta Psychiatrica Scandinavica. 1996;93(1):43–8. pmid:8919328
  25. 25. Ajdacic–Gross V, Bopp M, Gostynski M, Lauber C, Gutzwiller F, Rössler W. Age–period–cohort analysis of Swiss suicide data, 1881–2000. European Archives of Psychiatry and Clinical Neuroscience. 2006;256(4):207–14. pmid:16283596
  26. 26. Fienberg SE. Identification and estimation of age-period-cohort models in the analysis of discrete archival data. Sociological Methodology. 1979:1–67.
  27. 27. Durkheim E. Le suicide: étude de sociologie [Suicide. A sociological study] Paris: F. Alcan; 1897.
  28. 28. Merton RK. Social structure and anomie. American Sociological Review. 1938;3(5):672–82.
  29. 29. Ellwood DT. Teenage unemployment: Permanent scars or temporary blemishes? In: Freeman RB, Wise DA, editors. The youth labor market problem: Its nature, causes, and consequences: University of Chicago Press; 1982. p. 349–90.
  30. 30. Clark A, Georgellis Y, Sanfey P. Scarring: The psychological impact of past unemployment. Economica. 2001;68(270):221–41.
  31. 31. Gangl M. Welfare states and the scar effects of unemployment: A comparative analysis of the United States and West Germany. American Journal of Sociology. 2004;109(6):1319–64.
  32. 32. Dannefer D. Cumulative advantage/disadvantage and the life course: Cross-fertilizing age and social science theory. The Journals of Gerontology Series B: Psychological Sciences and Social Sciences. 2003;58(6):S327–S37.
  33. 33. DiPrete TA, Eirich GM. Cumulative advantage as a mechanism for inequality: A review of theoretical and empirical developments. Annual Review of Sociology. 2006:271–97.
  34. 34. Burton-Jeangros C, Widmer ED. Cumulative and compensatory effects over the life course. Schweizerische Zeitschrift für Soziologie. 2009;35(2):183–92.
  35. 35. Hobcraft J, Menken J, Preston S. Age, period, and cohort effects in demography: a review. In: Mason WM, Fienberg SE, editors. Cohort analysis in social research: Springer; 1985.
  36. 36. Bell A. Life-course and cohort trajectories of mental health in the UK, 1991–2008–A multilevel age–period–cohort analysis. Social Science & Medicine. 2014;120:21–30.
  37. 37. Mason KO, Mason WM, Winsborough HH, Poole WK. Some methodological issues in cohort analysis of archival data. American Sociological Review. 1973:242–58.
  38. 38. Glenn ND. Cohort analysts' futile quest: Statistical attempts to separate age, period and cohort effects. American Sociological Review. 1976;41(5):900–4.
  39. 39. Mason WM, Wolfinger NH. Cohort analysis. International Encyclopedia of the Social & Behavioral Sciences. Oxford: Pergamon; 2001. p. 2189–94.
  40. 40. Yang Y, Land KC. Age-period-cohort analysis. New models, methods, and empirical applications. Chapman, Hall, editors. Boka Raton, FL: CRC Press, Taylor & Francis Group; 2013.
  41. 41. Holford TR. Understanding the effects of age, period, and cohort on incidence and mortality rates. Annual Review of Public Health. 1991;12(1):425–57. pmid:2049144.
  42. 42. Tu Y-K, Chien K-L, Burley V, Gilthorpe MS. Unravelling the effects of age, period and cohort on metabolic syndrome components in a Taiwanese population using partial least squares regression. BMC Medical Research Methodology. 2011;11(1):1.
  43. 43. Chauvel L. The intensity and shape of inequality: The ABG method of distributional analysis. Review of Income and Wealth. 2014.
  44. 44. Caselli G, Vallin J, Wunsch G. Démographie: analyse et synthèse, tome 1 [Demography: Synthesis and analysis. I. Population dynamics] La dynamique des populations, Paris, INED. 2001.
  45. 45. Mason WM, Smith HL. Age-period-cohort analysis and the study of deaths from pulmonary tuberculosis: Springer; 1985.
  46. 46. Raftery AE. Social mobility measures for cross-national comparisons. Quality & Quantity. 1985;19(2):167–82.
  47. 47. Hill J, Gelman A. Data analysis using regression and multilevel/hierarchical models: Cambridge University Press; 2007.
  48. 48. Raftery AE. Bayesian model selection in social research. Sociological Methodology. 1995;25:111–64.
  49. 49. Bell A, Jones K. Explaining fixed effects: Random effects modeling of time-series cross-sectional and panel data. Political Science Research and Methods. 2015;3(01):133–53.
  50. 50. Gómez-Redondo R, Boe C. Decomposition analysis of Spanish life expectancy at birth: evolution and changes in the components by sex and age. Demographic Research. 2005;13:521–46.
  51. 51. Fernández-Navarro P, Barrigón M, Lopez-Castroman J, Sanchez-Alonso M, Páramo M, Serrano M, et al. Suicide mortality trends in Galicia, Spain and their relationship with economic indicators. Epidemiology and Psychiatric Sciences. 2015:1–10.
  52. 52. Valdes B, George K. Demographic analysis of AIDS mortality in Spain. Population. 2013;68(3):473–85.
  53. 53. Spijker J, Llorens AB. Mortality in Catalonia in the context of the third, fourth and future phases of the epidemiological transition theory. Demographic Research. 2009;20:129–68.
  54. 54. Chauvel L, Leist AK, Smith HL. Cohort factors impinging on suicide rates in the United States, 1990–2010. Annual Meeting of the Population Association of America, March 31—April 2, 2016, Washington, DC Full paper available at http://orbiluunilu/handle/10993/25339 [Internet]. 2016.
  55. 55. Galbraith JK, Crook A. The affluent society: Houghton Mifflin Boston; 1958.
  56. 56. Case A, Deaton A. Rising morbidity and mortality in midlife among white non-Hispanic Americans in the 21st century. Proceedings of the National Academy of Sciences. 2015;112(49):15078–83.
  57. 57. Wray M, Colen C, Pescosolido B. The sociology of suicide. Annual Review of Sociology. 2011;37:505–28.