The role of depression and use of alcohol and other drugs after partner suicide in the association between suicide bereavement and suicide: cohort study in the Danish population

Background Although suicide bereavement is associated with suicide and self-harm, evidence regarding mechanisms is lacking. We investigated whether depression and substance use (alcohol and/or other drugs) explain the association between partner suicide bereavement and suicide. Methods Linkage of nationwide, longitudinal data from Denmark for the period 1980–2016 facilitated a comparison of 22 668 individuals exposed to bereavement by a partner's suicide with 913 402 individuals bereaved by a partner's death due to other causes. Using causal mediation models, we estimated the degree to which depression and substance use (considered separately) mediated the association between suicide bereavement and suicide. Results Suicide-bereaved partners were found to have a higher risk of suicide (HRadj = 1.59, 95% CI 1.36–1.86) and of depression (ORadj 1.16, 95% CI 1.09–1.25) when compared to other-bereaved partners, but a lower risk of substance use (ORadj 0.83; 95% CI 0.78–0.88). An increased risk of suicide was found among any bereaved individuals with a depression diagnosis recorded post-bereavement (ORadj 3.92, 95% CI 3.55–4.34). Mediation analysis revealed that depression mediated 2% (1.68%; 95% CI 0.23%–3.14%; p = 0.024) of the association between suicide bereavement and suicide in partners when using bereaved controls. Conclusions Depression is a partial mediator of the association between suicide bereavement and suicide. Efforts to prevent and optimize the treatment of depression in suicide-bereaved people could reduce their suicide risk. Our findings might be conservative because we did not include cases of depression diagnosed in primary care. Further work is needed to understand this and other mediators.


Introduction
Consistent evidence demonstrates an association between suicide bereavement and subsequent risk of suicide (Pitman, Osborn, King, & Erlangsen, 2014, 2022).Apart from preventing exposure to suicide loss in significant others, intervention hinges upon identification of modifiable risk factors.However, we still lack insights into potential explanatory factors between suicide bereavement and suicide.These include psychiatric disorders (particularly depression), alcohol or other drug use, stigma, grief, loneliness, shared social and economic adversity, suicide suggestion, and genetic factors (inherited or due to homophily) (O'Connor & Nock, 2014;Pitman et al., 2014).It is important to understand these links to inform the content of postsuicide emotional support (postvention); an established policy adopted widely within suicide prevention strategies (Schlichthorst et al., 2022;World Health Organization, 2018).However, while evidence supports the effectiveness of postvention in reducing depression and anxiety there is no evidence it helps reduce risks of suicide (Andriessen et al., 2019;Linde, Treml, Steinig, Nagl, & Kersting, 2017;McDaid, Trowman, Golder, Hawton, & Sowden, 2008;Szumilas & Kutcher, 2011).
To date, the mechanisms underlying the elevated risk of suicide after suicide bereavement have not been investigated longitudinally to understand temporality.The two strongest candidate mediators are depression and substance (alcohol and/or other drugs) use, which are both prevalent after a suicide loss (Bolton et al., 2013;Erlangsen et al., 2017;Pitman et al., 2014;Spiwak et al., 2020) and each is associated with suicide (Fazel & Runeson, 2020).Qualitative studies in Australia, the UK, and the US provide accounts of depressive symptomatology after suicide bereavement (Entilli, Ross, de Leo, Cipolletta, & Kõlves, 2021;Pitman et al., 2018) linked to suicidality (Hunt, Young, & Hertlein, 2019), while longitudinal studies of suicide-bereaved partners and children confirm an increased probability of mood disorders (Appel et al., 2013;Erlangsen et al., 2017;Kuramoto et al., 2010;Pitman et al., 2014;Spiwak et al., 2020;Wilcox et al., 2010) and antidepressant treatment (Appel et al., 2016) after suicide loss.Qualitative research in the US and UK also documents excessive use of alcohol and other drugs after being bereaved by suicide (Eng et al., 2019;Hunt et al., 2019;Pitman, Stevenson, King, & Osborn, 2020) and other sudden unnatural causes (Drabwell et al., 2020;Pitman et al., 2020) as a means of coping with overwhelming thoughts and emotions.However, some individuals described reducing their use of alcohol or other drugs to better cope with emotions or to avoid substances perceived as contributing to the death of their relative or friend (Drabwell et al., 2020;Eng et al., 2019;Pitman et al., 2020).Longitudinal studies of suicide-bereaved offspring report an increased probability of alcohol or other drug use after suicide bereavement (Brent, Melhem, Donohoe, & Walker, 2009;Wilcox et al., 2010), possibly transient (Hamdan, Melhem, Porta, Song, & Brent, 2013), but longitudinal studies of suicidebereaved partners find a similar prevalence of alcohol or other drug use disorders to partners bereaved by other causes (Bolton et al., 2013;Erlangsen et al., 2017;Spiwak et al., 2020).Such inconsistencies may be explained by differences in measurement and samples.It is also possible that mechanisms differ by kinship group (Pitman et al., 2022).
Given the evidence supporting depression and alcohol and/or other drug use as candidate mediators, large longitudinal studies are required to investigate this, ideally by comparing suicide bereavement to other bereavements to account for the experience of bereavement per se.New methods of counterfactual mediation analysis (MacKinnon, Fairchild, & Fritz, 2007) are underused in suicide research.These offer a more robust means of assessing the relative magnitude of different pathways and mechanisms by which an exposure may affect an outcome because it makes the causal assumptions explicit, adjusts for confounders of each arm of mediational models, and takes into account interactions with exposure (VanderWeele, 2016a).The aim of this study was to estimate the relative proportions of the association between suicide bereavement and suicide mediated by depression and by substance use (alcohol and/or other drugs) when comparing suicide-bereaved partners to other-bereaved partners.

Study design and participants
We conducted a population-based cohort study, analyzing data from administrative national registers linked using a unique personal identification number assigned to all individuals in Denmark (Erlangsen & Fedyszyn, 2015).Population data from the Civil Registration System was linked with information on psychiatric and somatic hospital contacts from the Psychiatric Central Research Register (since 1970) (Mors, Perto, & Mortensen, 2011) and the National Patient Register (since 1977) (Schmidt et al., 2015) but linkage was not available to primary care data.
We included all Danish-born individuals aged 16 years and older who were living in Denmark between 1 January 1980 and 31 December 2016, comparing participants bereaved by the (a) suicide or (b) non-suicide death of a current/former partner (spouses, civil partners, and cohabitees) over that period.Cohort inception in 1980 provided us with over 10 years of previous data on psychiatric confounders as identified in the Psychiatric Central Research Register.
We focused on studying mechanisms in one kinship group (partners) because they offered the largest sample relative to parents, offspring, and siblings (Pitman et al., 2022).Exposure to partner bereavement was identified by linking data on individuals in the cohort to deceased current and former spouses/cohabitees, identifying causes of death from the Register of Causes of Death.Partners (opposite and same-sex married and registered partners and opposite sex cohabitees) were identified using an established linkage method based on data recorded in the Civil Registration System (Pedersen, 2011) (online Supplemental Methods S1).Individuals entered the cohort on the date of the index (first) bereavement and were followed up until the date of their suicide, death from a cause other than suicide, migration out of Denmark, second partner bereavement, or the end of follow-up (31 December 2016).

Outcome
Our outcome, suicide, was identified from the Register of Causes of Death (Helweg-Larsen, 2011) based on International Classification of Diseases (ICD)-8 and ICD-10 codes (online Supplemental Table S1).

Exposure
Relevant ICD-8 and ICD-10 codes were used to identify whether bereavement was due to suicide or other causes.The first partner bereavement over this period was identified as the index exposure, and individuals were censored at any subsequent partner bereavement.

Confounders
We chose eight confounders a priori, based on existing evidence (Pitman et al., 2014), which we used in each arm of the mediation model: sex; age; bereavement year; marital status (as distinct from cohabitation status, to capture the effect of divorce); household income level (quartiles); pre-bereavement history of any admission recording self-harm (in psychiatric or physical health settings); prebereavement psychiatric disorders recorded on any admission (including alcohol and/or other drug use); and pre-bereavement physical health conditions recorded on medical admission (online Supplemental Table S1).Bereavement year captured period effects and differing inception years for inpatient and outpatient data.All confounders were measured before the index bereavement.

Mediators
Depression was defined as a diagnosis recorded during psychiatric inpatient (since 1980) or outpatient (since 1995) contact.Substance use was defined broadly as an ICD-8/ICD-10 diagnosis of alcohol and/or other drug use disorders recorded during psychiatric inpatient (since 1980) or outpatient (since 1995) contact, including acute intoxication, as well as inpatient medical admission recording specific medical problems linked to alcohol use (Grissa, Rasmussen, Krag, Brunak, & Jensen, 2020).Both mediators were restricted to events recorded after the index bereavement.

Statistical analysis
The demographic and clinical characteristics of the sample were examined using complete case analysis.We used Cox proportional hazards regression to examine the associations between suicide bereavement and time to suicide in unadjusted and adjusted models.We assessed model assumptions including proportionality of hazards (online Supplemental Methods S2).In preliminary analyses we accounted for potential clustering effects, where multiple current/former partners might be bereaved by the same death, using multilevel Cox proportional hazards regression clustered on the individual who died.As the multilevel model provided comparable estimates to the single-level model, we opted for the latter.
Mediation by depression and substance use was assessed separately using the causal inference potential outcomes framework to decompose the total effect (TE) of suicide bereavement on suicide into four components: (a) controlled direct effect (CDE), (b) reference interaction (INT ref ), (c) mediated interaction (INT med ); and (d) pure indirect effect (PIE), and other derived values, including the total indirect effect (TIE), portion attributable to interaction (PAI), and the portion eliminated (PE) (Richiardi, Bellocco, & Zugna, 2013;Vanderweele, 2014;VanderWeele, 2016b) (online Supplemental Methods S4; Supplemental Box S1).From this output, we identified: the overall proportion of the association mediated (equivalent to 100 × [TIE/TE]); the proportion of the association attributable to interaction (effectively 100 × [PAI/TE]); and the proportion of the total effect of the exposure on the outcome that would be eliminated in the absence of the mediator (effectively 100 × PE/TE).The last of these is of policy relevance in conveying how much of the effect of the exposure can be prevented by intervening on the mediator.
We fitted the mediation model using the med4way command in Stata (Discacciati, Bellavia, Lee, Mazumdar, & Valeri, 2019;StataCorp, 2019).The outcome model (path C in Fig. 1) was fitted using Cox proportional hazards regression, while the mediation model (paths A and B) was fitted using logistic regression because med4way does not accommodate a time-to-event mediator and time-to-event outcome.
In our first sensitivity analysis we assessed the influence of competing risks (death from other causes; second bereavement; emigration) in our main association (Supplemental Methods 3).We then conducted an interaction test to investigate whether partner status (ex-versus current) modified this main association.In our second sensitivity analysis we repeated our mediation analyses restricted to records of inpatient admissions to ascertain the influence of having incorporated outpatient contact data within our main analyses.In our third sensitivity analysis we repeated our mediation analyses excluding those with widowed marital status to ascertain the influence of having included those with additional partner bereavements prior to our period of interest.
We repeated our mediation analyses restricted to records of inpatient admissions to ascertain the influence of having incorporated outpatient contact data within our main analyses.We also repeated our mediation analyses excluding those with widowed marital status to ascertain the influence of having included those with additional partner bereavements prior to our period of interest.
Data management was performed using SAS software 9.4 (SAS System for SunOP, 2003) and analyses were performed using Stata 17 software (StataCorp, 2021).

Sample characteristics
A total of 960 272 bereaved partners were identified of whom 24 202 (2.5%) were excluded due to missing data on income level (2.5%) and/or marital status (0.01%) (online Supplemental Table S2).Among the remaining 936 070 (97.5%) partners, 22 668 (2.4%) had been bereaved by suicide and 913 402 (97.6%) by other causes of death (Table 1).The majority of the analytic sample was female (67.3%).Individuals bereaved by a partner's suicide were more likely to be younger, female, to have a higher household income level and a past history of self-harm.

Main association
During almost 10 million person-years of follow-up, 3733 individuals died by suicide, of whom 182 (4.9%) had been bereaved by partner's suicide and 3551 (95.1%) bereaved by other causes.Individuals bereaved by a partner's suicide had a higher risk of subsequent suicide compared to those bereaved by other causes (HR adj : 1.59, 95% CI 1.36-1.86)(Table 2).The absolute risk of suicide after partner bereavement was 0.80% in individuals bereaved by a partner's suicide and 0.39% in individuals bereaved by other causes.
Suicide bereavement was associated with an increased risk of substance use in the unadjusted model (OR crude 2.12; 95% CI 2.01-2.25;Table 3) but this attenuated and changed direction when adjusted (OR adj 0.83; 95% CI 0.78-0.88).Post hoc adjustments (online Supplemental Table S3) showed that the attenuation was primarily due to the confounding effect of age, as the suicidebereaved group was younger than controls, and older people had a lower probability of substance use.Nevertheless, post-bereavement substance use was associated with an increased risk of suicide (OR crude 3.61; 95% CI 3.25-4.02;OR adj 1.59, 95% CI 1.41-1.79).
In a separate model for substance use, we found no evidence to support a mediating role in the association between suicide bereavement and suicide (Table 5).

Sensitivity analyses
When accounting for competing risks in our main association (online Supplemental Table S4), similar results were obtained to Psychological Medicine those from the main analysis (crude sub-distribution hazard ratio [sHR]: 1.99, 95% CI 1.71-2.31;sHR adj : 1.54, 95% CI 1.32-1.80;online Supplemental Table S5).No interaction was found with respect to current/ex-partner status ( p = 0.130; online Supplemental Table S6).When we restricted depression to inpatient records only for our main association, suicide bereavement was not associated with an increased risk of depression or substance use (online Supplemental Table S7), therefore we did not conduct mediation analyses using either measure.When we ran our mediation model for depression excluding those with widowed marital status, the proportion mediated rose from 1.68% to 1.74%, but with overlapping confidence intervals.

Main findings
Based on complete, nationwide data on all bereaved partners, we found evidence of an elevated risk of depression and of suicide after partner suicide bereavement compared with bereaved controls, but a reduced risk of substance use.Suicide-bereaved Dashed lines reflect potential confounding pathways.Path A represents the association between suicide bereavement and post-bereavement depression/substance use (modeled using logistic regression).Path B represents the association between post-bereavement depression/substance use (in all those bereaved) and suicide modeled using logistic regression.Path C represents the association between suicide bereavement and suicide (modeled using Cox proportional hazards regression).The same set of confounders was used for all pathways: sex, age, bereavement year, legal marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders, and any physical disorders.2017), as could the perception that the system failed the deceased, resulting in under-ascertainment of psychiatric disorder and substance use when compared to controls.Finally, interventions offered during outpatient and inpatient care may also have mitigated suicide risk.For these reasons, it is plausible that depression mediates a greater proportion of the association between suicide bereavement and suicide than our study using secondary care data suggests.Further work is needed to investigate mediators of this association using both primary and secondary care data, as well self-report variables (where available) to capture perceived stigma and reluctance to seek help.

Findings in the context of other studies
This is, to our knowledge, the first formal evaluation of potential mediators of suicide risk after suicide loss.Our findings extend existing evidence regarding the elevated risk of depression (Erlangsen et al., 2017) and suicide (Agerbo, 2003(Agerbo, , 2005;;Erlangsen et al., 2017) among suicide-bereaved partners compared with bereaved controls.In contrast, our novel finding of a reduced risk of substance use in the suicide-bereaved compared with bereaved controls (primarily due to the confounding effect of age) differs from previous findings of no differences (Bolton et al., 2013;Erlangsen et al., 2017).It is, however, consistent with the reduced risk of liver cirrhosis documented among suicide-bereaved spouses compared with bereaved controls (Erlangsen et al., 2017) and the accounts of people bereaved by suicide and other causes reporting a post-bereavement reduction in substance use (Drabwell et al., 2020;Eng et al., 2019;Pitman et al., 2020).This may be linked to the enhanced awareness described by suicide-bereaved individuals of their shared vulnerabilities to suicide, engendering a determination to safeguard their mental health after the loss (Pitman et al., 2017).We also noted that Danish partners bereaved by suicide had a significantly higher household income than partners bereaved by non-suicide causes, in keeping with a similar finding for bereaved partners in Canada (Spiwak et al., 2020) but in contrast with the converse finding for bereaved parents in Canada (Bolton et al., 2013).

Strengths and limitations
Strengths include analyzing longitudinal data on a large population-based cohort of almost a million people, followed for approximately 10 million person-years, with minimal loss to follow-up, low levels of missing data, and avoidance of selection biases or recall biases.The Danish registers have been evaluated as reliable with respect to psychiatric hospital contacts (Tøllefsen et al., 2015) and registration of suicide death (Helweg-Larsen, 2011;Tøllefsen et al., 2015) and unnatural deaths (Tøllefsen et al., 2015), with good validity of diagnostic codes and the timing of admission/discharge dates in the Danish National Patient Register (Schmidt et al., 2015) and of depression diagnoses and admission/outpatient dates in the Psychiatric Central Research Register (Mors et al., 2011).For our study, this was important in separating out mediators and confounders for each pathway.In utilizing the Danish Civil Registration household variable, we were able to identify cohabiting couples as well as those in legal unions (although not opposite-sex cohabiting couples), providing a more realistic representation of Danish family structure and a more comprehensive measure of exposure to partner suicide loss.The limitations of using routine Danish registry data include the under-ascertainment of mental disorder and self- b Household income quartiles represent total income within the household divided by the total number of adults living in the household, then categorized into quartiles based on national annual income averages.c The widowed category represented people who were bereaved by subsequent partner loss before the population registers had started, yet were bereaved while in a new partnership, so by default were registered as widowed by a former partner's death.d Data for these variables were from diagnoses (or self-harm) recorded on inpatient admissions (i.e.excluding outpatient data).
e Data for these variables were from diagnoses recorded on inpatient admissions and outpatient contacts.All p-values were <0.001 apart from depression (<0.01) and PTSD (non-significant) but these values are to be interpreted in the context of a large sample size, and therefore a high probability of statistically significant differences.
Table 2. Hazard ratios and 95% confidence intervals for the association between suicide bereavement and suicide compared with other bereavement (Fig. 1  odds ratio; CI, confidence interval.a Adjusted for sex, age, bereavement year, marital status, and household income level.b Final model: adjusted for all variables in adjustment 1, plus pre-bereavement history of self-harm, any psychiatric disorders (one binary variable), and any physical disorders (one binary variable).c N is solely for the outcome in each path; thus of those who are suicide-bereaved, 938 (4.1%) have depression and 1407 (6.2%) have substance use, and of those who are other-bereaved, 26744 (2.9%) have depression and 27597 (3.0%) have substance use.Of those who have depression, 510 (1.8%) die by suicide, and of those who do not have depression, 3223 (0.4%) die by suicide.Similarly, of those who have substance use, 383 (1.3%) die by suicide, and of those who do not have substance use, 3350 (0.4%) die by suicide.* p < 0.001.Psychological Medicine harm (due to lack of secondary care outpatient data pre-1995, primary care data, and data from people who do not seek or receive treatment) and not capturing other influences on suicide risk, such as personality traits, coping style, or social support.Our application of contemporary causal inference methods overcame the limitations of older methods of mediation by allowing for the adjustment of confounders in each arm of the mediation model, making causal assumptions explicit, allowing us to model interactions, and was particularly suitable for binary outcomes (Rijnhart, Valente, MacKinnon, Twisk, & Heymans, 2021).However, mediation assumes the absence of residual confounding between exposure-outcome, exposure-mediator, or mediator-outcome, and we acknowledge the possibility of residual confounding due to unmeasured pre-bereavement depression or self-harm (not captured in settings beyond psychiatric admission).Our mediation models could not assess the role of the other respective putative mediator as an intermediate confounder because med4way can only consider one mediator at a time, for ease of interpretability.Our use of this model also relied on a set of theoretical assumptions, and it is possible that not all were satisfied, introducing the potential for biased estimates of causal effects.
We also acknowledge the potential for collider bias where selecting bereaved individuals for comparison might result in distorted associations between variables (Holmberg & Andersen, 2022).Given the difficulties in aligning follow-up for a cohort comprising individuals bereaved in childhood, early adulthood, mid-life, and late life, and because mechanisms may differ by kinship, we opted only to assess bereaved partners.This diminished statistical power to detect mediation effects.As a putative mediator, our broad substance use definition (including medical complications) relied less on a clinician recording an alcohol and/or other drug use disorder diagnosis e during admission or on treatment-seeking for perceived substance use problems.However, our narrow definition of depression omitted untreated depression, pre-1995 cases of depression in outpatients, cases identified in primary care (as we lacked primary care) or during somatic admissions, or information from prescribing data (available from 1995).Ascertainment bias is also possible if suicide-bereaved individuals are more likely to be in contact with mental health professionals and therefore more likely to be admitted.

Clinical, policy, and research implications
The increased risk of depression after suicide loss and its contribution to suicide risk suggest a benefit of early identification of depression in suicide-bereaved partners (both ex-and current).Suicide risk mitigation might be achieved through appropriate bereavement support (Andriessen et al., 2019;Linde et al., 2017;McDaid et al., 2008;Szumilas & Kutcher, 2011), improving uptake through collaboration between bereavement services and primary care.Future research should use primary care and prescribing data to capture psychiatric disorders treated in primary and secondary care, investigating depression, substance use, and other putative mediators to advance our understanding of therapeutic targets.This will inform recommendations regarding appropriate support after suicide loss.The inclusion of postvention within a national suicide prevention strategy is not associated with changes in suicide mortality (Schlichthorst et al., 2022).However, a better understanding of mediators has the potential to improve the effectiveness of postvention and reduce suicides.

Conclusions
Suicide risk in suicide-bereaved partners was almost twice that for partners bereaved by other causes.Depression mediated a small proportion (2%) of this association, but this study likely underestimated the mediating role of depression by omitting depression treated in primary care.Risk of substance use was reduced after partner suicide loss, and substance use did not mediate the association between suicide bereavement and suicide.These findings, from a large, longitudinal representative dataset, identified early treatment of depression as a potential means of reducing the burden of suicide after suicide bereavement.
Supplementary material.The supplementary material for this article can be found at https://doi.org/10.1017/S0033291724000448 Funding statement.This work was supported by a grant from the American Foundation for Suicide Prevention (SRG-0-111-17).AP, Glyn Lewis and Gemma Lewis are also supported by the National Institute for Health Research (NIHR) University College London Hospital (UCLH) Biomedical Research Centre (BRC).Gemma Lewis is supported by a Sir Henry Dale Fellowship jointly funded by the Wellcome Trust and the Royal Society (grant number 223248/Z/21/Z).The funders were not involved in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
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.The study was approved by the Danish Data Protection Agency (Capital Region of Denmark; P-2020-774) and the UCL (University College London) Research Ethics Committee (approval number 14 075/001).Informed consent was waived.

Figure 1 .
Figure 1.Hypothesized model for the mediating effect of (A) depression and (B) substance use on the association between suicide bereavement and suicide.Legend: Panel A displays the hypothetical model for mediation by depression.Panel B displays the hypothetical model for mediation by substance use.Solid lines represent potential causal pathways.Dashed lines reflect potential confounding pathways.Path A represents the association between suicide bereavement and post-bereavement depression/substance use (modeled using logistic regression).Path B represents the association between post-bereavement depression/substance use (in all those bereaved) and suicide modeled using logistic regression.Path C represents the association between suicide bereavement and suicide (modeled using Cox proportional hazards regression).The same set of confounders was used for all pathways: sex, age, bereavement year, legal marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders, and any physical disorders.

Table 1 .
(Continued.) a Data are n (%), except age, which is summarized as median (IQR).All values were pre-bereavement unless otherwise specified.

Table 3 .
Odds ratios and 95% confidence intervals for each pathway (modeled separately) in the model for mediation of the association between suicide bereavement and suicide by depression (Fig.1, Panel A) and substance use(Fig. 1, Panel B)

Table 4 .
Mediation analyses describing the role of depression in the association between suicide bereavement and suicide Adjusted for sex, age, bereavement year, marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders (one binary variable), and any physical disorders (one binary variable).Note that in a mediation analysis, only adjusted estimates are presented because mediation can only be interpreted in the context of adjusting for all measured confounders (see online Supplemental Methods S4).
a b Where no stata label is given, this value was computed from other stata output values.

Table 5 .
Mediation analyses describing the role of substance use in the association between suicide bereavement and suicide Adjusted for sex, age, bereavement year, marital status, household income level, pre-bereavement history of self-harm, any psychiatric disorders (one binary variable), and any physical disorders (one binary variable).Note that in a mediation analysis, only adjusted estimates are presented because mediation can only be interpreted in the context of adjusting for all measured confounders (see online Supplemental Methods S4).b Where no stata label is given, this value was computed from other stata output values.
RER, relative excess risk; CI, confidence interval; see online Supplemental Box S1 for definitions.a 2280 Alexandra Pitman et al.