1,007
Views
0
CrossRef citations to date
0
Altmetric

Abstract

Suicidal ideation (SI) is a significant public health concern with increasing prevalence. Therapist-supported digital mental health interventions (DMHI) are an emergent modality to address common mental health problems like depression and anxiety, although less is known about SI. This study examined SI trajectories among 778 patients who participated in a therapist-supported DMHI using multilevel models during and up to 6-months post-treatment. Estimates of associated suicide attempts and deaths by suicide were calculated using published data linking PHQ-9-assessed SI to records of suicide attempts and deaths by suicide. The proportion of participants reporting no SI significantly increased between baseline and end-of-treatment (78.02% to 91.00%). Effect sizes of SI changes between baseline and end-of-treatment, 3-month, and 6-month follow-ups were 0.33 (95%CI = 0.27–0.38), 0.32 (95%CI = 0.27–0.38), and 0.32 (95%CI = 0.27–0.38), respectively. Results also indicated an estimated 30.49% reduction (95%CI = 25.15%-35.13%) in suicide attempts and death by suicide across treatment. This study provides preliminary evidence of the effectiveness of a therapist-supported DMHI in reducing SI.

Suicide continues to be a leading cause of preventable mortality in the United States (US). In 2019, it was estimated that approximately 5% of the US population reported suicidal ideation (SI) in the past year (Stone et al., Citation2021). Importantly, SI is a strong predictor of future suicide behaviors, such that an estimated 12.5% to 29% of those with SI end up attempting to die by suicide (Nock et al., Citation2008). Unfortunately, the US is currently experiencing a mental health crisis characterized by large-scale structural barriers (e.g., long waitlists, high costs, transportation issues, lack of trained providers) to delivering effective suicide prevention and treatment interventions (Moutier, Citation2021). As such, it is estimated that only 44% of people experiencing SI and 61% of those that have attempted suicide seek mental health treatment (Ahmedani et al., Citation2012; Kessler et al., Citation2005), pointing to a dire need for novel, scalable, and cost-effective solutions to reduce SI and associated health outcomes at the population level.

Digital mental health interventions (DMHI) offer one possible solution to address SI and associated mental health problems (Büscher et al., Citation2020; Torok et al., Citation2020). Prior research has consistently demonstrated that DMHIs incorporating cognitive-behavioral therapy (CBT) and other evidence-based modalities (e.g., mindfulness, behavioral activation, stress reduction) are more effective than treatment as usual, waitlist controls, and attention controls in reducing common mental health symptoms (Goldberg et al., Citation2022). DMHI interventions also provide a novel avenue for population-level scalability by lowering the barriers to treatment, improving cost-efficiency, reducing stigma, and ameliorating expected shortages of behavioral health providers (Andrilla et al., Citation2018; Askari et al., Citation2023; Kaur et al., Citation2023). Moreover, emergent research indicates that DMHIs may be useful in reducing SI (Büscher et al., Citation2020; Hetrick et al., Citation2017; Sander et al., Citation2020), yet there is comparatively less research about the association between therapist-supported DMHIs and reductions in SI. Therapist-supported DMHIs represent a unique modality for addressing SI, as they facilitate a remote continuous care approach whereby therapists first establish rapport and trust with patients that provide the foundation for a digital therapeutic alliance over the duration of the intervention (Forman-Hoffman et al., Citation2021; Krieger et al., Citation2023). As such, there remains a crucial need for real-world research that evaluates the potential effectiveness of therapist-supported DMHIs in reducing SI.

The primary goal of this study was to investigate the feasibility and preliminary effectiveness of a 12-week, therapist-supported DMHI in reducing SI among people experiencing clinically significant symptoms of depression and anxiety. We hypothesized that treatment week would be negatively associated with SI, with significant reductions in SI from baseline to end-of-treatment as well as at 3- and 6-months post-treatment. The secondary goal was to explore potential reductions in suicide attempts and mortality during and after the DMHI using estimates from a large electronic health records database (Rossom et al., Citation2017; Simon et al., Citation2013). We hypothesized that participation in the DMHI would be associated with significant reductions in estimated suicide attempts and mortality from baseline to end-of-treatment and both post-treatment follow-ups.

METHODS

Participants

The present study included patients treated with the Meru Health Program (MHP), an evidence-based therapist-supported DMHI (Goldin et al., Citation2019; Peiper et al., Citation2023). All data came from 778 adult participants in the United States who started the program on or after January 1, 2020 and had an end date on or before December 1, 2020. Participants entered the MHP via referrals from employee assistance programs or via their healthcare providers. Inclusion criteria of the MHP required patients to have at least mild levels of depression, anxiety, or burnout; own a smartphone; and not have an active substance use disorder, severe active suicidal ideation with a specific plan, severe active self-harm, or a history of psychosis or mania.

All enrolled patients provided informed consent to participate and have their de-identified data used for research purposes. Data collected as part of care were stored in Health Insurance Portability and Accountability Act-compliant electronic medical records that included protected health information. All data were encrypted in transit, and at rest. Institutional review board exemption for this analysis was obtained from the Pearl Institutional Review Board (#21-Meru-113) for analyses of previously collected and de-identified data.

Intervention Description

After the informed consent procedure, patients were enrolled in the intervention. The MHP is composed of several different evidence-based components or processes of evidence-based therapies delivered via a smartphone app. The program is self-guided, but incorporates a stepped and continuous care model that includes daily interaction with a dedicated, licensed clinical therapist and a medical doctor and psychiatrist available for consultation as needed. The program lasts 12 weeks and contains components of cognitive-behavioral therapy (CBT), behavioral activation, mindfulness-based stress reduction, sleep therapy, nutritional psychiatry, and heartrate variability biofeedback (Carney et al., Citation2017; Lehrer et al., Citation2013; Morgan, Citation2003; Sarris et al., Citation2015). The MHP sequentially delivers self-guided modules covering specific topics over a 12-week period: 1) mindfulness; 2) negative cognition; 3) mood and motivation; 4) rumination; 5) emotional regulation; 6) self-compassion; 7) introspection; 8) social functioning; 9) sleep habits; 10) healthy diet; 11) nutrition and mood; and 12) resilience.

Prior to the start of treatment, participants are trained on how to use the app, including how to participate in the anonymous group interaction and how to communicate with their assigned therapist via chat or phone/video calls. Each week of the program begins with an introductory video that gives the patient an overview of the topics covered that week. Patients are prompted to complete various practices, including CBT exercises and prompts from the app to record progress in a journal, as well as short audio-guided mindfulness meditation exercises on a daily or almost daily basis. Daily content and practices range from 5–15 minutes, except for the first day of each week, in which the weekly psychoeducation video lessons extend the content to a maximum of 25 minutes.

A licensed therapist provides support to participants via messaging (and less frequently, phone or video calls) as needed, and reviews practice logs using a provider dashboard and electronic medical records that detail individual progress (including participant engagement and patient-reported outcomes to date). In total, therapists allocate approximately 10–20 minutes (on average) per week per participant. Interaction between therapist and participant can be initiated via either party.

As a safety measure, therapists conduct a phone-based protocol assessment for any participants that show signs of mental deterioration during the intervention. In case of an emergency, such as having SI with intent to act or onset of psychotic symptoms, the intervention includes a written action plan for declining mental health, which all participants are required to review with their therapist prior to engaging with the intervention. In these situations, care coordinators and therapists will help connect the patient to immediate and local care outside of the program, however, a psychiatrist employed by Meru Health is also available for consultation in these situations.

Patients complete the intervention as a cohort of typically 8–15 patients. As such, cohorts can interact anonymously when sharing their practice experiences, such as providing and receiving support and feedback on their experiences as they navigate the intervention. Direct communication between participants, however, is not allowed. Instead, participants can post anonymous reflections on practices and lessons to the chat discussion board, to which their therapist can respond freely, and to which other group members can respond with pre-written empathy statements and/or emoticons.

Measures

SI and depression symptoms were collected during the intervention on a biweekly basis within the MHP app. For the 3- and 6-month post-treatment follow-ups, patients received emails with instructions on how to report their mental health symptoms in a secure web-based platform. Participants completed the PHQ-9 during and after the intervention. The PHQ-9 is a clinically validated and widely used instrument used to screen for depression (Kroenke et al., Citation2001). We used the first 8 items of the PHQ-9 to calculate the PHQ-8, which assessed purely depressive symptoms without the last item that assesses SI, the primary outcome of interest (Kroenke et al., Citation2009). We coded treatment baseline PHQ-8 symptoms using categories established for the PHQ-9 with the 9th item removed as None/Minimal (0–4), Mild (5–9), Moderate (10–14), Moderately Severe (15–19), and Severe (20–24) to assess baseline depression severity.

Self-reported measures collected during clinical intake included age, gender, number of lifetime major depressive episodes (none, first episode, or recurrent), lifetime suicide attempts (yes/no), lifetime prior psychiatric hospitalization (yes/no), lifetime exposure to a major traumatic event (yes/no), and current psychotropic medication use (yes/no). Intervention engagement was defined as the number of days active in the app each week and captured as a core app metric.

We estimated reductions in suicide attempts and death by suicide using published risk data aggregated from 297,290 adult outpatients who received mental healthcare or had a mental health diagnosis documented in primary care (Rossom et al., Citation2017). Based on these findings, there was a 0.40% risk of suicide attempt and 0.03% risk of suicide death for 0 (not at all), 1.60% and 0.12% risk for 1 (several days), 2.80% and 0.21% risk for 2 (more than half the days), and 4.00% and 0.30% risk for 3 (nearly every day).

Statistical Analysis

All statistical analyses were conducted with RStudio, Version 1.3.959. The study was preregistered on Open Science Framework (https://osf.io/eug8w/). Statistical significance was defined as a p-value of 0.05. Exploratory analyses including histograms as well as skew and kurtosis statistics were run for each variable to check for normality. Any variable that had a skew of +/- 2 was log transformed, except for the SI item (skew = 3.81). Descriptive statistics were calculated for each patient demographic and clinical variable. Outcome measures were analyzed using an intention-to-treat (ITT) analysis in which all participants with outcome measures at baseline were included, regardless of intervention engagement or attrition. Approximately 12.7% of data were missing from the final sample (note this study used ITT, so we included all treatment weeks 0–12 and 3- and 6-month follow-up data as missing for participants that dropped out). To assess whether data were missing completely at random (MCAR) we performed a parametric test using the naniar package (Tierney & Cook, Citation2018). Although the data were not MCAR (p < 0.001), we performed multiple imputation (10 imputations) using the mice package to account for missing data based on recent recommendations (van Buuren & Groothuis-Oudshoorn, Citation2011). Additional details about missing data procedures may be found in the supplemental materials.

First, we examined the impact of treatment week on SI by conducting a multilevel model (MLM) that nested treatment week (Level-1) within individuals (Level-2) across linear, Poisson, and negative binomial models to account for the skewed distribution of the SI item (skew = 3.81). We used the lme4 package for linear and negative binomial models and the glmmTMB package for the Poisson model (Bates et al., Citation2014; Brooks et al., Citation2017). Fixed effects including covariates were tested at the level of participants (Level-2). This statistical approach accounts for dependency within participants and introduces less bias due to missing data compared to traditional statistical analyses, such as repeated measures of analysis of variance (Hoffman & Stawski, Citation2009). We tested concurrent associations between treatment week and SI during the same week (Level-1).

Second, we used ggstatsplot package to conduct a nonparametric Friedman Test, an alternative to a repeated measures ANOVA for nonnormal data, to determine whether there were significant reductions in SI from baseline to end-of-treatment as well as at the 3- and 6-month post-treatment follow-ups (Patil, Citation2021). We then used the rstatix package to conduct nonparametric paired samples Wilcoxon Test contrasts that were corrected for multiple comparisons using the Holm Method to determine whether there was a significant reduction in SI from baseline to end-of-treatment as well as at the 3- and 6-month post-treatment follow-ups (Holm, Citation1979; Kassambara, Citation2020).

Third, we calculated frequency statistics on the response options of the SI item at baseline and end-of-treatment. We then extrapolated suicide attempt and mortality data from Rossom et al. (Citation2017) to calculate estimated suicide attempts and mortality at baseline, end-of-treatment, 3-months post-treatment, and 6-months post-treatment. The percent reduction in suicide attempts and mortality at post-treatment timepoints was calculated with the following equation: (baseline percent - post-treatment percent)/baseline percent*100. We used the equation 1/(baseline percent-post-treatment percent) to calculate the number needed to treat (NNT) for both suicide attempts and deaths by suicide outcomes.

RESULTS

Demographic and Clinical Characteristics

presents demographic and clinical characteristics for the 778 patients. Most patients were female (77.12%) with an average age of 40. The mean baseline PHQ-8 was 11.1, which falls into the moderate range of symptom severity. Just under half of the sample (44.22%) reported taking psychiatric medication, while 42.03% reported a past traumatic event. Approximately 4.5% of participants reported a lifetime suicide attempt and 5.01% reported a lifetime psychiatric hospitalization.

TABLE 1. Patient demographic and clinical characteristics (N = 778).

Changes in Suicidal Ideation

The proportion of participants reporting SI “not at all” significantly increased between baseline and end-of-treatment (78.02% to 91.00%, see ). Conversely, there was a significant decline in reported SI during treatment that was maintained at 3- and 6-month follow-ups (). Effect sizes of SI changes between baseline and end-of-treatment, baseline and 3-month, and baseline and 6-month follow-ups were Hedge’s g = 0.33 (95%CI = 0.27–0.38), Hedge’s g = 0.32 (95%CI = 0.27–0.38), and Hedge’s g = 0.32 (95%CI = 0.27–0.38), respectively, indicating moderate effect sizes. The significant declines in SI remained after adjusting for demographic and clinical characteristics across linear (β = −0.110, p < 0.001), Poisson (log mean = −0.158, p < 0.001), and negative binomial (log-mean = −0.158, p < 0.001) models (). Higher engagement was also significantly associated with declines in SI across all three models (p’s < 0.001). The trajectories of SI by depression symptom severity appeared to be steepest in the first few weeks of treatment and then tapered off toward the end-of-treatment and continued to be stable over the 3-and 6-month follow-up assessments for each of the severity groups (see for individual and pretreatment symptom severity group trajectories). Individual participant depression symptom trajectories are displayed in Supplementary Figure S1.

FIGURE 1. Severity of suicidal ideation during and after treatment (N = 778).

Notes: Suicidal ideation over the last 2 weeks is measured by the ninth item of the Patient Health Questionnaire-9 item (PHQ-9), with 0 = not at all, 1 = several days, 2 = more than half the days, and 3 = nearly every day. Non-significant Wilcoxon Test comparisons for changes in suicidal ideation are not shown. That is, there were no significant differences between end-of-treatment and 3-month follow-up, between end-of-treatment and 6-month follow-up, or between 3-month and 6-month follow-up.

FIGURE 1. Severity of suicidal ideation during and after treatment (N = 778).Notes: Suicidal ideation over the last 2 weeks is measured by the ninth item of the Patient Health Questionnaire-9 item (PHQ-9), with 0 = not at all, 1 = several days, 2 = more than half the days, and 3 = nearly every day. Non-significant Wilcoxon Test comparisons for changes in suicidal ideation are not shown. That is, there were no significant differences between end-of-treatment and 3-month follow-up, between end-of-treatment and 6-month follow-up, or between 3-month and 6-month follow-up.

FIGURE 2. Suicidal ideation during and after treatment by baseline depressive symptom severity.

Notes: Suicidal ideation over the last 2 weeks is measured by the ninth item of the Patient Health Questionnaire-9 item (PHQ-9), with 0 = not at all, 1 = several days, 2 = more than half the days, and 3 = nearly every day. Treatment week 0 is the pretreatment (baseline) assessment. Treatment week 12 is the end-of-treatment assessment. Treatment week 25 is the 3-month follow-up assessment. Treatment week 38 is the 6-month follow-up assessment. Symptom severity is measured by baseline PHQ-8 score as follows: less than mild = 0–4, mild = 5–9, moderate = 10–14, moderately severe = 15–19, severe = 20+. Light lines are individual trajectories and dark lines are pretreatment symptom severity group trajectories.

FIGURE 2. Suicidal ideation during and after treatment by baseline depressive symptom severity.Notes: Suicidal ideation over the last 2 weeks is measured by the ninth item of the Patient Health Questionnaire-9 item (PHQ-9), with 0 = not at all, 1 = several days, 2 = more than half the days, and 3 = nearly every day. Treatment week 0 is the pretreatment (baseline) assessment. Treatment week 12 is the end-of-treatment assessment. Treatment week 25 is the 3-month follow-up assessment. Treatment week 38 is the 6-month follow-up assessment. Symptom severity is measured by baseline PHQ-8 score as follows: less than mild = 0–4, mild = 5–9, moderate = 10–14, moderately severe = 15–19, severe = 20+. Light lines are individual trajectories and dark lines are pretreatment symptom severity group trajectories.

TABLE 2. Suicidal ideation during and after treatment.

TABLE 3. Adjusted associations of changes in suicidal ideation across meru health program treatment.

Correlates of Changes in Suicidal Ideation

Significant correlates of change in SI across all models included average days of weekly activity (-), having a history of recurrent MDEs versus no lifetime MDEs (+), and having either moderately severe or severe baseline depression symptoms versus less than mild depression symptoms (+). Other significant correlates included age (- in Poisson and negative binomial models only), having a single MDE versus no lifetime MDEs (+ in the Poisson model only), having a history of a psychiatric hospitalization (+ in the linear model only), and mild or moderate depression symptom severity versus less than mild depression symptoms (+ in Poisson and negative binomial models only).

Estimates of Change in Suicide Attempts and Deaths

Suicide attempt estimates fell from 0.75% to 0.52% from baseline to end-of-treatment and further declined to 0.50% and 0.49% at 3- and 6-month follow-up. For suicide mortality, estimates reduced from 0.033% to 0.003% between baseline and end-of-treatment that were maintained at the 3- and 6-month follow-ups (0.003% and 0.0003%, respectively). These represent a 30.49% reduction (95%CI = 25.15%–35.13%) in estimates of both suicide attempts and death by suicide between baseline and end-of-treatment. Extrapolating findings to number needed to treat (NNT) estimates for end-of-treatment, analyses revealed that one suicide attempt would be prevented for every 438 participants enrolled in the intervention and one death by suicide would be prevented for every 5841 participants enrolled.

DISCUSSION

The study investigated whether SI decreased during and following participation in a 12-week, therapist-supported DMHI for depression and anxiety symptoms. In addition, the study aimed to estimate reductions in suicide attempts and deaths by suicide among this sample. Overall, this study shows moderate evidence of preliminary feasibility and effectiveness of improving SI among patients who participated in a therapist-supported DMHI.

Our hypotheses were supported by our findings. First, treatment week was associated with significant decreases in SI across the MLM models independent of pretreatment depression symptom severity, demographic factors, and clinical characteristics. In addition to treatment week effects, these models revealed that a history of MDE, particularly recurrent history of MDE, was associated with higher SI across treatment and that higher pretreatment depression symptom severity was associated with greater SI as would be expected. Older age was also associated with lower SI, which coincides with findings that SI tends to have a negative association with age (Dennis et al., Citation2007). Moreover, the number of days patients engaged with therapy each week was associated with lower rates of SI, which is consistent with recent research indicating that twice vs once weekly therapy sessions is more effective at reducing depression symptoms and attrition (Bruijniks et al., Citation2020). This finding suggests DMHIs that prioritize frequent practice of psychotherapy skills may be particularly useful in lowering SI and other symptoms of depression. A history of suicide attempt and psychiatric hospitalization were not significantly associated with SI in the past two weeks, except for the latter in linear models, which did not fit the data as well as Poisson and negative binomial models (see Table S1). As previous research indicates that risk factors of SI differ from those for the transition from SI to attempt, additional studies will be necessary to better understand how DMHIs may impact different suicide pathways (e.g., ideation, planning, attempts) and complex etiological mechanisms (Nock et al., Citation2016)

Results also indicated that there was a significant overall reduction in SI from baseline to post-treatment, which was maintained at both 3- and 6-month follow-ups, indicating a robust and long-lasting association between treatment on SI trajectories across time. Importantly, while SI continued to decrease from post-treatment to 3- and 6-month follow-ups, there was no significant decrease, possibly due to a floor effect. Finally, there was a greater than 30% decrease in both estimated suicide attempts and deaths by suicide during the 12-week DMHI. The estimated NNT of 438 and 5,841 for attempts and deaths, respectively, suggests that therapist-supported DMHIs focusing on depression and anxiety symptoms may be a viable suicide prevention strategy through remote continuous monitoring and post-treatment follow-ups (Forman-Hoffman et al., Citation2021; Krieger et al., Citation2023).

While this study had a number of significant strengths including a large sample size of patients with a diverse history of psychiatric and clinical characteristics with an ITT study design, the study findings should be interpreted in light of several limitations. Because this study utilized a single arm intervention design with observational data, pragmatic randomized controlled trials that include relevant control conditions (e.g., treatment as usual, waitlist) are necessary to produce stronger evidence of the effectiveness of therapist-supported DMHIs in reducing SI. Nevertheless, the ITT design and use of multiple imputation procedures increased the study’s rigor compared to other studies that commonly employ completer-only analyses. Another limitation is the reliance on the 9th item of the PHQ-9 as our SI measure. Although the 9th item of the PHQ-9 has been shown to be a robust predictor of suicidal behavior in study populations (Louzon et al., Citation2016; McCusker, Citation2016; Rossom et al., Citation2017; Simon et al., Citation2016; Yarborough et al., Citation2021), some research indicates that this item may not be a sufficient screen for SI (Razykov et al., Citation2012). Furthermore, the predictive ability of the SI item may vary by demographic and clinical characteristics (Na et al., Citation2018), suggesting it may be limited in distinguishing individual-level suicide. Future studies should incorporate gold-standard assessments of SI (O’Connor et al., Citation2013; Posner et al., Citation2011; Ringer et al., Citation2018). Similarly, the reductions in suicide attempts and mortality found in this study were based upon estimates using recently published associations found in large electronic health record databases as opposed to actual data on attempts and mortality linked to the DMHI participants (Rossom et al., Citation2017). Studies that collect longitudinal follow-up data and leverage administrative data linkages are needed to better assess actual rates of suicide attempts and mortality (Kessler & Luedtke, Citation2021). This study used a single measure of intervention engagement to ensure statistical parsimony. Follow-up studies will more thoroughly examine multiple dimensions of engagement (e.g., patient-therapist interactions, cognitive-affective responses) and differential engagement patterns to understand potential effects on SI (Chien et al., Citation2020; Perski et al., Citation2017; Wu et al., Citation2021). Lastly, the DMHI in this study excludes patients with a current suicide plan, limiting generalizability to other DMHIs that may target people actively transitioning from ideation to action (Nock et al., Citation2016).

Suicide is a significant global health concern associated not only with the human cost, but also significant secondary societal and economic impacts. There is a pressing need for novel interventions that reduce barriers to treatment, decrease costs, and increase scalability in order to reduce these suicide outcomes. Thus, this study provides real-world evidence of the potential impact of therapist-supported DMHIs on reducing SI, depression, and anxiety, which has major implications for addressing the mental health crisis in the US.

AUTHOR NOTES

Benjamin W. Nelson, Verily Life Sciences, San Francisco, California, USA. Valerie L. Forman-Hoffman, Woebot Health, San Francisco, California, USA. Nicholas C. Peiper, Department of Epidemiology and Population Health, University of Louisville, Louisville, Kentucky, USA.

Supplemental material

Supplemental Material

Download MS Word (367 KB)

DISCLOSURE STATEMENT

Dr. Peiper reports research support from the National Institute of Mental Health (R44MH126836; PI: Peiper), Lilly Oncology, Gilead Sciences, and Norton Healthcare Foundation. Drs. Peiper and Forman-Hoffman report stock options from Meru Health. The authors were affiliated with Meru Health during the conduct of this study. Drs. Nelson, Forman-Hoffman, and Peiper are now affiliated with Verily Life Sciences, Woebot Health, and University of Louisville, respectively.

Additional information

Funding

Dr. Peiper reports research support from the National Institute of Mental Health (R44MH126836), Lilly Oncology, Gilead Sciences, and Norton Healthcare Foundation. Drs. Peiper and Forman-Hoffman report stock options from Meru Health. Drs. Forman-Hoffman and Nelson were employed by Meru Health during the conduct of this research. They are now with Woebot Health and Verily Life Sciences, respectively.

REFERENCES

  • Ahmedani, B. K., Perron, B., Ilgen, M., Abdon, A., Vaughn, M., & Epperson, M. (2012). Suicide thoughts and attempts and psychiatric treatment utilization: Informing prevention strategies. Psychiatric Services (Washington, D.C.), 63(2), 186–189. doi:10.1176/appi.ps.201100159
  • Andrilla, C. H. A., Patterson, D. G., Garberson, L. A., Coulthard, C., & Larson, E. H. (2018). Geographic variation in the supply of selected behavioral health providers. American Journal of Preventive Medicine, 54(6 Suppl 3), S199–S207. doi:10.1016/j.amepre.2018.01.004
  • Askari, M. S., Mauro, P. M., Kaur, N., & Keyes, K. M. (2023). Age, period, and cohort trends in perceived mental health treatment need and differences by mental health severity in the United States, 2008-2019. Community Mental Health Journal, 59(4), 631–640. doi:10.1007/s10597-022-01044-3
  • Bates, D., Mächler, M., Bolker, B., & Walker, S. (2014). Fitting linear mixed-effects models using lme4. In arXiv [stat.CO]. arXiv. http://arxiv.org/abs/1406.5823.
  • Brooks, M. E., Kristensen, K., van Benthem, K. J., Magnusson, A., Berg, C. W., Skaug, H. J., Mächler, M., & Bolker, B. M. (2017). Modeling zero-inflated count data with glmmTMB. bioRxiv, 132753. doi:10.1101/132753
  • Bruijniks, S. J. E., Lemmens, L. H. J. M., Hollon, S. D., Peeters, F. P. M. L., Cuijpers, P., Arntz, A., … Huibers, M. J. H. (2020). The effects of once- versus twice-weekly sessions on psychotherapy outcomes in depressed patients. The British Journal of Psychiatry : The Journal of Mental Science, 216(4), 222–230. doi:10.1192/bjp.2019.265
  • Büscher, R., Torok, M., Terhorst, Y., & Sander, L. (2020). Internet-Based Cognitive Behavioral Therapy to Reduce Suicidal Ideation: A Systematic Review and Meta-analysis. JAMA Network Open, 3(4), e203933. doi:10.1001/jamanetworkopen.2020.3933
  • Carney, C. E., Edinger, J. D., Kuchibhatla, M., Lachowski, A. M., Bogouslavsky, O., Krystal, A. D., & Shapiro, C. M. (2017). Cognitive Behavioral Insomnia Therapy for Those With Insomnia and Depression: A Randomized Controlled Clinical Trial. Sleep, 40(4), zsx019. doi:10.1093/sleep/zsx019
  • Chien, I., Enrique, A., Palacios, J., Regan, T., Keegan, D., Carter, D., … Belgrave, D. (2020). A Machine Learning Approach to Understanding Patterns of Engagement With Internet-Delivered Mental Health Interventions. JAMA Network Open, 3(7), e2010791. doi:10.1001/jamanetworkopen.2020.10791
  • Dennis, M., Baillon, S., Brugha, T., Lindesay, J., Stewart, R., & Meltzer, H. (2007). The spectrum of suicidal ideation in Great Britain: Comparisons across a 16-74 years age range. Psychological Medicine, 37(6), 795–805. doi:10.1017/S0033291707000013
  • Forman-Hoffman, V. L., Kvedar, J. C., & Ranta, K. (2021). Potential benefits of remote continuous care for depression. International Journal of Digital Health, 1(1), 15. doi:10.29337/ijdh.39
  • Goldberg, S. B., Lam, S. U., Simonsson, O., Torous, J., & Sun, S. (2022). Mobile phone-based interventions for mental health: A systematic meta-review of 14 meta-analyses of randomized controlled trials. PLOS Digital Health, 1(1), e0000002. doi:10.1371/journal.pdig.0000002
  • Goldin, P. R., Lindholm, R., Ranta, K., Hilgert, O., Helteenvuori, T., & Raevuori, A. (2019). Feasibility of a Therapist-Supported, Mobile Phone-Delivered Online Intervention for Depression: Longitudinal Observational Study. JMIR Formative Research, 3(1), e11509. doi:10.2196/11509
  • Hetrick, S. E., Yuen, H. P., Bailey, E., Cox, G. R., Templer, K., Rice, S. M., … Robinson, J. (2017). Internet-based cognitive behavioural therapy for young people with suicide-related behaviour (Reframe-IT): A randomised controlled trial. Evidence-Based Mental Health, 20(3), 76–82. doi:10.1136/eb-2017-102719
  • Hoffman, L., & Stawski, R. S. (2009). Persons as Contexts: Evaluating Between-Person and Within-Person Effects in Longitudinal Analysis. Research in Human Development, 6(2-3), 97–120. doi:10.1080/15427600902911189
  • Holm, S. (1979). A Simple Sequentially Rejective Multiple Test Procedure. Scandinavian Journal of Statistics, Theory and Applications, 6(2), 65–70. http://www.jstor.org/stable/4615733
  • Kassambara, A. (2020). rstatix: Pipe-friendly framework for basic statistical tests. version 0.7. 0. CRAN.
  • Kaur, N., Esie, P., Finsaas, M. C., Mauro, P. M., & Keyes, K. M. (2023). Trends in Racial-Ethnic Disparities in Adult Mental Health Treatment Use From 2005 to 2019. Psychiatric Services (Washington, D.C.), 74(5), 455–462. doi:10.1176/appi.ps.202100700
  • Kessler, R. C., Berglund, P., Borges, G., Nock, M., & Wang, P. S. (2005). Trends in suicide ideation, plans, gestures, and attempts in the United States, 1990-1992 to 2001-2003. JAMA, 293(20), 2487–2495. doi:10.1001/jama.293.20.2487
  • Kessler, R. C., & Luedtke, A. (2021). Pragmatic Precision Psychiatry-A New Direction for Optimizing Treatment Selection. JAMA Psychiatry, 78(12), 1384–1390. doi:10.1001/jamapsychiatry.2021.2500
  • Krieger, T., Bur, O. T., Weber, L., Wolf, M., Berger, T., Watzke, B., & Munder, T. (2023). Human contact in internet-based interventions for depression: A pre-registered replication and meta-analysis of randomized trials. Internet Interventions, 32, 100617. doi:10.1016/j.invent.2023.100617
  • Kroenke, K., Spitzer, R. L., & Williams, J. B. (2001). The PHQ-9: Validity of a brief depression severity measure. Journal of General Internal Medicine, 16(9), 606–613. doi:10.1046/j.1525-1497.2001.016009606.x
  • Kroenke, K., Strine, T. W., Spitzer, R. L., Williams, J. B. W., Berry, J. T., & Mokdad, A. H. (2009). The PHQ-8 as a measure of current depression in the general population. Journal of Affective Disorders, 114(1-3), 163–173. doi:10.1016/j.jad.2008.06.026
  • Lehrer, P., Vaschillo, B., Zucker, T., Graves, J., Katsamanis, M., Aviles, M., & Wamboldt, F. (2013). Protocol for Heart Rate Variability Biofeedback Training. Biofeedback, 41(3), 98–109. doi:10.5298/1081-5937-41.3.08
  • Louzon, S. A., Bossarte, R., McCarthy, J. F., & Katz, I. R. (2016). Does Suicidal Ideation as Measured by the PHQ-9 Predict Suicide Among VA Patients? Psychiatric Services (Washington, D.C.), 67(5), 517–522. doi:10.1176/appi.ps.201500149
  • McCusker, P. J. (2016). The Clinical Utility of PHQ-9 Item 9 for Suicide Prediction. Psychiatric Services (Washington, D.C.), 67(9), 1042–1042. doi:10.1176/appi.ps.670904
  • Morgan, D. (2003). Mindfulness-based cognitive therapy for depression: A new approach to preventing relapse. Psychotherapy Research : journal of the Society for Psychotherapy Research, 13(1), 123–125. doi:10.1080/713869628
  • Moutier, C. (2021). Suicide Prevention in the COVID-19 Era: Transforming Threat Into Opportunity. JAMA Psychiatry, 78(4), 433–438. doi:10.1001/jamapsychiatry.2020.3746
  • Na, P. J., Yaramala, S. R., Kim, J. A., Kim, H., Goes, F. S., Zandi, P. P., … Bobo, W. V. (2018). The PHQ-9 Item 9 based screening for suicide risk: A validation study of the Patient Health Questionnaire (PHQ)-9 Item 9 with the Columbia Suicide Severity Rating Scale (C-SSRS). Journal of Affective Disorders, 232, 34–40. doi:10.1016/j.jad.2018.02.045
  • Nock, M. K., Borges, G., Bromet, E. J., Cha, C. B., Kessler, R. C., & Lee, S. (2008). Suicide and suicidal behavior. Epidemiologic Reviews, 30(1), 133–154. doi:10.1093/epirev/mxn002
  • Nock, M. K., Kessler, R. C., & Franklin, J. C. (2016). Risk factors for suicide ideation differ from those for the transition to suicide attempt: The importance of creativity, rigor, and urgency in suicide research. Clinical Psychology: Science and Practice, 23(1), 31–34. doi:10.1111/cpsp.12133
  • O’Connor, E., Gaynes, B., Burda, B. U., Williams, C., & Whitlock, E. P. (2013). Screening for Suicide Risk in Primary Care: A Systematic Evidence Review for the U.S. Preventive Services Task Force. Agency for Healthcare Research and Quality (US). https://www.ncbi.nlm.nih.gov/pubmed/23678511.
  • Patil, I. (2021). Visualizations with statistical details: The “ggstatsplot” approach. Journal of Open Source Software, 6(61), 3167. doi:10.21105/joss.03167
  • Peiper, N. C., Nelson, B. W., Aschbacher, K., & Forman-Hoffman, V. L. (2023). Trajectories of depression symptoms in a therapist-supported digital mental health intervention: A repeated measures latent profile analysis. Social Psychiatry and Psychiatric Epidemiology, 58(8), 1237–1246. doi:10.1007/s00127-022-02402-y
  • Perski, O., Blandford, A., West, R., & Michie, S. (2017). Conceptualising engagement with digital behaviour change interventions: A systematic review using principles from critical interpretive synthesis. Translational Behavioral Medicine, 7(2), 254–267. doi:10.1007/s13142-016-0453-1
  • Posner, K., Brown, G. K., Stanley, B., Brent, D. A., Yershova, K. V., Oquendo, M. A., … Mann, J. J. (2011). The Columbia-Suicide Severity Rating Scale: Initial validity and internal consistency findings from three multisite studies with adolescents and adults. The American Journal of Psychiatry, 168(12), 1266–1277. doi:10.1176/appi.ajp.2011.10111704
  • Razykov, I., Ziegelstein, R. C., Whooley, M. A., & Thombs, B. D. (2012). The PHQ-9 versus the PHQ-8–is item 9 useful for assessing suicide risk in coronary artery disease patients? Data from the Heart and Soul Study. Journal of Psychosomatic Research, 73(3), 163–168. doi:10.1016/j.jpsychores.2012.06.001
  • Ringer, F. B., Soberay, K. A., Rogers, M. L., Hagan, C. R., Chu, C., Schneider, M., … Joiner, T. E. (2018). Initial validation of brief measures of suicide risk factors: Common data elements used by the Military Suicide Research Consortium. Psychological Assessment, 30(6), 767–778. doi:10.1037/pas0000519
  • Rossom, R. C., Coleman, K. J., Ahmedani, B. K., Beck, A., Johnson, E., Oliver, M., & Simon, G. E. (2017). Suicidal ideation reported on the PHQ9 and risk of suicidal behavior across age groups. Journal of Affective Disorders, 215, 77–84. doi:10.1016/j.jad.2017.03.037
  • Sander, L., Gerhardinger, K., Bailey, E., Robinson, J., Lin, J., Cuijpers, P., & Mühlmann, C. (2020). Suicide risk management in research on internet-based interventions for depression: A synthesis of the current state and recommendations for future research. Journal of Affective Disorders, 263, 676–683. doi:10.1016/j.jad.2019.11.045
  • Sarris, J., Logan, A. C., Akbaraly, T. N., Amminger, G. P., Balanzá-Martínez, V., Freeman, M. P., … Jacka, F. N. (2015). Nutritional medicine as mainstream in psychiatry. The Lancet. Psychiatry, 2(3), 271–274. doi:10.1016/S2215-0366(14)00051-0
  • Simon, G. E., Coleman, K. J., Rossom, R. C., Beck, A., Oliver, M., Johnson, E., … Rutter, C. (2016). Risk of suicide attempt and suicide death following completion of the Patient Health Questionnaire depression module in community practice. The Journal of Clinical Psychiatry, 77(2), 221–227. doi:10.4088/JCP.15m09776
  • Simon, G. E., Rutter, C. M., Peterson, D., Oliver, M., Whiteside, U., Operskalski, B., & Ludman, E. J. (2013). Does response on the PHQ-9 Depression Questionnaire predict subsequent suicide attempt or suicide death? Psychiatric Services (Washington, D.C.), 64(12), 1195–1202. doi:10.1176/appi.ps.201200587
  • Stone, D. M., Jones, C. M., & Mack, K. A. (2021). Changes in Suicide Rates - United States, 2018-2019. MMWR. Morbidity and Mortality Weekly Report, 70(8), 261–268. doi:10.15585/mmwr.mm7008a1
  • Tierney, N. J., & Cook, D. H. (2018). Expanding tidy data principles to facilitate missing data exploration, visualization and assessment of imputations. In arXiv [stat.CO]. arXiv doi:10.18637/jss.v000.i00
  • Torok, M., Han, J., Baker, S., Werner-Seidler, A., Wong, I., Larsen, M. E., & Christensen, H. (2020). Suicide prevention using self-guided digital interventions: A systematic review and meta-analysis of randomised controlled trials. The Lancet Digital Health, 2(1), e25–e36. doi:10.1016/s2589-7500(19)30199-2
  • Wu, A., Scult, M. A., Barnes, E. D., Betancourt, J. A., Falk, A., & Gunning, F. M. (2021). Smartphone apps for depression and anxiety: A systematic review and meta-analysis of techniques to increase engagement. Npj Digital Medicine, 4(1), 20. doi:10.1038/s41746-021-00386-8
  • van Buuren, S., & Groothuis-Oudshoorn, K. (2011). mice: Multivariate Imputation by Chained Equations in R. Journal of Statistical Software, 45(3), 1-67. doi:10.18637/jss.v045.i03
  • Yarborough, B. J. H., Stumbo, S. P., Ahmedani, B., Rossom, R., Coleman, K., Boggs, J. M., & Simon, G. E. (2021). Suicide Behavior Following PHQ-9 Screening Among Individuals With Substance Use Disorders. Journal of Addiction Medicine, 15(1), 55–60. doi:10.1097/ADM.0000000000000696