Abstract

This study used ecological momentary assessment (EMA) to explore the correlates of suicidal ideation (SI) instability in patients hospitalized for depression and SI. Thirty-nine adult inpatients were given smartphones with visual analogue scales to rate current depressed mood, anger/irritability, feeling socially connected, and SI three times a day throughout hospitalization. Affective Lability Scales (ALS) were also completed at baseline. SI instability was correlated with SI intensity, depressed mood instability, and social connection instability. Social connection instability was not associated with SI instability after controlling for depressed mood instability. ALS scores were not associated with EMA-derived SI instability. Participants with multiple past suicide attempts experienced greater SI instability. More research examining the clinical significance of SI instability is warranted.

INTRODUCTION

According to the DSM, major depression involves periods of low mood that are both stable (i.e., “present most of the day, nearly every day”) and persistent (i.e., “for at least 2 consecutive weeks”) (American Psychiatric Association, Citation2013, p. 163). In actuality, depressed individuals experience rapid and intense mood fluctuations over the course of hours/days (Bowen, Peters, Marwaha, Baetz, & Balbuena, Citation2017; Houben, Van Den Noortgate, & Kuppens, Citation2015; Köhling et al., Citation2016; Peters, Balbuena, Baetz, Marwaha, & Bowen, Citation2015; Thompson, Berenbaum, & Bredemeier, Citation2011; Thompson et al., Citation2012), including brief periods of high mood (Bylsma, Taylor-Clift, & Rottenberg, Citation2011; Peeters, Berkhof, Delespaul, Rottenberg, & Nicolson, Citation2006; Peeters, Nicolson, Berkhof, Delespaul, & deVries, Citation2003; Thompson et al., Citation2012), which are neither acknowledged in the DSM (American Psychiatric Association, Citation2013), nor captured by depression rating scales (Ben-Zeev, Young, & Madsen, Citation2009; Bowen, Balbuena, & Baetz, Citation2014; Solhan, Trull, Jahng, & Wood, Citation2009). The DSM-5 mixed features specifier allows for features of elevated mood (e.g., high energy, goal-directed activity) in depressive episodes, but the requirement that these be present for “the majority of days” discounts brief, recurrent mood elevations (American Psychiatric Association, Citation2013, p. 184). Mood instability is reported by up to 60% of people with depression when it is assessed directly by questionnaires (Peters et al., Citation2015; Solmi et al., Citation2016; Thompson et al., Citation2011), and is apparent when measured prospectively with ecological momentary assessment (EMA; Bowen et al., Citation2017; Bylsma et al., Citation2011; Houben et al., Citation2015; Köhling et al., Citation2016; Peeters et al., Citation2006; Solhan et al., Citation2009; Thompson et al., Citation2012; but see also Bos, de Jonge, & Cox, Citation2019; Koval, Pe, Meers, & Kuppens, Citation2013; Verdolini et al., Citation2019).

Mood instability is associated with suicidality and self-harm (Palmier-Claus, Taylor, Varese, & Pratt, Citation2012; Peters, Baetz, Marwaha, Balbuena, & Bowen, Citation2016; Peters, John, Bowen, Baetz, & Balbuena, Citation2018), even after controlling for depressed mood (Bowen et al., Citation2015; Peters, Balbuena, Marwaha, Baetz, & Bowen, Citation2016). Similar to mood, suicidal ideation (SI) can fluctuate considerably over a day (Forkmann et al., Citation2018; Hallensleben et al., Citation2018; Kleiman et al., Citation2017; Nock, Prinstein, & Sterba, Citation2009; Witte, Fitzpatrick, Joiner, & Schmidt, Citation2005; Witte, Fitzpatrick, Warren, Schatschneider, & Schmidt, Citation2006). Accordingly, there have been recent calls for EMA research to capture the dynamic nature of suicidality (Davidson, Anestis, & Gutierrez, Citation2017; Kleiman et al., Citation2017; Nock, Citation2016; Witte et al., Citation2006). So far, most EMA studies have examined factors that predict SI at subsequent time points (Ben-Zeev, Young, & Depp, Citation2012; Husky et al., Citation2017; Kleiman et al., Citation2017). It may also be pertinent to elucidate the factors associated with SI instability—the overall volatility and unpredictability of an individual’s desire to end their life. One study did find greater SI instability in individuals with multiple suicide attempts compared to those without (Witte et al., Citation2005, Citation2006), although contradictory results have been reported (Peters, Mereish, Solomon, Spirito, & Yen, Citation2019).

The main goal of this study was to use EMA to explore the effect of mood instability on SI instability in patients hospitalized for a major depressive episode. In addition to depressed mood, we also measured anger/irritability because it is associated with mood instability (Oliver & Simons, Citation2004), depression (Judd, Schettler, Coryell, Akiskal, & Fiedorowicz, Citation2013), and suicidality (albeit inconsistently, in clinical samples; Massimiliano, Perret, Turecki, & Geoffroy, Citation2018). Our second goal was to use EMA to measure instability in feeling socially connected/close to others, similar to the concept of belongingness from the interpersonal theory of suicide (Kleiman et al., Citation2017; Van Orden et al., Citation2010), and to examine potential associations with mood instability and SI instability. Past research has documented considerable within-person variability in feelings of loneliness (Kleiman et al., Citation2017), but it is unclear to what extent this correlates with mood and SI instability. It was hypothesized that instability in all three variables (i.e., depressed mood, anger/irritability, social connectedness), but not necessarily their average level of intensity, would be associated with SI instability, and that mood instability and social connectedness instability would be related. Additional goals were to compare the predictive value of questionnaire-rated mood instability on SI instability to EMA-derived mood instability; to compare SI instability between groups defined by the presence of past suicide attempts, in an attempt to replicate previous results (Witte et al., Citation2005, Citation2006); and finally, because mood instability has been associated with poorer treatment response (Stange et al., Citation2016) and symptom recurrence (Wichers et al., Citation2010), to determine if improvement in depression and SI throughout hospitalization are associated with mood and SI instability.

METHODS

Sample

The sample consisted of 39 patients (M age = 36.3 years, SD = 13.0; 69.2% female) hospitalized on a psychiatric unit with a chief diagnosis of a major depressive episode secondary to major depressive disorder or bipolar disorder (n = 12; 30.8%). To be included patients had to be 18-70 years old, be experiencing suicidal ideation at the time of admission, be admitted voluntarily, and not be withdrawing from any substances. All participants provided written informed consent and the study received local ethical approval. No compensation was provided.

Procedure

Recruitment began in December of 2016. Patients were referred to the study by their inpatient psychiatrist or social worker. A trained research assistant then met with potential participants, explained the study, obtained consent, administered the Mini-International Neuropsychiatric Interview (Sheehan et al., Citation1998) to confirm eligibility, showed participants how to use the phone, and had participants complete baseline questionnaire assessments (Time 1). Participants remained in the study until they were discharged. At this time questionnaire assessments were repeated (Time 2).

Materials and Measures

Beck Depression Inventory

Depressive symptoms were assessed at Time 1 (study entry) and Time 2 (study exit/discharge) with the Beck Depression Inventory (BDI; Beck, Steer, & Brown, Citation1996). The BDI has 21 items, scored 0 to 4, which are summed to yield a total score reflecting greater depression severity. The internal consistencies were good to excellent (α = .89 and .94). BDI change scores (Time 1–Time 2) were used to measure depression improvement throughout hospitalization.

Beck Scale for Suicide Ideation

SI at Time 1 and Time 2 was assessed with the Beck Scale for Suicide Ideation (BSS; Beck & Steer, Citation1991). The BSS has 21 items, scored 0–2, the first 19 of which are summed to yield a total score reflecting greater SI. The internal consistencies were excellent (α = .94 and .97). Responses to item 20 were used to assess past suicide attempts. BSS change scores (Time 1–Time 2) were used to measure SI improvement throughout hospitalization.

Affective Lability Scales

Mood instability at Time 1 was assessed with the short Affective Lability Scales (ALS; Oliver & Simons, Citation2004). The ALS consists of 18 items rated on a Likert scale from 1 to 4 (Oliver & Simons, Citation2004). There are six ALS subscales that represent rapid and intense fluctuations between different combinations of mood states: euthymia/anxiety, euthymia/depression, euthymia/elation, euthymia/anger, anxiety/depression, and depression/elation. The total ALS (the sum of all 18 items) and each subscale had adequate to excellent internal consistency (α = .75–.94).

EMA

Each participant was provided with an Android smartphone for the duration of the study. Phones were installed with Ethica (Hashemian, Qian, Stanley, & Osgood, Citation2012; Knowles, Stanley, & Osgood, Citation2014) to facilitate data collection. This application employed four questions, each rated on a visual analogue scale from 0 to 100, designed to reflect the following constructs: depressed mood (“How depressed do you feel right now?”), anger/irritability (“How irritable/angry do you feel right now?”), perceived social connectedness (“How close or connected to people do you feel right now?”), and suicidality (“How suicidal are you right now?”). Participants were prompted to answer each question three times a day at 09:00, 15:00, and 20:00, with a 60-minute window to answer each question, after which the response was set to missing if no response was given.

Analysis

Missing Data

Two participants did not complete any Time 1 questionnaires, and 14 participants were discharged before completing Time 2 questionnaires. Imputation was not attempted for these individuals, and analyses with questionnaire data had reduced sample sizes as a result (Time 1, n = 37; BDI/BSS change scores, n = 25). If questionnaires were completed but some items were missing then these were filled using random forest imputation (Stekhoven & Bühlmann, Citation2012). For each questionnaire, less than 1% of responses were missing.

Creation of EMA Variables

Responses to the EMA questions were used to create variables representing intensity and instability. Intensity was defined as the mean of all ratings for a question over the study period. Instability was defined as the mean squared successive difference (von Neumann, Kent, Bellinson, & Hart, Citation1941), which takes into account both the within-person variability and the temporal sequence between measurements (Solhan et al., Citation2009). Mean square successive difference scores were left skewed but were normalized after a square root transformation. Social connection instability was not normally distributed after transformation, but this corrected after removing one strong outlier (score = 92.8; range excluding this score = 2.32–44.8) that was larger than three times the IQR (19.8) added to the 75th percentile (33.0). Most intensity scores were also left skewed, so they were also normalized by square rooting. The final analysis used eight EMA variables: SI intensity (SIINT), SI instability (SIINST), depressed mood intensity (DEPINT), depressed mood instability (DEPINST), anger/irritability intensity (ANGINT), anger/irritability instability (ANGINST), social connection intensity (CONINT), and social connection instability (CONINST).

Exploratory Correlations

Pairwise correlations between all variables were calculated. Prior to calculating each correlation, scatter plots were used to assess for linearity and heteroscedasticity. Pearson’s r was calculated if assumptions were met. Nonparametric Kendall’s τB correlation coefficients were used to explore the associations for the following variables that were not normally distributed: ALS euthymia/anger, ALS anxiety/depression, BDI change scores, and BSS change scores.

Multiple Regression Models

Multiple linear regression models predicting SIINST were created to observe what variables remained associated with SIINST when multiple predictors were included in the same model. In general, we were interested in models that included different instability variables, pairs of EMA intensity and instability variables (e.g., DEPINT and DEPINST), and instability variables with their ALS counterparts (e.g., DEPINST and the ALS euthymia/depression subscale). We did not include any variables in the regression analysis that had nonsignificant relationships with SIINST in the correlation analyses. However, given the small sample size and possibility of type II error, we did add variables with small but nonsignificant p values if the effect size was not trivial (r ∼ .20 or higher; Dupont, Citation2009).

Group Comparisons

We conducted Wilcoxon rank-sum tests to compare SIINST between groups defined by past suicide attempts on the BSS: no past attempts (n = 9) vs. past attempt(s) (n = 28), one past attempt (n = 12) vs. multiple past attempts (n = 16), no past attempts (n = 9) vs. one past attempt (n = 12). We chose a nonparametric test due to the small and unbalanced group sizes, as well as concerns about data normality. Two participants had outlying scores, so the comparisons were repeated without their data to ensure the results did not differ. Differences in SIINST could be driven by differences in DEPINST, so we also compared DEPINST between any groups that had significant differences in SIINST.

RESULTS

EMA Usage and Compliance

The average number of days participants were enrolled in the study was 12.0 (SD = 5.61). VAS completion rates by study duration are presented in . The correlation between total VAS completion rates and study duration was not statistically significant (r = .07, p = .69). Total VAS completion rates for each participant were not significantly associated with Time 1 BDI scores (r = .26, p = .12), DEPINT (r = .17, p = .30), DEPINST (r = −.14, p = .40), Time 1 BSS scores (r = .19, p = .27), SIINT (r = .03, p = .88), or SIINST (r = −.10, p = .56).

FIGURE 1. Graph depicting EMA compliance as a function of VAS completion rates at each day of the study. EMA: ecological momentary assessment; VAS: visual analogue scales.

FIGURE 1. Graph depicting EMA compliance as a function of VAS completion rates at each day of the study. EMA: ecological momentary assessment; VAS: visual analogue scales.

Exploratory Correlations

Pairwise correlations between EMA variables are presented in and . SIINST was significantly positively correlated with SIINT, DEPINST, and CONINST. Bivariate correlations between EMA and questionnaire variables are presented in . ALS total and subscale scores were not significantly correlated with any EMA instability variable. Time 1 BDI scores (M = 42.1, SD = 10.8) were positively correlated with DEPINT, SIINT, and CONINST. Time 1 BSS scores (M = 20.4, SD = 10.1) were positively correlated with DEPINT, SIINT, and ANGINT; the correlation with SIINST was small and borderline significant (r = .28, p = .09). BDI and BSS change scores were not significantly correlated with any EMA variable.

FIGURE 2. Bivariate correlations between EMA variables. Variables were derived from ratings on visual analogue scales measuring suicidal ideation (SI), depressed mood (DEP), feeling socially connected/close to others (CON), and anger/irritability (ANG). Intensity variables were defined as the mean of all ratings. Instability variables were defined as the mean squared successive difference. Correlation not shown if r < .20 and p > .20, except for CONINT with DEPINT (r = −.22, p =.17) and ANGINST (r = −.21, p = .19). EMA: ecological momentary assessment. *p < .05; **p < .01; ***p < .001.

FIGURE 2. Bivariate correlations between EMA variables. Variables were derived from ratings on visual analogue scales measuring suicidal ideation (SI), depressed mood (DEP), feeling socially connected/close to others (CON), and anger/irritability (ANG). Intensity variables were defined as the mean of all ratings. Instability variables were defined as the mean squared successive difference. Correlation not shown if r < .20 and p > .20, except for CONINT with DEPINT (r = −.22, p =.17) and ANGINST (r = −.21, p = .19). EMA: ecological momentary assessment. *p < .05; **p < .01; ***p < .001.

TABLE 1. Pairwise correlations between EMA variables.

TABLE 2. Correlations between EMA indices (first column) and questionnaire scores (first row).

Multiple Regression Models

First we estimated a model with Time 1 BSS scores and SIINT predicting SIINST; the effect of SIINT remained significant (β = .61, SE = 1.02, p = .006), whereas the effect of Time 1 BSS scores became negative and nonsignificant (β = −.16, SE = 0.26, p = .47). Therefore, we did not include Time 1 BSS scores in any subsequent models. Based on the remaining results, we estimated models predicting SIINST with different combinations of DEPINST, CONINST, and SIINT (see ). Both DEPINST and CONINST remained significant after adding SIINT to the model. When DEPINST and CONINST were included in the same model the coefficient for CONINST became nonsignificant, whereas the coefficient for DEPINST remained significant.

TABLE 3. Multiple linear regression models predicting suicidal ideation instability.

We conducted additional analyses with ANGINST because the effect size was roughly “small” and it was one of the hypothesized correlates of SIINST. In a model with SIINT, the effect remained largely unchanged (β = .22, SE = 0.15, p = .12), although it became trivial in size with DEPINST in the model (β = −.03, SE = 0.17, p = .83); DEPINST had a significant effect (β = .54, SE = 0.18, p = .001).

We also estimated exploratory models with ALS euthymia/anxiety scores predicting SIINST given the size of the correlation (the largest of all ALS variables) and the relatively small p value (r = .27, p = .10). SIINT and CONINST were used as covariates as these had the largest associations with ALS euthymia/anxiety scores. The effect decreased after controlling for SIINT (β = .16, SE = 0.67, p = .30) or CONINST (β = .16, SE = 0.61, p = .29), and was trivial in size with both covariates in the model (β = .09, SE = 0.57, p = .54).

Group Comparisons

SIINST group means and comparisons are presented in . The only statistically significant difference was between participants with multiple past attempts compared to individuals with one past attempt. DEPINST scores did not significantly differ between these groups (M = 21.0 and 23.0, respectively; z = 0.33, p = .75).

TABLE 4. Mean suicidal ideation instability with comparisons by past suicide attempt status.

DISCUSSION

It should be first stated that given the large number of exploratory analyses, significant correlations may have occurred due to chance; the results could be useful for generating hypotheses, but replication is required before drawing firm conclusions.

Correlates of SI Instability

SI intensity was strongly associated with SI instability, demonstrating that patients who were more suicidal, on average, throughout hospitalization also tended to have greater SI instability. The potential significance of this finding is that rapid but temporary reductions in SI could be misinterpreted as improvement in patients who are likely to remain more severely suicidal during hospital admissions. The only other EMA variables clearly associated with SI instability were depressed mood instability and instability in feelings of social connectedness, although the former may be more proximally related to SI instability than the latter. In contrast to our expectations, anger/irritability instability was not significantly associated with SI instability, and the small nonsignificant association that was present appeared to be due to the correlation with depressed mood instability.

Among the questionnaire variables, Time 1 BSS scores had the largest association with SI instability, consistent with EMA-derived SI intensity also being correlated with SI instability. The finding that BDI and BSS change scores were unrelated to SI instability and depressed mood instability suggests that, in this sample, it would have been inappropriate to conflate instability and improvement when it comes to depression and suicidality. However, missing data makes this result difficult to interpret as we cannot rule out the possibility that participants who were unreachable at Time 2 had different outcomes.

Mood Instability Questionnaires vs. EMA Indices

One noteworthy finding is that ALS scores were not significantly associated with SI instability, although the correlation with the ALS euthymia/anxiety subscale was nontrivial in size and approached statistical significance, in keeping with previous research (Ben-Zeev et al., Citation2012; Berman, Citation2018; Forkmann et al., Citation2018; Husky et al., Citation2017). The ALS and EMA-derived mood instability indices purportedly measure the same phenomenon (Marwaha et al., Citation2014; Solhan et al., Citation2009), yet this finding suggests the ALS may not sufficiently discriminate individuals with unstable SI. ALS scores were also unrelated to the other EMA instability variables, further suggesting these may not be measuring the same underlying construct. A similar finding was reported by Links, Heisel, and Garland (Citation2003), who suggested the ALS may be a better index of affective intensity than affective instability. One ostensible difference is that the ALS measures shifts between discrete emotions (i.e., euthymia to anger, anxiety to depression) and not fluctuations within a particular emotional state (Oliver & Simons, Citation2004). A question for future research is whether shifts between emotions or within emotional states are similarly related to suicidality, and the extent to which these phenomena are independent.

SI Instability by Past Suicide Attempts

Consistent with a previous study of a student sample (Witte et al., Citation2005, Citation2006), we found that participants with multiple past suicide attempts experienced the highest degree of SI instability, although unlike this past study, we did not find higher SI instability in participants with a previous attempt compared to those without (Witte et al., Citation2005, Citation2006). This discrepancy could have been due to the sample differences, as most of the subgroup without past attempts would not have been depressed (Witte et al., Citation2005, Citation2006). A more recent longitudinal study of hospitalized adolescents found that SI instability did not predict future suicide attempts, nor was SI instability associated with the mean of all SI ratings over a 6-month period (Peters et al., Citation2018). SI in this study was measured weekly with questionnaires (Peters et al., Citation2018), so it may be misleading to compare “instability” variables with differing sampling frequencies (Ebner-Priemer & Sawitzki, Citation2007; Thompson et al., Citation2012). More research using diverse samples, additional measurement schemes and covariates, and longitudinal designs is needed to fully understand the clinical significance of SI instability, particularly in relation to future suicide attempts.

Feasibility of Smartphone EMA

The current study adds to the existing research supporting the feasibility of smartphone EMA among depressed inpatients (Ben-Zeev et al., Citation2012; Hallensleben et al., Citation2018; Husky et al., Citation2014; Kleiman et al., Citation2017). Previous studies have tended to limit sampling to one week (Ben-Zeev et al., Citation2012; Hallensleben et al., Citation2018; Husky et al., Citation2014), although we sampled participants for the duration of admission (up to 24 days). Therefore, it is notable that EMA compliance did not decrease as sampling durations increased. Furthermore, we did not find EMA compliance to be significantly associated with baseline depression/SI severity, or EMA-derived depression/SI intensity during hospitalization, further suggesting illness severity among depressed patients does not preclude the use of EMA.

Limitations

Major limitations were the small, heterogeneous sample (i.e., patients with unipolar or bipolar depression, and various comorbidities) and the lack of long-term follow-up that prevented us from considering outcomes such as readmission to hospital and future suicide attempts. It is possible that some of the smaller correlations in the .20–.29 range would have been statistically significant had a larger sample been recruited. We also did not compare patients admitted for a suicide attempt to the rest of the sample, nor were treatments monitored. Patients were not excluded if they had comorbid diagnoses that could have inflated estimates of mood/SI instability, such as personality disorders or PTSD (American Psychiatric Association, Citation2013; Marwaha et al., Citation2014; Santangelo et al., Citation2014), although doing so would have limited the generalizability of the sample. The fact that we did not measure anxiety/agitation is another limitation, as it tends to be associated with SI (Ben-Zeev et al., Citation2012; Berman, Citation2018; Forkmann et al., Citation2018; Husky et al., Citation2017). EMA compliance was generally satisfactory, although missing data on Time 2 questionnaires was common and further limits the results from these analyses. This study was also not designed to answer the question of whether depressed individuals experience greater mood instability than nondepressed individuals, which remains debated (Bos et al., Citation2019; Koval et al., Citation2013).

CONCLUSIONS

SI instability in depressed inpatients can be successfully captured with smartphone EMA. These preliminary results suggest that SI instability is associated with more intense suicidality and depressed mood instability during hospitalization. A single assessment with the ALS may not adequately distinguish depressed patients with greater SI instability from patients with more stable SI. SI instability also seems to be more pronounced in patients with multiple past suicide attempts, although the clinical significance of this finding in relation to future suicide attempts remains unknown. More research on SI instability using EMA is warranted.

Additional information

Notes on contributors

Evyn M. Peters

Evyn M. Peters, Department of Psychiatry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

Lisa Yue Dong

Lisa Yue Dong, Department of Psychiatry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

Teena Thomas

Teena Thomas, Department of Computer Science, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

Sara Khalaj

Sara Khalaj, Department of Psychiatry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

Lloyd Balbuena

Lloyd Balbuena, Department of Psychiatry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

Marilyn Baetz

Marilyn Baetz, Department of Psychiatry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

Nathaniel Osgood

Nathaniel Osgood, Department of Computer Science, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

Rudy Bowen

Rudy Bowen, Department of Psychiatry, University of Saskatchewan, Saskatoon, Saskatchewan, Canada.

REFERENCES