- Research article
- Open Open Peer Review
Spatial clusters of suicide in Australia
BMC Psychiatryvolume 12, Article number: 86 (2012)
Understanding the spatial distribution of suicide can inform the planning, implementation and evaluation of suicide prevention activity. This study explored spatial clusters of suicide in Australia, and investigated likely socio-demographic determinants of these clusters.
National suicide and population data at a statistical local area (SLA) level were obtained from the Australian Bureau of Statistics for the period of 1999 to 2003. Standardised mortality ratios (SMR) were calculated at the SLA level, and Geographic Information System (GIS) techniques were applied to investigate the geographical distribution of suicides and detect clusters of high risk in Australia.
Male suicide incidence was relatively high in the northeast of Australia, and parts of the east coast, central and southeast inland, compared with the national average. Among the total male population and males aged 15 to 34, Mornington Shire had the whole or a part of primary high risk cluster for suicide, followed by the Bathurst-Melville area, one of the secondary clusters in the north coastal area of the Northern Territory. Other secondary clusters changed with the selection of cluster radius and age group. For males aged 35 to 54 years, only one cluster in the east of the country was identified. There was only one significant female suicide cluster near Melbourne while other SLAs had very few female suicide cases and were not identified as clusters. Male suicide clusters had a higher proportion of Indigenous population and lower median socio-economic index for area (SEIFA) than the national average, but their shapes changed with selection of maximum cluster radii setting.
This study found high suicide risk clusters at the SLA level in Australia, which appeared to be associated with lower median socio-economic status and higher proportion of Indigenous population. Future suicide prevention programs should focus on these high risk areas.
Suicide remains an important public health problem in Australia, with approximately 2,000 suicide cases occurring nationally since the mid-1990s (an average annual incidence of approximately 9.7 per 100,000) [1, 2]. Suicide in Australia has also been shown to vary by area-based socio-economic strata , by measures of urban–rural residence , and by small-area geographic units [5, 6]. Geographic analyses to identify areas of high suicide risk have also been conducted in other countries, such as Belgium , Taiwan [8, 9], United Kingdom [10–15], United States [16, 17], and Brazil .
Spatial analysis methods have been used to identify high risk areas of suicide in previous studies [6, 8–11, 13, 17, 18]. However, some studies focused only on urban areas or discrete geographic areas (e.g., state level) [17, 18], and others used smoothed standardized mortality ratios (SMRs) to identify areas with high or low suicide incidences [9, 10]. Such smoothing methods may be inadequate in areas with sparse population and very few suicide cases [9, 10]. Studies in Brazil  and the United Kingdom  did not examine the pattern of suicide risk in rural and remote areas, which may substantially differ from urban areas. In Australia, a few reports and studies have examined the suicide pattern across the whole country [1, 2, 19], but the spatial pattern of suicide has not been investigated in the national wide to date. At the state level some studies on suicide patterns have also been implemented in New South Wales  the Northern Territory , and South Australia , but they did not use spatial analysis. A spatial analysis can identify areas with high suicide mortality, visualize the pattern of high risk areas and explore possible reasons for high risk, e.g. socioeconomic variety across different areas. Given the geographic and socio-demographic diversity across Australia, it is important to examine the suicide pattern in the whole of Australia. This study is based on our previous research which was conducted in Queensland .
Suicide data (1999–2003), including sex, age, country of birth, year and month of suicide and statistical local area (SLA) code were provided by Australian Bureau of Statistics (ABS). Access to recent suicide mortality data was unavailable as the related procedure is currently under review. The institutional ethics approval was granted by the Human Research Ethics Committee, Queensland University of Technology.
In 2001, Australia was divided into 1,346 SLAs in eight states and territories: New South Wales (NSW), Victoria (VIC), Queensland (QLD), South Australia (SA), Western Australia (WA), Tasmania (TAS), Northern Territory (NT) and Australian Capital Territory (ACT). 2001 Census data (CDATA 2001) included information on digital statistical boundaries, base maps, and SLA information (name, code, area (km²), longitude and latitude of the centroids, and population by age and gender). Some suicide cases (less than 1% of total suicides) were excluded from this study due to a lack of documentation on the SLA code. Socio-demographic data, such as Socio-Economic Indexes for Area (SEIFA), Indigenous status and unemployment rate at the SLA level were obtained from CDATA 2001. These data were used as potential explanatory factors for clusters. A higher SEIFA score indicates higher socio-economic status at the SLA level.
A series of statistical and Geographical Information System (GIS) methods were applied to analyse the data. GIS is a powerful tool to store, retrieve and display spatial data. The spatial locations and trends of disease outbreaks can be identified by GIS . GIS can also spatially explore disease aetiologies, such as socio-environmental factors associated with disease . Then GIS can help public health workers to design and implement effective disease control and prevention programs at a local level. Descriptive analysis was conducted to explore the characteristics of each variable, such as suicide cases, population and suicide age-adjusted standardised mortality (ASM) by gender . GIS and mapping approaches were applied to investigate the spatial distribution of suicide SMR by gender at the SLA level . A direct method, which was used in our previous Queensland study , was applied to calculate the ASM by gender for each SLA. Then the SMR was also calculated by using the mean of annual mortality of suicide in the whole Australia (1999–2003) as a reference, stratified by gender.
To identify high risk areas, we applied spatial cluster analysis to identify the randomly-distributed suicide cases and to explore primary clusters (i.e., those with the highest risk among all clusters) and secondary clusters (other high risk clusters with significance), using SaTScan [26, 27]. A Poisson regression model was performed to compute the mean relative risk (RR) of each cluster and likelihood ratio to identify the two types of clusters as circular windows. We also used maximum population size of cluster (covering less than 50%, 25% and 10% of total population) and maximum length of cluster radii (100 km and 400 km) to examine whether these factors could determine the pattern of clusters, as 100 km radii may cover most of urban areas and 400 km radii may cover most of rural or remote areas. Different population size and radii of cluster can be selected in SaTScan. As some documents indicated that most of suicide cases were aged between 15 and 54 [1, 2], the spatial clusters of suicide for this age group were specifically examined as two groups (15 to 34-year and 35 to 54-year). The median SEIFA and mean values of other socio-demographic factors of clusters were calculated and compared with the national average, in order to identify whether there were any differences in these factors between clusters and other areas.
11,586 suicide deaths were included in the analysis covering the whole study period, with 9,142 males and 2,444 females (male/female ratio: 3.74). 4,472 (38.6%) of total suicides were aged between 15 and 34 (3,577 males and 895 females) and 4,540 (39.2% of total) between 35 and 54 (3,579 males and 961 females). The metropolitan areas of capital cities covered only 0.5% of the total area examined but contained 63.4% of total population  and 59.4% of total suicide deaths. Rural and remote SLAs had sparse population density and many had no suicide cases during the study period. Table 1 indicates the distribution of area, population and suicide by SLA. All the variables differed substantially across SLAs.
Figure 1A indicates that the central south of Northern Territory (NT), north, inland and the eastern coast of Queensland (QLD), inland areas of New South Wales (NSW) and Western Australia (WA), and the eastern coast of Tasmania (TAS) had higher male SMRs than other areas. There were no male suicide cases during the study period in central and some southern areas of WA; central, north and southwest QLD; and northwest NSW. North and inland areas QLD, some areas in NSW and VIC, north and some western areas in WA, and eastern TAS had higher female suicide SMRs than other areas (Figure 1B). However, 47.6% of SLAs had no female suicide cases during the study period.
Figure 2 shows the spatial clusters of male suicide. Mornington Shine (A in Figure 2A and 2B) in the northwest of QLD was identified as the primary cluster. Bathurst-Melville area (B1, north of NT) and some suburbs of Adelaide in SA (E) were discovered as secondary clusters at both of the maximum radii of 100 km (Figure 2A) and 400 km (Figure 2B). Other secondary clusters changed with selection of different maximum lengths of cluster radii. All the detailed information of male suicide clusters at different cluster radii is provided in Table 2. For female suicide clusters, only one primary cluster was identified near Melbourne after adjustment for radius and age group (Table 2). We also tested the impact of different population size (i.e., clusters of less than 50%, 25%, 10% and 5% of total population) on the spatial pattern of suicide but did not find any significant difference.
The clusters of male suicide in different age groups were also presented. Mornington Shire was the primary cluster of male suicide for youths (15 to 34-year) at the maximum radii of 100 km (Figure 3A) and this cluster expanded to cover much larger areas of Northwest QLD (including Mornington Shire) at the maximum radii of 400 km (Figure 3B). Bathurst-Melville area was identified as a secondary cluster with different cluster radii selection (Figure 3A and 3B). Other secondary clusters changed with the selection of different maximum lengths of cluster radii (Figure 3A and 3B). For suicide aged between 35 and 54, only one cluster was found in the Central Coast of NSW at the maximum radii of 100 km (Figure 4A) and this cluster expanded to cover the whole Northeast NSW and parts of South QLD at the maximum radii of 400 km (Figure 4B). Table 3 provides detailed information of clusters at different age groups and cluster radii.
Table 4 shows the socio-demographic characteristics of high risk clusters for male suicide at all ages. Primary and secondary clusters had a higher proportion of Indigenous population than the national average level (about 39 times higher in the primary cluster; secondary clusters: 41% higher in the 100 km radii and 73% higher in the 400 km radii settings). The median SEIFA scores in the clusters were lower than the national median SEIFA score (Table 4), indicating that high risk clusters had lower socioeconomic status than the national average level. For female suicide, we did not find any significant difference in socio-demographic factors between cluster areas and other areas.
This study explored the geographical distribution of suicide and clusters of high risk in Australia. There was a higher risk of male suicide in the north of QLD, some areas in the east coast of QLD and TAS, inland areas in QLD, NSW and WA, and central south areas in NT (SMR > 2), compared with the national average. Female suicide incidence was significantly lower than male suicide incidence overall, and over 40% of all SLAs had no female suicide deaths in the study period. Only one cluster of female suicide was identified in this study.
In capital cities, the numbers of suicide cases in each year were relatively steady. Only some SLAs near Adelaide and Darwin were identified as clusters of high risk of male suicide. The Darwin metropolitan area has a smaller population than other capital cities, as well as a smaller average population at the SLA level. Thus the suicide mortality was high in some SLAs in Darwin. A report in South Australia indicated that the west of Adelaide, where the cluster in this study lies, had higher incidence of mental and behavioural disorders than that of the whole of Adelaide and Australia, which may be associated with suicidal behaviours . For female suicide, many SLAs with high SMR had only a small number of suicide cases. Thus the clusters of female suicide were not as obvious as male suicide. For example, Moreton Island near Brisbane had a SMR of 22.2 compared with the national female suicide incidence, but had only 1 female suicide. Thus these SLAs were not identified as high risk clusters by SaTScan. This phenomenon can also be found in some SLAs having a high male suicide SMR but with a very small population size (e.g., Fyshwick in Canberra).
In the spatial cluster analysis of total and 15 to 34-year male suicide, the Mornington Shire was identified as the whole or part of a primary cluster of high risk, as well as the secondary cluster of Bathurst-Melville area in both when setting the maximum radii at 100 km and 400 km. Our previous studies indicated that suicide was higher in the areas with larger proportion of Indigenous population than other areas and these areas usually had lower socioeconomic status [5, 6]. The Bathurst-Melville area also had a low SEIFA score and over 80% of local population were Indigenous. The findings in previous and current studies are similar which suggests that the same set of determinants of suicide clusters exist at both national and state levels. Additionally, social disruption and alcohol abuse may also contribute to the high suicide incidence, according to some studies [30, 31] and media reports [32–34] relating to the Mornington Shire and Bathurst-Melville. Due to very small numbers of 35 to 54-year suicides in Mornington Shire (2 cases) and Bathurst-Melville (1 case), these two areas were not identified as high risk clusters in this age group.
The sizes and positions of other secondary clusters changed with the selection of different radius and age structure. A cluster with a large radius limit could hide the information of smaller areas within the cluster, while the selection of clusters with a small radius limit may miss some significant high risk areas compared with a larger radius. SaTScan applies a circling approach to select all the geographical units in one place (e.g., northeast of NSW) as a particular cluster, in which these areas may be heterogeneous, especially when the radius setting is large [35, 36]. Thus some clusters tend to include both high risk areas and adjacent areas with low risk or even no suicide cases. Compared with the study in QLD , this study covered a much larger area, thus SaTScan could select more SLAs within one cluster (larger radius). This can explain that the cluster (400 km of maximum radius setting) in the east of QLD in this national study had 13 times the size (km2), 9 times of population and 7 times the number of suicide cases of that in the QLD study . There are some variations of significant clusters (especially secondary) between males at different age groups, due partly to the different distribution of age groups of population and suicides across SLAs (Table 1).
There are several strengths in this study. This is the first study to examine the spatial clusters of suicide at a national level in Australia. Clustered areas of high risk need to be identified to facilitate the assessment of factors associated with high suicide risk and to design effective public health interventions. This study explored the variations of spatial clusters at different settings, including various cluster radii and age groups. Finally, the method developed in this study may contribute to identifying high risk areas of other mental health problems or diseases and improving mental health promotion.
The limitations of this study should also be acknowledged. Firstly, the data are not current and covered a period when Australian suicide was declining after a peak in 1997. Thus it may not represent the current patterns of suicide, potentially limiting its use in current suicide prevention strategies. Secondly, detailed personal information (e.g., health status before death and suicide methods) was not available in this study. Thus it is difficult to assess how suicidal behaviours may be modified by individual-level factors in these areas. Finally, some risk factors such as drug and alcohol use, and potential modifying factors such as suicide prevention activities and the provision of healthcare services were not taken into account. These data are not routinely collected at the SLA level, but it is likely that such factors would affect spatial patterns of suicide.
Based on the findings of this study, some recommendations can be proposed. Firstly, even though most suicides occur in capital cities due to a large urban population, some rural and remote areas had high suicide risk, which is consistent with previous studies of rural suicide in Australia [3–6]. It is necessary to collect more detailed information (e.g., suicide methods, mental health status of suicide cases and general population) in high risk areas, to discern the causes of suicide and to help design specific suicide prevention activities in these areas. Previous studies have found that a series of suicide control and prevention activities targeting the general population, such as antidepressant use [37, 38], firearm and pesticide restrictions [39, 40], domestic gas detoxification [41, 42], primary care physician education [43, 44], and public education campaigns [45, 46], may be associated with reductions in suicidal behaviour . Investigating the extent to which such activities in local populations have been, or can be, implemented in high risk areas, is important in increasing the effectiveness of suicide prevention programs. The impact of climate [48, 49], socioeconomic factors [50–52] and natural disasters [53, 54] on suicide should also be addressed in future research, especially in identifying the variation of these impacts across different areas.
This study has described the spatial variation of suicide and clusters of high risk in Australia. The spatial and cluster analysis methods may have significant applications in mental health research and the development of effective suicide control and prevention strategies, especially after studying suicide patterns in more recent years. Exploring the spatiotemporal association between socio-environmental variables and suicide may also provide a better understanding of suicide cases at a local level and offer area-specific information to plan and implement suicide prevention activities.
3309.0 Suicides, Australia, 1921 to 1998. 2000, Australian Bureau of Statistics, Canberra, Australia
3309.0.55.001 Suicide: recent trends, Australia, 1993–2003. 2004, Australian Bureau of Statistics, Canberra, Australia
Page A, Morrell S, Taylor R, Carter G, Dudley M: Divergent trends in suicide by socio-economic status in Australia. Social Psychiatry and Psychiatric Epidemiology. 2006, 41: 911-917. 10.1007/s00127-006-0112-9.
Page A, Morrell S, Taylor R, Dudley M, Carter G: Further increases in rural suicide in young Australian adults: secular trends, 1979–2003. Social Science and Medicine. 2007, 65: 442-453. 10.1016/j.socscimed.2007.03.029.
Qi X, Tong S, Hu W: Preliminary spatiotemporal analysis of the association between socio-environmental factors and suicide. Environ Heal. 2009, 8: 46-10.1186/1476-069X-8-46.
Qi X, Tong S, Hu W: Spatial distribution of suicide in Queensland. Australia. BMC Psychiatry. 2010, 7: 10-
Hooghe M, Vanhoutte B: An ecological study of community-level correlates of suicide mortality rates in the Flemish region of Belgium, 1996–2005. Suicide and Life-Threatening Behavior. 2011, 41: 453-464. 10.1111/j.1943-278X.2011.00047.x.
Chang SS, Sterne JA, Wheeler BW, Lu TH, Lin JJ, Gunnell D: Geography of suicide in Taiwan: spatial patterning and socioeconomic correlates. Health and Place. 2011, 17: 641-645. 10.1016/j.healthplace.2011.01.003.
Chang SS, Gunnell D, Wheeler BW, Yip P, Sterne JA: The evolution of the epidemic of charcoal-burning suicide in Taiwan: a spatial and temporal analysis. PLoS Medicine. 2010, 7: 1-
Gunnell D, Wheeler B, Chang SS, Thomas B, Sterne JA, Dorling D: Changes in the geography of suicide in young men: England and Wales 1981–2005. Journal of Epidemiology and Community Health. 2010, in press
Middleton N, Sterne J, Gunnell DJ: An atlas of suicide mortality: England and Wales, 1988–1994. Health and Place. 2008, 14: 492-506. 10.1016/j.healthplace.2007.09.007.
Middleton N, Sterne J, Gunnell DJ: The geography of despair among 15-44-year-old men in England and Wales: putting suicide on the map. Journal of Epidemiology and Community Health. 2006, 60: 1040-1047. 10.1136/jech.2005.045302.
Saunderson TR, Langford IH: A study of the geographical distribution of suicide rates in England and Wales 1989–92 using empirical Bayes estimates. Social Science and Medicine. 1996, 43: 489-502. 10.1016/0277-9536(95)00427-0.
Rezaeian M, Dunn G, St Leger S, Appleby L: The ecological association between suicide rates and indices of deprivation in English local authorities. Social Psychiatry and Psychiatric Epidemiology. 2005, 40: 785-791. 10.1007/s00127-005-0960-x.
Rezaeian M, Dunn G, St Leger S, Appleby L: Do hot spots of deprivation predict the rates of suicide within London boroughs?. Health and Place. 2007, 13: 886-893. 10.1016/j.healthplace.2007.02.004.
Congdon P: The spatial pattern of suicide in the US in relation to deprivation, fragmentation and rurality. Urban Studies. 2011, 48: 2101-2122. 10.1177/0042098010380961.
Saman DM, Walsh S, Borowko A, Odoi A: Does place of residence affect risk of suicide? A spatial epidemiologic investigation in Kentucky from 1999 to 2008. BMC Publ Health. 2012, 12: 1-10.1186/1471-2458-12-1.
Minamisava R, Nouer SS, Neto OL, Melo LK, Andrade AL: Spatial clusters of violent deaths in a newly urbanized region of Brazil: highlighting the social disparities. Int J Heal Geogr. 2009, 8: 66-10.1186/1476-072X-8-66.
Large MM, Nielssen OB: Suicide in Australia: meta-analysis of rates and methods of suicide between 1988 and 2007. Medical Journal of Australia. 2010, 192: 432-437.
Suicide in New South Wales. 2004, We need to know more – The NSW Suicide Data Report, Sydney
Measey MA, Li SQ, Parker R, Wang Z: Suicide in the Northern Territory, 1981–2002. Medical Journal of Australia. 2006, 185: 315-319.
Lawrence D, Almeida OP, Hulse GK, Jablensky AV, Holman CD: Suicide and attempted suicide among older adults in Western Australia. Psychol Med. 2000, 30: 813-821. 10.1017/S0033291799002391.
Croner CM, Sperling J, Broome FR: Geographic information systems (GIS): new perspectives in understanding human health and environmental relationships. Statistics in Medicine. 1996, 15: 1961-1977. 10.1002/(SICI)1097-0258(19960930)15:18<1961::AID-SIM408>3.0.CO;2-L.
Norusis M: SPSS 16.0 guide to data analysis. 2008, Prentice Hall, Upper Saddle River, New Jersey, United States of America
MapInfo Professional Version 8.5 User Guide (Unabridged). 2006, MapInfo Corporation, New York, United States of America
Kulldorff M: SaTScan User Guide for version 9.0. 2010 http://www.satscan.org[31October2011][31October2011],
Kulldorff M: A spatial scan statistic. Communications in Statistics: Theory and Methods. 1997, 26: 1481-1496. 10.1080/03610929708831995.
3218.0 Regional Population Growth, Australia and New Zealand, 1991 to 2001. 2002, Australian Bureau of Statistics, Canberra, Australia
Population Health Profile of the Adelaide Western Division Of general Practice. Population Profile Series Number 86, Public Health Information Development Unit, University of Adelaide. 2005, Australian Institute of Health and Welfare, Adelaide, Australia
Parker R, Ben-Tovim DI: A study of factors affecting suicide in Aboriginal and 'other' populations in the Top End of the Northern Territory through an audit of coronial records. Australian and New Zealand Journal of Psychiatry. 2002, 36: 404-410. 10.1046/j.1440-1614.2001.01032.x.
Manne R: Pearson’s Gamble, Stanner’s Dream: The Past and Future of Remote Australia. The Monthly. 2007, http://www.themonthly.com.au/past-and-future-remote-australia-pearson-s-gamble-stanner-s-dream-robert-manne-590[11November2011][11November2011].,
Tiwi Islands deal with tragedy of youth suicide. 1999, Australian Broadcasting Corporation, Australia, http://www.abc.net.au/worldtoday/stories/s43678.htm [11 November 2011]
Franklin J: The Cultural Roots of Aboriginal Violence. Quadrant. 2008, 451: -http://www.quadrant.org.au/magazine/issue/2008/451/the-cultural-roots-of-aboriginal-violence [11 November 2011]
Mornington Island's suicide problem. 2000, Australian Broadcasting Corporation, Australia, http://www.abc.net.au/pm/stories/s105887.htm [11 November 2011]
Chen J, Roth RE, Naito AT, Lengerich EJ, Maceachren AM: Geovisual analytics to enhance spatial scan statistic interpretation: an analysis of U.S. cervical cancer mortality. Int J Heal Geogr. 2008, 7: 57-10.1186/1476-072X-7-57.
Glaz J, Balakrishnan N: Scan Statistics and Applications. 1999, Birkhäuser, Boston, United States of America
Olfson M, Shaffer D, Marcus SC, Greenberg T: Relationship between antidepressant medication treatment and suicide in adolescents. Arch Gen Psychiatry. 2003, 60: 978-982.
Hall WD, Mant A, Mitchell PB, Rendle VA, Hickie IB, McManus P: Association between antidepressant prescribing and suicide in Australia, 1991–2000: trend analysis. British Medical Journal. 2003, 326: 1008-10.1136/bmj.326.7397.1008.
Lester D, Leenaars A: Suicide rates in Canada before and after tightening firearm control laws. Psychological Reports. 1993, 72: 787-790. 10.2466/pr0.19188.8.131.527.
Ludwig J, Cook PJ: Homicide and suicide rates associated with implementation of the Brady Handgun Violence Prevention Act. JAMA. 2000, 284: 585-591. 10.1001/jama.284.5.585.
Lester D: The effect of the detoxification of domestic gas in Switzerland on the suicide rate. Acta Psychiatr Scand. 1990, 82: 383-384. 10.1111/j.1600-0447.1990.tb01406.x.
Wiedenmann A, Weyerer S: The impact of availability, attraction and lethality of suicide methods on suicide rates in Germany. Acta Psychiatr Scand. 1993, 88: 364-368. 10.1111/j.1600-0447.1993.tb03474.x.
Rihmer Z, Belso N, Kalmar S: Antidepressants and suicide prevention in Hungary. Acta Psychiatr Scand. 2001, 103: 238-239.
Marusic A, Roskar S, Dernovsek M, et al: An attempt of suicide prevention: the Slovene Gotland Study. 2004, In Program and abstracts of the 10th European Symposium on Suicide and Suicidal Behavior, Copenhagen, Denmark, ICS A/S; 2004
Jorm AF, Christensen H, Griffiths KM: The impact of beyond blue: the national depression initiative on the Australian public’s recognition of depression and beliefs about treatments. Australian and New Zealand Journal of Psychiatry. 2005, 39: 248-254.
Paykel ES, Hart D, Priest RG: Changes in public attitudes to depression during the Defeat Depression Campaign. Br J Psychiatry. 1998, 173: 519-522. 10.1192/bjp.173.6.519.
Mann JJ, Apter A, Bertolote J, Beautrais A, Currier D, Haas A, Hegerl U, Lonnqvist J, Malone K, Marusic A, Mehlum L, Patton G, Phillips M, Rutz W, Rihmer Z, Schmidtke A, Shaffer D, Silverman M, Takahashi Y, Varnik A, Wasserman D, Yip P, Hendin H: Suicide prevention strategies: a systematic review. JAMA. 2005, 294: 2064-2074. 10.1001/jama.294.16.2064.
Lambert G, Reid C, Kaye D, Jennings G, Esler M: Increased suicide rate in the middle-aged and its association with hours of sunlight. Am J Psychiatry. 2003, 160: 793-795. 10.1176/appi.ajp.160.4.793.
Nicholls N, Butler CD, Hanigan I: Inter-annual rainfall variations and suicide in New South Wales, Australia, 1964–2001. Int J Biometeorol. 2006, 50: 139-143. 10.1007/s00484-005-0002-y.
Huisman M, Oldehinkel AJ: Income inequality, social capital and self-inflicted injury and violence-related mortality. Journal of Epidemiology and Community Health. 2009, 63: 31-37.
Chan WS, Yip PS, Wong PW, Chen EY: Suicide and unemployment: what are the missing links?. Archives of Suicide Research. 2007, 11: 327-335. 10.1080/13811110701541905.
Fragar LJ, Gray EJ, Franklin RJ, Petrauskas V: A picture of health: a preliminary report of the health of country Australians. 1997, The Australian Agricultural Health Unit, P.O. Box 256, Moree, NSW 2400, Australia
Fuller J, Kelly B, Sartore G, Fragar L, Tonna A, Pollard G, Hazell T: Use of social network analysis to describe service links for farmers' mental health. Aust J Rural Heal. 2007, 15: 99-106. 10.1111/j.1440-1584.2007.00861.x.
Sartore G: Drought and its effect on mental health–how GPs can help. Australian Family Physician. 2007, 36: 990-993.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-244X/12/86/prepub
XQ was funded by the QUT Postgraduate Research Award Scholarship, and ST was supported by a NHMRC Research Fellowship (#553043). We also thank Dr. Lyle Turner from Queensland University of Technology for examining the written languages of the manuscript.
We declare that we have no competing interests.
XQ designed the study, implemented all statistical analyses and drafted the manuscript. ST conceptualised the idea and revised the study protocol, especially the research design and data analysis. WH provided advice on statistical analyses and interpretation of the results. AP helped interpret the results and assisted in drafting the manuscript. All the authors contributed to the preparation of the final manuscript and approved the submission.