Hostname: page-component-7c8c6479df-995ml Total loading time: 0 Render date: 2024-03-29T08:11:07.353Z Has data issue: false hasContentIssue false

Understanding the characteristics and mechanisms underlying suicide clusters in Australian youth: a comparison of cluster detection methods

Published online by Cambridge University Press:  06 August 2020

N.T.M. Hill*
Affiliation:
Orygen, Centre for Youth Mental Health, University of Melbourne, Melbourne, Australia Telethon Kids Institute, Perth, Australia
L.S. Too
Affiliation:
Centre for Mental Health, Melbourne School of Population and Global Health, University of Melbourne, Melbourne, Australia
M.J. Spittal
Affiliation:
Centre for Mental Health, Melbourne School of Population and Global Health, University of Melbourne, Melbourne, Australia
J. Robinson
Affiliation:
Orygen, Centre for Youth Mental Health, University of Melbourne, Melbourne, Australia
*
Author for correspondence: Nicole T.M. Hill, E-mail: nicole.hill@telethonkids.org.au
Rights & Permissions [Opens in a new window]

Abstract

Aims

There is currently no gold-standard definition or method for identifying suicide clusters, resulting in considerable heterogeneity in the types of suicide clusters that are detected. This study sought to identify the characteristics, mechanisms and parameters of suicide clusters using three cluster detection methods. Specifically, the study aimed to: (1) determine the overlap in suicide clusters among each method, (2) compare the spatial and temporal parameters associated with different suicide clusters and (3) identify the demographic characteristics and rates of exposure to suicide among cluster and non-cluster members.

Methods

Suicide data were obtained from the National Coronial Information System. N = 3027 Australians, aged 10–24 who died by suicide in 2006–2015 were included. Suicide clusters were determined using: (1) poisson scan statistics, (2) a systematic search of coronial inquests and (3) descriptive network analysis. These methods were chosen to operationalise three different definitions of suicide clusters, namely clusters that are: (1) statistically significant, (2) perceived to be significant and (3) characterised by social links among three or more suicide descendants. For each method, the demographic characteristics and rates of exposure to suicide were identified, in addition to the maximum duration of suicide clusters, the geospatial overlap between suicide clusters, and the overlap of individual cluster members.

Results

Eight suicide clusters (69 suicides) were identified from the scan statistic, seven (40 suicides) from coronial inquests; and 11 (37 suicides) from the descriptive network analysis. Of the eight clusters detected using the scan statistic, two overlapped with clusters detected using the descriptive network analysis and one with clusters identified from coronial inquests. Of the seven clusters from coronial inquests, four overlapped with clusters from the descriptive network analysis and one with clusters from the scan statistic. Overall, 9.2% (12 suicides) of individuals were identified by more than one method. Prior exposure to suicide was 10.1% (N = 7) in clusters from the scan statistic, 32.5% (N = 13) in clusters from coronial inquest and 56.8% (N = 21) in clusters from the descriptive network analysis.

Conclusion

Each method identified markedly different suicide clusters. Evidence of social links between cluster members typically involved clusters detected using the descriptive network analysis. However, these data were limited to the availability information collected as part of the police and coroner investigation. Communities tasked with detecting and responding to suicide clusters may benefit from using the spatial and temporal parameters revealed in descriptive studies to inform analyses of suicide clusters using inferential methods.

Type
Original Articles
Creative Commons
Creative Common License - CCCreative Common License - BYCreative Common License - NCCreative Common License - ND
This is an Open Access article, distributed under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives licence (http://creativecommons.org/licenses/by-nc-nd/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is unaltered and is properly cited. The written permission of Cambridge University Press must be obtained for commercial re-use or in order to create a derivative work.
Copyright
Copyright © The Author(s), 2020. Published by Cambridge University Press

Introduction

Suicide clusters involve multiple suicides that occur closer in time or place than would normally be expected using statistical inference or community expectation (Hawton et al., Reference Hawton, Hill, Gould, John, Lascelles and Robinson2019). However, multiple suicides that are connected socially are also taken very seriously and are included in recent public health guidelines for responding to suicide clusters (Public Health England, Reference Public Health England2015; Palmer et al., Reference Palmer, Inder, Shave and Bushnell2018; Public Health England, 2019). Young people aged 10–24 years are two to four times more likely to be involved in a suicide cluster, compared to adults in the general population (Gould et al., Reference Gould, Wallenstein, Kleinman, O'Carroll and Mercy1990b; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016). Whilst early studies suggest between 1 and 13% of suicides in young people occur as part of a suicide cluster (Gould et al., Reference Gould, Wallenstein, Kleinman, O'Carroll and Mercy1990b), recent estimates range between 5.2 and 57% (Robertson et al., Reference Robertson, Skegg, Poore, Williams and Taylor2012; Cheung et al., Reference Cheung, Spittal, Williamson, Tung and Pirkis2013; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016). However, these differences appear to depend largely on the type of parameters and methods used to detect suicide clusters.

To date, several high income countries have implemented surveillance systems for the detection of suicide clusters (Clinical Advisory Services Aotearoa, 2016; Griffin et al., Reference Griffin, Dillion, Arensman, Corcoran, Williamson and Perry2017; Public Health England, 2020) as well as public health guidelines for the response and prevention of suicide clusters in the community (Centers for Disease Control, 1988; Health Service Executive (HSE), 2011; Centre for Health Policy Programs and Economics, 2012; Public Health England, Reference Public Health England2015; Palmer et al., Reference Palmer, Inder, Shave and Bushnell2018; Public Health England, 2019). Yet, despite significant public health investment, there is currently no gold-standard definition or method for detecting suicide clusters (Niedzwiedz et al., Reference Niedzwiedz, Haw, Hawton and Platt2014). For example, existing public health guidelines define a suicide cluster on the basis of both statistical inference and community expectation and note that suicide clusters may or may not involve social links between cluster members (Centers for Disease Control, 1988; Health Service Executive (HSE), 2011; Centre for Health Policy Programs and Economics, 2012; Public Health England, Reference Public Health England2015; Palmer et al., Reference Palmer, Inder, Shave and Bushnell2018; Public Health England, 2019).

The absence of a gold-standard definition of suicide clusters has resulted in considerable methodological heterogeneity in the way suicide clusters are operationalised and detected. For example, early studies of suicide clusters described the relationships between suicide descendants and found that many cluster members shared social links as friends or acquaintences (Bechtold, Reference Bechtold1988; Davidson et al., Reference Davidson, Rosenberg, Mercy, Franklin and Simmons1989; Wilkie et al., Reference Wilkie, Macdonald and Hildahl1998; Poijula et al., Reference Poijula, Wahlberg and Dyregrov2001; Wissow et al., Reference Wissow, Walkup, Barlow, Reid and Kane2001). This led to the common cited hypothesis that suicide contagion, more accurately known as the social transmission of suicidal behaviour (whereby exposure to suicide facilitates suicidal behaviour in others), is a key mechanism underlying the development of suicide clusters (Hawton et al., Reference Hawton, Hill, Gould, John, Lascelles and Robinson2019).

In the past two decades, there has been a methodological shift from descriptive studies of suicide clusters to inferential methods such as Poisson scan statistics. These studies used population data from suicide registries combined with geoinformation systems to determine whether suicides are greater than statistically expected within a particular time and place (Gould et al., Reference Gould, Wallenstein and Kleinman1990a; Williamson et al., Reference Williamson, Tung, Pirkis, Derek Cheung and Spittal2014; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016; Sy et al., Reference Sy, Shaman, Kandula, Pei, Gould and Keyes2019). Whilst early descriptive studies of suicide clusters were reported in small community settings such as a schools (Poijula et al., Reference Poijula, Wahlberg and Dyregrov2001), inpatient units (Taiminen et al., Reference Taiminen, Kallio-Soukainen, Nokso-Koivisto, Kaljonen and Helenius1998) and remote indigenous communities (Bechtold, Reference Bechtold1988; Wilkie et al., Reference Wilkie, Macdonald and Hildahl1998; Wissow et al., Reference Wissow, Walkup, Barlow, Reid and Kane2001), the shift towards inferential studies of suicide clusters has established the presence of suicide clusters in large nationwide studies (Gould et al., Reference Gould, Wallenstein and Kleinman1990a, Reference Gould, Wallenstein, Kleinman, O'Carroll and Mercyb; Cheung et al., Reference Cheung, Spittal, Pirkis and Yip2012; Jones et al., Reference Jones, Gunnell, Platt, Scourfield, Lloyd, Huxley, John, Kamran, Wells and Dennis2013; Williamson et al., Reference Williamson, Tung, Pirkis, Derek Cheung and Spittal2014; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016; Sy et al., Reference Sy, Shaman, Kandula, Pei, Gould and Keyes2019).

Despite being methodologically very different, evidence from both descriptive and inferential studies of suicide clusters are often treated as synonymous. For example, population-based studies of suicide clusters suggest that the close temporal and spatial proximity of suicides within a suicide cluster are a proxy measure of social contiguity among cluster members (Gould et al., Reference Gould, Wallenstein and Kleinman1990a, Reference Gould, Wallenstein, Kleinman, O'Carroll and Mercyb; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016; Sy et al., Reference Sy, Shaman, Kandula, Pei, Gould and Keyes2019). However, the ecological design of these studies means that it is unclear whether individuals were actually exposed to the suicide of another or shared social links with cluster members. Significant differences in the size and duration of suicide clusters are also evident between methods. Descriptive studies of suicide clusters range between 2 and 11 suicides over a maximum 5-year period (Niedzwiedz et al., Reference Niedzwiedz, Haw, Hawton and Platt2014). In contrast, suicide clusters reported using the scan statistic range between 3 and 1500 suicides over a maximum to 2.5-year period (Cheung et al., Reference Cheung, Spittal, Williamson, Tung and Pirkis2013; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016; Sy et al., Reference Sy, Shaman, Kandula, Pei, Gould and Keyes2019). In some instances, inferential studies have identified suicide clusters that span multiple states and territories (Cheung et al., Reference Cheung, Spittal, Williamson, Tung and Pirkis2013; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016; Sy et al., Reference Sy, Shaman, Kandula, Pei, Gould and Keyes2019).

Distinguishing between the outcomes of different cluster detection methods has important implications for informing the way that communities detect suicide clusters. For example, previous studies of suicide clusters using the scan statistic have selected specific spatial and temporal parameters based on the size and duration of clusters reported in descriptive studies (Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016; Sy et al., Reference Sy, Shaman, Kandula, Pei, Gould and Keyes2019). However, since descriptive studies have typically involved small community settings, it is unclear whether the spatio-temporal parameters used to guide inferential studies of suicide clusters are generalisable to broader populations and settings. Furthermore, distinguishing between suicide clusters which involve exposure to suicide, or links between cluster members, has the potential to improve a community's response to suicide clusters. Arguably, suicide clusters which involve social links between suicide descendants, may warrant different interventions and preventative approaches compared to those which involve no apparent links between cluster members.

This study sought to identify the characteristics and mechanisms of suicide clusters using three different, but common cluster detection methods, namely, the scan statistic, coronial inquests into suicide clusters and descriptive network analysis. These methods were chosen in order to operationalise the different definitions of suicide clusters, including those which are statistically significant, those which are significant to communities and those which involve social links among cluster members. Using the same nationwide data we aimed to: (1) determine the overlap between suicide clusters for each cluster detection method, (2) compare the spatial and temporal parameters associated with suicide clusters and (3) identify the demographic characteristics and prior of exposure to suicide among cluster and non-cluster members.

Methods

Data source and case ascertainment

The study was approved by the Justice Department Human Research Ethics Committee (CF/15/13188). Youth suicides were identified in the National Coronial Information System (NCIS). The NCIS is an online database that records external causes of death in Australia. Each death in the NCIS is accompanied by an individual case number which is linked to a set of core demographic variables including age, sex, indigenous status, employment status, marital status, date of death and location of death (Saar et al., Reference Saar, Bugeja and Ranson2017). Additional data sources in the NCIS include narrative text from police reports, coroner's reports (including inquest findings), autopsy reports and toxicology findings.

The following cases were included in the analysis: (1) the case was closed and the cause of death was determined as suicide (recorded as intentional self-harm by the coroner), (2) the person was aged between 10 and 24 years at the time of their death and (3) the death occurred in Australia between 1 January 2006 and 31 December 2015 (N = 3365). The year 2015 was selected as the study endpoint because over 90% of suicides from this time period were closed and therefore fully investigated by the coroner (National Coroners Information System, 2019). A case was excluded if: (1) the death occurred outside Australia, (2) the cause of death was undetermined or (3) the case did not include at least one coroner or police report (N = 93). Finally, suicides which occurred in the state of South Australia (N = 245) were excluded due to the limited availability of narrative texts in police and coroner reports. As a result, a total of 3027 cases were included in the analyses.

Methods for detecting suicide clusters

The scan statistic

The scan statistic was used to identify suicide clusters that were statistically significant using inferential methods. The space-time scan statistic was used to identify the presence of spatio-temporal suicide clusters in Australia using the software SaTScan version 9.4.4. (Kulldorff, Reference Kulldorff2006). Each suicide was acquired at the individual level and aggregated by the month of death and SA2. The 2011 Australian standard population (the midpoint of our study period) was used to estimate age-adjusted rates of suicide in young people aged 10–24 years. Population estimates were recorded for each SA2 from the 2011 Estimated Resident Population recorded by the Australian Bureau of Statistics.

The Poisson discrete scan statistic was used to detect high relative rates of suicide incidence and was conducted separately in each state and territory. The likelihood of each possible cluster was assessed using Monte Carlo simulations. Consistent with previous studies (Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016), the minimum time window was set at seven days to a maximum of 730 days (two years) and the maximum spatial window at 10% of the population at risk with a maximum radius of 100 km (Jones et al., Reference Jones, Gunnell, Platt, Scourfield, Lloyd, Huxley, John, Kamran, Wells and Dennis2013). This resulted in a set of cylindrical scanning windows where the base represents the area of the suicide cluster and the height represents the duration of the cluster. Clusters were included if their p value was <0.10 to account for the statistically rare incidence of suicides (Too et al., Reference Too, Pirkis, Milner and Spittal2017). Clusters are referred to as ‘possible clusters’ in the 0.10 < p < 0.05 range and ‘clusters’ if p < 0.05.

Coroner inquests into suicide clusters

Coroner inquests into suicide clusters were used to identify suicide clusters that were perceived to be significant by communities that had experienced multiple suicides. A custom query was built to identify all cases where the cause of death was intentional self-harm and the term ‘cluster’ or ‘contagion’ was included in the coroner reports. No limitations were placed on the distance or duration between suicides included in a coronial inquest. A case was identified as a cluster member if it was included in a coronial inquest into a suicide cluster in Australia during the study period. A case was excluded if: (1) the term ‘cluster’ did not refer to a cluster of suicidal behaviour (e.g. ‘Cluster B personality traits’) or (2) the suicide did not meet the eligibility criteria for included cases in the NCIS outlined above. A total of 55 case records were identified in the initial search and were manually inspected by the first author (N.T.M.H) for cluster membership. A total of seven independent coronial inquests into suicide clusters were identified, encompassing 40 individuals who died by suicide. Further details on the case ascertainment of cluster members is provided in S1 Methods in the supplement.

Descriptive network analysis

Descriptive network analysis was used to identify suicide clusters that comprised social links between three or more suicide descendants (Larkin and Beautrais, Reference Larkin and Beautrais2012; Public Health England, Reference Public Health England2015). The narrative text from police and coroner's reports was examined for evidence of social links between suicide descendants. Cases were linked if: (1) the young person was known to have been exposed to the suicide of another young person (e.g. they knew a friend who died by suicide; or they knew a peer from their school who died by suicide), (2) the police or coroner report referred to the first and last name of the person who previously died by suicide (the index case), (3) the police or coroner report referred to the date of death (e.g. month and year) and described at least one other characteristic that could identify the index case based on information included in the case records (e.g. the name of the school) and (4) There was evidence of social links between at least three or more individuals who died by suicide. No limitations were placed on the distance or duration between linked cases. A case was not linked and included in the network analysis if: (1) the police or coroner's report did not provide sufficient information on the index suicide death(s) or did not sufficiently match an eligible case in the NCIS, or (2) the suicide did not meet the eligibility criteria for included cases in the NCIS outlined above. Further details on the case ascertainment of linked cases are provided in S1 Methods in the supplement.

Comparison of cluster detection methods

Cluster membership was analysed as a binary variable representing the presence or absence of the outcome (cluster and non-cluster membership). Location of death was based on Statistical Area Level 2 (SA2) recorded in the NCIS. SA2s are general-purpose areas which represent communities that interact socially and economically. They represent a population range of 3000–25 000 persons and have an average population of ~10 000 persons (Australian Bureau of Statistics, 2018).

The comparison of cluster detection methods included: (1) the duration of the suicide cluster (the period of time between the index and last suicide death in the cluster), (2) the proportion of geospatial overlap between cluster detection methods (e.g. the number of suicide clusters that were detected by more than one method), (3) the spatial congruence between overlapping clusters (e.g. the proportion of SA2 spatial units that were identified within overlapping clusters) and (4) the number of overlapping cluster members (e.g. the number of individual cluster members that were identified by more than one cluster method).

The characteristics of cluster members were identified using the core demographic variables in the NCIS (age, sex, indigenous status, employment status, date of death, location and manner of death). Evidence of prior exposure to suicide was recorded for each case based on information recorded in the police and/or coroner's reports (e.g. family history of suicide or the death of a friend, peer or colleague). Differences between cluster and non-cluster members were compared using Pearson's chi-square test of independence. Fisher's exact probability test was employed when over 20% of cells had expected counts <5. All analyses, with the exception of the scan statistic, were conducted using R version 3.6.2. The geocode function of the ggmap package was used to identify the geocoordinates of each SA2 and to visualise the location of suicide clusters for each cluster detection method.

Results

Characteristics of suicide clusters

The scan statistic

Eight spatial–temporal clusters of high relative risk were detected in six Australian states and territories (Fig. 1). A total of 69 young people were identified as cluster members, accounting for 2.3% of suicides. The number of individual's involved in a single suicide cluster ranged between 3 and 20 and the duration of each suicide cluster ranged between 6 and 725 days (Table 1). Cluster members were more likely than non-cluster members to be from Aboriginal and Torres Strait Islander origins and to reside in a remote or very remote location. Evidence of prior exposure to suicide was reported among 10.1% (7/69) of cluster members compared to 6.6% (195/2958) of non-cluster members but this difference was non-significant (p = 0.38). The remaining demographic characteristics were comparable among cluster and non-cluster members (online Supplementary Table S4).

Fig. 1. Geographical locations of suicide clusters identified in Australian youth in 2006-2015 using the scan statistic. Geographic regions are based on Statistical Area 2 geographies (red). The circles represent a group of SA2s involved in a single cluster. All suicide clusters that were detected using the scan statistic were limited to a 100 km radius but are visualised using entire SA2 boundaries.

Table 1. Characteristics of suicide clusters detected using the scan statistic, coronial inquests and descriptive network analysis of social linked clusters

SA2, Australian Statistical Geography Standard – Statistical Area Level 2; S.D., standard deviation.

a The number of SA2s (per cluster).

b Complete (100% overlap among individuals and SA2s), partial (<100 overlap between cluster members and SA2s), none (no overlap between cluster members and SA2s).

c The number and proportion of SA2s that overlap (per cluster).

d The number and proportion of individuals that overlap (per cluster).

e Mean (S.D.).

Coronial inquests

Seven coronial inquests into suicide clusters were detected in three Australian states and territories (Fig. 2). A total of 40 young people were identified as cluster members, accounting for 1.3% of suicides. Six out of seven suicide clusters were described as being greater than expected on the basis of statistical inference and one suicide cluster involved social links between cluster members in an inpatient unit. The number of individuals involved in a single suicide cluster ranged between 3 and 8 and the duration of each suicide cluster ranged between 121 and 914 days (Table 1). Cluster members were more likely to be aged 18 years or less (70%, 28/40) compared to non-cluster members (29.7%, 887/2987, p < 0.01) and were more likely to be from Aboriginal and Torres Strait Island origins (55 v. 14.6% p < 0.01), and residing in a remote location at the time of death (50 v. 7.20%, p < 0.01). Evidence of prior exposure to suicide was reported among 32.5% (14/40) of cluster members compared to 6.43% (192/2990) of non-cluster members and this difference was statistically significant (p < 0.001). The remaining demographic characteristics were comparable among cluster and non-cluster members (online Supplementary Table S5).

Fig. 2. Geographical locations of suicide clusters identified in Australian youth in 2006–2015 using information from coronial inquests into suicide clusters. Geographic regions are based on Statistical Area 2 geographies (orange). The circles represent a group of SA2s involved in a single coronial inquest.

Descriptive network analysis

Eleven suicide clusters comprising seven triads and four tetrads were detected in four Australian states and territories (Fig. 3). A total of 37 young people were identified as cluster members, accounting for 1.3% of suicides. The number of individuals involved in a single suicide cluster ranged between 3 and 4 and the duration of each suicide cluster ranged between 31 and 1626 days (Table 1). Cluster members were more likely to be aged 18 years or less (62.2%, 23/37) compared to non-cluster members (29.83%, 892/2990, p < 0.01). Evidence of prior exposure to suicide was reported among 64.9% (24/37) of cluster members compared to 6.09% (182/2990) of non-cluster members and this difference was statistically significant (p < 0.001). The remaining demographic characteristics were comparable among cluster and non-cluster members (online Supplementary Table S6).

Fig. 3. Geographical locations of suicide clusters identified in Australian youth in 2006–2015 using descriptive network analysis of social linked suicide clusters. Geographic regions are based on Statistical Area 2 geographies (purple). The numbered nodes represent the sequence in which the suicide occurred beginning from the index death (1). The arrows represent the direction of exposure to suicide, and the connection between individuals in a suicide cluster.

Comparison of cluster detection methods

Overall, <50% of suicide clusters were identified by more than one cluster method (Table 1, Fig. 4a–c). The majority of overlapping clusters involved those that were detected from coronial inquests and the descriptive network analysis (four out of eight clusters). However, spatial congruence between overlapping clusters was low (<100% of SA2s). The duration of suicide clusters ranged between 6 and 725 days using the scan statistic, 121–914 days using information from coronial inquests, and 31–1626 days using descriptive network analysis. Overall, one-third (12/133) of cluster members were identified by more than one cluster detection method, however the number of young people that were identified by more than one method ranged from 1 to 3 individuals.

Fig. 4. (a) Overlapping clusters that occurred across Australia; (b) overlapping clusters that occurred in smaller regions across Queensland (QLD); (c) overlapping clusters that occurred smaller regions across Victoria (VIC). All suicide clusters that were detected using the scan statistic occurred within a 100km radius but are visualised using SA2 boundaries

Discussion

This study identified suicide clusters in Australian youth aged 10–24 who died by suicide in 2006–2015. Comparative analysis of the three cluster detection methods showed considerable heterogeneity in the location, duration and size of suicide clusters, as well as both the number and demographic characteristics of cluster members. Together these findings demonstrate that the ways in which suicide clusters are defined and operationalised can result in markedly different suicide clusters. This was corroborated by our comparison of cluster members which showed that 90.8% of individuals who were involved in a single suicide cluster were not identified in overlapping clusters, using alternate cluster detection methods.

Although inferential approaches for the detection of suicide clusters have many empirical advantages (e.g. the facilitation of evidence based-decision making), the accuracy and precision of suicide clusters detected using the scan statistic depends on pre-determined spatial and temporal parameters, for which there is currently no-gold standard. Results from both the coronial inquest into suicide clusters and the descriptive network analysis revealed that the maximum duration of suicide clusters ranged between 2.5 and 4.5 years, respectively. Together, these findings suggest that the predetermined 2-year parameter used in the scan statistic analyses did not correspond with the maximum duration of suicide clusters that were perceived as significant by local communities, nor suicide clusters which involved social links between cluster members. Since both the coronial inquest into suicide clusters and the descriptive network analyses were conducted on a nationwide scale and therefore overcome some of the limitations of earlier descriptive studies (e.g. small study effects and limited generalisability) results from the present study have the potential to better inform the parameters used in inferential studies of suicide clusters.

In the present study, six out of seven suicide clusters involved in coronial inquests verified by the coroner as being greater than statistically expected. Despite this, there was limited geospatial overlap between coronial inquests and suicide clusters that were identified using the scan statistic. It is likely, that these differences were influenced by broad area aggregation and subsequent differences in the geospatial boundaries used for each method (Nelson and Brewer, Reference Nelson and Brewer2017). Suicide clusters that were detected using the scan statistic, for example, were aggregated by SA2s, however suicide clusters that were detected using coronial inquests did not correspond with precise SA2 boundaries. These findings suggest that spatial parameters comprising broad area aggregations may have limited translational benefit for local communities and stakeholders that are tasked with responding to suicide clusters. The use of point-level geocoded data has been shown to effectively distinguish between the suicide profiles of small community settings and has the potential to overcome spatial biases that arise due to broad area aggregation in studies of suicide clusters (Torok et al., Reference Torok, Shand, Phillips, Meteoro, Martin and Larsen2019). Future research which investigates relevance of different spatial parameters among stakeholders who are responsible for the coordination of targeted suicide prevention and cluster response activities at the local and regional level is therefore warrented.

Exisiting guidelines for the detection response and prevention of suicide clusters recommend that community stakeholders use both inferential methods and descriptive methods to identify the links between cluster members. Although the descriptive network analyses were conducted using nation-wide data, information on prior exposure to suicide was not consistently recorded in the NCIS, resulting in linkage among only 60% of exposed cases. Routine collection of exposure to suicide as part of the police and coroner investigation as well as the inclusion of exposure to suicide in the core-data set of suicide cluster surveillence systems has the potential to improve the sensitivity and specificity of suicide clusters that are characterised by social links between cluster members.

Lastly, results from the present study provide some insight into the mechanisms underlying suicide clusters. Lifetime prevalence of exposure to suicide among cluster members that were detected using the scan statistic was 10.1% compared to over one third of cluster members involved in 32.5% of clusters identified through coronial inquests, and 52.6% of cluster members identified using descriptive network analyses. These findings are contrary to previous studies that suggest that the close spatial and temporal parameters used in the scan statistic is consistent with a social transmission hypothesis (Gould et al., Reference Gould, Wallenstein and Kleinman1990a; Robinson et al., Reference Robinson, Too, Pirkis and Spittal2016; Helbich et al., Reference Helbich, Plener, Hartung and Blüml2017; Sy et al., Reference Sy, Shaman, Kandula, Pei, Gould and Keyes2019). Whilst the prevention of the social transmission of suicidal behaviour is a core component of existing cluster response guidelines in Australia (Centre for Health Policy Programs and Economics, 2012) and internationally (Centers for Disease Control, 1988; Health Service Executive (HSE), 2011; Public Health England, Reference Public Health England2015; Palmer et al., Reference Palmer, Inder, Shave and Bushnell2018; Public Health England, 2019), results from the present study suggest that not all clusters involve social links between cluster members.

Limitations

Although the NCIS is a comprehensive database on coroner determined causes of death, information included in the NCIS is retrospective and may not be immediately available to communities that are experiencing a suicidal crisis. Furthermore, data included in our analysis was dependent on the narrative text of police and coroner's reports. Consequently, information on prior exposure to suicide may not have been recorded if it was not included in the police or coroner investigation, or if informants involved in the investigation did not have knowledge of the young person's exposure to suicide. Lastly, the limited geospatial congruence between suicide clusters that were detected using the scan statistic are an artefact of the spatial and temporal parameters used in the analysis (Cheung et al., Reference Cheung, Spittal, Williamson, Tung and Pirkis2013; Nelson and Brewer, Reference Nelson and Brewer2017), specifically aggregation of suicides using SA2s. Future research should therefore seek to acquire point data in order to detect suicide clusters using inferential methods with greater accuracy and precision.

Conclusions

The present study showed that the same data – when analysed in different ways – may lead to the detection of different suicide clusters. Whilst some suicide clusters involved social links between suicide descendants, not all cluster members had a history of prior exposure to suicide. The collection of exposure to suicide as part of routine police and coroner investigations has the potential to improve the ways communities respond to suicide clusters by distinguishing between suicide clusters that may or may not be driven by social links between cluster members. Finally, results from the present study have the potential to inform the parameters used in inferential studies of suicide clusters, particularly among local suicide prevention teams that are tasked with responding to suicide clusters.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S2045796020000645.

Data

The NCIS, the database used in the current study, has no public access due to the sensitivity of the suicide data.

Acknowledgements

None.

Financial support

NH is a PhD student and was supported by the Australian Rotary Health PhD Partnership Scholarship. JR is supported by a National Health and Medical Research Council (NHMRC) Career Development Fellowship. LST was supported by a NHMRC Early Career Fellowship. MS is a recipient of an Australian Research Council Future Fellowship (project number FT180100075) funded by the Australian Government.

Conflict of interest

None.

Ethical standards

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. This study received full ethical approvals from the Justice Department Human Research Ethics Committee (CF/15/13188).

References

Australian Bureau of Statistics (2018) Australian Statistical Geography Standard (ASGS). Canberra: Australian Bureau of Statistics.[Online]. Available at https://www.abs.gov.au/websitedbs/D3310114.nsf/home/Australian+Statistical+Geography+Standard+(ASGS).Google Scholar
Bechtold, DW (1988) Cluster suicide in American Indian adolescents. American Indian and Alaska Native Mental Health Research 1, 2635.CrossRefGoogle ScholarPubMed
Centers for Disease Control (1988) CDC Recommendations for a community plan for the prevention and containment of suicide clusters. Morbidity and Mortality Weekly Report 37, 12.Google Scholar
Centre for Health Policy Programs and Economics (2012) Developing A Community Plan for Preventing and Responding to Suicide Clusters. Melbourne: Melbourne School of Population Health, The University of Melbourne.Google Scholar
Cheung, YTD, Spittal, MJ, Pirkis, J and Yip, PSF (2012) Spatial analysis of suicide mortality in Australia: investigation of metropolitan-rural-remote differentials of suicide risk across states/territories. Social Science & Medicine 75, 14601468.CrossRefGoogle ScholarPubMed
Cheung, YT, Spittal, MJ, Williamson, MK, Tung, SJ and Pirkis, J (2013) Application of scan statistics to detect suicide clusters in Australia. PLoS One 8, e54168.CrossRefGoogle ScholarPubMed
Clinical Advisory Services Aotearoa (2016) Coronial Suspected Suicide Data Sharing Service (CDS). Aotearoa: Clinical Advisory Services, Aotearoa.Google Scholar
Davidson, LE, Rosenberg, ML, Mercy, JA, Franklin, J and Simmons, JT (1989) An epidemiologic study of risk factors in two teenage suicide clusters. The Journal of the American Medical Association 262, 26872692.CrossRefGoogle ScholarPubMed
Gould, MS, Wallenstein, S and Kleinman, M (1990 a) Time-space clustering of teenage suicide. American Journal of Epidemiology 131, 7178. Gould MS, Wallenstein S and Kleinman M (1990b) Time-space clustering of teenage suicide. American Journal of Epidemiology 131, 71–78.CrossRefGoogle ScholarPubMed
Gould, MS, Wallenstein, S, Kleinman, MH, O'Carroll, P and Mercy, J (1990 b) Suicide clusters: an examination of age-specific effects. American Journal of Public Health 80, 211212.CrossRefGoogle ScholarPubMed
Griffin, E, Dillion, CB, Arensman, E, Corcoran, P, Williamson, E and Perry, IJ (2017) National self-harm registry Ireland: annual report 2016.Google Scholar
Hawton, K, Hill, NTM, Gould, M, John, A, Lascelles, K and Robinson, J (2019) Clustering of suicides in children and adolescents. The Lancet Child & Adolescent Health 4, 5867.CrossRefGoogle ScholarPubMed
Health Service Executive (HSE) (2011) Responding to Murder Suicide and Suicide Clusters Guidance Document. Ireland: Health Service Executive (HSE).Google Scholar
Helbich, M, Plener, PL, Hartung, S and Blüml, V (2017) Spatiotemporal suicide risk in Germany: a longitudinal study 2007–11. Scientific Reports 7, 18.CrossRefGoogle ScholarPubMed
Jones, P, Gunnell, D, Platt, S, Scourfield, J, Lloyd, K, Huxley, P, John, A, Kamran, B, Wells, C and Dennis, M (2013) Identifying probable suicide clusters in wales using national mortality data. PLoS One 8, e71713.CrossRefGoogle ScholarPubMed
Kulldorff, M (2006) Satscan: Software for the spatial and space-time scan statistics.Google Scholar
Larkin, GL and Beautrais, AL (2012) Geospatial Mapping of Suicide Clusters. Te Pou o Te Whakaaro Nui.Google Scholar
National Coroners Information System (2019) Operational Statistics: Case Closure [Online]. VIC, Australia: National Coronial Information System. Available at https://www.ncis.org.au/about-the-data/operational-statistics/ [Accessed 28 October 2019].Google Scholar
Nelson, JK and Brewer, CA (2017) Evaluating data stability in aggregation structures across spatial scales: revisiting the modifiable areal unit problem. Cartography and Geographic Information Science 44, 3550.CrossRefGoogle Scholar
Niedzwiedz, C, Haw, C, Hawton, K and Platt, S (2014) The definition and epidemiology of clusters of suicidal behavior: a systematic review. Suicide and Life Threatening Behavior 44, 569581.CrossRefGoogle ScholarPubMed
Palmer, S, Inder, M, Shave, R and Bushnell, J (2018) Postvention Guidelines for the Management of Suicide Clusters. Aotearoa: Clinical Advisory Services.Google Scholar
Poijula, S, Wahlberg, KE and Dyregrov, A (2001) Adolescent suicide and suicide contagion in three secondary schools. International Journal Emergency Mental Health 3, 163168.Google ScholarPubMed
Public Health England, (2015) Identifying and responding to suicide clusters and contagion: A practice resource. London, UK: Department of Health.Google Scholar
Public Health England (2019) Identifying and Responding to Suicide Clusters: A Practice Resource. London, UK: Public Health England.Google Scholar
Public Health England (2020) Thames Valley Real Time Suicide Surveillance System [Online]. UK: Public Health England. Available at https://www.nspa.org.uk/wp-content/uploads/2018/07/Thames-Valley-Real-time-Suicide-Surveillance-pack.pdf [Accessed 26 June 2020].Google Scholar
Robertson, L, Skegg, K, Poore, M, Williams, S and Taylor, B (2012) An adolescent suicide cluster and the possible role of electronic communication technology. Crisis 33, 239245.CrossRefGoogle ScholarPubMed
Robinson, J, Too, LS, Pirkis, J and Spittal, MJ (2016) Spatial suicide clusters in Australia between 2010 and 2012: a comparison of cluster and non-cluster among young people and adults. BMC Psychiatry 16, 417426.CrossRefGoogle ScholarPubMed
Saar, E, Bugeja, L and Ranson, DL (2017) National coronial information system: epidemiology and the coroner in Australia. Academic Forensic Pathology 7, 582590.CrossRefGoogle ScholarPubMed
Sy, KTL, Shaman, J, Kandula, S, Pei, S, Gould, M and Keyes, KM (2019) Spatiotemporal clustering of suicides in the US from 1999 to 2016: a spatial epidemiological approach. Social Psychiatry and Psychiatric Epidemiology 12, 14711482.CrossRefGoogle Scholar
Taiminen, TJ, Kallio-Soukainen, K, Nokso-Koivisto, H, Kaljonen, A and Helenius, H (1998) Contagion of deliberate self-harm among adolescent inpatients. Journal of the American Academy of Child & Adolescent Psychiatry 37, 211217.CrossRefGoogle ScholarPubMed
Too, LS, Pirkis, J, Milner, A and Spittal, MJ (2017) Clusters of suicides and suicide attempts: detection, proximity and correlates. Epidemiology & Psychiatric Science 26, 491500.CrossRefGoogle ScholarPubMed
Torok, M, Shand, F, Phillips, M, Meteoro, N, Martin, D and Larsen, M (2019) Data-informed targets for suicide prevention: a small-area analysis of high-risk suicide regions in Australia. Social Psychiatry and Psychiatric Epidemiology 54, 12091218.CrossRefGoogle ScholarPubMed
Wilkie, C, Macdonald, S and Hildahl, K (1998) Community case study: suicide cluster in a small Manitoba community. The Canadian Journal of Psychiatry 43, 823828.CrossRefGoogle Scholar
Williamson, MK, Tung, SJ, Pirkis, J, Derek Cheung, YT and Spittal, MJ (2014) Predictors of suicides occurring within suicide clusters in Australia, 2004–2008. Social Science and Medicine 118, 135142.Google Scholar
Wissow, LS, Walkup, J, Barlow, A, Reid, R and Kane, S (2001) Cluster and regional influences on suicide in a Southwestern American Indian tribe. Social Science and Medicine 53, 11151124.CrossRefGoogle Scholar
Figure 0

Fig. 1. Geographical locations of suicide clusters identified in Australian youth in 2006-2015 using the scan statistic. Geographic regions are based on Statistical Area 2 geographies (red). The circles represent a group of SA2s involved in a single cluster. All suicide clusters that were detected using the scan statistic were limited to a 100 km radius but are visualised using entire SA2 boundaries.

Figure 1

Table 1. Characteristics of suicide clusters detected using the scan statistic, coronial inquests and descriptive network analysis of social linked clusters

Figure 2

Fig. 2. Geographical locations of suicide clusters identified in Australian youth in 2006–2015 using information from coronial inquests into suicide clusters. Geographic regions are based on Statistical Area 2 geographies (orange). The circles represent a group of SA2s involved in a single coronial inquest.

Figure 3

Fig. 3. Geographical locations of suicide clusters identified in Australian youth in 2006–2015 using descriptive network analysis of social linked suicide clusters. Geographic regions are based on Statistical Area 2 geographies (purple). The numbered nodes represent the sequence in which the suicide occurred beginning from the index death (1). The arrows represent the direction of exposure to suicide, and the connection between individuals in a suicide cluster.

Figure 4

Fig. 4. (a) Overlapping clusters that occurred across Australia; (b) overlapping clusters that occurred in smaller regions across Queensland (QLD); (c) overlapping clusters that occurred smaller regions across Victoria (VIC). All suicide clusters that were detected using the scan statistic occurred within a 100km radius but are visualised using SA2 boundaries

Supplementary material: File

Hill et al. supplementary material

Hill et al. supplementary material

Download Hill et al. supplementary material(File)
File 2.3 MB