Camden and Islington Mental Health and Social Care Trust and Department of Psychiatry and Behavioural Sciences, Royal Free and University College London Medical Schools, London
Department of Biostatistics and Computing Institute of Psychiatry, London
National Confidential Inquiry into Suicide and Homicide by People with Mental Illness, School of Psychiatry and Behavioural Sciences, University of Manchester, UK
Correspondence: Dr Nigel McKenzie, Highgate Mental Health Centre, Dartmouth Park Hill, London N19 5JG, UK. E-mail: n.mckenzie{at}ucl.ac.uk
![]() |
ABSTRACT |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Aims To determine whether space time and spacetimemethod clustering occur in a national case register of those who had recent contact with mental health services and had died by suicide and to estimate the suicide imitation rate in this population.
Method Knox tests were used for spacetime and spacetimemethod clustering. Model simulations were used to estimate effect size.
Results Highly significant spacetime and spacetimemethod clustering was found in a sample of 2741 people who died bysuicide over 4 yearswho hadhadrecent contact with one of 105 mental health trusts. Model simulations with an imitation rate of 10.1% (CI 417) reproduced the observed spacetimemethod clustering.
Conclusions This study provides indirect evidence that imitative suicide occurs among people with mental illnesses and may account for about 10% of suicides by current and recent patients.
![]() |
INTRODUCTION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
METHOD |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Since 1996 information on deaths with a verdict of suicide or an open verdict in a coroners court has been forwarded to the NCI, who then submitted identifying details to the main hospitals or trusts providing mental health services in the victims health district. Hospital records were checked to identify those who had had contact with mental health services in the 12 months preceding their death by suicide. A questionnaire was sent to the responsible medical officer (RMO) requesting further information about the suicide and care provided in the period before death. For the purposes of this study, date of death, method of suicide and coded identities for trust and person completing the questionnaire were used to test for clustering of suicide in time, space and by method. Prior ethical approval was obtained.
To investigate clustering, all possible pairs of suicides were considered and, following Knox (1964), the number of pairs close in space and time (or space, time and method) according to chosen criteria was taken as the test statistic. Knox showed that under certain assumptions this statistic follows a Poisson distribution under the null hypothesis of independence of suicide location and time. A permutational approach suggested by Mantel (1967) enables the distribution of the test statistic to be derived empirically, avoiding such assumptions. The spatial labels of the suicides are randomly permuted while holding the time labels fixed (or vice versa). The number of close pairs is calculated for each permutation. A one-sided P value of the test is given by:
P=(1+number of permutations where value of test statisticsobserved value)/(1+number of permutations).
Similarly, to test the null hypothesis of independence of suicide location, time and method the labels of two of the variables can be independently permuted to derive the distribution of the spacetimemethod test statistics under the null hypothesis.
The Knox procedure required selection of criteria for closeness in space, time and method.
Closeness in space
The selection of a criterion for closeness in space required taking into
account the model of suggestion as a cause for clustering:
closeness in space should define an appropriate communication
unit whose members become aware of the suicide of one of their number
and may go on to imitate the suicidal behaviour. It was assumed that patients
meet and interact socially primarily at the level of a geographical sector
served by a single community mental health team and ward team under the
clinical leadership of a consultant psychiatrist (the RMO). Some contact would
be expected between patients of adjacent sectors within a single trust,
allowing news of a suicide to spread within a trust. Data on sectors were not
collected as such, although where the RMO completed the NCI questionnaire it
was possible to use the identity of the RMO as a proxy for sector. This had
certain limitations: the RMO did not always complete the questionnaire,
leading to potential gaps in the data, and it was evident from descriptive
analysis of the data that there was a fairly high rate of turnover of RMOs, so
that the same RMO was not necessarily covering the same sector for the whole
period of data collection. The trust was therefore our primary choice as a
variable for categorising communication units, and pairs of cases were defined
as close in space if they occurred in the same trust. We repeated the analysis
defining suicide cases as close in space if their suicide was recorded by the
same RMO in the same trust.
A further consideration was mergers between trusts. Trusts were set up in the mid-1990s by the then government as part of the creation of the internal market in healthcare. They were typically based on the services provided by one or two local hospitals. The current government made changes to the commissioning of healthcare and encouraged trusts to merge into larger units. A number of the merged trusts comprised geographically dispersed community teams and in-patient units much larger than the ideal communication unit referred to previously. It was unlikely that news of a suicide would spread through all the different constituent sites. Hence it was decided to include only trusts that did not merge before the end of the study period. This also reduced the possibility that changes in management structure could have given rise to gaps in identifying cases that had been in contact with mental health services.
Closeness in time
There was no a priori principle on which to base the criterion of
closeness in time. It might be expected that news of a suicide would take some
time to disseminate through the patient population and the recollection of the
suicide would remain with a patient for some time and might influence suicidal
behaviour some time after the index event. This might happen, for instance, if
a patient later experienced a period of low mood and suicidal ideation and the
previous suicide seemed to offer a way out. It seemed plausible that a suicide
could influence behaviour for several months and it could even be that the
1-year anniversary of a suicide might influence another patient to take their
own life. The test procedure was repeated for a range of plausible threshold
values for closeness from 30 to 360 days. As the different
threshold test statistics are highly correlated, the significance level should
not be greatly affected by multiple testing.
Closeness in method
Suicides were defined as close in method if the method employed was the
same using the classifications given in
Table 1. The percentages of
suicides according to method in the sample studied are also shown. Cases for
which the suicide method did not fall into one of the broad categories or was
not known were excluded from the assessment of spacetimemethod
clustering.
|
Choice of study period
Since systematic gaps in the data could also give rise to spacetime
clustering, steps were taken to ensure that the data were as complete as
possible. The NCI assessed the accuracy of detecting a previous contact with
mental health services and found a 97% detection rate
(Appleby et al, 1999).
By comparing the accumulated sample at two points 1 year apart, mid-2001 and
mid-2002, a period from February 1996 to February 2000 was identified when the
annual number of suicides had built up to a fairly constant level, indicating
that data collection was approximately complete. As additional safeguards, to
ensure full reporting: (a) only those trusts were selected that had a first
case on or before the first day of the study period and a last case on or
after the last day; (b) where trusts subsequently merged the merged trust did
not have a first case before February 2000. The optimum study period was
chosen so that it maximised the number of suicides to be included within these
limits.
Estimation of effect size
If significant clustering were found it would be important to estimate the
effect size. The non-parametric test did not automatically provide estimates
of parameters that could lead to an estimate of numbers of imitative suicides.
However, the test statistic (observed number of close pairs) and its empirical
distribution under the null hypothesis provide some information about suicide
imitation parameters.
We defined an excess pair statistic as the difference between the observed number of close pairs and the number expected under independence. This is affected by the delay time (between index case and imitative suicide) and rate of imitative suicide. Assuming that imitative suicides occur in the same space unit (and by the same method) as the index case we expect that the excess pairs will reach a maximum when the threshold used to define closeness in time approaches the true maximum delay in imitative suicide, T: with increasing time threshold the observed number of close pairs, and hence the excess pairs statistic, gradually includes more imitative suicides close to their index cases until T is reached. However, as the time threshold increases beyond T, more and more pairs involving imitative cases are also included in the expected number of pairs under independence and hence excluded from the excess pairs statistic. The combined effect of these two opposing mechanisms should result in a maximum value for excess pairs at time threshold T. It can be shown (see data supplement 1 to the online version of this paper) that under certain restrictive assumptions the excess pairs statistic at threshold T provides an estimate of the number of imitative suicides and the relative excess (number of excess pairs divided by the sample size) an estimate of the suicide imitation rate.
To obtain an unbiased estimate of the suicide imitation rate and to quantify its precision we used simulation models. This approach (see data supplement 2 to the online version of this paper) entails simulating values of the test statistics from a suicide model with a given imitation rate to generate a distribution under the model. Such distributions are generated for a range of possible suicide rates and then the suicide rate is estimated by the rate of the model that fits the observed value of the test statistic most closely. Since each computer simulation took an appreciable time to complete, we limited the number of simulations to 200 for each possible suicide rate. An attractive feature of the chosen procedure for simulating is that it maintains the marginal distribution of suicide times and locations and can be thought of as a generalisation of the Mantel permutation procedure.
![]() |
RESULTS |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Spacetime clustering and spacetimemethod clustering were tested for separately and the results are shown in Tables 2 and 3. Each table shows the total number of possible distinct pairs of suicides and the observed and expected numbers of close pairs for increasing thresholds of closeness in time. Significant spacetime clustering (Table 2) and spacetimemethod clustering (Table 3) were found for time thresholds from 30 to 360 days.
|
|
The relative excess pairs close in space and time (Table 2) provides an estimate of the suicide imitation rate and increases from 3.8% at 30 days to reach a maximum value of 13.7% at a 210-day time threshold. (The pattern of steady increase to a maximum value followed by decrease remained when values of relative excess pairs were calculated for delay times <30 days and >360 days.) Assuming that imitative suicide is the sole reason for spacetime clustering and such suicides occur in the same trust as the index cases, the maximum delay between an index case and an imitative case can be estimated as in the region of 69 months. A model simulation with a maximum imitation delay of 7 months gave an estimation of 13.3% (95% CI 322) for imitative suicides as a percentage of all suicides that copy the act of suicide of an index case but not necessarily the method of the index case.
The relative excess pairs close in space, time and method (Table 3) reaches a maximum value of 13.0% at a 300-day time threshold. Assuming a true maximum delay of 10 months, the model simulation including method gave an estimation of 10.1% (95% CI 417) for imitative suicides as a percentage of all suicides that copy the act and method of suicide of an index case.
The clustering analysis was repeated using RMO as the space variable. The optimum study period was determined as 845 days, during which 328 RMOs reported 888 cases of suicide. Spacetime clustering was again highly significant for time thresholds from 60 to 360 days. Spacetimemethod clustering did not reach significance, perhaps because reduced numbers limited the power to detect clustering. The relative excess pairs statistic reached a maximum value of 10.2% at a time threshold of 8 months.
![]() |
DISCUSSION |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Imitation as cause of clustering
The observed clustering might have been caused by several factors operating
singly or together. The first of these is imitation of suicidal behaviour. If
this were the sole cause of the clustering, a model used to simulate the
effect of imitation gave a possible effect size of about 10% (95% CI
417) of suicides imitating the method of and being close in time to an
index case in the same trust. Imitations appear to build up in number steeply
initially and then level off over a 7- to 10-month time scale.
Strengths of study
A strength of the study is the much larger numbers of cases and locations
analysed than in previous studies, leading to greater statistical power to
detect clustering. The methodology also has the advantage of being sensitive
only to spacetime or spacetimemethod interactions and so
is not confounded by local differences in rates or method of suicide that do
not change over the study period, or changes over time affecting all locations
equally, such as seasonal variations
(Preti, 2000;
Hakko et al,
2002).
Other possible causes of observed clustering
A weakness of the study, shared by other studies of clustering, is that the
evidence for imitative suicide is indirect and other causes for the observed
clustering cannot be ruled out.
Quality of care or socio-economic conditions
A change in local factors, such as the quality of care or socio-economic
conditions, that alters the suicide rate in some trusts but not others can
result in timespace clustering. It is less plausible, however, that
this mechanism on its own could also account for the observed
spacetimemethod clustering of suicides. The time scale of about
9 months suggested by the analysis, with clustering also observed at time
thresholds down to 30 days, seems too short for differential changes in the
quality of care in trusts or other local factors affecting the suicide rate to
have occurred.
Missing data
Systematic gaps in data collection can also give rise to apparent
clustering. This possibility was minimised by including only trusts that
identified a first case on or before the start of the study period and a last
case on or after the last day, thereby ensuring as far as possible that the
trusts had systems in place for reporting during the whole of the study
period. Trusts that merged during the study period were excluded, thereby
eliminating possible gaps in reporting caused by changes in management
structure after a merger. In addition the NCI conducted an audit of the
accuracy of reporting by trusts and found a 97% identification rate of cases
(Appleby et al,
1999).
Coronerscourts
Variations over time between coroners courts in identifying suicides
and cause of death could also cause apparent clustering. It seems unlikely,
however, that there could have been sufficient variation between
coroners courts in identifying cases on the timescale suggested by the
data.
Findings from previous studies
Support for imitation as an explanation of the observed clustering of
suicides among people in contact with mental health services is given by
studies which have explored imitation of suicidal behaviour in the general
population. It seems likely that imitation would occur to an equal or greater
degree among people with mental illnesses. Various mechanisms have been
proposed: low mood and low self-esteem may render an individual less able to
resist copying a behaviour that seems to offer a way out. Of three previous
quantitative studies of clustering of suicides among those with mental
illnesses only one found significant clustering
(Haw, 1994) although two found
clinical evidence suggesting that imitation had occurred
(Modestin & Würmle,
1989; Taiminen & Helenius,
1994). The latter studies may have had sample sizes that were too
small to detect clustering that was present.
Conclusion
If imitation is implicated as a causal factor in a significant percentage
of suicides, it will be important to consider how best to reduce its impact as
part of a drive to cut the national suicide rate among people with mental
illnesses (Department of Health,
2002). Suggestions for prevention of suicide
epidemics were made by Rissmiller & Rissmiller
(1990) but more research is
required to identify effective strategies, and parallel efforts should be made
to raise mental health professionals awareness of this phenomenon.
![]() |
DATA SUPPLEMENT 1 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Assumption A4 that an index case gives rise to only one imitative case is an approximation. However, for the observed imitation rate of about 10%, the probability of an index case giving rise to two or more imitations decreases rapidly with the number of imitative cases.
A suicide can be classified as either imitative or spontaneous, where spontaneous suicides occur purely by chance.
By assumptions A1 to A3, pairs of cases can be close in time and space either: (1) as a result of imitation (such pairs consist of an imitative suicide and its index case) or (2) by chance.
The expected number of close pairs P is given by:
![]() |
The observed number of pairs that are close in space and time, O, is an unbiased estimate of the theoretical quantity P.
The mean number of close pairs of the permutation distribution, E, resulting from the Mantel procedure carried out with the whole sample size, n, is an estimate of the number of close pairs expected by chance under the independence model. Provided imitative suicides are rare, E can be considered an estimate of P(2), that is the number of close pairs expected by chance under the suicide imitation model. However, in this case the estimator suffers from upwards bias since the Mantel procedure counts all close pairs, including those from permutations that allocate imitative cases close in space and time to their respective index cases.
Owing to assumption A4, the number of close pairs due duetoimitationis, to
imitation is, in fact, the number of imitative suicides in the sample,
s. Thus P(1) = s and a downwards-biased
estimate of s is given by the excess pairs statistic
![]() |
![]() |
![]() |
DATA SUPPLEMENT 2 |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
With assumptions A1 to A3 from Appendix 1, the following algorithm generates a set of suicide locations and times for n cases that are consistent with the specified model and maintain the observed marginal distributions of suicide time and location.
The set of n cases was chosen as the set of observed cases within the selected study period. A run-in period had to be specified during which suicide times, locations and methods remained fixed to ensure that there were index cases available for imitation at the start. The chosen run-in period was T days prior to the start of the study period.
The procedure was extended to accommodate the additional assumption that an imitative suicide employs the same method as the index case. A set M of methods is defined containing each method as often as it has been observed. A spontaneous case is allocated a method at random without replacement. If, following the procedure above, the selected identity specifies an imitative suicide and a spatial location has been selected, the method is chosen to be that previously allocated to this space unit. Should the method no longer be available in the set of methods, the choice of space unit (and associated index case) is discarded and a new space unit chosen until a match is found.
![]() |
Clinical Implications and Limitations |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
LIMITATIONS
![]() |
REFERENCES |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Appleby, L., Shaw, J., Amos, T., et al
(1999) Suicide within 12 months of contact with mental health
services: national clinical survey. BMJ,
318, 1235
1239.
Department of Health (2002) National Suicide Prevention Strategy for England. London: Department of Health.
Gould, M. S. (2001) Suicide and the media.
Annals of the New York Academy of Sciences,
932, 200
224.
Gould, M. S., Wallenstein, S., Davidson, L. (1989) Suicide clusters: a critical review. Suicide and Life-Threatening Behavior, 19, 17 29.[Medline]
Hakko, H., Rasanen, P., Tiihonen, J., et al (2002) Use of statistical techniques in studies of suicide seasonality, 1970 to 1997. Suicide and Life-Threatening Behavior, 32, 191 208.[CrossRef][Medline]
Haw, C. M. (1994) A cluster of suicides at a London psychiatric unit. Suicide and Life-Threatening Behavior, 24, 256 266.[Medline]
Joiner, T. E. (1999) The clustering and contagion of suicide. Current Directions in Psychological Science, 8, 89 92.[CrossRef]
Knox, G. (1964) Epidemiology of childhood leukaemia in Northumberland and Durham. British Journal of Preventive and Social Medicine, 18, 25 34.[Medline]
Mantel, N. (1967) The detection of disease clustering and a generalized regression approach. Cancer Research, 27, 200 220.
Modestin, J. & Würmle, O. (1989) Role of modelling in in-patient suicide: a lack of supporting evidence. British Journal of Psychiatry, 155, 511 514.[Abstract]
Phillips, D. P. (1974) The influence of suggestion on suicide: substantive and theoretical implications of the Werther effect. American Sociological Review, 39, 340 354.
Preti, A. (2000) Seasonal variation and meteotropism in suicide. Current Opinion in Psychiatry, 13, 655 660.[CrossRef]
Rissmiller, D. J. & Rissmiller, F. (1990) Inpatient suicide epidemics and suggestions for prevention. Hospital and Community Psychiatry, 41, 922 924.[Medline]
Stack, S. (2000) Media impacts on suicide: a quantitative review of 293 findings. Social Science Quarterly, 81, 957 971.
Stack, S. (2003) Media coverage as a risk
factor in suicide. Journal of Epidemiology and Community
Health, 57, 238
240.
Taiminen, T. J. & Helenius, H. (1994) Suicide clustering in a psychiatric hospital with a history of a suicide epidemic: a quantitative study. American Journal of Psychiatry, 151, 1087 1088.[Abstract]
Taiminen, T. J., Salmenpera, T. & Lehtinen, K. (1992) A suicide epidemic in a psychiatric hospital. Suicide and Life-Threatening Behavior, 22, 350 363.[Medline]
Zemishlany, Z., Weinberger, A., Ben-Bassat, M., et al (1987) An epidemic of suicide attempts by burning in a psychiatric hospital. British Journal of Psychiatry, 150, 704 706.[Abstract]
Received for publication March 25, 2004. Revision received February 17, 2005. Accepted for publication February 19, 2005.
Related articles in BJP:
HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Psychiatric Bulletin | Advances in Psychiatric Treatment | All RCPsych Journals |