Introduction
Subsyndromal depression (SSD) or subthreshold depression is a common affective condition that has received increasing attention in the literature during the last decade. The concept of SSD can be described as depressive symptoms below the threshold of what is called a syndromal or a major depressive episode (Meeks et al., Reference Meeks, Vahia, Lavretsky, Kulkarni and Jeste2011). Compared to syndromal depression, SSD is generally associated with lower degrees of impairment and suffering (Ayuso-Mateos et al., Reference Ayuso-Mateos, Nuevo, Verdes, Naidoo and Chatterji2010), but the relevance of SSD is substantial because of a higher prevalence in the population (Meeks et al., Reference Meeks, Vahia, Lavretsky, Kulkarni and Jeste2011), and consequently, for example, the population-attributable fraction of mortality is comparable between the two (Cuijpers et al., Reference Cuijpers, Vogelzangs, Twisk, Kleiboer, Li and Penninx2013).
As affective conditions vary greatly in intensity over time, longitudinal studies of depressive symptoms have certain advantages for their understanding. SSD compared with non-depression (ND) in old individuals is longitudinally associated with impaired activities of daily living (ADLs; Barry et al., Reference Barry, Allore, Bruce and Gill2009) and cognitive function (Han et al., Reference Han, McCusker, Cole, Abrahamowicz and Capek2008), reduced self-perceived health (Callahan et al., Reference Callahan, Hui, Nienaber, Musick and Tierney1994), worse psychiatric outcomes with higher risk of major depression (Beekman et al., Reference Beekman2002), and higher mortality (Cuijpers et al., Reference Cuijpers, Vogelzangs, Twisk, Kleiboer, Li and Penninx2013). However, few studies have investigated longitudinal outcomes in very old individuals (age 80+). Very old individuals differ in general from the young old (age 60–80) with respect to the degree of multimorbidity, frailty, and social dependence.
SSD in the very old has been investigated with respect to natural history (Scafato et al., Reference Scafato2012), psychiatric comorbidity (Laborde-Lahoz et al., Reference Laborde-Lahoz, El-Gabalawy, Kinley, Kirwin, Sareen and Pietrzak2014), growth trajectories and cohort variations (Yang, Reference Yang2007), cognitive functions (Han et al., Reference Han, McCusker, Cole, Abrahamowicz and Capek2008), and mortality (Rapp et al., Reference Rapp, Gerstorf, Helmchen and Smith2008; Stek et al., Reference Stek, Vinkers, Gussekloo, Beekman, van der Mast and Westendorp2005). However, to our knowledge only two papers have described outcomes of ADLs (Bacon et al., Reference Bacon, Heeren, Keysor, Stuver, Cauley and Fredman2016; Hajek et al., Reference Hajek2017), and none have investigated the development of loneliness in very old individuals with SSD. Moreover, the previous studies on mortality in very old individuals with SSD and depression have shown slightly different results. In the Leiden 85+ study (Stek et al., Reference Stek, Vinkers, Gussekloo, Beekman, van der Mast and Westendorp2005), SSD and depression had an effect on mortality only due to loneliness. In the BASE study (Rapp et al., Reference Rapp, Gerstorf, Helmchen and Smith2008), the effect of SSD and depression on mortality was significant for young old but not for very old individuals, when adjusting for sex and comorbidity.
The aim of the present study was to investigate longitudinal outcomes of SSD in very old individuals. Based on previous studies on SSD in the young old and on syndromal depression in very old age, we hypothesized that in very old individuals, SSD, compared to ND, would be associated with lower function of ADLs (Hajek et al., Reference Hajek2017), worse self-perceived health (Callahan et al., Reference Callahan, Hui, Nienaber, Musick and Tierney1994), lower cognitive functions (Han et al., Reference Han, McCusker, Cole, Abrahamowicz and Capek2008), a higher degree of loneliness (Houtjes et al., Reference Houtjes, van Meijel, van de Ven, Deeg, van Tilburg and Beekman2014), depressive symptoms (Beekman et al., Reference Beekman2002), and higher mortality (Cuijpers et al., Reference Cuijpers, Vogelzangs, Twisk, Kleiboer, Li and Penninx2013).
Methods
Participants and procedure
Data were collected from the prospective observational cohort study “Elderly in Linköping Screening Assessment” (ELSA 85), with a population-based design covering 85-year-old individuals in the south of Sweden (Nagga et al., Reference Nagga, Dong, Marcusson, Skoglund and Wressle2012). All individuals born in 1922 and living in Linköping municipality in Sweden (n = 650) in 2007 were invited to participate (Figure 1), and no exclusion criteria were used.
Of the 650 eligible individuals, 496 (76%) provided written informed consent and answered a postal questionnaire (T0). The next phase (T1) of baseline measures included a home visit by an occupational therapist who interviewed 371 participants (75% of 496) and performed assessments using the depression screening instrument, Geriatric Depression Scale (GDS-15; Sheikh and Yeasavage, Reference Sheikh, Yeasavage and Brink1986). The baseline phase also included a visit to the local memory clinic where history taking, a medical examination, and administration of cognitive tests occurred. Three additional measure waves were completed: after 1 year (T2), after 5 years (T3), and after 8 years at age 93 (T4), at which questionnaires including the GDS-15 and other assessments were readministered. In parallel, information about mortality status and date of death for deceased participants was obtained from the national tax agency from baseline and up to 9 years later, when only n = 133 (36%) of the T1 sample were still alive. The present study was approved by the Research Ethics Committee of Linköping University, Sweden (141-06, 2012/332-31 and 2014/455-31).
Measures
Categories of depressive symptoms: Definitions of SSD vary in the literature (Cuijpers et al., Reference Cuijpers, Vogelzangs, Twisk, Kleiboer, Li and Penninx2013; Meeks et al., Reference Meeks, Vahia, Lavretsky, Kulkarni and Jeste2011). In this study, we chose to define SSD and syndromal depression on the basis of the results from GDS-15 (Sheikh and Yeasavage, Reference Sheikh, Yeasavage and Brink1986). GDS-15 is a valid and reliable self-report screening instrument assessing depressive symptoms in older people, and it contains 15 questions with yes/no answers, where 1 point is awarded for each question (de Craen et al., Reference de Craen, Heeren and Gussekloo2003). In this study, 0–2p were interpreted as having ND, 3–5p were interpreted as having SSD, and >5p were interpreted as having syndromal depression (Stek et al., Reference Stek, Vinkers, Gussekloo, van der Mast, Beekman and Westendorp2006).
Sociodemographic factors: The questionnaire included questions about sex and number of years of education.
Functional abilities: The questionnaires included the self-report instrumental activity measure (Grimby et al., Reference Grimby, Andren, Daving and Wright1998) assessing instrumental ADLs (I-ADLs) through perceived difficulty in the performance of eight different activities on a four-level scale (too difficult (1), great difficulties, some difficulties, and no difficulties (4)): locomotion outdoors, preparing a simple meal, cooking, using public transportation, small-scale shopping, large-scale shopping, cleaning, and washing. These eight values were converted into a single summary score (8–32p) of I-ADLs. Basic ADLs were assessed through questions about the perceived need for assistance in four different activities (bathing/showering, dressing, toilet visits, and eating) on a three-level scale (need for no (3), little, or much assistance (1)). These four values were converted into a single summary score (4–12p) of basic ADLs.
Self-perceived health status: At home visit, an EuroQol-visual analog scale (EQ-VAS; EuroQol Group, 1990; Rabin and de Charro, Reference Rabin and de Charro2001) was used. This is a generic instrument that assesses self-perceived health status in terms of EQ-VAS in which the end points are labeled “best imaginable health state” (100) and “worst imaginable health state” (0).
Cognitive functions: The overall cognitive function was assessed with the Mini-Mental State Examination (MMSE; Folstein et al., Reference Folstein, Folstein and McHugh1975), cognitive speed was assessed with a composite measure of the sum of the standardized values from Trail Making Test part A (Reitan, Reference Reitan1992) and part one of the Victoria version of the Stroop test (Regard, Reference Regard1981), and executive functions were assessed with a similar composite measure of the sum of the standardized values from part three of the Stroop test and the Parallel Serial Operations test (Nordlund et al., Reference Nordlund, Rolstad, Hellstrom, Sjogren, Hansen and Wallin2005).
The feeling of loneliness was measured with a single-item, author-constructed question based on a four-level scale of frequency (often (4), sometimes, seldom, or never (1); Park et al., Reference Park, Jang, Lee, Haley and Chiriboga2013).
Somatic multimorbidity or multiple coexisting chronic diseases were operationalized as the number of chronic diseases from an author-constructed predetermined list of 12 disease categories (Fortin et al., Reference Fortin, Stewart, Poitras, Almirall and Maddocks2012): arrhythmic heart disease, chronic heart failure or myocardial infarction, other vascular disease, hypertension or hyperlipidemia, diabetes mellitus, thyroid disease, respiratory disease, joint disease, central nervous system disease, gastrointestinal disease, urinary incontinence, and malignancy.
Statistical analyses
The SPSS statistical package (version 23.0) was used for data analyses. Differences in characteristics between participants and nonparticipants were investigated using chi square test and two-sidedt tests for each measure wave. Missing values (described more in the results section) were excluded listwise for the analyses. Linear mixed models were used for the longitudinal morbidity analyses because of the large amount of missing data, thereby calculating parameter estimates, 95% confidence intervals (95% CIs) and p values. The missing data were not missing at random. Effect sizes of the results were estimated as standardized differences by dividing the change in absolute score by the baseline standard deviation. Mortality was analyzed using univariate and multivariate Cox regressions with baseline data as independent variables, calculating hazard ratios (HRs), 95% CIs, and p values. The rather large proportions of nonresponse for the two variables of cognitive speed and executive functions were handled through multiple imputations before the Cox regressions, using five automatic imputations in all missing values of the independent variables of the regressions, according to the method of fully conditional specification in SPSS (Liu and De, Reference Liu and De2015). All the hypothesized independent variables and the sociodemographic variables of sex and educational level were selected for the multivariable Cox regression, using the enter method. The level of significance was set to α = 0.05.
Results
Sample characteristics and missing data description
Descriptive data of the participants at baseline are presented in Table 1. The distribution of the variables at baseline for each category of depressive symptoms has been described in more detail in a previous study (Ludvigsson et al., Reference Ludvigsson, Marcusson, Wressle and Milberg2016). The participation rate gradually diminished (Figure 1) with each measure wave, with only n = 55 (15% of n = 371 at baseline) responding to GDS-15 at the 8-year follow-up.
GDS-15, Geriatric Depression Scale; SSD, subsyndromal depression; ND, non-depression; ADLs, activities of daily living; I-ADLs, instrumental activities of daily living; MMSE, Mini-Mental State Examination; SD, standard deviation; EQ, EuroQol; VAS, visual analog scale.
a A more thorough comparison of the baseline values for the three categories of depressive symptoms, including tests with logistic regression, is presented in Ludvigsson et al. (Reference Ludvigsson, Marcusson, Wressle and Milberg2016).
b Adapted housing, i.e. sheltered housing, nursing home or dementia care.
Losses to follow-up were considered largely to be missing not at random, and factors significantly associated with attrition between each measure wave are presented in Table 2. The attrition was generally associated with adapted housing, lower levels of basic ADLs, I-ADLs, cognitive functions, and higher scores on GDS-15. Between T0 and T1, the attrition was also associated with female sex and lower educational level. The item nonresponse proportions ranged from 0% on most variables to higher proportions on the composite measures: 17% (64/371) missing on cognitive speed and 36% (132/371) on executive function, numbers that were reduced by means of multiple imputations before the Cox regression analyses.
GDS-15, Geriatric Depression Scale; SSD, subsyndromal depression; ND, non-depression; ADLs, activities of daily living; I-ADLs, instrumental activities of daily living; MMSE, Mini-Mental State Examination; SD, standard deviation; EQ-VAS, EuroQol-visual analog scale.
Values in bold letters show significance at the level of α = 0.05.
a More women and fewer men among those lost to follow-up, a recurrent tendency in all the measure waves, though the association was significant only between T0 and T1.
b More individuals with adapted housing and fewer with independent housing among those lost to follow-up.
c No comparison possible, as the first measure wave of this variable was at T1.
d More with SSD and fewer with ND among those lost to follow-up.
Longitudinal morbidity in SSD compared with ND
Table 3 shows the effect of SSD compared with ND over 8 years of follow-up on a number of outcome variables. Compared with ND, SSD was associated with higher morbidity in most of the hypothesized variables: functional abilities, self-perceived health status, social function, and mental health function. For example, in all measure waves, the presence of SSD instead of ND lowered the basic ADL scores of 0.58 points (standardized difference of 0.57) and lowered the I-ADLs scores of 5.4 points as a mean (standardized difference of 0.90). These effects on morbidity were still significant when adjusting for the time aspect and other potential confounders.
GDS-15, Geriatric Depression Scale; SSD, subsyndromal depression; ND, non-depression; ADLs, activities of daily living; I-ADLs,instrumental activities of daily living; MMSE, Mini-Mental State Examination; SD, standard deviation; EQ-VAS, EuroQol-visual analog scale.
Values in bold letters show significance at the level of α = 0.05.
a Diagonal covariance structure.
b Adjusted for sex, somatic multimorbidity and time. Variables other than basic ADLs and I-ADLs also adjusted for I-ADLs.
c A time by group interaction term was used as a fixed variable for the analyses, but this variable did not change the results significantly and was therefore left out for the final model.
d Ordinal data for measure waves T1–T4, except for the variable depressive symptoms for which the longitudinal outcome corresponds only to measure waves T2 to T4, as the T1 value was the independent variable in the analysis.
The effect of time or age on morbidity, on the other hand, was significant only for the variables of basic ADLs, MMSE, and I-ADLs where, for example, the 8-year movement from T1 to T4 meant a decrease of I-ADL score of 5.4 points as a mean (standardized difference of 0.90). The effect of SSD on cognitive functions was significant in the initial unadjusted analyses, but after adjusting for relevant factors the effect was no longer significant. Figure 2 shows the longitudinal development of the categories SSD and ND, and the participants with SSD compared with ND at age 85 showed a higher subsequent combined proportion of SSD and syndromal depression at each wave of follow-up, corresponding to a relatively greater morbidity in terms of depressive symptoms.
Mortality in SSD compared with ND
Table 4 shows mortality over 9 years of follow-up as a function of different independent variables at age 85. All the presented variables except sex and educational level had a significant association with mortality, according to univariate regressions. The variable SSD compared with ND showed the largest effect size of the independent variables, with SSD increasing the risk of death 1.6 times.
SSD, subsyndromal depression; ND, non-depression; ADLs, activities of daily living; I-ADLs, instrumental activities of daily living; MMSE, Mini-Mental State Examination; SD, standard deviation; EQ-VAS, EuroQol-visual analog scale; HR, hazard ratio.
Values in bold letters show significance at the level of α = 0.05.
a Final model built on n371, whereas 238 (64%) died, 133 (36%) were censored. Five automatic imputations were done in every independent variable with missing data, according to the method fully conditional specification in SPSS. Number of imputed values was 64 in the variable cognitive speed, 132 in the variable executive function, and lower numbers in the rest of the variables. Cases available for the multiple Cox regression were 163/371 (44%) before multiple imputation, and 371/371 (100%) after multiple imputation.
b The model was tested for proportional hazards, p = 0.889.
c Without imputations in this variable HR was 1.60 (1.21–2.11) and p = 0.001 in the univariate regression.
d When using an alternative ordinal variable (including all the three categories ND, SSD, and Depr; without any imputed values) instead of the dichotomous variables (SSD vs ND) and (Depr vs ND), the results of the multivariate regression were HR = 1.10 (0.783–1.55) and p = 0.578 for the categorical step ND-SSD, and HR = 1.03 (0.508–2.09) and p = 0.936 for the categorical step ND-Depr in the multivariate regression.
e Without imputations in this variable, HR was 1.35 (0.794–2.30) and p = 0.268 in the univariate regression.
f The results for the variable loneliness are presented for each categorical step: 1–2 comparison never-seldom, 1–3 comparison never-sometimes, 1–4 comparison never-often.
When combining all the independent variables into a multivariable Cox regression, only three of the variables showed significant association with mortality, namely, male sex, lower cognitive speed, and lower I-ADLs.
Discussion
Longitudinal morbidity in very old individuals with SSD
We have examined the longitudinal outcomes of SSD in very old individuals, in a population-based sample of 85-year-old individuals who were followed for 8 years. SSD compared with ND was associated with all the hypothesized variables of morbidity, even though such association with cognitive functions was no longer significant when adjusting for potential confounders. SSD compared with ND at baseline was further associated with 1.6 times higher risk of mortality according to our univariate Cox regressions, but against our hypotheses the association was not significant after controlling for other relevant variables in the multivariate regression. As presented in the following discussion, the findings were mainly in line with previous studies on depression in very old individuals as well as with previous studies on SSD in young old individuals (Blazer, Reference Blazer2003; Meeks et al., Reference Meeks, Vahia, Lavretsky, Kulkarni and Jeste2011).
The presence of SSD compared with ND was associated with lower levels of basic ADLs and I-ADLs over time, and the effect size was larger for I-ADLs compared to the former. These effect sizes were comparable to previous studies on younger individuals, though such effect sizes depend greatly on the type of measure (Schillerstrom et al., Reference Schillerstrom, Royall and Palmer2008). The association between depressive symptoms and ADLs over time has been described as complex though bidirectional, and a recent thorough analysis indicated that depressive symptoms drive ADL impairment more than vice versa (Hajek et al., Reference Hajek2017).
Also, the relation between cognitive functions and depressive symptoms has been described as bidirectional in both demented and nondemented study samples (Koenig et al., Reference Koenig2015). In this study, we found that SSD was associated with lower overall cognitive function as well as with lower cognitive speed and executive function over time but when adjusting for sex, somatic multimorbidity and I-ADLs, the association was not significant. Possibly the absence of a significant association in the latter was a result of selective losses to follow-up, as lower cognitive function was significantly associated with nonparticipation in consecutive measure waves. Other possible explanations would be that a real association was obscured by the selective survival of cognitively healthier individuals before the age of 85 or by competing risks from other causes in the oldest ages.
SSD was associated with lower self-perceived health and about equally strongly associated with frequency of loneliness over 8 years of follow-up, which was in line with the previous studies (Houtjes et al., Reference Houtjes, van Meijel, van de Ven, Deeg, van Tilburg and Beekman2014; Park et al., Reference Park, Jang, Lee, Haley and Chiriboga2013). Regarding psychiatric outcomes, SSD was associated with higher levels of depressive symptoms over time, which was significant when adjusting for potential confounders. For those with SSD at baseline, the combined risk of SSD or depression at each consecutive measure wave was elevated compared with participants with ND at baseline, even when considering the higher risk of loss to follow-up. This finding is in accordance with the previous findings for SSD in younger ages, and with a dimensional view on affective disorders in which SSD and depression are adjacent states with a common origin (Beekman et al., Reference Beekman2002; Judd, Reference Judd2012).
Our study did not show any significant association between time or age and depressive symptoms, and whether such an association really exists has been repeatedly investigated—several previous studies have shown that there is no such association when controlling for confounders and cohort effects (Blazer, Reference Blazer2000; Yang, Reference Yang2007).
Mortality in very old individuals with SSD
Turning to the effect of SSD on mortality, the univariate Cox regression showed a rather strong association in terms of HR, which subsided when controlling for sex, somatic multimorbidity, I-ADLs and other factors of the multivariate regression. Most probably the interactions between depressive symptoms and different mediators of mortality are reciprocal in a complex system of mediators with sometimes synergistic effects on mortality (Schulz et al., Reference Schulz, Drayer and Rollman2002).
Previous studies have further shown a weaker or a nonexistent association between depressive symptoms and mortality in very old individuals, and a stronger association in younger ages (Cuijpers et al., Reference Cuijpers, Vogelzangs, Twisk, Kleiboer, Li and Penninx2013; Rapp et al., Reference Rapp, Gerstorf, Helmchen and Smith2008). In contrast to the two previous studies by Rapp et al. and Stek et al., respectively, that have investigated the association between SSD and mortality in very old ages, we found that this association was significant when adjusting only for sex and comorbidity (Table S1 published as supplementary material online) (Rapp et al., Reference Rapp, Gerstorf, Helmchen and Smith2008; Stek et al., Reference Stek, Vinkers, Gussekloo, Beekman, van der Mast and Westendorp2005). However, in our study the association subsided when adding the covariate I-ADLs, and, therefore, our conclusion is (similar to the conclusion of Rapp et al.) that an independent association has been observed between SSD and mortality in lower ages, but not in the very old ages. Possible explanations for the difference between age-groups are selective survival, where individuals with depressive symptoms die earlier in life while more resilient individuals live longer, and also the dilution effects on the association as comorbidity and frailty increase with age and increase mortality independently (Vaughan et al., Reference Vaughan, Corbin and Goveas2015). Maybe such dilution effects with increasing ages can explain why the association was not significant between SSD at baseline and mortality over the 9 years of follow-up in our study of very old individuals, when considering interacting factors.
The significance of SSD for morbidity and mortality in the aged population
The magnitudes of the association between SSD and the outcome variables were generally small in this study, for example, with percentage effects of the measure range 7 for basic ADLs and 23 for I-ADLs. However, the effects were clearly at significant levels; and in comparison with syndromal depression, the smaller effects of SSD on diverse outcomes at an individual level might mean a comparable effect for the population attributable fraction of the same conditions, due to a higher prevalence of SSD than of syndromal depression. The association between several longitudinal outcomes of morbidity and SSD compared with ND in our study indicates that the relatively mild condition of SSD has substantial influences on cognitive, functional, and psychosocial aspects of health, not only in younger adults but also in very old individuals. Interventions to prevent or treat SSD among very old individuals might correspondingly reduce these different aspects of morbidity.
Future studies could extend our understanding of possible causal aspects of the relationships between depressive symptoms and morbidity in old age by using tighter measure waves while also controlling for mediating factors clearly defined a priori. A more thorough analysis of single depressive symptoms could also improve the possibilities to disentangle moderators or relevant phenotypes in the heterogeneity of SSD and syndromal depression. In addition, there is a need for prevention studies of SSD in old individuals, to assess whether interventions might reduce morbidity, and to improve the quality of life (Gilbody et al., Reference Gilbody2017). Regardless of whether SSD in old individuals is considered as pathology or normality (Ludvigsson et al., Reference Ludvigsson, Milberg, Marcusson and Wressle2015), such interventions could possibly improve life and decrease the need of health care for many people, as SSD is common and the number of very old individuals is steadily increasing around the world.
Limitations
This study has several limitations. First, only 57% of the eligible population or 75% of the baseline participants responded to the GDS measurement, and systematical differences were observed between the responders and nonresponders. Therefore, the generalizability is limited, especially for certain subgroups of very old individuals, namely, females, those with adapted housing, higher educational level and lower ADL functions. Second, the subsequent attrition between the measure waves was substantial as in many longitudinal studies, with more attrition among individuals with lower levels of functions and higher levels of depressive symptoms. Even if the problem of missing data for some calculations was reduced with the multiple imputations method, the attrition probably had the consequence of falsely weakened association, though the exact proportions and directions of this selection bias were not estimated. The relatively small baseline sample sizes also limit the power of the results. Furthermore, SSD was defined by the screening instrument of GDS-15 instead of the Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition classification registry, which is of relevance for the validity. The long time range of the study as well as the relatively large number of measure waves are the strengths of this study, while the rather long and irregular time intervals in between the waves limited the possibilities to make causal inferences.
Conclusion
This 8-year longitudinal study showed that morbidity in a population sample of 85-year-old individuals was significantly elevated in those with SSD compared with ND, even after controlling for relevant covariates. Mortality was also elevated in those with SSD according to univariate regressions, but this association vanished in multiple regressions. This contributes to previous diverging study results on the mortality of depressive symptoms in very old individuals and confirms previous results, showing that even mild depressive symptoms are associated with significant levels of morbidity in the population. The increasing proportions of very old individuals necessitate a better understanding of the consequences of the whole range of depressive symptoms in old individuals as well as studies on prevention and treatment in the subject area. Regardless of whether SSD is considered as pathology or normality, such interventions could possibly improve quality of life and reduce morbidity and the need for health care among many old people.
Conflicts of interest
The authors have no conflicts of interest to declare.
Funding
This project was supported by grants from Region Östergötland, Sweden. The funding source had no involvement in the study design, data collection, analysis, interpretation of the data, or in writing the manuscript.
Description of authors’ roles
All authors contributed to the design of the study. Mikael Ludvigsson, Ewa Wressle, and Jan Marcusson collected the data. Mikael Ludvigsson managed the literature search and analyses. Mikael Ludvigsson and Anna Milberg undertook the statistical analyses and wrote the first draft of the manuscript. All authors contributed to and have approved the final manuscript.
Acknowledgments
We thank Mr Mats Fredriksson at Forum Östergötland, Linköping University, for help with the statistical analyses. This project was supported by grants from Region Östergötland, Sweden. The funding source had no involvement in the study design, data collection, analysis, interpretation of the data, or in writing the manuscript.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S1041610219001480