Corrigendum: A Theory-Based Longitudinal Investigation Examining Predictors of Self-Harm in Adolescents With and Without Bereavement Experiences

[This corrects the article DOI: 10.3389/fpsyg.2020.01153.].

In the original article, there was a mistake in the data matching and resulting statistics/text as published. A data processing error occurred whilst combining participant data across both time points, which led to a proportion of the baseline data being incorrectly matched to the follow-up data. Upon rectifying this error and re-running the analyses, we found that the main results and conclusions remain unaffected, such that bereavement experiences (by suicide or non-suicide causes) continue to show no association with self-harm group membership at baseline nor at follow-up. We did, however, find that the specific descriptive and inferential statistics have changed slightly throughout the paper as a result of the corrected matching. Our cross-sectional findings remain largely unaffected, although one variable (maladaptive coping) which was previously approaching significance has now reached statistical significance. Furthermore, two additional variables have emerged as significant in our longitudinal analyses (endorsing stigmatising beliefs, and glorifying/normalising beliefs about suicide predict prospective self-harm), and one variable (family self-harm) is no longer a significant predictor of prospective self-harm. The overall conclusions of the paper are not affected, and the additional predictors which have emerged do not have a major influence on the interpretation of our findings. The corrected paragraphs and tables appear below, as well as additional comments in the Results and Discussion section and the insertion of Table 5b, which address the changes in the results noted above.
Corrections have been made to Abstract, Methods, and Results sections, as shown below: [ABSTRACT] Methods: A 6-month prospective questionnaire study was conducted with 185 Scottish adolescents aged 11-17 (115 adolescents also completed the questionnaire at follow-up). The questionnaire included measures on experiences with bereavement and lifetime engagement in self-harm, as well as measures of defeat, entrapment, social support, coping, and other psychological variables.
Results: At baseline, 11% of young people reported exposure to a suicide death, and 62% to a non-suicide death. In addition, 21% of pupils reported ever engaging in self-harm, while 24% had experienced self-harm ideation without engaging in it. Cross-sectional multivariate logistic regressions showed that maladaptive coping, family social support, glorifying/normalising beliefs about suicide, and family selfharm were significantly associated with self-harm group membership (control, ideation, or enactment groups). At followup, 9% of pupils reported exposure to a suicide death and 11% to a non-suicide death for the first time. A total of 29% of the sample reported self-harm at T2 (8% of participants for the first time), and 23% reported self-harm ideation without engaging in it. Multivariate analyses found that stigmatising beliefs about suicide, glorifying/normalising beliefs about suicide, and self-harm ideation at baseline were the only variables to predict self-harm group membership prospectively. Bereavement experiences, whether by suicide or non-suicide, did not predict self-harm group status at baseline nor at follow-up.
Corrections have also been made to Materials and Methods, Participants, Paragraphs 1-3. The corrected text is shown below: [

MATERIALS AND METHODS] Participants
A total of 185 pupils (aged 11-17, M = 13.16, SD = 1.49) were recruited at T1 from nine secondary schools across Scotland. This sample consisted of individuals retained after removing participants with >50% missing data (n = 2) or who did not provide data on any of the SSHTB outcome measures (n = 22), and including participants from T2 who only provided data once and not at baseline (i.e., were absent at the first time point, or baseline data was removed due to missingness but T2 questionnaire was complete; n = 15). Of the T1 sample, 85 stated they were male, 97 female, 2 other, and 1 did not respond. Approximately half were in Year 1 (predominantly aged 12-13) of secondary school (n = 91, 49.2%) and described their ethnicity as White (n = 164, 88.65%), consistent with the last Scottish Census (96.1%; National Records of Scotland, 2011). The percentage of pupils entitled to free school meals, as a proxy measure of Socioeconomic Status (SES), ranged from 4.74 to 20.99% between schools (M = 14.29, SD = 5.19), slightly lower than previous Scottish studies (e.g., mean of 17.8% in Russell et al., 2018), though comparable to the national average of 14.4% (Scottish Government, 2018).
One hundred and fifteen individuals (aged 12-18, M = 13.65, SD = 1.52; 46 male, 67 female, 2 other) provided data for T2, which could be matched to corresponding baseline data. This sample was retained after removing data from respondents whose T2 participant identifier codes could not be confidently matched to their baseline data (n = 31), who had >50% missing data (n = 3), or who did not respond to any of the outcome measures (n = 5).
Corrections have also been made to Measures, Motivational Phase Variables, sections Entrapment, Coping, Self-Esteem, Social Support, and Stigma. The corrected sections are shown below: Entrapment Entrapment was measured using the Entrapment Scale (Gilbert and Allan, 1998), which evaluates perceptions of being unable to escape from one's current situation or circumstances. Sixteen self-report items, e.g., "I have a strong desire to escape from things in my life, " are rated on a 5-point scale from 0 (never) to 4 (always), reflecting how frequently they have been experienced. The scale showed high internal consistency at α = 0.94. Coping Coping was assessed using the Brief COPE (Carver, 1997), measuring the degree to which a person uses a specific strategy to deal with difficult or stressful situations. The 28-item measure covers various strategies, e.g., "I have been using alcohol or other drugs to make myself feel better, " which are evaluated using a 4-point scale from 1 (I haven't been doing this at all) to 4 (I've been doing this a lot). Several scoring methods have been proposed. As per Moore et al. (2011) andBlomgren et al. (2016), we differentiated between adaptive coping (16 items covering active coping, planning, positive reframing, humour, acceptance, religion, use of emotional support, and use of instrumental support) and maladaptive coping (12 items on self-distraction, denial, substance use, behavioural disengagement, venting, and self-blame). Cronbach's α was high for both the adaptive and maladaptive subscales, at 0.83 and 0.74, respectively.
Self-esteem The Rosenberg Self-Esteem Scale (Rosenberg, 1965) was used to measure self-esteem by asking about self-worth and positive or negative feelings about oneself. Ten items are answered on a 4-point scale from 0 (Strongly Agree) to 3 (Strongly Disagree), with higher scores indicating greater self-esteem, e.g., "I take a positive attitude toward myself." Internal consistency was high at α = 0.91.

Social support
The Multidimensional Scale of Perceived Social Support (MSPSS; Zimet et al., 1988) was used to assess the perceived adequacy of social support that an individual receives from family, friends, and significant others. Each of the three categories is assessed with four items, given on a 7-point scale ranging from 1 (very strongly disagree) to 7 (very strongly agree), e.g., "I can talk about my problems with my family." A total overall score or three subscale scores for the different sources of support can be calculated by summing the relevant items. Subscales were used here to differentiate the influence of different sources of support; internal consistency was high (family α = 0.88, friends α = 0.90, and significant others α = 0.88).

Stigma
Attitudes toward people who die by suicide were measured using the Short Form of the Stigma of Suicide Scale (SOSS; Batterham et al., 2013b), which asks participants to rate how much they agree or disagree with words describing people who take their own lives. Sixteen items, e.g., "irresponsible, " "lonely, " "noble, " are rated on a 5-point scale from Strongly Disagree to Strongly Agree. Subscales of stigma, isolation/depression, and glorification/normalisation can be calculated by summing the relevant items for each subscale; internal consistency for the respective subscales was α = 0.83, α = 0.85, and α = 0.75.
Corrections have also been made to Materials and Methods, Data Analytic Plan, Paragraph 1. The corrected paragraph is shown below:

Data Analytic Plan
Missing data was dealt with using multiple imputation, as Little's MCAR test was non-significant, χ 2 (5) = 5.09, p = 0.405, and data was deemed to be most likely missing completely at random. A total of m = 68 imputations were generated based on 68% of cases having incomplete data (as suggested by White et al., 2010). Analyses were conducted using SPSS Version 27, which supports pooled analyses based on imputed datasets for several statistical tests; however, some analyses are not supported by this function. In such cases, parameter estimates were manually averaged across the 68 imputed datasets, an approach also taken by Jones et al. (2014) when dealing with imputed data in SPSS. Microsoft Excel 2013 was used to manually pool parameter estimates where necessary.
Additionally, corrections have been made to the Results section.
Corrections have been made to the Results, Prevalence of Bereavement and Self-Harm at Baseline (T1) and Prevalence of Bereavement and Self-Harm at Follow-Up (T2) sections. The corrected paragraphs are shown below: [RESULTS] Prevalence of Bereavement and Self-Harm at Baseline (T1) A comparison of those who took part at baseline only and those who participated at both time points revealed no significant differences on most of the demographic or studied variables, apart from SES, family social support, and lifetime self-harm (those followed-up came from schools with a lower proportion of pupils receiving free school meals, reported lower family social support, and were more likely to report self-harm at baseline). Descriptive statistics of continuous study variables for all participants across all self-harm groups are shown in Table 1. At baseline (n = 185), 136 (73.51%) young people reported that someone among their immediate family and/or someone else close had died; 21 (11.35%) of which knew someone who had died by suicide (making up the suicide exposed group), while the remaining 115 (62.16%) people were exposed to a non-suicide death. 38 (20.54%) pupils reported having ever engaged in selfharm behaviours during their lifetime (enactment group), while a further 44 (23.78%) reported past self-harm ideation with no history of behaviours (ideation group). Thus, the control group at baseline consisted of 103 (55.68%) individuals with no history of self-harm or suicidal thoughts or behaviours.

Prevalence of Bereavement and Self-Harm at Follow-Up (T2)
At follow-up, 81 participants of the T2 sample of n = 115 reported that someone among their immediate family and/or someone else close had died. Seventeen (14.78%) individuals overall reported knowing someone who had died by suicide (suicide exposed group), of which 10 (8.70%) were reported for the first time since T1. A further 66 (57.39%) individuals were exposed to a non-suicide death, with 13 (11.30%) reported for the first time since baseline. It is worth noting that two individuals responded 'no' to the death of an immediate family member or anyone close, but 'yes' to experiencing a suicide death of family or friends.
At follow-up, 33 (28.70%) adolescents reported ever engaging in self-harm, with 9 (7.83%) of these for the first time between Time 1 and Time 2. A further 26 (22.61%) individuals reported having experienced self-harm ideation (with no actions) at follow-up. The control group at T2 therefore comprised of 56 (58.70%) individuals who reported no history of self-harm ideation or behaviours at follow-up.
Corrections have also been made to Results, Cross-Sectional Associations Between Motivational and Volitional Phase Variables and Self-Harm at Baseline (T1), sections Motivational Phase Variables and Volitional Phase Variables. The corrected sections are shown below:

Volitional Phase Variables
A similar logistic regression analysis was conducted to examine volitional phase variables and their association with selfharm group status at baseline. Univariate analyses showed that ideation and enactment groups did not differ from controls on any variable. The ideation group differed from the enactment group only on family self-harm (OR = 0.17, 95% CI = 0.04-0.70, p = 0.014), where those who self-harmed were more likely to report this experience ( Table 4). Neither experiencing a suicide nor a non-suicide death were associated with self-harm group membership. A multivariate analysis was not conducted as only one variable emerged as a significant predictor in this analysis.
Corrections have also been made to Results, Longitudinal Associations Between Motivational and Volitional Phase Variables and Self-Harm at Follow-Up (T2), sections Motivational Phase Variables and Volitional Phase Variables. The corrected sections are shown below: Longitudinal Associations Between Motivational and Volitional Phase Variables and Self-Harm at Follow-Up (T2) Motivational Phase Variables A hierarchical multinomial logistic regression examined whether motivational phase variables were associated with lifetime self-harm group 6-months later. In univariate analyses, participants in the ideation group were significantly more likely than controls to have reported self-harm ideation at baseline (Table 5a; OR = 0.08, 95% CI = 0.02-0.34, p = 0.001). Those in the enactment group differed from controls on family social support (OR = 0.53, 95% CI = 0.33-0.87, p = 0.011), stigmatising beliefs about suicide (OR = 0.27, 95% CI = 0.09-0.78, p = 0.016), glorifying/normalising beliefs about suicide (OR = 0.31, 95% CI = 0.12-0.76, p = 0.010), and self-harm ideation at baseline (OR = 0.05, 95% CI = 0.01-0.26, p = 0.001). The enactment group did not differ from the ideation group on any motivational phase variable. Significant univariate predictors associated with self-harm were entered into a multivariate analysis (Table 5b), which found that three factors continued to be associated with self-harm group membership: the ideation group were less likely to report glorifying/normalising beliefs about suicide (OR = 0.19, 95% CI = 0.06-0.66, p = 0.009), and more likely to report selfharm ideation at baseline (OR = 0.03, 95% CI = 0.00-0.20, p < 0.001) compared to controls. In addition, the enactment group were less likely to hold stigmatising beliefs about suicide (OR = 0.18, 95% CI = 0.05-0.73, p = 0.016), less likely to hold glorifying/normalising beliefs about suicide (OR = 0.11, 95% CI = 0.03-0.46, p = 0.002), and more likely to have reported selfharm ideation at baseline (OR = 0.02, 95% CI = 0.00-0.22, p = 0.002), compared to controls.

Volitional Phase Variables
Another analysis was conducted to examine volitional phase variables and their association with self-harm group status prospectively. In univariate tests, none of the variables, including experiencing a suicide or a non-suicide death, emerged as significant predictors of self-harm group membership ( Table 6). A multivariate analysis was therefore not required.
Overall, cross-sectional analyses showed that maladaptive coping, family social support and endorsing glorifying/normalising beliefs about suicide (motivational phase variables) and family self-harm (volitional phase variable) were significant predictors of self-harm group status. Longitudinally, endorsing stigmatising beliefs about suicide, endorsing glorifying/normalising beliefs about suicide, and self-harm ideation at baseline (motivational phase variables) predicted self-harm group at follow-up.
Lastly, corrections were also made to the Discussion section.
Corrections have been made to Discussion, Paragraph 2. The corrected paragraph is shown below:

[DISCUSSION]
Results partially supported the hypotheses both crosssectionally and longitudinally. Although several variables (defeat, entrapment, maladaptive coping, family social support, and endorsing glorifying/normalising beliefs about suicide) predicted self-harm group membership in univariate analyses, only maladaptive coping, social support from family members and endorsing glorifying/normalising beliefs about suicide remained significant multivariate predictors within the motivational phase of the model. Family self-harm was the only predictor among the volitional phase variables to predict self-harm group cross-sectionally. Results of longitudinal analyses showed that endorsing stigmatising beliefs about suicide, endorsing glorifying/normalising beliefs about suicide, and self-harm ideation (motivational phase test) at baseline predicted selfharm group membership 6-months later, with none of the volitional phase variables emerging as significant predictors of prospective self-harm.
Corrections have also been made to Discussion, IMV Model Psychological Variables as Predictors of Self-Harm, All Paragraphs IMV Model Psychological Variables as Predictors of Self-Harm As predicted, maladaptive coping was found to be associated with self-harm group at baseline, and both ideation and enactment group participants differed from controls but not each other. This is in keeping with the large body of research that associates low levels of coping skills with suicidal thoughts and feelings (Gooding et al., 2015). Research particularly highlights that avoidant or emotion-focused (rather than problem-focused) coping styles are associated with self-harm in adolescents (Guerreiro et al., 2013), including the use of alcohol and drugs, behavioural disengagement and self-blame strategies, which were measured by the maladaptive coping subscale used in this study. The finding that only the maladaptive subscale was significant may reflect that self-harm might represent a coping style in itself; Laye-Gindhu and Schonert-Reichl (2005) suggest that selfharm is an emotion-focused strategy that serves to regulate affect. Indeed, research on motivations for self-harm among adolescents has found that getting relief from a terrible state of mind was the strongest predictor of self-harm (Rasmussen et al., 2016), which may account for the significant association found here.
Our cross-sectional findings reflect previous research (Kleiman and Liu, 2013;O'Connor and Nock, 2014) and theory (O'Connor, 2011;O'Connor and Kirtley, 2018) showing that levels of social support are significantly associated with suicide risk. In a recent large-scale study, Wan et al. (2019) found that lower social support was significantly associated with self-reported NSSI, suicidal ideation and suicide attempts among young people aged 10-20 years old. Our finding that only family social support was associated with self-harm group membership is consistent with Cheng and Chan (2007); using a translated version of the MSPSS, they found that the impact of family social support was stronger than that of friends in predicting suicidality among adolescents. Similarly, Tabaac et al. (2016) reported that social support from family and significant others was associated with suicidal ideation, but only family social support was associated with suicide attempts. They suggest that family members may represent a closer and more permanent source of support than other social groups, particularly for adolescents dealing with stressful life events.
As predicted, endorsing stigmatising beliefs about suicide significantly differentiated controls from enactment groups in longitudinal analyses, and glorifying/normalising beliefs differentiated controls from both ideation and enactment groups prospectively. It was found that higher levels of such beliefs were associated with being in the control group. The finding of a significant association at baseline in the motivational phase variable of glorifying/normalising beliefs about suicide was contrary to IMV model predictions, as ideation and enactment groups were not expected to differ. The ideation group were more likely to endorse glorifying or normalising beliefs about suicide than the enactment group at baseline. Previous research using the same SOSS measure (Batterham et al., 2013a) also showed that suicidal ideation was associated with greater glorification of suicide, as well as less stigma toward suicide, whereas suicide attempts were not associated with any attitude subscale (stigma, isolation/depression, or glorification/normalisation). One possible explanation for both cross-sectional and longitudinal findings is that individuals who self-harm are more likely to have been exposed to similar behaviours in others (Dhingra et al., 2015;Mars et al., 2019b), and increased exposure has been shown to reduce stigma (e.g., in relation to mental disorders; Jorm and Wright, 2008); in this study, experiencing self-harm of family members was indeed associated with self-harm group status crosssectionally, which may account for the lack of an association with glorifying/normalising beliefs among the enactment group. Interestingly, self-harm group status was not associated with suicide bereavement. Given the small numbers of young people bereaved by suicide it was beyond the scope of this research to compare different bereavement groups. However, Bartik et al. (2015) found that those bereaved by suicide were more likely than the general population to perceive suicide as stigmatising and in glorifying or normalising terms, and less likely to attribute it to isolation and depression. Future research should therefore endeavour to assess how attitudes impact help-seeking among those who are suicide bereaved, to better understand the relationship between attitudes and self-harm.
The finding that baseline self-harm ideation predicted selfharm group at follow-up is consistent with past research (Ribeiro et al., 2016) and the theoretical assertion that ideation/intention is a proximal predictor of engagement in behaviours (O'Connor, 2011;O'Connor and Kirtley, 2018). We also found that selfharm among family members could predict self-harm group membership at baseline. In a UK population-based cohort study, Mars et al. (2019a) showed that exposure to family self-harm was a predictor of future suicide attempts among adolescents who reported suicidal thoughts (but not those who engaged in NSSI). O'Connor et al. (2009a) found that adolescents who engaged in repeat self-harm over a 6-month period were also significantly more likely to have family and friends who self-harmed than those who did not report self-harm; however, only family selfharm remained a significant predictor in multivariate analyses. These findings may be explained by familial transmission of suicidal behaviour (O'Connor et al., 2009a;Pitman et al., 2014), possibly through increased risk from shared environmental stressors or genetic factors, or transmission of psychopathology and impulsive aggression (Brent et al., 2002;Melhem et al., 2007). On the other hand, the finding that self-harm of friends did not predict self-harm group status here may also be attributed to a lack of statistical power, as numerous studies have suggested a role for social modelling of self-harm among non-family members. Self-harm among peers significantly predicted future suicidal behaviour in four large-scale studies across various countries (De Leo and Heller, 2008), where the sample sizes ranged from n = 731 to 11,572, depending on the study time point. This effect is observed in studies specifically with adolescents (Hawton et al., 2002;Doyle et al., 2015). Given the small sample in this study, further work is required to test this using a larger dataset.
Overall, some support for the IMV model was found. That several factors did not predict self-harm group membership cross-sectionally nor longitudinally may likely be the result of limited statistical power. The baseline self-harm groups consisted of 38 people in the enactment group, 44 in the ideation group, and 103 controls. At follow up, there were only 33 individuals in the self-harm enactment group, 26 in the ideation group, and 56 controls. While the sample sizes were deemed adequate for the analyses chosen, they may not have been sufficient to detect group differences, if these existed, where cell sizes were small (e.g., one third of cells had values of less than 5 in the suicide death and family self-harm variables at T2). Risk factors for self-harm can vary significantly over time and even within a day (Kashyap et al., 2015), so estimating future outcomes from measures taken 6 months earlier is particularly challenging, especially when using a small sample. It is also possible that the IMV model does not appropriately model the relationship between certain variables, or may not be applicable to young people in a Scottish context. Given the absence of an association between various established risk factors (including defeat and entrapment) and self-harm in multivariate analyses, additional research is needed to determine whether these findings hold with a larger sample, and ultimately whether the model requires further refinement. Future research should also examine the difference between internal and external entrapment; we refrained from exploring this due to the small sample size.
Corrections have also been made to Discussion, Implications, Paragraph 1. The corrected paragraph is shown below: Implications These findings offer an important contribution to the limited literature on adolescent bereavement experiences and their relation to self-harm. Results highlight that self-harm ideation and behaviours are prevalent among Scottish youth, and a large proportion of adolescents have also been bereaved or exposed to the death of someone close to them. Given the potential consequences of bereavement, and particularly suicide bereavement with its association to adverse outcomes, understanding the extent and nature of this experience among adolescents is essential. As a test of a theoretical model, support is promising for some aspects of the IMV model, in particular identifying stigmatising and glorifying/normalising beliefs about suicide, and self-harm ideation as predictors of future behaviours. At the same time, evidence which was not wholly consistent with regards to the role of other variables within the IMV model, such as the impact of experiences of suicide loss, requires further investigation with larger samples to assess their placement within the model. In addition to guiding future research and theory refinement, our findings have implications on targets for clinical interventions and postvention. Efforts aimed at enhancing healthy coping skills, increasing family cohesion and social support, and addressing beliefs and attitudes about suicide (such as viewing suicide as stigmatised or glorified), and targeting self-harm ideation before it becomes severe, may be especially effective.
Lastly, corrections have been made to Tables 1, 2, 3, 4, 5 and 6, and a new Table 5b has been added.
Copyright © 2021 del Carpio, Rasmussen and Paul. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.