Introduction
Mental disorders in old age are associated with higher general health service use (Luppa et al. Reference Luppa, Sikorski, Motzek, Konnopka, Konig and Riedel-Heller2012). Yet, the use of mental health services is often poorly matched to needs and only one-third of all older adults with mental disorders report ever seeking help for their complaints such as psychiatric medication or psychotherapy (Wittchen & Pittrow, Reference Wittchen and Pittrow2002; Cepoiu et al. Reference Cepoiu, McCusker, Cole, Sewitch, Belzile and Ciampi2008; Busch et al. Reference Busch, Neuner, Aichberger, Hapke, Riedel-Heller and Luppa2013). Reasons for this misallocation of health care resources are unclear: under-diagnosis of mental health problems in old age is one likely barrier to adequate treatment (Mitchell et al. Reference Mitchell, Rao and Vaze2010), but negative attitudes towards psychiatric medication or psychotherapy could also play a role. Looking specifically at attitudes towards psychotherapy and pharmacotherapy, a study among depressed primary care patients in the USA indicated a preference for counselling or psychotherapy over antidepressant medication (Gum et al. Reference Gum, Arean, Hunkeler, Tang, Katon, Hitchcock, Steffens, Dickens and Unutzer2006), which corresponds to population studies that consistently show stronger preferences for psychotherapy over pharmacotherapy (Schomerus et al. Reference Schomerus, Schwahn, Holzinger, Corrigan, Grabe, Carta and Angermeyer2012). However, only very few older persons with depression receive psychotherapy, most patients of older age receiving treatment are taking antidepressant medication (Morichi et al. Reference Morichi, Dell'Aquila, Trotta, Belluigi, Lattanzio and Cherubini2015). On the one hand, this is simply because guidelines for geriatric depression recommend antidepressant medication as an initial treatment (Alexopoulos, Reference Alexopoulos2011), probably still reflecting the false belief that older persons do not benefit from psychotherapy. On the other hand, access to psychotherapy is limited (Unützer et al. Reference Unützer, Katon, Callahan, Williams, Hunkeler, Harpole, Hoffing, Della Penna, Noel, Lin, Tang and Oishi2003) and depression in older age is often accompanied by cognitive deficits, which might result in under-diagnosis of mood disorders (Morimoto & Alexopoulos, Reference Morimoto and Alexopoulos2013; Koethe et al. Reference Koethe, Mattern and Herpertz2014).
When looking at age-related differences in attitudes, three effects need to be distinguished. Age-effects refer to regular changes in attitudes over the lifespan, regardless of the year a person is born. Besides these changes over the lifespan, cohort effects might also cause variation in attitudes between older and younger persons. Cohort effects refer to attitudes that develop similarly in a certain birth cohort, for example, due to a shared cultural environment in childhood. The third time related effect on attitudes, changes that affect the general population regardless of age or birth cohort, is referred to as period effects. They reflect general attitude changes within a specific time frame, caused for example by economic or societal changes, but potentially also by population wide antistigma campaigns. Put briefly, age effects characterise the effect of ageing, while cohort effects characterise the effects of upbringing and period effects characterise societal changes.
So far, few studies examined attitudes of older persons towards mental health care and the existing cross-sectional studies based on population-based and clinical samples revealed inconsistent results (Leaf et al. Reference Leaf, Bruce, Tischler and Holzer1987; Kessing et al. Reference Kessing, Hansen, Demyttenaere and Bech2005). A study in Denmark found more negative attitudes towards antidepressant medication among older compared with younger adults being treated for affective disorders (Kessing et al. Reference Kessing, Hansen, Demyttenaere and Bech2005). A population-based study in the USA of attitudes towards mental health services in 1980 found that both younger and older participants are more sceptical towards mental health care seeking compared with those aged 25–64 years (Leaf et al. Reference Leaf, Bruce, Tischler and Holzer1987). In contrast, analyses of the population-based National Comorbidity Survey 1990–1992 in the USA found a positive correlation between age and willingness to seek mental health treatment (Mojtabai, Reference Mojtabai2007). So far, it is thus unclear whether attitudes towards using mental health care differ between younger and older persons.
Regarding period effects, attitudes of the general population towards psychiatric treatment have changed over the last decades. A meta-analysis of time-trend studies from several Western industrialised countries showed that in more recent surveys, more people recommended drug treatment of mental disorders than did so in earlier surveys (Schomerus et al. Reference Schomerus, Schwahn, Holzinger, Corrigan, Grabe, Carta and Angermeyer2012). Accordingly, prescription rates particularly for antidepressant medication are rising. Between 1999 and 2010, antidepressant use in the USA increased from 6 to 10% of the population (Olfson & Marcus, Reference Olfson and Marcus2009; Mojtabai & Olfson, Reference Mojtabai and Olfson2014). Similarly, in OECD countries, consumption of antidepressants rose from 3 to 6 defined daily doses (DDD) per 100 inhabitants between 2000 and 2011 (OECD. Health at a Glance 2013: http://www.oecd.org/health/healthtaglance (accessed 26.5.2014)). A time trend study over 21 years from Germany showed that between 1990 and 2011, the approval rates of psychotherapy for depression rose from 57 to 71%, and for schizophrenia from 62 to 82% (Angermeyer et al. Reference Angermeyer, Matschinger and Schomerus2013a , Reference Angermeyer, Matschinger and Schomerus b ). A Swedish study recently reported an overall increase in more positive attitudes toward psychotherapy and medication between 1976 and 2014 (Mirnezami et al. Reference Mirnezami, Jacobsson and Edin-Liljegren2015). While suggesting a strong period-effect, these studies did not investigate simultaneously the effect of age, cohort and time period on treatment approval. Hence, it is unclear to what extent age of the respondents, changes in birth cohorts, or general societal changes are responsible for the development of preferences for psychotropic medications or psychotherapy over time.
Aims of the study
In this paper, we examine age related differences in attitudes towards mental health treatment by disentangling age, cohort and period effects. In particular, we use age, period and cohort analyses based on three pooled cross-sectional studies of the German population to examine the approval of either psychiatric medication or psychotherapy for a person with either depression or schizophrenia.
Methods
Surveys
Our study is based on data from three methodologically identical population-based cross-sectional surveys among German citizens aged 18 years and older conducted in 1990, 2001 and 2011 (Angermeyer & Matschinger, Reference Angermeyer and Matschinger1997; Angermeyer et al. Reference Angermeyer, Beck and Matschinger2003, Reference Angermeyer, Matschinger and Schomerus2013a , Reference Angermeyer, Matschinger and Schomerus b ; Schomerus et al. Reference Schomerus, Van der Auwera, Matschinger, Baumeister and Angermeyer2015). Aim of these surveys was to assess various mental health-related attitudes of the general population cross-sectionally, and to enable time-trend analyses of attitudes. In all surveys, the samples were drawn using an identical random sampling procedure with the three stages sample points, households and individuals within the target households. Target households within the sample points were determined according to the random route procedure (Gabler & Hoffmeyer-Zlotnik, Reference Gabler and Hoffmeyer-Zlotnik1997). Target individuals within households were selected using random digits. Informed consent was considered to have been given when individuals agreed to complete the interview. Fieldwork was carried out by GETAS (Hamburg) in 1990, and by USUMA (Berlin) in 2001 and 2011. Data in 1990 were gathered in two waves, wave 1 in the old German states (April 1990) and wave 2 in the new German states (November 1990) shortly after re-unification.
Interview
All interviews were conducted face-to-face using pencil and paper. In all surveys, the fully structured interview was identical regarding wording and the sequence of questions. At the beginning of the interview respondents were presented with a vignette of a diagnostically unlabelled case history depicting a person with mental illness. While several mental disorders had been studied in single surveys, vignettes on schizophrenia and major depression had been used in three consecutive surveys and could thus be subject of our age-period-cohort (APC) analysis. The symptoms described in both vignettes fulfilled the criteria of DSM-IV for schizophrenia or depression and had been validated by five experts in psychopathology blinded to the diagnosis (vignettes see supplement). While in 1990 and 2011 the sex of the individual presented in the vignettes varied at random, in 2001 only responses to a male vignette were elicited. For the present study, we thus only use responses elicited with the male vignette. We asked respondents treatment recommendation for the person described in the vignette. Presenting them with a catalogue of treatment options, we asked: ‘Although you might not be expert, you may have an idea what should best be recommended to someone with a problem described in the vignette. Please tell me for each option to what extent you would recommend it or would not recommend it.’ Eliciting treatment recommendations for a third person is an established method to assess mental health literacy (Reavley & Jorm, Reference Reavley and Jorm2012). Recently it has been shown that treatment recommendations and personal treatment preferences are closely related (Rüsch et al. Reference Rüsch, Müller, Ajdacic-Gross, Rodgers, Corrigan and Rössler2014). Answers were recorded using a five-point Likert scale with the anchors 1 = ‘I would strongly recommend it’ and 5 = ‘I would strongly advise against it’, plus a ‘don't know’ category. For our present analysis, we use answers to the items ‘using psychotherapy’ and ‘taking psychotropic drugs’.
Sample description
Our combined study sample from all surveys consisted of 9046 participants, 4505 subjects receiving the depression vignette and 4541 participants receiving the schizophrenia vignette. Since in both samples from 1990 the population from the old (West-)German states was underrepresented, we weighted this sample to achieve a proportion of 80:20 for old and new German states for the schizophrenia and depression sample, respectively. This is representative of the overall German population and identical to the old/new-proportions in the samples from 2001 and 2011. Participants with missing values in one of the analysed variables were excluded. To aid interpretation of results, we also excluded participants answering with ‘don't know’ (vignette: depression + recommendation: psychotherapy N = 356; vignette: depression + recommendation: psychotropic medication N = 443; vignette: schizophrenia + recommendation: psychotherapy N = 299; vignette: schizophrenia + recommendation: psychotropic medication N = 457). Our final study sample thus included N = 4120 (psychotherapy) and N = 4034 (psychotropic medication) for the depression vignette and N = 4212 (psychotherapy) and N = 4061 (psychotropic medication) for the schizophrenia vignette.
Statistical analysis
Disentangling age, cohort and period effects is methodologically challenging because of the perfect collinearity of the three time variables (period = cohort + age). To extract the ‘pure age related effect’ that is not confounded by cohort or period effects, more advanced statistical methods than simple ordinary least-squares regression are required. For each of the two disorders, we performed an APC analysis to test and separate associations of age at examination, birth cohort and time period with our two items of treatment recommendation (psychotherapy and psychotropic medication). This collinearity problem is referred to as the identification problem because only two of the three variables can be put into a regression model simultaneously. Instead, we used partial least-squares regressions (PLS), which allow using perfectly collinear variables in one regression model. In PLS components are extracted that are weighted combinations of the original variables. These components are chosen according to their relationship to the outcome variable. They seek to maximise the covariance between the outcome and those weighted components. As a result, the first component has a greater covariance with the outcome than the second and so on. In a three-dimensional space with three covariates as input variables a maximum of three PLS components can be extracted which are orthogonal to each other. The regression coefficients for the three variables can directly be obtained from their weights in each component (Tu et al. Reference Tu, Krämer and Lee2012; Jiang et al. Reference Jiang, Gilthorpe, Shiely, Harrington, Perry, Kelleher and Tu2013). For the mathematical equations see supplemental material. Bastien et al. (Reference Bastien, Esposito-Vinzi and Tenenhaus2005) extended the linear PLS to generalised linear models (GLM), including binary and ordinal regression. We used the R-package plsRglm that is available on the CRAN archive and maintained by Frédéric Bertrand. First, we used a GLM PLS with logit link and binomial distribution to model binary treatment recommendation for psychotherapy and psychotropic medication, respectively. We dichotomised the five-point Likert scales (1: ‘I would strongly recommend’ & ‘I would recommend’ 0: ‘I would neither recommend or advise against it’ & ‘I would advise against it’ & ‘I would strongly advise against it’). We included age at interview and birth cohort as continuous variables and two dummy variables for time period (i.e. 2011 and 2001 v. 1990) as covariates. We performed 1000 bootstrap replications to obtain 95% confidence intervals (CIs). To additionally examine non-linear relationships for age and birth cohort, we used restricted cubic splines (Harrell, Reference Harrell2001) in PLS regression with five knots for age at examination and birth cohort. Wald tests of coefficients of the second, third and fourth spline transformation equal to zero were used to obtain p values for nonlinearity. Wald tests indicated that age and cohort were linearly related to all outcomes. To assess the model fit of linear and non-linear PLS, we compared the pseudo R 2, which is a measure of the variance of the outcome variable that can be explained by the covariates, analogous to ordinary least-squares regression. Because PLS penalises variables according to their variance when extracting the components, covariates were scaled to have unit variance and zero mean. To validate our results, we conducted two additional sensitivity analyses detailed in the supplement. First, we included those participants that had answered with the ‘don't know’ category by coding them as denial of treatment (i.e. ‘I would strongly advise against it’) in the binary logistic regression. Second, we performed an ordered logistic regression included in the plsRglm package with the five response categories of the Likert scale (setting respondents with ‘don't know’ to missing). To test the proportional odds assumption underlying the ordered logit model, we used a graphical method proposed by Harrell (Reference Harrell2001) because other statistical tests have the tendency to reject the null hypothesis. Analyses were performed using the statistical software R (version 3.0.2).
Results
In the pooled cross-sectional data age at examination ranged from 18 to 93 and birth cohort from 1900 to 1993. Table 1 shows the number of individuals in each age group and stratified by type of vignette for 1990, 2001 and 2011 (see Supplementary Table S1 for the number of individuals per birth cohort). For more details on the distribution of age and birth cohort with respect to their recommendation see Supplementary Fig. S1. The mean age was 47 for both samples and mean birth cohort was 1952. The depression and schizophrenia samples had 44 and 46% males, respectively.
Table 1. Number of subjects of each age group at each time period
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170110052235-11746-mediumThumb-S2045796015001134_tab1.jpg?pub-status=live)
Model fit
In all PLS analyses, a two-component model showed better model-fit than a one-component model and a similar model fit to the three-component model (see Supplementary Table S2). We thus report results for the two-components logistic PLS models only (Table 2). In depression, the two components PLS model explained 4.07% of the variance (pseudo R 2) for recommending psychotherapy and 3.63% for psychotropic medication (one-component model: psychotherapy 3.74%; psychotropic medication 2.98%). In total, 69.4% of all participants would recommend psychotherapy and 36.3% psychotropic medication in case of depression irrespective of age, birth cohort or period. For schizophrenia, the model with two PLS components explained 4.85% of the variance for the treatment recommendation psychotherapy and 9.23% for psychotropic medication (one-component model: psychotherapy 4.36%; psychotropic medication 8.24%). Overall 81.3% recommend psychotherapy and 50.4% psychotropic medication as treatment for schizophrenia.
Table 2. Treatment recommendations towards a person with schizophrenia or depression: age, period and cohort effects.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170110052235-54744-mediumThumb-S2045796015001134_tab2.jpg?pub-status=live)
Logistic PLS regression with two components. Unstandardised coefficients.
a CI after 1000 non-parametric Bootstrap replicates.
Age effects
For the depression vignette, age was significantly associated with an increase of both treatment recommendations with similar effect sizes (psychotherapy OR per decade = 1.09; psychotropic medication OR per decades = 1.10) independent from cohort and period effects (Table 2). Over the lifespan of 60 years the odds for recommending psychotherapy increases to OR = 1.69 and to OR = 1.82 for psychotropic medication (Fig. 1). In schizophrenia, both recommendations also showed a significant age effect (Table 2). Likewise, we observed an OR of 1.06 per decade for the recommendation of psychotherapy and even higher, OR = 1.16 per decade, for recommending psychotropic medication. Over the lifespan (60 years) this cumulates to an increase of ORs to 1.44 and 2.55 for recommending psychotherapy and psychotropic medication, respectively (Fig. 1).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20170110052235-43455-mediumThumb-S2045796015001134_fig1g.jpg?pub-status=live)
Fig. 1. Cumulative effect of age, cohort and period on treatment recommendations for schizophrenia and depression.
Cohort effects
In depression, an effect of birth cohort on treatment recommendation could only be observed for psychotherapy (OR = 1.07 per decade) with greater approval in persons born in more recent birth cohorts. Over a time period from cohorts 1930 to 1990 this increases OR to 1.53.
For schizophrenia, both treatment recommendations also showed a significant cohort effect with similar effect sizes (ORs per decade of 1.08 and 1.06 for psychotherapy and medication, respectively). Participants born in later years indicated stronger recommendation for both treatment options, irrespective of age and period effects. Comparing early and most recent birth cohorts over 60 years (1930–1990), this effect cumulated to an OR = 1.62 and OR = 1.46 for psychotherapy and psychotropic medication, respectively (Fig. 1).
Taken together, age and cohort effects indicate higher attitudes both over the lifespan and in coming years, when more recent cohorts grow older.
Period effects
For time period, we calculated the effects for the three distinct time points. In depression, time period had a strong effect on treatment recommendations. Participants from more recent surveys indicated higher approval for the treatment options than the reference sample from 1990. The OR for choosing psychotherapy was 2.361 in 2001 and 2.25 in 2011, when 1990 served as reference category. For psychotropic medication, ORs were 2.43 in 2001 and 1.52 in 2011 (with 1990 as reference category). Using 2001 as reference category yielded no significant change between 2001 and 2011, neither for psychotherapy nor for psychotropic medication.
Similar to depression a significant period effect could be observed in schizophrenia for both treatment recommendations, particularly for psychotropic medication. People from the more recent surveys were more likely to recommend both treatment options. The OR for recommending psychotherapy were 2.97 in 2001 and 2.36 in 2011, compared with 1990. The decline between 2001 and 2011 was not significant when 2001 served as the reference group. For psychotropic medication the OR in 2001 was 3.59 and 3.34 in 2011, which is a relatively stable approval for this treatment between 2001 and 2011. Again, the decrease in treatment recommendation between 2001 and 2011 was not statistically significant.
Using the full sample by including the ‘don't know category’ as potential denial of treatment did not change the overall associations (Supplementary Table S3). Similarly, ordinal logistic regression on the five-point Likert scale showed overall similar results to those presented in the main analysis (Supplementary Table S4) indicating a robust model.
In summary, we found improving attitudes over time. Recommendation increases with growing age, in more recent birth cohorts, and generally improved over the last 20 years (Fig. 1). Differences in the APC pattern for treatment recommendations between schizophrenia and depression were primarily observed for recommending psychotropic medication.
Discussion
The first positive conclusion we can draw from our analyses is that no treatment recommendation showed a significant decrease irrespective of vignette, growing age, birth cohort or time period. The by far strongest effects were observed for time period, especially between 1990 and 2001. Obviously there has been a general rise in acceptance for professional mental health treatment, a trend that has been described and discussed previously (Angermeyer & Matschinger, Reference Angermeyer and Matschinger2005; Mojtabai, Reference Mojtabai2007). In our study, this trend peaked in 2001, with no significant further increase in treatment recommendations between 2001 and 2011. Disentangling this trend from any age and cohort effects, our study corroborates the notion that psychiatric treatment has generally gained acceptance over the last 20 years (Angermeyer et al. Reference Angermeyer, Matschinger and Schomerus2013b ). These changes have been attributed to changes within psychiatric services following psychiatric reform, to increasing mental health literacy of the general population, but also to a more biological understanding of mental disorders (Angermeyer & Matschinger, Reference Angermeyer and Matschinger2005; Mojtabai, Reference Mojtabai2007; Angermeyer et al. Reference Angermeyer, Matschinger and Schomerus2013b ). Since our study is the first to look at age and cohort effects in treatment recommendations for mental disorders, we will focus our discussion on this novel part of our results. Our findings show that there is no reason to suspect that particularly psychotherapy receives less support among older compared with younger persons. Especially for depression, psychotherapy as well as psychotropic medication show a similar increase in popularity with growing age. In our sample, 71.5% of the participants aged 70 years and above recom mended psychotherapy and 47.1% recommended medication. In schizophrenia, the recommendation for psychotropic medication showed a much stronger increase than psychotherapy but both treatment options are highly recommended in old age (psychotherapy 78.6% and psychotropic medication 60.2% in people of age ≥70 in our sample). The higher recommendation rates in schizophrenia could reflect the greater perceived severity of schizophrenia as depicted in our vignettes, highlighting a stronger necessity for treatment. Accordingly, the percentage of the subjects in the survey 2011 considering the person described as suffering from a ‘mental illness in a medical sense’ was 83% in schizophrenia compared with 65% in depression (Schomerus et al. Reference Schomerus, Matschinger and Angermeyer2013). If schizophrenia is more frequently seen as a medical problem, a higher rate of professional treatment recommendations seems plausible.
Depression among older persons is very common, affecting between 10 and 15% of older persons living in the community (Beekman et al. Reference Beekman, Copeland and Prince1999). Facing the challenge of an ageing society and a rising number of geriatric patients in the future, adequate treatment options are required. Although antidepressant use especially in the elderly has increased (Olfson & Marcus, Reference Olfson and Marcus2009), only one-third of older patients achieve remission (McGovern et al. Reference McGovern, Kiosses, Raue, Wilkins and Alexopoulos2014). Pharmacotherapy with antidepressants in older persons is also associated with higher all-cause mortality (Coupland et al. Reference Coupland, Dhiman, Morriss, Arthur, Barton and Hippisley-Cox2011). Older persons are more likely to exhibit adverse drug effects due to multi-morbidity and the high number of prescribed medications (Dines et al. Reference Dines, Hu and Sajatovic2014). Moreover, meta-analyses show similar effect sizes for pharmacotherapy and psychotherapy in older persons with depression (Pinquart et al. Reference Pinquart, Duberstein and Lyness2007). Provision of psychotherapy should thus complement pharmacotherapy or could even be a safer alternative to medication. It is becoming increasingly recognised that older people represent a suitable target group for psychotherapy due to their better emotional regulation skills, increased wisdom and resilience mechanisms accumulated over many years of experience (Laidlaw, Reference Laidlaw2013). Our findings support the willingness of this patient group to undertake psychotherapy in case of depression.
The particular merits of disentangling age, cohort and period effects are illustrated by the fact that both older age and a more recent birth cohort are associated with stronger approval of mental health treatment. This finding might level differences between older and younger respondents in cross-sectional assessments, but our analyses indicate that prospectively, the popularity of professional treatment among older persons will further increase when younger cohorts grow older. In depression, this is particularly true for psychotherapy, where younger cohorts that are even more open minded towards psychotherapy will become the new generation of older people. The strong period effects found in this study, yielding continuously higher recommendation rates for psychotherapy in more recent surveys compared with 1990, also supports this trend. Conversely, attitudes towards pharmacotherapy of depression seem to have peaked: they do not improve in younger cohorts, and between 2001 and 2011 there seems to be an overall (although insignificant) decline in approval. Generally, the strong period effects found in this study show that attitudes towards mental health treatment are changeable, which is an encouraging finding for future awareness and prevention programs.
Treatment recommendations for a person with schizophrenia also increased with age. Regarding specific treatment options for older people, this seems less important. First, approval rates for both medication and psychotherapy are generally high in schizophrenia, so age-related effects are less relevant. Second, schizophrenia mainly starts at younger age and is thus of less personal relevance for older persons. Nevertheless, our concordant findings on schizophrenia and depression suggest that greater openness particularly towards psychotherapy is not restricted to a specific diagnosis, but is reflecting a general change of mental health related attitudes. This is especially visible in the strong period effect of increased overall approval of treatment options between 1990 and 2001.
Our study has several limitations that need to be considered. First, we elicited treatment recommendations for another person, and not personal treatment preferences. However, there is evidence that recommendations for others and personal preferences are closely related (Rüsch et al. Reference Rüsch, Müller, Ajdacic-Gross, Rodgers, Corrigan and Rössler2014). Examining recommendations is also relevant in itself, because advice from a friend or family member is one of the main sources of information in case of help seeking for psychiatric disorders (Rüsch et al. Reference Rüsch, Müller, Ajdacic-Gross, Rodgers, Corrigan and Rössler2014). But it must also be noted that at the end the recommendation of the doctor has a great influence on treatment decision, too. Thus, it may be interesting to analyse the changes in treatment preferences of primary care physicians.
Other factors besides the three-time variables examined in this study might determine recommendations regarding mental treatment options. These may be sex, education or prior contact to people with mental disorders. This study aimed to examine and compare age relative to cohort and period effects. Future studies will help to examine whether this effects are attributable to other potential confounding, effect-modifying factors and mediators.
In summary, our study shows that both medication and psychotherapy are becoming more popular over the lifespan. Increasing the provision of psychotherapy for older persons in need would be perfectly in line with treatment preferences elicited among the general population.
Supplementary materials
The supplementary materials referred to in this article can be found at http://dx.doi.org/10.1017/S2045796015001134