Dynamic simulation modelling is increasingly being recognised as a valuable decision-support tool to help guide investments and actions to address complex public health issues such as suicide. In particular, participatory system dynamics (SD) modelling provides a useful tool for asking high-level ‘what if’ questions, and testing the likely impacts of different combinations of policies and interventions at an aggregate level before they are implemented in the real world.
We developed an SD model for suicide prevention in Australia, and investigated the hypothesised impacts over the next 10 years (2015–2025) of a combination of current intervention strategies proposed for population interventions in Australia: 1) general practitioner (GP) training, 2) coordinated aftercare in those who have attempted suicide, 3) school-based mental health literacy programs, 4) brief-contact interventions in hospital settings, and 5) psychosocial treatment approaches.
Findings suggest that the largest reductions in suicide were associated with GP training (6%) and coordinated aftercare approaches (4%), with total reductions of 12% for all interventions combined. This paper highlights the value of dynamic modelling methods for managing complexity and uncertainty, and demonstrates their potential use as a decision-support tool for policy makers and program planners for community suicide prevention actions.
Suicide remains a leading cause of death in young people, and contributes significant social, economic and health system costs to Australia. The extent to which specific suicide prevention policies and programs over the past two decades have had an impact on population-level suicide rates in Australia is unclear.1 This is despite government and nongovernment investment, senate committee inquiries, implementation frameworks, and efforts to improve coordination of programs and services.2-5
The ability to design effective policy for suicide prevention in populations is challenged by the complex aetiology of suicide, limited generalisable evidence for interventions that prevent suicide in populations, and the quality and timeliness of surveillance data on suicidal behaviour (and other intermediate indicators) to measure the effect of suicide prevention initiatives.
There are also often incongruous priorities between researchers, policy makers and service providers. The time, scope and rigour required for conducting research to implement and evaluate combined multilevel approaches to suicide prevention using traditional epidemiological approaches often do not match with the priorities and timelines of service providers, policy makers and local communities, who usually require local or context-specific information.6,7
Although there is some evidence for effective prevention strategies in particular contexts8, it is unlikely that a single intervention will have population-level impacts on suicide rates, given the interplay of aetiological factors associated with suicidal behaviour within a complex health and social services system. Interventions that may have promising effects on preventing suicide in one context may not be generalisable to other contexts, and their likely effects over time and in combination are unknown. These issues can be addressed by acknowledging this complexity, and applying systems modelling tools and approaches to support decision making and guide prevention initiatives.6
This article addresses some of these challenges by combining objective data on suicidal behaviour with evidence reviews of suicide prevention strategies to develop a system dynamics (SD) model. This is a potential tool for policy makers to provide timely and population-specific information on the incidence of suicidal behaviour. The model demonstrates a potential platform for integrating diverse evidence sources into an analytic tool that can allow policy makers to explore the likely impacts of different policy and intervention scenarios over the short and longer terms in a robust, risk-free and low-cost way. That is, it can be used to conduct virtual experiments where real-world studies may not be feasible.
The use of system science modelling methods, such as SD (and also agent-based modelling and discrete event simulation), is becoming more prominent in the health sector because of the potential of these models to address policy making challenges for complex problems.7 These methods have a long history in engineering, economics, defence and ecology, and have more recently been useful in appropriately characterising the complex nature of the health sector – including evaluating operational aspects of healthcare capacity and demand, patient flows, and disease screening7, as well as a range of disease risk factors and outcomes.6,9,10 An SD approach has not previously been applied to the topic of suicide; however, the complex aetiology of suicidal behaviour and the range of potentially interacting prevention strategies suggest that this approach may be useful to inform policy responses.
We developed an SD model to reflect a series of suicidal states in the Australian population, ranging from no suicidal thoughts or behaviours, to suicidal ideation, to the planning of suicidal acts, to attempted suicide and completed suicide (Figure 1; see also Supplementary Figure 1 for model structure, available from: www.saxinstitute.org.au/tables-suicide-prevention-tool). This progression from ideation to behaviour, modified by level of intent and lethality of method, is broadly consistent with current conceptualisations of suicidal behaviour, although debate remains about the nomenclature of the range of suicidal behaviour that occurs in populations.11 Conceptualising the pathway of suicidal behaviour in this manner also enabled representative data sources to be used to validate the model on the two key outcomes of interest: cases of hospital-treated attempted suicide and cases of suicide.
Figure 1. Summary structure of the system dynamics model of Australian suicide (click to enlarge)
The model was developed using historical time series of hospital-treated attempted suicide and cases of suicide from 2000 to 2014, for males and females separately. Models were stratified by sex to reflect the three- to four-fold difference in the incidence of suicide between males and females. Nationally representative estimates of suicidal ideation and suicide planning for males and females were obtained from the most recent National Survey of Mental Health and Wellbeing12, and the level of suicidal intent for suicide attempts was obtained from the Australian National Survey of Self Injury13 (Supplementary Tables 1–4, available from: www.saxinstitute.org.au/tables-suicide-prevention-tool). The most common methods of suicidal behaviour for suicide attempts and suicides were also incorporated, as a reflection of the differing methods used by males and females, and also the differing levels of suicidal intent and lethality. These methods comprised ‘self-poisoning’ (International Statistical Classification of Diseases and Related Health Problems [ICD-10] codes X60–X69), ‘hanging’ (X70), ‘drowning’ (X71), ‘firearms’ (X72–X74), ‘cutting’ (X78), ‘jumping from a high place’ (X80) and ‘other’ (X75–X77, X79, X81–X84).
A recent systematic review of peer-reviewed intervention studies was used to extract intervention effects for different health service–related suicide prevention initiatives14 (Supplementary Table 5, available from: www.saxinstitute.org.au/tables-suicide-prevention-tool). The strategies described in this review were selected because they are currently the subject of a large-scale community trial (www.lifespan.org.au), and have been the recent focus of debate among population health researchers15,16 and of government funding strategies. The population preventive fractions reported in this review were used as inputs to intervention mechanisms that modify net rates of movement of the population from one state to another. The resulting estimated number of suicides averted (as described below) incorporates these hypothesised intervention mechanism effects, which include sources of real-world inertia and delay (such as time to scale up interventions, service provider adoption, patient uptake of interventions), and their impact on feedback flows in the model.
Three approaches were considered to assess the hypothesised impact of the interventions described above over a forecast period of 2015–2025. First, each of the service interventions was considered separately to forecast its impact on attempted suicide and suicide. Second, all the service interventions were considered in combination. Third, hypothetical changes in the lethality of each of the selected methods of suicidal behaviour were considered. For illustrative purposes, these model runs assumed a 20% absolute reduction in lethality for each of the methods described above – that is, to assess the hypothesised impact on subsequent suicide if the lethality of (for example) hanging could be reduced by 20%. Finally, the combined effects of the health service interventions plus the combined effects of reductions in lethality of method were considered.
The interactive model can be accessed online (saxinstitute.scem.uws.edu.au/netsims/suicide-1/suicide_v403/index.html) to investigate alternative scenarios.
The model estimated a historical time series (2000–2014) for males and females in the Australian population (Figure 2A) that was within 5% of the observed suicide data for an equivalent period. The model forecast an increase in suicide frequency from 2015 to 2025 from 2244 to 2861 for males and 713 to 935 for females.
Comparative runs of the model investigating the separate impacts of interventions suggest that general practitioner (GP) training in the identification and management of patients presenting with risk factors for suicide (such as mental illness) or suicidal ideation, and coordinated assertive aftercare for those who present to hospital services following a suicide attempt result in the largest proportion of suicides prevented over the forecast period. For GP training, approximately 6% of suicide cases were prevented in males and females (Figure 2B; Table 1). The other interventions had negligible impacts on suicide trends when simulated individually. In combination, the estimated proportion of prevented suicides for all service interventions was 12% for both males and females (Figure 2B; Table 1).
Simulated reductions in lethality of method also resulted in declines in suicide over the forecast period (2015–2025). For both males and females, the largest proportion of prevented suicides was for hanging (7% in males; 5% in females) (Figure 2C; Table 1). In combination, the estimated proportion of prevented suicides for a reduction in lethality of method of 20% for all of the selected methods was 10% for males and females (Figure 2C; Table 1).
In models combining the effects of all health service interventions plus reductions in lethality of method, the estimated proportion of prevented suicides was 22% for both males and females (Figure 2D; Table 1).
Figure 2. Model outputs for Australian suicide (2001–2025) and comparative model runs for the forecast period 2015–2025 (click to enlarge)
Table 1. Summary of estimated reductions in suicide from 2015 to 2025 using combinations of selected interventions and reductions in lethality of method
Males | Females | ||||||
Measure | Type | Cumulative casesa | Cases averted | % reduction | Cumulative casesa | Cases averted | % reduction |
Total suicides (no intervention) | na | 28 791 | 0 | 0.0 | 9 289 | 0 | 0.0 |
Service interventions | General practitioner training | 27 108 | 1 684 | 5.8 | 8 741 | 548 | 5.9 |
Coordinated aftercare | 27 550 | 1 241 | 4.3 | 8 877 | 413 | 4.4 | |
Brief-contact intervention | 28 698 | 93 | 0.3 | 9 255 | 35 | 0.4 | |
School programs | 28 431 | 361 | 1.3 | 9 171 | 119 | 1.3 | |
Psychosocial programs | 28 140 | 652 | 2.3 | 9 068 | 222 | 2.4 | |
Combined service interventions | 25 358 | 3 434 | 11.9 | 8 153 | 1 137 | 12.2 | |
20% reduction in method lethality | Poisoning | 28 418 | 374 | 1.3 | 8 951 | 338 | 3.6 |
Hanging | 26 777 | 2 015 | 7.0 | 8 793 | 497 | 5.3 | |
Gases and vapours | 28 483 | 309 | 1.1 | 9 222 | 68 | 0.7 | |
Firearms | 28 514 | 277 | 1.0 | 9 274 | 15 | 0.2 | |
Combined methods above | 25 812 | 2 979 | 10.3 | 8 371 | 918 | 9.9 | |
Combined service interventions and reduction in lethality | na | 22 572 | 6 219 | 21.6 | 7 243 | 2 046 | 22.0 |
This paper describes the application of an SD model of suicide prevention in Australian males and females, which also incorporates current best evidence relating to strategies for prevention.14 In comparisons of the simulated effects of prevention strategies, the largest estimated impacts on declines in suicide were for GP training17 and coordinated aftercare.18 Simulated reductions in lethality of method by 20% for selected methods of suicidal behaviour (that is, assuming that interventions were available to reduce lethality of method) also suggested that the largest declines would be achieved with reductions in the lethality of hanging.
The National Mental Health Commission2 has called for a 50% reduction in the number of suicides and suicide attempts over the next decade. Based on the interventions currently being considered as the best approaches to prevent suicide, findings suggest that this 50% reduction will not be achieved by these strategies alone. Additionally, the methods that are most amenable to intervention (such as self-poisoning, jumping and firearms) appear to have the smallest impact on suicide declines, whereas hypothesised changes to lethality of hanging – a method that is not amenable to intervention – had the largest estimated impact on suicide declines over the forecast period.
Although the SD model is able to account for behavioural and population dynamics that affect population-level changes over time in intentional self-harm, aggregate models such as these are an oversimplification that do not consider the complex individual mechanisms and trajectories associated with suicidal behaviour. Agent-based models may be more appropriate in capturing these mechanisms. However, the application of SD applies a level of mathematical rigour that has not previously been applied to trends in suicide, and is consistent with observed data inputs at the aggregate population level.
The model is also calibrated on routinely collected data on mortality and hospital admissions. There is some underenumeration of suicide cases due to the misclassification of suicides to ICD codes relating to unintentional injury and events of ‘undetermined intent’; however, this does not preclude the use of these data, particularly in studies of time-series analyses. Additionally, hospital admissions data are likely to underestimate the total population level of attempted suicide, in that the data on suicide attempts capture only those cases that are serious enough to warrant medical intervention.
The interventions incorporated into the model represent a limited set of potentially important interventions. The extent to which the underlying relative risk estimates are applicable to the Australian and other country contexts warrants consideration, as do the representativeness and precision of the prevalence estimates used to calculate the population preventive fractions. However, in the context of using SD models to inform the prioritisation of evidence based interventions, SD models provide hypothesised effects that can subsequently be assessed using empirical data. Other evidence on potentially effective interventions from observational studies (for example, relating to educational achievement19, unemployment20 or early intervention for mental illness21) may be additional areas for consideration in subsequent dynamic models of suicidal behaviour. SD models can be calibrated for any given geographic context or subpopulation of interest, and are strengthened through participatory model-building exercises involving local stakeholders – a central principle that ensures that the resulting decision-support tool is directly relevant to decision makers.6
This article demonstrates how an SD model can be used as a potential decision-support tool for policy makers to ask ‘what if’ questions, and to test the potential impacts of intervention and policy scenarios over time before implementing solutions in the real world. Recent advances in modelling software and transparent participatory approaches are enabling suicide prevention researchers, policy makers, practitioners, consumer representatives and modellers to collaboratively develop such models and draw on diverse evidence sources. This offers promise in providing tools that are better aligned to the decision-support needs of policy makers and facilitate consensus building for action. This modelling tool is proposed as an adjunct to standard program design and evaluation approaches, and does not preclude the use of prospective evaluation studies based on primary data collection to address deficiencies in local data and assess joint effects of intervention combinations. Rather, SD models provide an additional way to engage stakeholders, and prioritise resources and interventions in the design and evaluation of suicide prevention initiatives.
This article is based on research funded by a Western Sydney University Research Partnership Grant.
Copyright:
© 2017 Page et al. This work is licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International Licence, which allows others to redistribute, adapt and share this work non-commercially provided they attribute the work and any adapted version of it is distributed under the same Creative Commons licence terms.