Skip to main content

Suicide risk among veterans using VHA justice-involved services: a latent class analysis



Justice-involved Veterans experience notable risk for psychosocial stressors (e.g., homelessness) and psychiatric multimorbidity, which can result in complex clinical presentations. However, research examining how such factors coalesce to impact risk for suicide remains limited.


We conducted a latent class analysis of 180,454 Veterans accessing Veterans Health Administration (VHA) justice-related services from 2005 to 2018.


A four-model class membership solution was identified. Among these classes, risk for suicide was highest among Veterans with greater psychiatric burden, with risk most notable among those with high VA service use. Veterans seeking healthcare primarily focused on substance use disorders or with low psychiatric burden and service use had a lower risk for suicide.


Psychiatric multimorbidity is salient as it relates to suicide among Veterans accessing VHA justice-related services. Further evaluation of existing VHA services for this population and methods of augmenting and enhancing care for justice-involved Veterans with histories of co-occurring psychiatric conditions may be beneficial in facilitating suicide prevention efforts.

Peer Review reports


According to the most recent Bureau of Justice Statistics report, Veterans comprise approximately 8% of individuals currently in prison within the United States [1]. Moreover, a sizable proportion of Veterans are currently on parole or probation, suggesting the percentage of Veterans with any interaction with the criminal justice system is likely much higher [2].

While rates of criminal justice system involvement among Veterans do not appear to be elevated relative to non-Veteran adults, justice-involved Veterans differ from justice-involved non-Veterans in several ways. In particular, justice-involved Veterans have higher documented rates of violent offenses compared to justice-involved non-Veteran adults [1, 3,4,5]. For example, in a recent report by Department of Justice [1], they noted that male Veterans in state (69.3%) and federal (24.7%) prisons had significantly higher rates of violent crime offenses relative to male non-Veterans (state: 56.7%; federal: 12.6%). Such offenses result in lengthier sentencing, societal stigma, and difficulty transitioning back into the community following incarceration [6,7,8].

Justice-involved Veterans also often experience psychiatric conditions and symptoms (e.g., substance use, serious mental illness, posttraumatic stress disorder). Data suggest that justice-involved Veterans have higher rates of several conditions and symptoms than both non-Veterans involved with the justice system and Veterans without a history of criminal justice involvement [9,10,11,12,13]. In addition, justice-involved Veterans experience psychosocial stressors (e.g., housing and occupational difficulties) with high propensity, which can further impede functioning [14]. Decreased psychosocial functioning (e.g., unemployment), coupled with stigma related to legal involvement, can impact perceptions of oneself, the world, and the future, which may subsequently influence risk for suicide [6, 15].

To meet the needs of justice-involved Veterans, the Veterans Health Administration (VHA) has developed a number of programs specifically tailored for this population [16]. For example, Veterans Justice Outreach (VJO) and Health Care for Reentry Veterans programs focus on connecting and engaging justice-involved Veterans in health and social services. Although these programs often prove beneficial to recovery, the Veteran population that accesses VJO and Health Care for Reentry Veterans services often have a high concentration of suicide risk factors [17]. Potential factors among Veterans accessing these services that may contribute to elevated suicide risk include psychiatric multimorbidity and stressors (e.g., housing instability). Such risk factors may result in clinical complexity and difficulties engaging Veterans in healthcare, thus impacting ability to participate and benefit from intervention. In fact, when psychiatric conditions and housing instability were accounted for in a prior examination of suicide risk among Veterans accessing VHA justice-related services, the effect of justice involvement substantially attenuated and, at times, was no longer significant [17].

Given the clinical complexity and presence of multiple drivers of suicide risk among members of this population, multivariate analytic approaches capable of identifying groups at risk is warranted. Specifically, latent class analytic approaches have demonstrated the ability to identify phenotypes of individuals based on a number of factors (e.g., psychiatric diagnosis, trauma exposure, service use). As such, our team combined multiple national datasets to identify Veterans accessing VHA justice-related services who did and did not die by suicide. A latent class analytic approach was used to understand how health comorbidity, psychiatric conditions, and VHA service use coalesced to generate latent classes among the VHA-using justice-involved Veteran population. Such classes were then examined as they related to suicide mortality to inform which Veterans were at greatest risk for suicide.


Data source and participants

The current study was conducted using data collected via a retrospective chart review of Veterans’ VHA electronic health records (EHR). We identified justice-involved Veterans based on any use of VHA justice-related services (i.e., Veterans Justice Outreach; Healthcare for Re-entry Veterans)Footnote 1 between 1/2005 and 12/2018.Footnote 2

For this cohort, we obtained data regarding sociodemographics (i.e., age, race, ethnicity, rurality), service connection, psychiatric diagnoses (i.e., posttraumatic stress disorder [PTSD], adjustment disorder, anxiety disorder, depressive disorder, bipolar disorder, alcohol use disorder, substance use disorder, personality disorder; schizophrenia, and other psychotic disorder), documented history of traumatic brain injury (TBI), military sexual trauma (MST; from the universal VA MST screen), and medical comorbidity (i.e., Charlson Comorbidity Index [18]). Diagnoses were identified using ICD codes (VHA EHR).Footnote 3 In addition, we accessed data regarding use of VHA mental health care, overall VHA care, and homeless services.Footnote 4 Mental health care, as well as VHA use (not including mental health or justice-related services), were identified using stop and bed codes and categorized based on tertile values for the distribution of frequency of unique visits over the study period: mental health care - low (≤ 30, moderate (> 30 and ≤ 141), or high (> 141); and overall use - low (≤ 65), moderate (> 65 and ≤ 185), or high (> 185) use. Due to high levels of homeless service engagement among justice-involved Veterans, VHA homeless service use was also included [14], which we examined as a dichotomous variable (yes/no).

VHA EHR data were combined with data from the U.S. Veterans Eligibility Trends and Statistics (USVETS) and VA/DoD Mortality Data Repository (MDR) National Death Index data. Sociodemographic data from USVETS was cross-referenced with VHA EHR data to enhance the validity of the dataset. MDR was used to determine mortality status within our cohort (who was alive versus had died, including by suicide or other causes).

In total, we identified 197,737 justice-involved Veterans with available data for the current examination. Given the focus on suicide mortality, we removed 17,283 Veterans classified by MDR as having died by a non-suicide-related cause. This resulted in a final cohort for the latent class analysis of 180,454 Veterans, 822 (0.5%) of whom died by suicide from 2005 to 2018. The current study was approved by the Colorado Multiple Institutional Review Board.

Analytic Plan

Class enumeration

Class enumeration was performed using Mplus v.8.5 [19]. We used latent class analysis to identify subgroups of justice-involved Veterans based on psychiatric diagnosis, documented TBI, MST screening results, medical comorbidity, and VHA service use (i.e., overall, mental health, and homeless service use). Latent class analysis is a mixture modelling approach which uses categorical observed variables to identify unobserved, mutually exclusive subgroups within a sample [20]. Response patterns for observed variables are used to classify individuals into mutually exclusive latent classes [21]. We fit models with one through eight latent classes, relying on information criteria (e.g., AIC; BIC), the Lo-Mendell-Rubin (LMR) adjusted likelihood ratio test, and clinical interpretability to select a parsimonious, statistically sound, and theoretically plausible model [20, 22]. Conditional item probabilities were used to characterize and label the classes. Such probabilities represent the likelihood of experiencing the observed trait (e.g., PTSD) assuming membership in the unobserved (i.e., latent) class. Upon selecting the model with the optimal class solution, we conducted a descriptive analysis of covariates among the total sample and stratified by class membership. A Veteran’s class membership was considered the class for which the Veteran had the highest likelihood of membership relative to other classes (i.e., posterior membership probabilities).

Latent class membership as a predictor of suicide

We then examined the indicator of Veterans’ class membership, identified using posterior membership probabilities, as a predictor of death by suicide. Unadjusted logistic regression was used to model suicide (yes/no) as a function of latent class membership to determine if the risk for suicide differs across classes (note that given the rare outcome, the odds ratio is a reasonable estimate of the relative risk and we therefore interpret results as such). Subsequently, we fit an adjusted logistic regression model, where suicide was modeled as a function of latent class membership, age (i.e., 18–39 years, 40–49 years, 50–59 years, 60–69 years, ≥ 70 years), sex, race (i.e., White, Black, non-Black underrepresented racial group [comprised of American Indian/Alaska Native, Asian, Native Hawaiian/Pacific islander, and “Other” given relatively small sample sizes for each]), ethnicity, rurality, and service connection (yes/no). Statistical significance of all parameter estimates was evaluated at α = 0.05. Additionally, due to the size of the sample, both 95% and 99% CI are reported.


Table 1 includes sample characteristics of the cohort. Overall, the cohort was largely male (93.9%), White (64.4%), and non-Hispanic (91.4%). Psychiatric diagnoses were highly prevalent within our sample, with substance use disorder (65.2%), alcohol use disorder (63.1%), and depressive disorders (63.6%) the most frequent. In addition, documented medical comorbidity was low for most of the sample (76.6%). More than half of justice-involved Veterans in our sample had a history of VHA homeless service use (62.0%).

Table 1 Sample Characteristics of Veterans Using VHA Justice-Related Services (N = 180,454)

Class enumeration

The LMR adjusted likelihood ratio test was not statistically significant for the three-class solution (p = 0.333), indicating that adding a third class did not improve model fit to a statistically significant degree compared to the two-class solution. This provided statistical evidence in favor of the two-class solution. However, examination of information criteria (e.g., Bayesian Information Criterion) indicated that the addition of classes improved model fit, with incremental improvement beginning to diminish upon adding a fourth class (see Supplementary Table 1). Thus, we considered classes three through five as candidate class solutions. Upon carefully inspecting the conditional item probabilities for classes three through five, we chose the four-class solution because it was the most parsimonious model and included classes that were both qualitatively distinct and clinically relevant. In addition, conditional probabilities for the four-class solution were characterized by superior class separation and within-class homogeneity relative to the three- and five-class solutions. Further, average probability for most likely class membership was high across classes, ranging from 0.76 to 0.93. See Table 2 for conditional item probabilities for the four-class solution.

Table 2 Conditional Item Probabilities Across Latent Classes for the Four-Class Solution

The High psychiatric burdenFootnote 5/high service user (30.27%) group was characterized by a high probability of PTSD, adjustment disorder, depression, anxiety, alcohol use disorder, and other substance use disorder. Veterans in this class were likely to receive high levels of mental health (i.e., > 141 visits) and overall VHA (i.e., > 185 visits) services and were highly likely to receive VHA homeless services. This class was distinguished from the High psychiatric burden/moderate service user group, which was also prevalent (29.79%), and carried a similar psychiatric burden but were most likely to receive moderate levels of mental health services (i.e., > 30 and  141 visits) and VHA (i.e., > 65 and  185 visits) services. Additionally, this class had a moderate probability (0.502) of receiving VHA homeless services. The Low psychiatric burden/low service user (25.53%) group was characterized by a low likelihood of all psychiatric conditions, and members were likely to receive low levels of mental health (i.e.,  30 visits) and VHA (i.e.,  65 visits) services. Members of this class were unlikely to receive homeless services. Finally, the Primarily substance use disorder (13.3%) group had a high probability of alcohol use disorder and substance use disorder, but not other psychiatric conditions, and were most likely to receive moderate levels of mental health and VHA services. Additionally, they were highly likely to receive homeless services.

Latent class membership as predictor of suicide

A crude logistic regression model was conducted to determine the relationship between latent class membership and suicide (see Table 3). The Low psychiatric burden/low service user group was chosen as the reference given prior research noting the association between psychiatric diagnosis and suicide risk (see Supplementary Table 2 for comparisons with other latent classes as the reference group).

Table 3 Crude and Adjusted Models Examining Class Membership as it Relates to Suicide Risk

Overall, the crude model was significant (χ2(3) = 22.82, p < 0.001, Nagelkerke R2 = 0.002). The groups experiencing the greatest risk for suicide were the High psychiatric burden/moderate service user (1.29 times greater risk for suicide relative to the Low psychiatric burden/low service user group, p = 0.0057) and the High psychiatric burden/high service user groups (1.22 times greater risk for suicide relative to the Low psychiatric burden/low service user group, p = 0.038). Interestingly, in the crude model, Veterans in the Primarily substance use-related group had lower (i.e., 0.77 times) risk for suicide relative to Veterans with low psychiatric burden and service use (p = 0.049).

The model was subsequently adjusted for demographic factors (i.e., age, sex, race, ethnicity, rurality) and service connection (Table 3). Again, the overall model was significant (χ2(12) = 128.95, p < 0.001, Nagelkerke R2 = 0.036). Notably, with the inclusion of these factors, risk appeared to increase largely within the high psychiatric burden and high service use group. Specifically, Veterans in the High psychiatric burden/high service user group had 1.88 times greater risk for suicide compared to the Low psychiatric burden/low service user group (p = 0.001). Similarly, Veterans in the High psychiatric burden/moderate service user group had 1.47 times greater risk for suicide relative to the Low psychiatric burden/low service user group (p = 0.046). However, differing from the crude model, there was no significant difference in risk for suicide between the Primarily substance use-related and Low psychiatric burden/low service user groups (p > 0.05).


Building upon prior work noting risk for suicide among justice-involved Veterans [15, 17], this study provides the first latent class analysis of factors associated with risk for suicide among Veterans using VHA justice-involved services. Four latent classes were identified, with risk most pronounced in adjusted models among justice-involved Veterans experiencing high psychiatric burden and moderate to high VHA service use.

Such risk is consistent with the notion that psychiatric clinical complexity may indeed be a primary factor in driving suicide risk among justice-involved Veterans. Given the psychiatric multimorbidity and stressors that justice-involved Veterans experience, these Veterans likely present to time-limited care settings with myriad psychosocial concerns. As such, it may be difficult for VHA providers working with these Veterans to balance provision of services, including suicide prevention care [14]. For example, intervention may focus on a current stressor (e.g., upcoming eviction) rather than providing an evidence-based mental health treatment (e.g., cognitive behavioral therapy for depression). While addressing the most pressing stressor is often a critical first component of treatment, especially if it is the Veteran’s preference, additional assessment of psychiatric need could help address underlying diagnoses and inform future treatment planning [14]. Given this, identifying optimal service delivery that can reduce patterns of psychosocial stress (e.g., housing or unemployment) as well as psychiatric symptoms among this population of Veterans is needed.

Despite this, the relationship between service delivery and suicide risk in this cohort of Veterans was complex. Specifically, no significant differences in risk for suicide were found between moderate and high service groups of Veterans experiencing high psychiatric burden (see Supplementary Table 2). A potential explanation for this is that type of VHA service, rather than frequency of use, may be more relevant to understanding suicide risk. Indeed, the VHA services used in this sample could have varied substantially in potential to address and prevent suicide risk (e.g., evidence-based psychotherapy for a psychiatric diagnosis). Moreover, we generated a total number of appointments accessed in the study time frame, but this number was broad and not indicative of treatment adherence nor engagement. Therefore, an initial step to understanding the utility of existing programing with respect to suicide risk is likely a mixed-method investigation of extant VHA services accessed by justice-involved Veterans. In particular, further understanding of specific VHA services which can reduce subsequent risk for suicide, as well as Veteran reported barriers and facilitators to engaging in those services, is likely necessary to ensure that these Veterans receive optimal suicide prevention efforts.

Interestingly, risk did not substantially differ between justice-involved Veterans in the Low psychiatric burden/low service use and Primarily substance use-related groups. Several potential explanations exist for this. First and foremost, it is possible that existing VHA substance use programing may be beneficial for this subset of Veterans in decreasing risk. Indeed, Ilgen and colleagues [23] noted that Veterans who completed substance use treatment had decreased subsequent risk for suicide. In addition, it is also possible that other psychiatric diagnoses germane to the high psychiatric burden groups (e.g., depression, PTSD) may have been a larger driver of suicide risk than substance use-related diagnoses. Nonetheless, additional research remains necessary to elucidate drivers of risk within each of our identified groups.

VHA homeless service use was associated with group membership in both the high psychiatric burden and primarily substance use-related groups. Such an association makes sense given prior research documenting the relationship between mental health conditions and first-time as well as chronic homelessness [24]. However, given differing suicide risk between these groups, additional research is likely needed to determine how housing instability may differentially impact those with more psychiatric burden.

Finally, medical comorbidity (Charlson Comorbidity Index [18]), TBI, and MST did not emerge as factors distinguishing latent classes. This finding is surprising given that these factors are associated with psychiatric conditions, service use, and suicide [25,26,27,28]. Such results may in part be related to medical and other comorbidities not being identified among members of this cohort. Additional research is warranted and may benefit from efforts to examine severity of TBI (e.g., mild, moderate, severe) and type of MST (e.g., sexual harassment or sexual assault), as well as whether these factors differentially impact justice-involved Veterans’ suicide risk based on gender (e.g., association of MST in women versus men).

Despite the strength of findings, results should be interpreted within the context of limitations. First, this examination was specific to Veterans accessing VHA justice-related services. Because of this, findings may not generalize to Veterans who do not qualify for comprehensive VHA care or who interact with the criminal justice system and who do not access VHA justice-related services. Moreover, a sizable portion of our sample was moderately to highly engaged in VHA services, including mental health care. As such, many of these Veterans may have been routinely accessing services which focused on addressing psychiatric conditions as well as which assessed the presence of suicide risk factors. In a similar vein, Veterans accessing mental health services would be more likely to have a psychiatric diagnosis conferred on their EHR.

We relied on the VHA EHR, rather than conducting diagnostic interviews. As such, we are unable to determine the accuracy of diagnoses. In addition, we examined substance use disorder diagnoses as a unidimensional construct. This approach limits results given differing risk for suicide based on type of substance (e.g., opioid vs. cocaine [29]), as well as factors which may have moderated suicide risk (e.g., early vs. sustained remission).

An additional consideration is that we were unable to employ contemporary approaches to predicting distal outcomes while accounting for measurement error in latent class membership. Specifically, while we attempted to use the 3-step BCH method to examine the relation between latent class membership and death by suicide [30], we were unable to achieve model convergence, likely due to the low base rate of suicide mortality. Given this, we instead examined the relationship between class membership and suicide mortality using binary logistic regression, which may have minimized understanding of the contribution of important covariates (i.e., age, sex, race, ethnicity, rurality, service connection). It’s worth noting that the latent classes associated with the interpreted model exhibited good measurement properties, including acceptable average posterior probabilities among those assigned to each class (0.76-0.93) which exceeded the acceptable range of 0.70 [20], limiting concerns regarding measurement error biasing the results. Nonetheless, future research, inclusive of a larger suicide mortality sample, should consider examining the effect of latent class membership upon suicide among justice-involved Veterans while utilizing approaches capable of accounting for latent class measurement error.

Finally, we used EHR codes to identify use of VHA justice-related and homeless services, rather than using Homelessness Operations and Management and Evaluation System (HOMES) data [31]. Given that HOMES is an administrative database that tracks use of VA homeless services separate from the EHR, some justice-involved Veterans may not have been captured in our cohort.

Justice-involved Veterans experience several risk factors for suicide (e.g., psychiatric diagnoses, homelessness). We provided the first latent class analysis of suicide mortality among members of this population. We found that psychiatric burden and moderate to high service use are factors associated with risk. Investigation of current and novel programing for this population may be a potential method of enhancing suicide prevention efforts for this at-risk population of Veterans.

Data availability

All data are available upon request contingent upon IRB and VA guidelines and approval.


  1. Justice-involved services were determined based on outpatient stop codes (591, 592) mirroring the method from Palframan and colleagues [17].

  2. Timeframe was determined based on data completeness (U.S Veterans Eligibility Trends and Statistics complete from 2005 onward) and availability (National Death Index data available through 2018) at the time of analysis.

  3. Health conditions were confirmed based on presence of 1 inpatient or 2 outpatient documented encounters.

  4. Homeless services were determined based on ICD codes (ICD-9: V60.0; ICD-10: Z59.0), outpatient stop codes (501, 504, 507, 508, 511, 522, 528, 529, 530, 590), and inpatient specialty codes (28, 29, 37, and 39).

  5. “Psychiatric burden” is indicative of risk for psychiatric diagnoses (i.e., posttraumatic stress disorder [PTSD], adjustment disorder, anxiety disorder, depressive disorder, bipolar disorder, alcohol use disorder, substance use disorder, personality disorder; schizophrenia, and other psychotic disorder).


  1. Department of Justice. Survey of Prison Inmates., 2016: Veterans in Prison. Retrieved from; 2021.

  2. Glaze LE. Correctional populations in the United States, 2010. Retrieved from; 2011.

  3. Marshall AD, Panuzio J, Taft CT. Intimate partner violence among military veterans and active duty servicemen. Clin Psychol Rev. 2005;25:862–76.

    Article  PubMed  Google Scholar 

  4. Mumola CJ. Veterans in prison or Jail. Bureau of Justice Statistics; 2000.

  5. Noona ME, Mumola CJ. Veterans in State and Federal Prison, 2004. Department of Justice; 2007.

  6. Desai A, Holliday R, Borges LM, Blue-Howells J, Clark S, Stimmel M, Wortzel HS. Facilitating successful reentry among justice-involved veterans: the role of veteran and offender identity. J Psychiatr Pract. 2021;27:52–60.

    Article  PubMed  Google Scholar 

  7. Tewksbury R, Lees M. Perceptions of sex offender registration: collateral consequences and community experiences. Sociol Spectr. 2006;26(3):309–34.

    Article  Google Scholar 

  8. Travis J, Western B, Redburn FS. The growth of incarceration in the United States: exploring causes and consequences. National Academies Press; 2014.

  9. Blodgett JC, Avoundjian T, Finlay AK, Rosenthal J, Asch SM, Maisel NC, Midboe AM. Prevalence of mental health disorders among justice-involved veterans. Epidemiol Rev. 2015;37:163–76.

    Article  PubMed  Google Scholar 

  10. Finlay AK, Owens MD, Taylor E, Nash A, Capdarest-Arest N, et al. A scoping review of military veterans involved in the criminal justice system and their health and healthcare. Health Justice. 2019;7:6.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Holliday R, Forster JE, Desai A, Miller C, Monteith LL, Schneiderman AI, Hoffmire CA. Association of lifetime homelessness and justice involvement with psychiatric symptoms, suicidal ideation, and suicide attempt among post-9/11 veterans. J Psychiatr Res. 2021;144:455–61.

    Article  PubMed  Google Scholar 

  12. Holliday R, Hoffmire CA, Martin WB, Hoff RA, Monteith LL. Associations between justice involvement and PTSD and depressive symptoms, suicidal ideation, and suicide attempt among post-9/11 veterans. Psychol Trauma. 2021;13:730–9.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Taylor EN, Timko C, Nash A, Owens MD, Harris AHS, Finlay AK. Posttraumatic stress disorder and justice involvement among military veterans: a systematic review and meta-analysis. J Trauma Stress. 2020;33:804–12.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Holliday R, Desai A, Stimmel M, Liu S, Monteith LL, Stewart KE. Meeting the health and social service needs of veterans who interact with the criminal justice system and experience homelessness: a holistic conceptualization and recommendations for tailoring care. Curr Treat Options Psychiatry. 2022;9:174–85.

    Article  Google Scholar 

  15. Holliday R, Martin WB, Monteith LL, Clark SC, LePage JP. Suicide among justice-involved veterans: a brief overview of extant research, theoretical conceptualization, and recommendations for future research. J Soc Distress Homeless. 2020;30:41–9.

    Article  Google Scholar 

  16. Holliday R, Desai A, Gerard G, Liu S, Stimmel M. Understanding the intersection of homelessness and justice involvement: enhancing veteran suicide prevention through VA programming. Fed Pract. 2022;39:8–11.

    PubMed  PubMed Central  Google Scholar 

  17. Palframan KM, Blue-Howells J, Clark SC, McCarthy JF. Veterans Justice Programs: assessing population risks for suicide deaths and attempts. Suicide Life Threat Behav. 2020;50:792–804.

    Article  PubMed  Google Scholar 

  18. Stavem K, Hoel H, Skjaker SA, Haagensen R. Charlson comorbidity index derived from chart review or administrative data: Agreement and prediction of mortality in intensive care patients. Clin Epidemiol. 2017;9:311–20.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Muthén LK, Muthén BO. Mplus User’s Guide.Eighth Edition. Muthén & Muthén;1998–2017.

  20. Masyn KE. Latent class analysis and finite mixture modeling. In: Little TD, editor. The Oxford Handbook of quantitative methods in psychology (volume 2). Oxford Handbooks Online; 2013. pp. 1–120.

  21. Collins LM, Lanza ST. Latent class and latent transition analysis: with applications in the social behavioral, and Health Sciences. Wiley; 2010.

  22. Muthén B. Statistical and substantive checking in growth mixture modeling: comment on Bauer and Curran (2003). 2003.

  23. Ilgen MA, Harris AHS, Moos RH, Tiet QQ. Predictors of a suicide attempt one year after entry into substance use disorder treatment. Alcohol Clin Exp Res. 2007;31:635–42.

    PubMed  Google Scholar 

  24. Creech SK, Johnson E, Borgia M, Bourgault C, Redihan S, O’Toole TP. Identifying mental and physical health correlates of homelessness among first-time and chronically homeless veterans. J Community Psychol. 2015;43:619–27.

    Article  Google Scholar 

  25. Bahraini NH, Simpson GK, Brenner LA, Hoffberg AS, Schneider AL. Suicidal ideation and behaviours after traumatic brain injury: a systematic review. Brain Impair. 2013;14:92–112.

    Article  Google Scholar 

  26. Hammond FM, Corrigan JD, Ketchum JM, Malex JF, Dams-O’Conner K, Hart T, et al. Prevalence of medical and psychiatric comorbidities following traumatic brain injury. J Head Trauma Rehabil. 2019;34:E1–10.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Kimerling R, Gima K, Smith MW, Street A, Frayne S. The Veterans Health Administration and military sexual trauma. Am J Public Health. 2007;97:2160–6.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Kimerling R, Makin-Byrd K, Louzon S, Ignacio RV, McCarthy. JF. Military sexual trauma and suicide mortality. Am J Prev Med. 2016;50:684 – 91.

  29. Bohnert KM, Ilgen MA, Louzon S, McCarthy JF, Katz IR. Substance use disorders and the risk of suicide mortality among men and women in the US Veterans Health Administration. Addiction. 2017;112(7):1193–201.

    Article  PubMed  Google Scholar 

  30. Nylund-Gibson K, Grimm RP, Masyn KE. Prediction from latent classes: a demonstration of different approaches to include distal outcomes in mixture models. Struct Equ Modeling. 2019;26(6):967–85.

    Article  Google Scholar 

  31. Tsai J, Kasprow WJ, Rosnheck RA. Latent homeless risk profiles of a national sample of homeless veterans and their relation to program referral and admission patterns. Am J Public Health. 2013;103:219–47.

    Article  Google Scholar 

Download references


The views expressed are those of the authors and do not necessarily represent the views or policy of the VA or the United States Government.


This work was funded by an intramural grant provided by the National Center on Homelessness among Veterans (PI: Holliday). Additionally, this work was supported in part by the Department of Veterans Affairs (VA); VA National Center on Homelessness among Veterans; VHA Homeless Programs Office; VA Office of Mental Health and Suicide Prevention; and VA Rocky Mountain Mental Illness Research, Education and Clinical Center for Suicide Prevention.

Author information

Authors and Affiliations



R. H., L. M., and L. B. conceptualized the framework for this project and were Investigators on the funded grant. R.H. wrote the main manuscript text. A. K. and (A) S. facilitated with conducting statistical analyses. J. F., M. S., S. C., S. L., L. M., and L. (B) provided editorial feedback.

Corresponding author

Correspondence to Ryan Holliday.

Ethics declarations

Ethics approval and consent to participate

All methods were carried out in accordance with extant guidelines and regulations. Study procedures were approved by the Colorado Multiple Institutional Review Board. As the study utilized archival electronic medical record data, a waiver of informed consent was requested and approved for this study.

Consent for publication

Not applicable.

Competing interests

Dr. Brenner reports grants from the VA, Department of Defense (DoD), National Institutes of Health (NIH), and the State of Colorado, editorial renumeration from Wolters Kluwer, and royalties from the American Psychological Association and Oxford University Press. In addition, she consults with sports leagues via her university affiliation. Drs. Holliday, Kinney, Forster, and Monteith report grant funding from VA. Drs. Holliday and Forster also disclose grant funding from DoD and NIH. Dr. Monteith reports grant funding from DoD. The authors have no additional conflicts of interest to disclose.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Holliday, R., Kinney, A.R., Smith, A.A. et al. Suicide risk among veterans using VHA justice-involved services: a latent class analysis. BMC Psychiatry 23, 235 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Justice involvement
  • Veteran
  • Suicide
  • Latent class analysis