Skip to main content

Trajectories of work disability among individuals with anxiety-, mood/affective-, or stress-related disorders in a primary healthcare setting

Abstract

Background

Anxiety-, mood/affective-, or stress-related disorders affect up to one-third of individuals during their lives and often impact their ability to work. This study aimed to delineate trajectories of work disability (WD) among individuals diagnosed with anxiety-, mood/affective-, or stress-related disorder in primary healthcare and to examine associations between trajectory group membership and sociodemographic, clinical, and clinical-related factors.

Methods

The study population included working-age individuals, aged 22–62 years, living in Stockholm County, Sweden, who experienced a new episode of any anxiety-, mood/affective, or stress-related disorder in primary healthcare in 2017 (N = 11,304). Data were obtained from Swedish national and regional registers and were linked using pseudonymised unique personal identification numbers. The primary outcome was days with WD (sum of sickness absence and disability pension days) during the three years before and three years after a diagnosis of anxiety-, mood/affective-, or stress-related disorders in primary healthcare. A zero-inflated Poisson group-based trajectory model was used to identify groups of individuals with similar patterns of WD over the study period, with a multinomial logistic regression used to examine associations of sociodemographic, clinical, and clinical-related factors with trajectory group membership.

Results

Four distinct trajectory groups were found, high increasing (5.1%), with high levels, from 16 to 80 days of WD in six-monthly intervals during follow-up, peak (11.1%), with a peak in WD, up to 32 days of WD, around the time of the diagnosis, low increasing (12.8%), with an increase in days of WD from 4 to 22 during the study period, and constant low (71.1%), with almost no WD over the study period. In multinomial regression models, diagnostic category, psychotropic medication use, a diagnosis of a psychiatric disorder within secondary healthcare, age at diagnosis, and occupation were associated with WD trajectory groups.

Conclusions

Around two-thirds of individuals treated for a new episode of any anxiety-, mood/affective-, or stress-related disorder in primary healthcare have an excellent prognosis regarding WD. Several sociodemographic and clinical characteristics were associated with group membership; these factors could identify individuals at risk of long-term welfare dependency and who might benefit from interventions to promote a return to work.

Peer Review reports

Background

Anxiety-, mood/affective-, or stress-related disorders affect around one-third of individuals during their working lives [1]. These disorders often affect the ability to work and are among the most common causes of long-term sickness absence and the granting of disability pension in Sweden [2,3,4,5]. These disorders, particularly stress-related disorders, have increased substantially during the last two decades [4, 6]. Possible reasons for this might be that work demands have increased simultaneously, and help-seeking for these disorders has changed [7]. Inability to work because of anxiety-, mood/affective-, or stress-related disorders can be very costly, both from a social and an employer perspective, but not least for the individuals who risk permanent exclusion from the labour market and risk being dependent on welfare benefits. A previous population-based study reported that as many as two-thirds of young adults (aged 20 to 35 years) treated for their first episode of a common mental disorder (CMD, which included anxiety, mood/affective, or stress-related disorders) in specialised healthcare settings had poor prognoses concerning later labour market participation [5]. This may be expected, as those treated for CMDs in secondary healthcare can be expected to have more severe clinical presentations than those treated within primary care and, therefore, are more likely to have long periods outside the labour market. However, given that most individuals who develop anxiety-, mood/affective-, or stress-related disorders are treated solely within primary care settings (without referral to secondary care), it is crucial to determine the long-term work trajectories for people treated for anxiety-, mood/affective-, or stress-related disorders in this setting as this represents the vast majority of people with these disorders [8].

Identifying risk and protective factors for long-term labour market marginalisation among individuals with anxiety-, mood/affective-, or stress-related disorders is essential to developing strategies to promote recovery. Systematic psychological therapies (e.g., cognitive behavioural therapy and psychodynamic psychotherapy) and pharmacological therapies (e.g., antidepressants) have been previously found to be associated with some improvements in symptoms and functioning among people with CMDs [9,10,11,12,13]. However, previous studies investigating associations between treatment factors and work disability among individuals with CMD have typically utilised small samples, making it difficult to detect differences between different subtypes of CMDs [14].

To address these knowledge gaps, we used data from national and regional registers obtained from Swedish healthcare and other government authorities to examine trajectories of work disability among individuals with anxiety-, mood/affective-, or stress-related disorders. These high-quality registers capture periods of work disability and provide information on sociodemographic and work-related factors (e.g., gender, age, educational level and profession) as well as clinical-related factors (e.g., type of disorder, comorbid somatic and psychiatric conditions and psychological and pharmacological therapies) that may be associated with trajectories of work disability [2, 5].

Days out of work due to illness is a significant public health issue as more extended periods of work disability tend to increase labour market marginalisation as well as declining mental illness. Finding individuals at most risk of ending up outside the labour market is a prioritised issue to avoid permanent exclusion from the labour market and worsening of mental health [15, 16]. We, therefore, aimed to delineate trajectories of work disability among individuals diagnosed with a new episode of any anxiety-, mood/affective-, or stress-related disorder in primary healthcare. Our primary aim was to assess the heterogeneity in the population and determine how belonging to these trajectory groups was associated with sociodemographic, work-related, and clinical factors.

Methods

Setting

In Sweden, primary healthcare has been commissioned by the government to prevent and treat individuals with mild to moderate mental health problems, including anxiety-, mood/affective-, or stress-related disorders. Patients with severe psychiatric disorders or elevated risk for suicide are immediately referred to secondary psychiatric clinics (to which they may also self-refer). However, most patients are referred from primary to secondary healthcare only after undergoing one or several treatment attempts in the primary care setting.

Data sources

Data were sourced from various Swedish national and regional registers, linked individually using pseudonymised unique personal identification numbers. The VAL database is the data storage solution for healthcare visits within Region Stockholm, Sweden, which includes individual-level data on primary care visits in Stockholm County. The Longitudinal Integration Database for Health Insurance and Labour Market Studies (LISA) contains data on sociodemographic and work-related factors. The Micro-Data for Analyses of Social Insurance (MiDAS) included data on work disability. The National Patient Register (NPR) included data on treatment in inpatient and specialised outpatient healthcare. The Prescribed Drug Register (PDR) included information on drug purchases. An overview of the registers utilised is provided in Supplementary Table 1.

Study population

In the VAL database, we first identified all individuals with any recorded diagnosis of anxiety-, mood/affective-, or stress-related disorders during 2017 (N = 124,963) [17]. The study population included working-aged individuals in Stockholm experiencing a new episode of anxiety- (International Statistical Classifications of Diseases – version 10 (ICD-10) [18]: F40-F42), mood/affective- (ICD-10: F32-F39) or stress-related disorders (ICD-10: F40-F42) in primary healthcare during the year 2017. This year was chosen to have data for healthcare visits three years before and after the inclusion date, as reliable diagnostic data has only been available since 2014. Individuals were required to be registered residents in Stockholm County and of working age (22–62 years) during each calendar year between 2014 and 2019, inclusive (N = 84,275). The reason for choosing 22 years as the lower age range was to ensure that all participants in the study were eligible to be granted disability pension three years before inclusion, as it is not possible to receive disability pension before the age of 19 in Sweden. To identify individuals experiencing a new episode of anxiety-, mood/affective-, or stress-related disorder, that is, a visit to primary healthcare due to any of these disorders in 2017, without having a visit to either primary or secondary healthcare due to anxiety-, mood/affective-, or stress-related disorders three years before, all individuals with an inpatient, specialised outpatient, or primary care contact for any anxiety-, mood/affective-, or stress-related disorders diagnosis in the three years before cohort entry were excluded from the study (N = 32,868). Further exclusions were made for individuals who, at any point in the three years before cohort entry, had spells of sickness absence attributed to any anxiety-, mood/affective-, or stress-related disorder, purchased psychotropic medications indicated for anxiety-, mood/affective-, or stress-related disorders (Anatomical Therapeutic Chemical Classification, ATC: N06A, N05B, N05C, R06AD01) or had visits in primary care where psychological therapy (including both systematic and non-systematic psychological treatment) had been administered (N = 11,635). Finally, individuals with severe mental illnesses, including schizophrenia-spectrum disorders (ICD-10: F20-F29), bipolar disorder (ICD-10: F30-F31), and organic mental disorders (ICD-10: F00-F09), and those who had been granted full-time disability pension before cohort entry, were excluded as these conditions are associated with high levels of disability and therefore may confound the association between anxiety-, mood/affective-, and stress-related disorders and our primary outcome. As anxiety-, mood/affective-, and stress-related disorders are highly co-morbid with many different psychological disorders, we accounted for other mental comorbidities (ICD-10: F10-F19 and F44-F99) by including these as covariates in the analyses. The final study population consisted of 11,304 individuals.

Outcome measures

The primary outcome of this study was days of work disability, defined as the sum of net days, that is, one combined day of work disability. Here, one day, two half days or four quarters of a day equals one net day with sickness absence and disability pension benefits. Work disability was measured at six-monthly intervals over the study period, which spanned three years before and three years after cohort entry (i.e., the date of the first diagnosis of a mood/affective-, anxiety-, or stress-related disorder in 2017).

All individuals aged 19–64 in Sweden can be granted disability pension if their work capacity is reduced due to illness. Individuals aged 19–29 may be granted time-restricted disability pension if they have an impaired work capacity or have not completed compulsory education. Individuals aged 30–64 years may be granted permanent disability pension if they have an expected lifelong duration of reduced work capacity. Individuals aged 16 and above with a certain income level from work can receive benefits concerning sickness absence. Employers are typically responsible for payment of sickness absence benefits during the first 14 days of a sick leave spell, after which the Social Insurance Agency covers payments. For this reason, only information on sickness absence spells longer than 14 days is readily available in national registers and was used to compute the outcome measure.

When computing the outcome measure, the extent to which sickness absence and disability pension were granted were taken into account; for example, two days of part-time sickness absence at a 50% extent was combined to equal one net day of work disability. In the study, disability pension and sickness absence were summed up to encompass days on work disability.

Sociodemographic and work-related factors

Sociodemographic factors included gender (male vs. female), country of birth (Sweden vs. elsewhere), years of education (0–9 vs. 10–12 vs. > 12 years), family composition (cohabiting with no children living at home vs. cohabiting with children living at home vs. single with no children living at home vs. single with children living at home), place of residence (large cities vs. other), and days of unemployment (none vs. any). These factors were measured on the 31st of December 2016, the year before the diagnosis of an anxiety-, mood/affective-, or stress-related disorder. The occupation group was measured in November 2016 based on employment status. It was categorised according to the Swedish Standard Classification of Occupations 2012 (SSYK 2012) as (1) occupations with university requirements (SSYK 2012): 1111–3522), (2) administrative roles and customer service, (SSYK 2012: 4111–5419), (3) non-managerial workers (SSYK 2012: 6111–9629, 0110–0310), (4) not gainfully employed and (5) no information on the occupation. Age was estimated during the cohort entry year (2017) and categorised as 22–29 vs. 30–39 vs. 40–49 vs. 50–59 vs. 60–62 years.

.

Clinically related factors

Type of diagnostic category included anxiety disorders (ICD-10: F40-F42), mood/affective disorders (ICD-10: F32-F39), and stress-related disorders (ICD-10: F43), that is, reaction to severe stress, adjustment disorders, acute stress and post-traumatic stress disorders (PTSD) from primary healthcare. The treatment variables included the receipt of systematic psychological treatment (defined as psychodynamic treatment (PDT), cognitive psychological therapy (CT), cognitive behavioural therapy (CBT), mentalization-based therapy (MBT), Eye Movement Desensitization and Reprocessing (EMDR), systemic therapy (ST), dialectical-behavioural therapy (DBT), interpersonal therapy (IPT), and other psychological treatment), administered by psychotherapists in primary care (treatment codes (KVÅ codes) in the registers: DU008- DU011, DU013-DU014, DU020-DU022) and purchases of psychotropic medications including antidepressants (Anatomical Therapeutic Chemical code (ATC): N06A), anxiolytics (ATC: N05B), hypnotics and sedatives (ATC: N05C) and alimemazin (R06AD01). These exposures were measured over the three years following the new episode of any anxiety-, mood/affective-, or stress-related disorder. Moreover, diagnoses of visits from inpatient/specialised outpatient healthcare were measured for all psychiatric disorders (ICD-10: F00-F99).

Statistical methods

Zero-inflated Poisson group-based trajectory models (ZIP GBTM) were used to identify groups of individuals with similar patterns of work disability over the study period [19]. This method is a type of finite mixture model appropriate for analysing evolving longitudinal processes. Instead of assuming the existence and forms of hypothesised trajectories of work disability, the method can identify distinct and meaningful patterns of work disability from the data itself [20]. These models assume the presence of unobserved latent groups, where each group follows a unique trajectory of work disability over time. GBTMs can accommodate various data types; in this instance, a zero-inflated Poisson (ZIP) distribution was determined to provide the best fit based on the distribution of work disability days (this outcome is highly skewed, with a large proportion of individuals having zero days of work disability). A ZIP model provides a more appropriate fit for data with these characteristics than other models such as the standard Poisson or a censored normal model). Proportions of the population belonging to each latent group and the probability of group membership can also be estimated from the model.

Due to the explorative nature of the GBTM method, different models with varying numbers of latent groups and other sets of parameter constraints were compared to identify the best-fitting model. Nagin [20] suggests a practical two-stage approach for model selection. The first stage involves determining the optimal number of trajectory groups, and the second stage involves identifying the optimal polynomial form of the trajectories. We fitted models that included two to six trajectory groups with varying polynomial shapes (from linear to quartic). Several diagnostic metrics were considered when selecting the optimal model; these included the Akaike information criterion (AIC), the Bayesian information criterion (BIC), the average posterior probability of assignment (APPA), the mean squared error (MSE), relative entropy, and the and the odds of correct classification (OCC).

After identifying the optimal trajectories of work disability, individuals were assigned trajectory groups based on estimated group membership probabilities. A multinomial logistic regression was then conducted to examine the associations between sociodemographic, work-related and clinical covariates with trajectory group assignment. Although we refer to these as exposures, they are measured at a single time point. Still, the outcome measure incorporates data regarding work disability from before and after diagnosis. Each covariate was excluded from the saturated model to assess the relative contribution of the different covariates; the reduced and saturated models were then compared using likelihood-based metrics, including likelihood ratio tests and Nagelkerke R2. Marginal effect plots are presented for the exposures and the covariates with the highest relative importance in the models. Statistical analyses were performed in R (version 4.2.2) using the ‘latrend’, ‘nnet’ and ‘effects’ packages.

Results

Characteristics

The characteristics of the 11,304 individuals treated for a new episode of any anxiety-, mood/affective-, or stress-related disorder in primary care during 2017 are presented in Table 1. Most individuals (53.8%) were aged 30 to 49 years at cohort entry and were born in Sweden (72.7%); just over half (50.1%) were female. The most common family situation category was single without children (47.1%), and nearly half (49.0%) had more than twelve years of education. With regards to work-related factors, less than 1 in 10 (8.3%) were unemployed in the year before cohort entry; the most common occupations were those that required a university education (40.8%), followed by administration and customer service occupations (25.2%).

Table 1 Baseline characteristics, overall and across trajectory groups for different trajectory groups for individuals diagnosed with CMDs in Stockholm County in 2017 (n = 11,304)

About 3 out of 10 individuals received at least one session of systematic psychological treatment during the three years after the diagnosis of any anxiety-, mood/affective-, or stress-related disorder, whereas the majority had their session of systematic psychological treatment within one year after the diagnosis. At least half of the population purchased psychotropic medications during the three-year follow-up (55.2%). In the study, 75% of individuals were diagnosed with just one of the three categories: anxiety-, mood/affective-, or stress-related disorders, while the remainder (25%) had combinations of these conditions. Over half of the study population (52.2%) had a diagnosis of a stress-related disorder. Within this group, 33.2% were diagnosed solely with a stress-related disorder, and 19.0% had a stress-related disorder alongside an anxiety- or mood/affective disorder. Anxiety disorders were found in just under half (46.9%) of the population, where 27.5% experienced solely an anxiety disorder and 19.4% had an anxiety disorder combined with a mood/affective- or stress-related disorder. Mood/affective disorders were found in about one-third (29.6%) of the population, with 14.3% having only a mood/affective disorder and 15.3% experiencing this alongside an anxiety- or stress-related disorder. The most common combination of disorders found in the study was between stress-related and anxiety disorders, which affected 9.8% of the population.

Trajectory groups

The model with four trajectory groups produced the lowest AIC, BIC and MSE metrics, while the remaining metrics were similar across all models. When identifying the optimal polynomial form of the trajectories, AIC and BIC decreased marginally as complexity was added to the model; however, the four-group quadratic model was selected as the most parsimonious model as it produced the highest APPA and relative entropy metrics and with the lowest MSE. The average levels of work disability across the study period, estimated from the overall model, are shown in Fig. 1 for the four WD trajectory groups. The trajectory groups were: high increasing (with 5.1% of the included population), characterised by a high starting level of average work disability (16 six-monthly intervals ), which continued to increase during follow-up; peak (11.1%), characterised by lower levels of average work disability at the start and end of the study period but with a peak around the time of diagnosis of the new episode of anxiety-, mood/affective-, or stress-related disorders; low increasing (12.8%), characterised by a lower starting level of average work disability which increased over the entire study period; and constant low (71.1%), with low levels of average work disability during the whole study period.

Fig. 1
figure 1

Trajectory groups of work disability around the year 2017 (t0) from 2014 (-3) to 2020 [3] divided into six-month intervals

Characteristics and predictors across trajectory groups

In the multinomial regression model, the most important predictors of trajectory group membership were found to be the type of diagnostic category, i.e., anxiety disorders, mood/affective disorders, and stress-related disorders (4.8% difference in R2 when introduced into the model), diagnosis of any psychiatric disorder in secondary healthcare (4.0%) following the new episode of anxiety-, mood/affective-, or stress-related disorders, profession (1.2%), age at cohort entry (1.0%) and purchases of psychotropic medication (1.0%). Marginal effect plots, showing the estimated probability of trajectory group membership across different levels of these covariates, along with psychological treatment, are shown in Fig. 2.

Fig. 2
figure 2

Predicted probability plots showing the probability of trajectory group membership

The fully adjusted multinomial logistic regression with estimated odds ratios (OR) and 95% confidence intervals (CI) are shown in Table 2, where the high-increasing, peak and low-increasing trajectory groups were compared with the constant low group which consisted of a higher proportion of 22–29-year-olds (20.5%) and fewer individuals aged over 50 (25.2%) compared to the other classes with higher levels of work disability (11.6 − 14.0% 22–29 years old; 32.2 − 39.5% over 50, Table 1). The Constant low group also comprised a high proportion of individuals with a high educational level (51.2%).

Table 2 Fully adjusted odds ratios (ORs) of trajectory group membership associated with sociodemographic, work- and health-related factors, compared to the reference group (constant low trajectory of work disability) in individuals with anxiety-, mood/affective-, or stress-related disorders treated in primary health care settings in Stockholm, Sweden

Additional crude multinomial logistic regressions were performed for all covariates. In this crude analysis, a larger number of significant associations were found between treatment variables and work disability; however, after adjusting for the other covariates, these effects decreased, indicating that a portion of the effects are accounted for in the variation of the different variables.

Sociodemographic and work-related factors

Differences in trajectories of work disability were found between professions. The fully adjusted odds ratios (ORs) reveal that non-managerial workers (OR: 1.72; 95% CI: 1.25–2.38) and workers within the administration and customer services (OR: 1.37; 95% CI:1.05–1.79) had a higher probability of belonging to the high-increasing trajectory as well as to the low-increasing trajectory (OR: 1.50; 95% CI: 1.11–1.54 and OR: 1.30; 95% CI: 1.21–1.84, respectively) than individuals in occupations with university requirements. With regards to age, individuals in the age group 60–62 years had a considerably higher probability of belonging to the high-increasing trajectory with high levels of work disability (OR: 6.28; 95% CI: 4.08–9.68), but also to the low-increasing trajectory of work disability (OR: 2.21; 95% CI: 1.63–3.01) compared to the youngest age group of 22–29 years. Generally, higher ages were associated with an increasing probability of belonging to the high-increasing trajectory of work disability (OR range 1.68–6.28).

Clinical-related factors

The type of diagnostic category explained the most significant variance in trajectory group membership according to the R2 value. Compared to individuals with anxiety disorders alone, those with stress-related disorders, either in isolation (OR: 2.29; 95% CI: 1.71–3.06) or in combination with anxiety disorders (OR: 2.83; 95% CI: 2.00–4.01) or mood/affective disorders (OR: 5.36; 95% CI: 3.74–7.68), had the highest probability of being in the high-increasing trajectory of work disability. A similar trend was also seen for belonging to the low-increasing trajectory group, albeit with lower odds ratios. Receiving a diagnosis for a psychiatric disorder in inpatient/specialised outpatient healthcare within three years following the new episode of any anxiety-, mood/affective-, or stress-related disorder was strongly associated with an increased probability of belonging to the high increasing (OR 6.79; 95% CI: 5.49–8.40), peak (OR 3.68; 95% CI: 3.12–4.33) as well as the low increasing (OR 2.14; 95% CI: 1.82–2.53) trajectories.

Purchases of psychotropic medication after the diagnosis of any anxiety-, mood/affective-, or stress-related disorder were associated with higher odds of following a high-increasing (OR: 2.34; 95% CI: 1.83–2.98), peak (OR: 1.83; 95% CI: 1.53–2.19) or low-increasing (OR: 1.51; 95% CI: 1.28–1.78) trajectory of work disability, compared to the constant low group. Furthermore, there were tendencies for individuals who received psychological therapy in the year following the episode of the anxiety-, mood/affective-, or stress-related disorder to be more likely to belong to the high-increasing group than the constant low group, as compared to those who did not receive psychological therapy (OR: 1.38; 95% CI: 1.09–1.75). Individuals who received systematic psychological treatment in all three years of the follow-up period had even higher odds of following the high-increasing trajectory of work disability (OR: 1.63; 95% CI: 1.16–2.30), compared to the constant low group in contrast to individuals who did not receive therapy.

Discussion

In this sizeable register study of individuals experiencing a new episode of anxiety-, mood/affective-, or stress-related disorders treated in primary healthcare settings, we identified four distinct trajectory groups: high-increasing (5.1% of the population), characterised by high levels of work disability during follow-up, peak (11.1%), with a peak around the time of diagnosis, low-increasing (12.8%), with increases in work disability throughout the study period and constant-low (71.1%), with low levels of average work disability over the study period. Stress-related disorders, either alone or in combination with mood/affective disorders or anxiety disorders, a diagnosis of any psychiatric disorder in secondary healthcare following the diagnosis of an anxiety-, mood/affective-, or stress-related disorder, purchases of psychotropic medication, having a non-managerial, administrative, or customer service occupation, and higher age at diagnosis were most strongly associated with trajectory group membership.

Trajectory groups

Nearly three-quarters of all participants in this study experienced a constant low trajectory of work disability, with none or very little work disability, during the three-year follow-up period. In a previous study by our group, we examined trajectories of work disability in a population of individuals treated for CMDs in secondary healthcare (who can be assumed to have a higher level of illness severity than the current population) [5]; in comparison, individuals in the present study had a markedly better connection to the labour market. Consistent with our prior work, we observed that individuals who were previously treated for psychiatric disorders in secondary care had a much higher risk of following a high work disability trajectory. Those belonging to the trajectory group with low increasing work disability had a relatively modest increasing trend of work disability, with about 20 days of WD in the last six-month interval. Whilst long-term work disability is considered a poor outcome, periods of short-term sickness absence may be beneficial for maintaining a sustainable work-life [21]. The curve of days of work disability for those following the low-increasing trajectory did, however, increase during the follow-up, which might lead to high levels of work disability some years after the diagnosis. Here, measures might be taken early to prevent long future spells of work disability.

The increasing-high trajectory group, who experienced a rapidly increasing pattern of work disability some years before the diagnosis of the anxiety-, mood/affective-, or stress-related disorder, continued to grow in the years after the diagnosis and accounted for an average of 80 days with work disability three years after the diagnosis. From a welfare perspective, this group is arguably the most challenging as such individuals are at risk of being dependent on welfare benefits over the long term. Reducing days of work disability within this group has the potential to both decrease welfare costs for societies and increase well-being among individuals, thereby enhancing a sustainable working life.

The Peak trajectory of work disability, indicating a recovery in the time after the diagnosis, is, from a welfare perspective, the most optimal group. Those following this trajectory have an increase in work disability around the time of diagnosis but have a relatively rapid recovery and have a much-reduced work disability two years after the diagnosis. This group might, therefore, not require extra support to reduce days on work disability.

Predictors of trajectory groups

Sociodemographic and work-related factors

We found that profession was an important factor in explaining the differences between trajectory groups of work disability among individuals diagnosed with anxiety-, mood/affective-, or stress-related disorders in primary health care. Non-managerial workers had a high probability of belonging both to the high-increasing trajectory and the low-increasing trajectory of work disability compared to individuals with occupations with university requirements. Earlier studies describe this and have suggested that such findings may be explained by less flexible working patterns, meaning that individuals cannot decide when or where to complete their work activities [22]. Also, those working within the administration and customer support had a high probability of belonging to either of the increasing trajectories of work disability. Therefore, these groups might be prioritised by the stakeholders within rehabilitation to reduce the risk of having long periods of work disability.

In general, age was considered an essential determinant of work disability. Therefore, our finding that age was among the strongest predictors of trajectory group membership is not unexpected. Nevertheless, the exceptionally high probability of sickness absence among those aged 60–62 is a novel finding. The retirement age has risen worldwide, including in Sweden [23], and is expected to increase further in parallel with global trends in life expectancy. Policies, not least Agenda 2030, stipulate a sustainable work-life throughout life [24]. Whilst these findings require further replication, they suggest that enabling persons over 60 to work with anxiety-, mood/affective-, or stress-related disorders is crucial to reaching sustainable work among individuals near retirement.

Clinical-related factors

Those with stress-related disorders, either alone or in combination with mood/affective disorders and anxiety disorders, were associated with the highest probability of belonging to the high increasing trajectory. Such findings are important given that studies have shown that stress-related disorders appear to be increasing [25]. One must, however, keep in mind that other reactions to severe stress (ICD-10: F48.8) in Sweden allow for much more extended periods of sickness absence compared to depressive disorders or anxiety disorders, so the number of individuals diagnosed with exhaustion disorder might, therefore, be a consequence of the administrative regulations within the Social insurance scheme rather than a marker of severity of the disorder [26]. Job demands have increased much during the past decades, and if high job control and support do not mitigate the effect of high job demands, the health consequences can lead to work disability [7, 27,28,29]. Also, the length of the sickness period at baseline seems to affect the ability to have a sustainable working life in the future. In a previous study from our research group, we found that individuals diagnosed with stress-related disorders in isolation were less likely to receive psychological treatment and psychotropic medication [30] and, as a consequence, may experience a higher degree of work disability.

With regards to clinical-related factors, we observed that individuals who had purchased psychotropic medications had a higher probability of belonging to the high-increasing work disability group. Such findings are likely to reflect confounding by indication, whereby individuals who are the most unwell (and therefore most likely to receive work disability payments) are the most likely to receive medication and systematic psychological treatment. We found, most probably due to confounding by indication, only limited evidence to suggest that receipt of neither a systematic psychological treatment nor medication was associated with the trajectory groups of work disability. These findings are, to some extent, also consistent with previous studies [10, 11]. Using the current study design, we cannot determine whether the most unwell individuals experienced an improvement in work disability following pharmacological treatment. An implication of this is that studies using statistical methods to account for treatment allocations in real-world practice are needed to fully assess the effect of confounding by indication and thereby estimate the impact of treatment on work disability.

Strengths and limitations

The strengths of this study include the use of high-quality register data, including individual-level information on several sociodemographic, work and clinical-related covariates [31, 32]. Also, the population-based design minimised the follow-up loss, including all adults between 22 and 62 who lived in Stockholm County and were diagnosed with anxiety-, mood/affective-, or stress-related disorders, and the study period of six years enabled us to assess the development of sickness absence both before and after a diagnosis.

There are also limitations worth mentioning. Data on sickness absence only includes information on spells longer than 14 days, as the employers cover the first two weeks and do not register these days. This might slightly underestimate the risk of work disability among those with anxiety-, mood/affective-, or stress-related disorders. Given the study design, this is not a significant problem as we are interested in more prolonged sickness absence spells. Moreover, we only know about receiving psychological therapy in primary care. Although VAL includes many private clinics, we may lose some individuals treated in private clinics. As our study included only individuals residing in Region Stockholm, a highly affluent region in Sweden, our findings may not be generalised to individuals living in rural parts of a country. Due to the risk of confounding by indication, our findings regarding associations between treatment receipt and work disability must be interpreted cautiously. Another crucial limitation regarding the validity of results pertains to the lumping together of very different types of disorders under rubrics, such as “stress-related disorders” (including both post-traumatic disorder and adjustment disorder) and “anxiety disorders” (including a broad range of very different disorders such as obsessive-compulsive disorder (OCD) and panic disorder). According to the Swedish National Board of Health and Welfare, these disorders have significantly differential recommendations for prescribing sick leave [26]. Our study findings cannot be generalised to all people with anxiety-, mood/affective-, or stress-related disorders, as our sample is based on individuals who seek healthcare and who received a diagnosis in primary care. In Europe, for example, refugees have been reported to have limited access to mental healthcare compared to their native-born peers [33]. In a review, about 50–60% of those who fulfil the criteria for being diagnosed with a CMD did not get any professional help or treatment [34]. Our analysis strategy is designed to capture changes in the outcome variable over time. For this, a meaningful event, in this case, a diagnosis of anxiety-, mood/affective-, or stress-related disorders, is required to measure changes occurring before and after this date (T0 in our analysis). Had we not implemented this restriction, individuals entering the cohort could be at markedly different stages of illness, which would likely influence their level of work disability. As we had seven different combinations of subgroups of anxiety-, mood/affective-, or stress-related disorders, we could not examine separate trajectories for all these groups. Moreover, sex was not strongly associated with membership in any trajectory groups in our study; we therefore chose not to stratify on sex. Assessing sex differences within treatment for anxiety-, mood/affective-, or stress-related disorders would, however, be the scope of future research. Finally, the results should be interpreted in light of the fact that the exposures were measured at the time of the diagnosis, which means that the outcome partly preceded the exposures.

Conclusions

Whilst the majority of individuals treated for anxiety-, mood/affective-, or stress-related disorders in primary care seem to have a good prognosis concerning work disability, many individuals do experience increasing levels of work disability in the three years that follow. Those having a diagnosis of stress-related disorders, being aged 60 years and above, receiving psychotropic medication, having a diagnosis of a psychiatric disorder from an inpatient and specialised outpatient healthcare, and having a non-managerial occupation or being employed within the administration and customer service were associated with a higher likelihood of extended work disability. As such, these individuals could benefit from additional support early in the disease spell.

Data availability

The data used in this study cannot be publicly available due to privacy regulations. According to the General Data Protection Regulation, the Swedish law SFS 2018:218, the Swedish Data Protection Act, the Swedish Ethical Review Act and the Public Access to Treatment of common mental disorders in primary care. Readers may contact Kristina Alexanderson (kristina.alexanderson@ki.se) regarding the data.

References

  1. Steel Z, Marnane C, Iranpour C, Chey T, Jackson JW, Patel V, Silove D. The global prevalence of common mental disorders: a systematic review and meta-analysis 1980–2013. Int J Epidemiol. 2014;43(2):476–93.

    Article  PubMed  PubMed Central  Google Scholar 

  2. FORTE. Psykisk ohälsa, arbetsliv och sjukfråvaro - en kunskapsöversikt (Mental illness, work life and sickness absence - A knowledge overview). [In Swedish]. Swedish Research Council for health, Working Life and Welfare; 2015.

  3. Thisted CN, Nielsen CV, Bjerrum M. Work participation among employees with Common Mental disorders: a Meta-synthesis. J Occup Rehabil. 2018;28(3):452–64.

    Article  PubMed  Google Scholar 

  4. OECD. Sick on the job? Myths and realities about mental health and work. 2012.

  5. Helgesson M, Tinghog P, Wang M, Rahman S, Saboonchi F, Mittendorfer-Rutz E. Trajectories of work disability and unemployment among young adults with common mental disorders. BMC Public Health. 2018;18(1):1228.

    Article  PubMed  PubMed Central  Google Scholar 

  6. Lidwall U, Olsson-Bohlin C. Sjukskrivning för reaktioner på svår stress ökar mest (sick leave due to reactions on severe stress increases most) [In Swedish]. 2016.

  7. Marklund S, Gustafsson K, Aronsson G, Leineweber C, Helgesson M. Working conditions and compensated sickness absence among nurses and care assistants in Sweden during two decades: a cross-sectional biennial survey study. BMJ open. 2019;9(11):e030096.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Castro-Rodríguez JI, Olariu E, Garnier-Lacueva C, Martín-López LM, Pérez-Solà V, Alonso J, Forero CG. Diagnostic accuracy and adequacy of treatment of depressive and anxiety disorders: a comparison of primary care and specialized care patients. J Affect Disord. 2015;172:462–71.

    Article  PubMed  Google Scholar 

  9. Cipriani A, Furukawa TA, Salanti G, Chaimani A, Atkinson LZ, Ogawa Y, et al. Comparative efficacy and acceptability of 21 antidepressant drugs for the acute treatment of adults with major depressive disorder: a systematic review and network meta-analysis. Lancet. 2018;391(10128):1357–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Nigatu YT, Liu Y, Uppal M, McKinney S, Rao S, Gillis K, Wang J. Interventions for enhancing return to work in individuals with a common mental illness: systematic review and meta-analysis of randomized controlled trials. Psychol Med. 2016;46(16):3263–74.

    Article  CAS  PubMed  Google Scholar 

  11. Salomonsson S, Hedman-Lagerlöf E, Öst LG. Sickness absence: a systematic review and meta-analysis of psychological treatments for individuals on sick leave due to common mental disorders. Psychol Med. 2018;48(12):1954–65.

    Article  PubMed  Google Scholar 

  12. Nieuwenhuijsen K, Verbeek JH, Neumeyer-Gromen A, Verhoeven AC, Bültmann U, Faber B. Interventions to improve return to work in depressed people. Cochrane Database Syst Rev. 2020;10(10):Cd006237.

    PubMed  Google Scholar 

  13. Bortolotti B, Menchetti M, Bellini F, Montaguti MB, Berardi D. Psychological interventions for major depression in primary care: a meta-analytic review of randomized controlled trials. Gen Hosp Psychiatry. 2008;30(4):293–302.

    Article  PubMed  Google Scholar 

  14. Finnes A, Enebrink P, Ghaderi A, Dahl J, Nager A, Öst LG. Psychological treatments for return to work in individuals on sickness absence due to common mental disorders or musculoskeletal disorders: a systematic review and meta-analysis of randomized-controlled trials. Int Arch Occup Environ Health. 2019;92(3):273–93.

    Article  PubMed  Google Scholar 

  15. McKee-Ryan F, Song Z, Wanberg CR, Kinicki AJ. Psychological and physical well-being during unemployment: a meta-analytic study. J Appl Psychol. 2005;90(1):53–76.

    Article  PubMed  Google Scholar 

  16. Helgesson M, Johansson B, Nordqvist T, Lundberg I, Vingard E. Sickness absence at a young age and later sickness absence, disability pension, death, unemployment and income in native swedes and immigrants. Eur J Public Health. 2015;25(4):688–92.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Stockholm R, Region, Stockholm. Frågor och svar för forskare om hälsodatabaser (Questions and Answers for Researchers about Health Databases), In Swedish 2022 [ https://www.regionstockholm.se/om-regionstockholm/forskning-och-innovation/centrum-for-halsodata/fragor-och-svar-for-forskare-om-halsodatabaser/

  18. WHO. International Statistical Classification of Diseases and related health problems, 10 revision (ICD 10). WHO; 2010.

  19. Nagin D. Group-based modeling of development. Harvard University Press; 2005.

  20. Nagin DS, Tremblay RE. What has been learned from group-based trajectory modeling? Examples from physical aggression and other problem behaviors. Annals Am Acad Political Social Sci. 2005;602(1):82–117.

    Article  Google Scholar 

  21. Hansen CD, Andersen JH. Sick at work—a risk factor for long-term sickness absence at a later date? J Epidemiol Commun Health. 2009;63(5):397–402.

    Article  CAS  Google Scholar 

  22. Helgesson M, Gustafsson K, Leineweber C. Psychosocial and organisational work factors as predictors of sickness absence among professionally active adults with common mental disorders. BMC Psychiatry. 2023;23(1):543.

    Article  PubMed  PubMed Central  Google Scholar 

  23. von Nordheim F, Kvist J. Regulating the retirement age—lessons from nordic pension policy approaches. Regul Gov. 2023;17(3):644–57.

    Article  Google Scholar 

  24. Ávila-Gutiérrez MJ, Suarez-Fernandez de Miranda S, Aguayo-González F. Occupational Safety and Health 5.0—A Model for Multilevel Strategic Deployment Aligned with the Sustainable Development Goals of Agenda 2030. Sustainability. 2022;14(11):6741.

  25. Engström L-G, Janson S. Stress-related sickness absence and return to labour market in Sweden. Disabil Rehabil. 2007;29(5):411–6.

    Article  PubMed  Google Scholar 

  26. Welfare). STSNBoHa. Försäkringsmedicinskt beslutsstöd (Insurance Medical Decision Support) Stockholm2023 [ https://roi.socialstyrelsen.se/fmb

  27. Aronsson G, Marklund S, Leineweber C, Helgesson M. The changing nature of work - job strain, job support and sickness absence among care workers and in other occupations in Sweden 1991–2013. SSM Popul Health. 2021;15:100893.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Fagerlind Ståhl A-C, Ståhl C, Smith P. Longitudinal association between psychological demands and burnout for employees experiencing a high versus a low degree of job resources. BMC Public Health. 2018;18(1):915.

    Article  PubMed  PubMed Central  Google Scholar 

  29. (Försäkringskassan) SSA. Sjukfrånvaro i psykiatriska diagnoser - En registerstudie av Sveriges arbetande befolkning i åldern 20–69 år (Sickness absence in psychiatric diagnoses - A register study of Sweden’s working population aged 20–69). Stockholm; 2021. Report No.: 1654–8574.

  30. Cullen AE, Lindsäter E, Rahman S, Taipale H, Tanskanen A, Mittendorfer-Rutz E, Helgesson M. Patient factors associated with receipt of psychological and pharmacological treatments among individuals with common mental disorders in a Swedish primary care setting. BJPsych Open. 2023;9(2):e40.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Ludvigsson JF, Andersson E, Ekbom A, Feychting M, Kim J-L, Reuterwall C et al. External review and validation of the Swedish national inpatient register. BMC Public Health. 2011;11.

  32. Ludvigsson JF, Almqvist C, Bonamy AK, Ljung R, Michaelsson K, Neovius M, et al. Registers of the Swedish total population and their use in medical research. Eur J Epidemiol. 2016;31(2):125–36.

    Article  PubMed  Google Scholar 

  33. Lebano A, Hamed S, Bradby H, Gil-Salmerón A, Durá-Ferrandis E, Garcés-Ferrer J, et al. Migrants’ and refugees’ health status and healthcare in Europe: a scoping literature review. BMC Public Health. 2020;20(1):1039.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Kohn R, Saxena S, Levav I, Saraceno B. The treatment gap in mental health care. Bull World Health Organ. 2004;82(11):858–66.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

AFA Insurance, grant number 200061, supported this work. We used data from the REWHARD consortium funded by the Swedish Research Council (grant no. 2017 − 00624).

Open access funding provided by Karolinska Institute.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualisation: M.H., A.E.C., E.P.; methodology: M.H., A.E.C., E.P.; software: E.P; validation: M.H., A.E.C., E.P.; formal analysis: E.P, A.E.C., M.H; resources: M.H.; data curation: E.P.; writing original draft preparation: M.H., A.E.C., E.P; writing—review and editing: M.H., A.E.C., E.P; project administration: M.H.; funding acquisition: M.H. All authors have read and agreed to the published version of the manuscript.

Corresponding author

Correspondence to Magnus Helgesson.

Ethics declarations

Ethics approval and consent to participate

The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. An ethical permit was granted by the Regional Ethic Board in Stockholm. Informed consent is generally required from all subjects or their legal guardian(s) in studies, including human subjects. As this is a register study, the Regional Ethics Board waived the requirement to obtain personal consent from participants of this study.

Consent for publication

Not applicable.

Competing interests

Dr Cullen has received consultancy fees from Stratenym Inc and Symmetron Ltd. Drs Taipale and Mittendorfer-Rutz have received grants from Janssen paid to their employer. Dr Taipale reports personal fees from Gedeon Richter, Janssen, Lundbeck and Otsuka.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Helgesson, M., Pettersson, E., Lindsäter, E. et al. Trajectories of work disability among individuals with anxiety-, mood/affective-, or stress-related disorders in a primary healthcare setting. BMC Psychiatry 24, 623 (2024). https://doi.org/10.1186/s12888-024-06068-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12888-024-06068-5

Keywords