A Markov chain model for studying suicide dynamics: an illustration of the Rose theorem
© Yip et al.; licensee BioMed Central Ltd. 2014
Received: 23 November 2013
Accepted: 12 June 2014
Published: 19 June 2014
High-risk strategies would only have a modest effect on suicide prevention within a population. It is best to incorporate both high-risk and population-based strategies to prevent suicide. This study aims to compare the effectiveness of suicide prevention between high-risk and population-based strategies.
A Markov chain illness and death model is proposed to determine suicide dynamic in a population and examine its effectiveness for reducing the number of suicides by modifying certain parameters of the model. Assuming a population with replacement, the suicide risk of the population was estimated by determining the final state of the Markov model.
The model shows that targeting the whole population for suicide prevention is more effective than reducing risk in the high-risk tail of the distribution of psychological distress (i.e. the mentally ill).
The results of this model reinforce the essence of the Rose theorem that lowering the suicidal risk in the population at large may be more effective than reducing the high risk in a small population.
KeywordsAn illness and death model Markov chain model Suicide Rose theorem
Suicide has become a major public health issue throughout the world. About a million people kill themselves every year, and more than half of these cases occur in Asia [1, 2]. It is estimated that about 2.5% of loss of disability-adjusted life year (DALY) is due to suicide and deliberate self-harm . Suicide is especially a major concern in Asia due to its large population size with a relatively high suicide rate and limited resources in suicide prevention [3, 4]. In response to the growing concern, different national strategies have been implemented in a number of countries to reduce the number of suicides and the results are somewhat mixed . The programs in Australia, Finland, Norway and Sweden had little or no impact on reducing suicide rates among youth and the general population. However, some more promising results have recently observed in South Korea and Taiwan. There are continuously debates in identifying cost-effective approaches in preventing suicides. Suicide has traditionally been viewed as a mental health issue that is addressed primarily through clinical intervention, especially by providing services for the treatment of depression and other mental illnesses (bipolar disorder, schizophrenia). However, it has been suggested that the role of mental illness in suicide risk is not as significant as expected, especially in the East [1, 6–8]. Approximately two-thirds of all people who commit suicide did not receive any specialist psychiatric care in the year before their death . The World Health Organization (WHO) and many national suicide prevention strategies (for example, those in the USA, Austria, Australia, Ireland, New Zealand, and the UK) have proposed a public health approach for suicide prevention, rather than treating it as a medical problem only [3, 10–13]. The public health approach involves three layers of intervention: universal, selective and indicated. This public health approach acknowledges the importance of both high-risk and population-based strategies of suicide prevention, and requires a multi-sector effort to tackle the problem at multiple levels: in the community (universal strategies), among specific population subgroups (selective strategies), and among those at a particularly high risk of suicidal behavior (indicated strategies). The public health approach is particularly apt for suicide prevention in Asia where mental health services have not been well developed and awareness and services for depression and mental illness in the community remains inadequate. The Rose Theorem states that a large number of people exposed to a low risk may generate more cases than a small number of people who are exposed to a high risk .
Here we adopt an illness and death model [15–17] to model the transition of suicide risk in the population. We make use of the Hong Kong mental morbidity and suicide data to illustrate the applicability of this model. Some empirical results to demonstrating the effectiveness of suicide prevention effort by modifying some parameters of the Markov model will be provided. The underlying framework of this study is not only applicable to suicide research, but could also be extended for other disease prevention purposes and public health challenges.
A suicide dynamic model
where w n , x n , y n , z n are respectively the population that is healthy, mentally ill, dead people due to other causes, and dead people due to suicide, in the n-th year.
Where R = (1 + p 31 + p 41)(p 12 + p 32 + p 42) + p 21(1 + p 32 + p 42), and Q 0 is the initial total population size (for example, say around 7 million for Hong Kong in 2011). In other words, Eq. (3) is the unique equilibrium of the dynamic system of people in various stages in the Markov model in Figure 2. Furthermore, given the explicit formulas in Eq. (3), the coefficients of assessing the effects on the number of suicides can also be computed by simply taking partial derivatives.
And similarly defined for the other parameters, p ij .
Based on the morbidity data of mental illness supplied by the Hospital Authority of the Hong Kong Government (2012), there were about 150,000 individuals who have received psychiatric service in the government funded hospitals. And there were about 20,000 new psychiatric cases out of the remaining 6.85 million healthy persons annually. Based on the Coroner court files it is estimated that about 30% of the population have received psychiatric treatment and the suicide rate for mentally ill patients was estimated to be 250 per 100,000, while the suicide rate among the non-psychiatric person is about 8.34 per 100,000, i.e., a 30 times differential between psychiatric and non-psychiatric person, which is quite consistent with figures in other countries . The recovery rate from mental illness was estimated at 6.67 per 100 persons and the overall death rate is about 6.0 per 1000 . Hence, the estimated values of the coefficients of p ij for the Markov model in Figure 2 are given as follows:
p 21 = 0.00286 (incidence rate of mentally ill),
p 12 = 0.0667 (incidence rate of recovery among the mentally ill),
p 31 = 0.006 (death rate from causes other than suicide among the healthy),
p 32 = 0.012 (death rate from causes other than suicide among the mentally ill),
p 41 = 0.0000834 (suicide rate among the healthy),
p 42 = 0.0025 (suicide rate among the mentally ill),
and Q 0 = 7 million in 2012.
From Eq. (3): then one gets
w = 0.960 × Q 0 (healthy)
x = 0.0338 × Q 0 (mentally ill)
y = 0.00616 × Q 0 (death from causes other than suicide)
and z = 0.000165 × Q 0 (suicide).
In other words, the Markov dynamic system suggested that the percentage of mentally ill rises from the initial stage of 2.14% and stabilizes at about 3.4% of total population in an equilibrium stage, and the risk of suicides stabilizes at about 16.5 per 100,000 and the death rate from other causes 6.16 per 1,000, which are quite consistent to the recent figures .
The magnitude of the difference is about 30 times (67.2/2.3) between the non-mentally ill and the mentally ill groups. If the cost of making reduction of risk for health and mentally ill is the same, it would lead to the same reduction in suicide number for producing 30 times effects in the mentally ill population with reference to the non-mentally ill group. This is due to the relatively smaller population size of mentally ill in comparing to the healthy population.
Likewise, the suicide rate for different sets of parameters can be computed. Here, we give the number of suicide (z) and its changes (Δz) against variations of the suicide risks of healthy (p 41 ) and mentally ill population (p 42 ), with other parameters held constant as above.
Sensitivity analyses of number of suicides respective to the parameters of p 41 and p 42
p 42 /p 41
Δp 42 b
p 42 fixed at 0.0025
p 41 fixed at 0.0000834
p 42 /p 41 fixed at 30
The analysis of our model suggested that the effect of reducing a relatively small suicidal risk for a large population is more effective than in reducing suicide risk among the mentally ill. The result is consistent with Lewis et al.  who have shown that high-risk (indicative) strategies would only have a modest effect on suicide prevention within a population, even if effective interventions could be developed, and that the UK government’s target for suicide reduction was more likely to be achieved using population-wide strategies aimed at lowering risk among the whole population. In contrast to the high-risk approach, a population based approach shall be more radical and with larger influence on the general population. As envisioned by Geoffrey Rose, the population strategy of prevention seeks to shift the entire population’s distribution of “suicide tendency” toward a more favorable direction. This includes promoting the mental health of individuals who are not currently identified as being above a cut-point or a threshold for mental illness or distress.
Indeed our model identifies, theoretically, a reduction of suicide incidence in the healthy population would be more “effective” but in practice, this would only be true if it is at least as practicable to reduce the incidence in the healthy population as compared to the mentally ill. The cost to preventing suicide among mentally ill by intensive clinical care and management would be high but could be more effective in reducing the number of suicide. However, the impact on reduction of number of suicide would still be small as suicide is still a rare event despite the excessive suicidal risk among the mentally ill. It needs to be at least 30 times more effective than in the general population in order to be able to save the same number of suicides. Indeed those who suffered from mentally ill in our community needs better clinical care and treatment but the impact on the number of suicides would still be limited. Furthermore, there are some promising universal programs, for example, restriction of means , responsible media reporting [24, 25] and some community based programs  have demonstrated their effectiveness in reducing suicide numbers in the community and the cost is not that expensive. However, it does rely on participation and co-operation of stakeholders in the community. It is important to develop and evaluate the effectiveness of some population based suicide prevention strategies. Potential effective population approaches for suicide prevention have been so far understudied. It seems that community based support and participation in preventing suicides is not only desirable but also essential [27–29].
There are a few limitations of this study. First, we assume there is a one to one replacement of a death by a living individual. It is a mathematical assumption that allows us to compute the state of equilibrium. It facilitates the computation rather than occurring in real life necessarily. Nevertheless, the essence of the results of the analysis would not be much affected. Second, the results would also depend on the values of p ij , which could be very different in various countries. Hence the results themselves are country specific and may differ between contexts. However, the overall effect should be similar so long as the relationship between the suicide rates among the healthy versus the mentally ill is in the range we have assumed. Also, our calculations do not reflect the differential costing of making a unit change of the death rate from mental ill and healthy individuals. However, it is likely that the cost to reducing a unit of suicide risk among the mentally ill may be more expensive than that of healthy ones but their effect to reducing number of suicides is still uncertain. Our model can be expanded to a formal cost-utility analysis if the specific intervention programs for each strategy (high risk approach versus population approach) could be cost.
Suicide rate is unlikely to decline sharply through unidimensional effort on preventing the mentally ill patients from suicide, given their high suicidal risk. A more integrated approach of mental health service is needed to be effective to reducing the number of suicide among the mentally ill [28, 29]. Furthermore, clinical care for mentally ill is labour intensive and on individual base, its effect on reducing number of suicides has yet to be established. The sensitivity analysis of this model suggests that lowering the suicidal risk in the population at large would be more effective than reducing the high risk in a small population. It underscores the importance of population-based strategies to suicide prevention.
The research is supported by a Strategic Public Policy Research (HKU-SPPR-12) and Chief Executive Community Project funded by the Hong Kong Jockey Club Charity Trust. The useful comments by the two reviewers are gratefully acknowledged.
- Chen YY, Wu KC, Yousuf S, Yip PSF: Suicide in Asia: opportunities and challenges. Epidemiol Rev. 2012, 34: 129-144.View ArticlePubMedGoogle Scholar
- Suicide in Asia: causes and prevention. Edited by: Yip PSF. 2008, Hong Kong: Hong Kong University PressGoogle Scholar
- World Health Organization: Public health action for the prevention of suicide: a framework. http://apps.who.int/iris/bitstream/10665/75166/1/9789241503570_eng.pdf,
- Yip PSF, Liu KY, Law CK: Years of life lost from suicide in china, 1990–2000. Crisis. 2008, 29: 131-136.View ArticlePubMedGoogle Scholar
- De Leo D, Milner A, Xiangdong W: Suicidal behavior in the western pacific region: characteristics and trends. Suicide Life Threat Behav. 2009, 39: 72-81.View ArticlePubMedGoogle Scholar
- Yip PSF, Liu KY, Hu J, Song XM: Suicide rates in china during a decade of rapid social changes. Soc Psychiatry Psychiatr Epidemiol. 2005, 40: 792-798.View ArticlePubMedGoogle Scholar
- Phillips MR, Li X, Zhang Y: Suicide rates in china, 1995–99. Lancet. 2002, 359: 835-840.View ArticlePubMedGoogle Scholar
- Chen EY, Chan WS, Wong PW, Chan SS, Chan CL, Law YW, Beh PS, Chan KK, Cheng JW, Liu KY, Yip PSF: Suicide in hong Kong: a case–control psychological autopsy study. Psychol Med. 2006, 36: 815-825.View ArticlePubMedGoogle Scholar
- Law YW: PhD thesis. Silent suicides: studies on the non-contact group of suicide. 2012, The University of Hong Kong, Department of Social Work and Social AdministrationView ArticleGoogle Scholar
- Hoven CW, Wasserman D, Wasserman C, Mandell DJ: Awareness in nine countries: a public health approach to suicide prevention. Leg Med (Tokyo). 2009, 11 (Suppl 1): S13-S17.View ArticleGoogle Scholar
- Ueda K, Matsumoto Y: National strategy for suicide prevention in Japan. Lancet. 2003, 361: 882-View ArticlePubMedGoogle Scholar
- Yip PSF, Law YW: Multidisciplinary and interdisciplinary approaches to suicide prevention. Suicide From A Global Perspective: Public Health Approaches. Edited by: Shrivastava A, Kimbrell M, Lester D. 2012, New York: Nova Science Publishers, 63-74.Google Scholar
- Yip PSF: A public health approach to suicide prevention. Hong Kong J Psychiatry. 2005, 15: 29-31.Google Scholar
- Rose G: The Strategy of Preventive Medicine. 1992, Oxford: Oxford University PressGoogle Scholar
- Chiang CL: An Introduction to Stochastic Processes and Their Applications. 1980, New York: KriegerGoogle Scholar
- Lagakos SW: A stochastic model of censored survival data in the presence of an auxiliary variable. Biometrics. 1976, 32: 551-559.View ArticlePubMedGoogle Scholar
- Yip PSF: An inference procedure for a capture and recapture experiment with time-dependent capture probabilities. Biometrics. 1989, 45: 471-479.View ArticleGoogle Scholar
- Yip PSF, Lam KF: An inference procedure for an illness and death model. Stoch Anal Appl. 1997, 15: 125-135.View ArticleGoogle Scholar
- Yim PH, Yip PSF, Li RH, Dunn EL, Yeung WS, Miao YK: Suicide after discharge from psychiatric inpatient care: a case–control study in hong Kong. Aust N Z J Psychiatry. 2004, 38: 65-72.View ArticlePubMedGoogle Scholar
- Appleby L, Shaw J, Amos T, McDonnell R, Harris C, McCann K, Kiernan K, Davies S, Bickley H, Parsons R: Suicide within 12 months of contact with mental health services: national clinical survey. BMJ. 1999, 318: 1235-1239.View ArticlePubMedPubMed CentralGoogle Scholar
- Hong Kong Census and Statistics Department: Hong Kong Annual Digest of Statistics 2012 Edition. 2012, Hong KongGoogle Scholar
- Lewis G, Hawton K, Jones P: Strategies for preventing suicide. Br J Psychiatry. 1997, 171: 351-354.View ArticlePubMedGoogle Scholar
- Yip PSF, Caine E, Yousuf S, Chang SS, Wu KC, Chen YY: Means restriction for suicide prevention. Lancet. 2012, 379: 2393-2399.View ArticlePubMedGoogle Scholar
- Fu KW, Chan YY, Yip PSF: Newspaper reporting of suicides in Hong Kong, Taiwan and Guangzhou: compliance with WHO media guidelines and epidemiological comparisons. J Epidemiol Community Health. 2011, 65: 928-933.View ArticlePubMedGoogle Scholar
- Niederkrotenthaler T, Fu KW, Yip PS, Fong DY, Stack S, Cheng Q, Pirkis J: Changes in suicide rates following media reports on celebrity suicide: a meta-analysis. J Epidemiol Community Health. 2012, 66: 1037-1042.View ArticlePubMedGoogle Scholar
- Wong PW, Fu KW, Chan KY, Chan WS, Liu PM, Law YW, Yip PSF: Effectiveness of a universal school-based programme for preventing depression in Chinese adolescents: a quasi-experimental pilot study. J Affect Disord. 2012, 142: 106-114.View ArticlePubMedGoogle Scholar
- Knox KL, Litts DA, Talcott GW, Feig JC, Caine ED: Risk of suicide and related adverse outcomes after exposure to a suicide prevention programme in the US Air force: cohort study. BMJ. 2003, 327: 1376-View ArticlePubMedPubMed CentralGoogle Scholar
- While D, Bickley H, Roscoe A, Windfuhr K, Rahman S, Shaw J, Appleby L, Kapur N: Implementation of mental health service recommendations in England and wales and suicide rates, 1997–2006: a cross-sectional and before-and-after observational study. Lancet. 2012, 379: 1005-1012.View ArticlePubMedGoogle Scholar
- Conwell Y, Farley-Toombs C: Evidence supports the obvious: suicides need not happen. Lancet. 2012, 379: 981-982.View ArticlePubMedPubMed CentralGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/14/625/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.