- Research article
- Open Access
- Open Peer Review
Changes in monthly unemployment rates may predict changes in the number of psychiatric presentations to emergency services in South Australia
BMC Emergency Medicinevolume 15, Article number: 16 (2015)
To determine the extent to which variations in monthly Mental Health Emergency Department (MHED) presentations in South Australian Public Hospitals are associated with the Australian Bureau of Statistics (ABS) monthly unemployment rates.
Times series modelling of relationships between monthly MHED presentations to South Australian Public Hospitals derived from the Integrated South Australian Activity Collection (ISAAC) data base and the ABS monthly unemployment rates in South Australia between January 2004–June 2011.
Time series modelling using monthly unemployment rates from ABS as a predictor variable explains 69 % of the variation in monthly MHED presentations across public hospitals in South Australia. Thirty-two percent of the variation in current month’s male MHED presentations can be predicted by using the 2 months’ prior male unemployment rate. Over 63 % of the variation in monthly female MHED presentations can be predicted by either male or female prior monthly unemployment rates.
The findings of this study highlight that even with the relatively favourable economic conditions, small shifts in monthly unemployment rates can predict variations in monthly MHED presentations, particularly for women. Monthly ABS unemployment rates may be a useful metric for predicting demand for emergency mental health services.
With Europe experiencing a major economic downturn accompanied by marked increases in unemployment, there have been concerns raised about the impact on population mental health and suicide rates [1–6].
Families and communities experience poorer mental health from higher rates of unemployment [7, 8]. Generally, the negative impact of unemployment on mental health is more evident among men than women . At a sociological level an increase in unemployment rate is a marker of social fragmentation and thereby contributes to mental health issues . This association has been long established, and in 1897 Durkheim  argued that employment protected an individual against suicide by better integrating the individual with their society. Australia, unlike most Western nations, did not go into recession following the Global Financial crisis in 2007. However, there was an economic slow-down with an accompanying rise in unemployment from a low of 4.0 % in early 2008 to a peak of 5.8 % in late 2009 before declining to 5.0 % in April 2012 . Emergency Departments (ED) in Australia have significant numbers of patients presenting with mental health problems as their primary problem, with this group accounting for between 4.5 and 5.6 % of all ED presentations . There is evidence most people who receive Emergency Department Mental Health (MHED) care are unemployed  and ED might be their preferred pathway to crisis care during economic slowdowns.
While time-trend analysis has been used to analyse associations between unemployment and suicide , no study has used this method to analyse possible links between unemployment and MHED presentations. We were interested to determine whether in a country with relatively low levels of unemployment, the number of MHED presentations was associated with unemployment rates and whether this varied by gender. In this study, we analysed the effects of changes in rates of unemployment in South Australia on numbers of MHED presentations using a time-trend analysis between 2004 and 2011.
Seasonally adjusted monthly unemployment rates for South Australia from January 2004 to December 2011 were sourced from ABS (Catalogue 6202.0–Labour Force, Australia, May 2012) . The total number of people presenting to South Australian tertiary hospital EDs during the same period were sourced from the SA Health database, ISAAC, (Integrated South Australian Activity Collection). In South Australia, there are no separate psychiatric emergency departments and all mental health/psychiatric emergency presentations occur at general hospital emergency departments. All presentations with ICD codes in the range F00–F99 were categorized as psychiatric ED presentations. Gender was determined from demographic variables of each presentation. The rest were categorized as physical illness presentations.
The number of mental health presentations to ED was organized into 3 outcome series–all mental health presentations, female mental health presentations and male mental health presentations. Each of these time series was seasonally adjusted to remove seasonal variations of ED presentations using the ratio-to-moving-average Seasonal Decomposition method. All subsequent analysis was then performed on seasonally adjusted series.
As data was extracted from de identified administrative data sets, ethics approval was not sought.
The analysis involved a combination of cross correlation and ARIMA modelling. The first step of analysis involved a series of cross correlation computations between predictor (total unemployed, male unemployed and female unemployed) and outcome (MHED total presentations, male MHED presentations, female MHED presentations) time series variables using “arima” and “ccf” from the package “stats” in R 2.13 (http://www.r-project.org/). The residuals of both series were used as input for the cross-correlation analysis. Cross correlation coefficients between output series and input series with a lag k at time t were considered significant if the correlation coefficients were outside the standard error margin of 0.2.
We then used the ARIMA function in the SPSS (Chicago, IL, USA) Forecasting Module to derive a series of ARIMA models with unemployment as the predictor variable and mental health ED presentations as the outcome variables, by specifying auto regression (AR) with the lags that were shown to be significant in cross correlation analysis, with either undifferenced or first differenced Integration (I) components and Moving Average (MA) components of either  or [1,2]. We selected the best model of the above.
Automatic detection of outliers was made and the outliers were modelled accordingly, thus trimming was not performed. The stationary R-squared measurement determined was used as an indicator of goodness of fit for the models. The Ljung-Box Q provided the diagnostic statistics to check the presence of structure in the observed series that was not accounted for by the model. Significance values less than 0.01 were considered indicating the presence of structure in the observed series, which was not accounted for by the model; therefore, we ignored the model if it had a significant value less than 0.01. The Bayesian Information Criteria (BIC) is also reported as a comparative goodness of fit measure for each of the models derived. We also used SPSS version 19.0 Expert Modeler (Chicago, IL, USA) to automatically determine the best-fitting univariate ARIMA models for each of the outcomes time series variables that had shown significant cross correlation with predictor variables.
We compared the metrics of the models using unemployment as predictor variables with the best-fitting univariate ARIMA model derived by the Expert Modeller.
Unemployment and physical illness presentations to ED
We initially looked at associations between unemployment rate and physical illness presentations to ED (see Fig. 1). A cross correlation of prewhitened data showed no significant correlation between these two time series for varying lags of unemployment.
Unemployment and mental health presentations to ED
Non-mental health ED presentations ranged from 22,450 to 33,697 per month (Mean 28,282.8; SE 234.40). MHED presentations ranged from 978 to 1421 per month (Mean 1226.5; SE 10.11). Table 1 shows the MHED diagnostic groups with affective disorders being the most frequent grouping.
We initially used the SPSS Expert modeller to determine the best univariate ARIMA model for predicting Mental Health ED presentations. The best univariate model was ARIMA (0,1,3) with stationary-R-squared of 0.36 and Ljung-Box p = 0.70, and lowest BIC = 7.81. We then examined any association between total MHED presentations and unemployment rate (Fig. 2). A cross correlation of prewhitened series showed significant correlation between the unemployment and presentations, at current month (Correlation coeff = 0.22 (SE = 0.2)) and a lag of 2 months (Correlation coeff = 0.36 (SE = 0.2)).
We then derived transfer function ARIMA models with MHED presentations as the dependent variable and unemployment as a predictor variable. We tested several different models by specifying lags of 0 and 2 months as the AR, and varying the combinations of I = 0 or 1 and MA = 0 or [1,2]. The best model was ARIMA (AR = 0, I = 0, MA = [1,2]) that reached the highest stationary R-Squared of 0.69, Ljung-Box p = 0.19 and also the lowest BIC of 8.20. This analysis demonstrated a significant improvement in explaining variations of Mental Health ED presentations was achieved by introducing unemployment as a predictor variable, the stationary R-Squared improving from 0.36 to 0.69.
Effect of gender
Male mental health presentations to ED
We then looked at the effect of gender. The best univariate model for predicting male MHED presentations was an ARIMA (0,1,1) that had stationary R-squared of 0.27 and Ljung-Box p = 0.16, BIC = 7.26.
We repeated cross-correlations analysis and observed significant associations between male unemployment rates and MHED presentations of males, most significantly with 2 months prior male unemployment rates (corr coeff = 0.36 (SE = 0.2)). The best ARIMA model using male unemployment rate with lag 2 as the predictor variable was an ARIMA (AR = 2, I = 1, MA = [1,2]) with stationary R-Squared of 0.32 and Ljung-Box p = 0.13, BIC = 7.42. This analysis shows that by introducing a male unemployment predictor variable with a lag of 2 months, for male MHED presentations, the stationary R-Squared was improved from 0.27 to 0.32.
A cross correlation analysis between female unemployment and male MHED presentations showed no significant association.
Female mental health presentations to ED
The best univariate model for predicting female MHED presentations was a simple univariate model that had stationary R-squared of 0.17 and Ljung-Box p = 0.05, BIC = 7.27. The cross-correlation analysis between female unemployment rates and female MHED presentations showed significant correlations with both current (Corr Coeff = 0.23; SE = 0.2) and previous months (Corr Coeff = 0.29, SE = 0.2) female unemployment rates. The best ARIMA model for predicting female MHED presentations using current month’s female unemployment rate as the independent predictor variable was (ARIMA (0,0,0)) with a stationary R-Squared of 0.63 and Ljung-Box p = 0.37, and lowest BIC = 7.00. By using the previous month’s female unemployment rate as predictor variable the best model observed was ARIMA (1,1,0) with stationary R-Squared of 0.50 and Ljung-Box p = 0.48, BIC = 7.20.
We also examined whether there was an impact of male unemployment on female MHED presentations. A cross correlation showed significant correlations between the female presentations and lag 2 of male unemployment (CorrCoeff = 0.24, SE = 0.2). The best ARIMA models for predicting female MHED presentations using 2 months’ prior male unemployment rate as the independent predictor variable was an ARIMA (2,0,0) model with a stationary R-Squared of 0.64 and Ljung-Box p = 0.95, BIC = 7.10 (Fig. 3).
This analysis demonstrated that by introducing either the current or previous month’s female unemployment rate or 2 months’ prior male unemployment rate as a predictor variable, a significant improvement was achieved in the ability to predict female MHED presentations.
Our findings suggest that unemployment might be an important social determinant increasing the number of visits to the medical emergency room for reasons of mental health. Unemployment rates were predictive of MHED presentations for both men and women in South Australia. Most of the association between unemployment rates and MHED presentations, however, was attributable to changes in ED use by women. While there is no clear evidence from life-events literature  that women are more sensitive than men to adverse life-events like unemployment, it is well established that women are more likely to seek help ; women have higher rates of depression ; and there is an over-representation of women with nonfatal suicidal behaviour presenting to EDs . Women show more ability to empathize and this tendency might increase their liability for emotional symptoms in the presence of acute disappointments such as a family crisis arising from unemployment .
Before we address the possible policy implications of this study we need to consider several of its limitations. Firstly, correlation tests assess associations and do not establish causality. Thus, the association between unemployment and ED mental health presentations might be an indirect one, caused by other unidentified factors. However the study was carried out in an economy with close to full employment, so other factors related to economies in recession would not have confounded the associations found. Furthermore, while statistically significant, whether or not the effect of the association is clinically significant is not clear from our study. A further limitation relates to the fact that not all diagnoses used in the analysis were made by psychiatrists, some were made by ED physicians. Finally, the study was not designed to explore how presentation of different psychiatric diagnostic categories varied over time.
With an economy not in recession, these findings indicate that relatively small increases in unemployment rates might predict observed increases in emergency mental health demand. Unemployment data might therefore be useful for predicting the requirements for MHED service provision. Emergency responses may require more flexible health care and social service responses and funding arrangements. If further research confirms our finding then an argument for hospitals to plan ahead for increasing demand for MHED services when unemployment rates rise, akin to hospitals planning for increased admissions during wintertime could be developed .
The findings from this study highlight that even with the relatively favourable economic conditions in Australia, small shifts in unemployment are associated with MHED presentations, particularly for women. Further studies should examine the effect of unemployment for spouses and on family systems since our study suggests that fluctuations in male unemployment can impact on MHED presentations for women.
Barr B, Taylor-Robinson D, Scott-Samuel A, McKee M, Stuckler D. Suicides associated with the 2008–10 economic recession in England: time trend analysis. BMJ. 2012;345, e5142.
Stuckler D, Basu S, Suhrcke M, Coutts A, McKee M. Effects of the 2008 recession on health: a first look at European data. Lancet. 2011;378(9786):124–5.
Kentikelenis A, Karanikolos M, Papanicolas I, Basu S, McKee M, Stuckler D. Health effects of financial crisis: omens of a Greek tragedy. Lancet. 2011;378(9801):1457–8.
Stuckler D, Basu S, Suhrcke M, Coutts A, McKee M. The public health effect of economic crises and alternative policy responses in Europe: an empirical analysis. Lancet. 2009;374(9686):315–23.
Chang SS, Stuckler D, Yip P, Gunnell D. Impact of 2008 global economic crisis on suicide: time trend study in 54 countries. BMJ. 2013;347:f5239.
Haw C, Hawton K, Gunnell D, Platt S. Economic recession and suicidal behaviour: possible mechanisms and ameliorating factors. Int J Soc Psychiatry. 2015;61(1):73–81.
Bloomer E, Allen J, Donkin A, Findlay G, Gamsu M. The impact of the economic downturn and policy changes on health inequalities in London. London: UCL Institute of Health Equity; 2012.
Rook K, Dooley D, Catalano R. Stress transmission: the effects of husbands’ job stressors on the emotional health of their wives. J Marriage Fam. 1991;53(1):165–77.
Paul KI, Moser K. Unemployment impairs mental health: meta-analyses. J Vocat Behav. 2009;74(3):264–82.
Blakely TA, Collings SC, Atkinson J. Unemployment and suicide. Evidence for a causal association? J Epidemiol Community Health. 2003;57(8):594–600.
Durkheim E. Suicide New York. Free Press, 1897.
Australian Bureau of Statistics. 1350.0–Australian economic indicators. 2012.
Knott J, Pleban A. Mental health presentations to the ED department. Melbourne, Victoria, Australia: Victorian Government Department of Human Services; 2005.
Australian Bureau of Statistics. 6202.0–labour force, Australia. 2012. http://www.abs.gov.au/ausstats/abs@.nsf/mf/6202.0.
Kendler KS, Thornton LM, Prescott CA. Gender differences in the rates of exposure to stressful life events and sensitivity to their depressogenic effects. Am J Psychiatry. 2001;158(4):587–93.
Parker G, Brotchie H. Gender differences in depression. Int Rev Psychiatry. 2010;22(5):429–36.
Kessler RC, Berglund P, Demler O, Jin R, Merikangas KR, Walters EE. Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the National Comorbidity Survey Replication. Arch Gen Psychiatry. 2005;62(6):593–602.
Schrijvers DL, Bollen J, Sabbe BG. The gender paradox in suicidal behavior and its impact on the suicidal process. J Affect Disord. 2012;138(1–2):19–26.
Goldberg D. The aetiology of depression. Psychol Med. 2006;36(10):1341–7.
Fullerton KJ, Crawford VL. The winter bed crisis–quantifying seasonal effects on hospital bed usage. QJM. 1999;92(4):199–206.
Kiernan M, Toro PA, Rappaport J. Economic predictors of mental health service utilization: a time-series analysis. Am J Community Psychol. 1989;17(6):801–20.
Jin RL, Shah CP, Svoboda TJ. The impact of unemployment on health: a review of the evidence. CMAJ. 1995;153(5):529–40.
We would like to acknowledge Megan Winsall for her assistance in the formatting of this paper and Mental Health Informatics Unit, Country Health SA, SA Health for facilitating access to data from ISAAC database. The project did not receive external sources of funding.
The authors declare that they have no competing interests.
NB, TB, GS, RA and JS contributed to the conception of the study, analysis and writing. GT contributed to the time series statistical analysis method. CP and SA contributed to the writing of the manuscript. All authors read and approved the final manuscript.