vlog

[Skip to Navigation]
Sign In
Figure 1. Cohort Selection

DMDC indicates Defense Manpower Data Center; USVETS, US Veterans Eligibility Trends and Statistics.

Figure 2. Crude Hazard Rates of Suicide Among Theater, Nontheater, and All Veterans of Vietnam War Era and US Male Population by Number of Years Since Entry to Follow-Up or Calendar Year

US male population suicide rates are based on calendar year of death (obtained from CDC Wonder). Veteran suicide rates were calculated for one-year intervals from 1979 to 2019. The first year interval includes 5.5 months to 1.5 years (1 year). All subsequent rates are for time intervals that begin at midpoint and end at midpoint of the subsequent year (eg, 2 years = 1.6 years − 2.5 years). Accordingly, the separate time intervals for veterans, 1 to 41 years, cross calendar years. Crude rates for both US males and veterans are per 100 000 persons.

Table 1. Demographic and Military Service Characteristics by Theater and Nontheater Status
Table 2. Results of Cox Proportional Hazards Regression Models of Suicide Mortality Risks Associated With Theater Status
1.
Centers for Disease Control and Prevention. Preventing suicide. 2022. Accessed July 7, 2023.
2.
Curtin SC, Hedegaard H, Ahmad FB. Provisional numbers and rates of suicide by month and demographic characteristics: United States, 2020: Vital Statistics Rapid Release, no 16. National Center for Health Statistics. November 2021. Accessed July 12, 2023.
3.
US Department of Veterans Affairs, Office of Mental Health and Suicide Prevention. Suicide among veterans and other Americans 2001-2014. 2016. Accessed July 20, 2023.
4.
US Department of Veterans Affairs. 2022 National Veteran Suicide Prevention annual report. Accessed July 12, 2023.
5.
Anderson R. Vietnam legacy: veterans’ suicide toll may top war casualties. Seattle Times. March 18, 1981.
6.
Langone J. The war that has no ending. پDZ. 1985;(June):44-54.
7.
Pollock DA, Rhodes P, Boyle CA, Decoufle P, McGee DL. Estimating the number of suicides among Vietnam veterans. Am J Psychiatry. 1990;147(6):772-776. doi:
8.
Centers for Disease Control Vietnam Experience Study. Postservice mortality among Vietnam veterans. Ѵ. 1987;257(6):790-795. doi:
9.
Boehmer TK, Flanders WD, McGeehin MA, Boyle C, Barrett DH. Postservice mortality in Vietnam veterans: 30-year follow-up. Arch Intern Med. 2004;164(17):1908-1916. doi:
10.
Watanabe KK, Kang HK, Thomas TL. Mortality among Vietnam veterans: with methodological considerations. J Occup Med. 1991;33(7):780-785. doi:
11.
Watanabe KK, Kang HK. Military service in Vietnam and the risk of death from trauma and selected cancers. Ann Epidemiol. 1995;5(5):407-412. doi:
12.
Anderson H, Hanrahan LP, Jensen M, Laurin D, Yick W, Weigman P. Wisconsin Vietnam Veteran Mortality. Wisconsin Division of Health; 1986.
13.
Kogan MD, Clapp RW. Mortality Among Vietnam Veterans in Massachusetts, 1972-1983. Massachusetts Department of Public Health, Division of Health Statistics; 1985.
14.
Lawrence CE, Reilly AA, Quickenton P, Greenwald P, Page WF, Kuntz AJ. Mortality patterns of New York State Vietnam Veterans. Am J Public Health. 1985;75(3):277-279. doi:
15.
Forsberg CW, Estrada SA, Baraff A, et al. Risk factors for suicide in the Vietnam-era twin registry. Suicide Life Threat Behav. 2022;52(4):631-641. doi:
16.
Kang HK, Cypel Y, Kilbourne AM, et al. HealthViEWS: mortality study of female US Vietnam era veterans, 1965-2010. Am J Epidemiol. 2014;179(6):721-730. doi:
17.
Bullman TA, Kang HK. Posttraumatic stress disorder and the risk of traumatic deaths among Vietnam veterans. J Nerv Ment Dis. 1994;182(11):604-610. doi:
18.
Bullman TA, Kang HK. The risk of suicide among wounded Vietnam veterans. Am J Public Health. 1996;86(5):662-667. doi:
19.
US Department of Veterans Affairs, Office of Mental Health and Suicide Prevention. 2020 National Veteran Suicide Prevention annual report. Accessed July 20, 2023.
20.
US Department of Veterans Affairs, Office of Mental Health and Suicide Prevention. 2021 National Veteran Suicide Prevention annual report. Accessed July 26, 2023.
21.
VHA Data Portal. United States Veterans Eligibility Trends and Statistics (USVETS). Accessed July 19, 2023.
22.
VA Information Resource Center (VIReC). United States Veterans Eligibility Trends and Statistics (USVETS): a new data source with socioeconomic variables. 2019. Accessed July 18, 2023.
23.
Title 10–Armed Forces: this title was enacted by act Aug. 10, 1956, ch. 1041, § 1, 70A Stat. 1. Accessed July 25, 2023.
24.
Defense Manpower Data Center. Vietnam File. July 1999. Accessed July 14, 2023.
25.
US Department of Veterans Affairs. Using data to prevent veteran mortality. 2022. Accessed July 10, 2023.
26.
Hoffmire CABS, Barth SK, Bossarte RM. Reevaluating suicide mortality for veterans with data from the VA-DoD mortality data repository, 2000-2010. Psychiatr Serv. 2020;71(6):612-615. doi:
27.
World Health Organization. Manual of the International Statistical Classification of Diseases: Inquiries, and Causes of Death, Ninth Revision. World Health Organization; 1977.
28.
World Health Organization. International Statistical Classification of Diseases and Related Health Problems, Tenth Revision. World Health Organization; 1992.
29.
SAS Institute SAS/STAT. 9.3 Users Guide. SAS Institute Inc; 2011.
30.
Stone DM, Mack KA, Qualters J. Notes from the field: recent changes in suicide rates, by race and ethnicity and age group—United States, 2021. MMWR Morb Mortal Wkly Rep. 2023;72(6):160-162. doi:
31.
Allison P. Survival Analysis Using SAS: A Practical Guide. 2nd ed. SAS Institute Inc; 2010.
32.
Kleinbaum DG, Klein M. Survival Analysis: A Self-Learning Text. 3rd ed. Springer; 2012. doi:
33.
Breslow NE, Day NE. Statistical methods in cancer research: volume II—the design and analysis of cohort studies. IARC Sci Publ. 1987;(82):1-406.
34.
Centers for Disease Control and Prevention. CDC WONDER: compressed mortality file. May 12, 2023. Accessed July 21, 2023.
35.
Kim HJ, Fay MP, Feuer EJ, Midthune DN. Permutation tests for joinpoint regression with applications to cancer rates. Stat Med. 2000;19(3):335-351. doi:
36.
Therneau TM, Grambsch PM. Testing proportional hazards. In: Modeling Survival Data: Extending the Cox Model. Statistics for Biology and Health. Springer; 2000. doi:
37.
Morral AR, Schell TL, Smart R. Comparison of suicide rates among US veteran and nonveteran populations. Ѵ Netw Open. 2023;6(7):e2324191. doi:
38.
Bakst SS, Braun T, Zucker I, Amitai Z, Shohat T. The accuracy of suicide statistics: are true suicide deaths misclassified? Soc Psychiatry Psychiatr Epidemiol. 2016;51(1):115-123. doi:
39.
Snowdon J, Choi NG. Undercounting of suicides: where suicide data lie hidden. Glob Public Health. 2020;15(12):1894-1901. doi:
1 Comment for this article
EXPAND ALL
Corrections for Suicide Risk Among US Veterans With Military Service During the Vietnam War
Tim Bullman, MA | Health Outcomes Military Exposures, Veterans Health Administration, Department of Veterans Affairs (VA), Washington, DC
In the Original Investigation titled “Suicide Risk Among U.S. Veterans with Military Service During the Vietnam War”[1], the frequency distribution of branch of military service for theater veterans is incorrect.

While doing more in-depth research of Marine Vietnam veterans for another study, we discovered that the number of Vietnam, ie, theater, Marines we reported in our paper was less than that reported in various sources. The underreporting in our original manuscript is related to the business rules used by the study to determine branch of service. Reprioritizing the business rules guiding incorporation of sources from which
branch was selected has resulted in a frequency distribution of branch of military service among theater veterans that more accurately reflects what is reported in additional reputable historical sources.
We have therefore revised branch of service for all theater veterans in Table 1. Analyses that were branch specific, Cox model results in Table 2, were recalculated based on new numbers of military branch for theater veterans. While the hazard ratios (HR)s and confidence intervals for branch specific risks did change, the results in terms of no branch having an increased risk of suicide associated with Vietnam deployment did not change. The changes in branch for theater veterans did not result in any changes in interpretation of data or conclusions.

The correction for branch of service for theater veterans in Tables 1 and 2 are as follows:

Table 1

Demographic/

Military Service Characteristics

Theater

(n=2,465,343)

Non-Theater

(n=7,122,976)



No. (%)

No. (%)

Branch of Service





Army

1,388,214 (56.3)

2,957,728 (41.5)

Marines

346,367 (14.1)

376,560 (5.3)

Air Force

486,450 (19.7)

1,251,967 (17.6)

Navy

241,102 (9.8)

1,586,229 (22.3)

Coast Guard



2902 (0.1)

71,906 (1.0)

Otherb

7 (0.0)

9,354 (0.1)

Unknown

301 (0.0)

869,232 (12.2)

Table 2:

Outcome/Cohort



No.

H.R.

Vietnam

Service

(95% C.I.)

Branch of Service





Suicide/Army Only





Theater

9123

0.86

(0.84-0.88)

Non-Theater

20544

Referent

Suicide/Marines Only





Theater

2556

1.05

(1.00-1.12

Non-Theater

3281

Referent

Suicide/Navy Only





Theater

1522

0.84

(0.79-0.88)

Non-Theater

11158

Referent

Suicide/Air Force Only





Theater

2732

0.83

(0.79-0.87)

Non-Theater

8069

Referent

From Results:

1. “Theater veterans were more likely to have been enlisted, n=2,228,986 (90.4%) and served in the Army, n=1,338,214 (56.3%) versus non-theater veterans, 4,249,372 (59.7%) enlisted and 2,957,728 (41.5%) Army.”

2. “Assessing branch specific suicide risks, only Marines theater veterans did not have a statistically significant decreased risk of suicide (HR=1.05; 95%, C.I., 1.00-1.12).” There are no additional errors or changes required for this article and I have requested that the article be corrected. All authors agree with these corrections.

Reference
1. Bullman TA, Akhtar FZ, Morley SW, et al. Suicide Risk Among US Veterans With Military Service During the Vietnam War. JAMA Netw Open. 2023;6(12):e2347616. doi:10.1001/jamanetworkopen.2023.47616.



CONFLICT OF INTEREST: None Reported
READ MORE
Original Investigation
Public Health
December28, 2023

Suicide Risk Among US Veterans With Military Service During the Vietnam War

Author Affiliations
  • 1Health Outcomes Military Exposures, Veterans Health Administration, Department of Veterans Affairs (VA), Washington, DC
  • 2Veterans Integrated Service Network (VISN) 2, Center of Excellence for Suicide Prevention, Department of VA, Washington, DC
  • 3VA Palo Alto Health Care System, Menlo Park, California
  • 4Stanford University School of Medicine, Stanford, California
  • 5Center for Excellence for Suicide Prevention, VA Finger Lakes Healthcare System, Department of VA, Canandaigua, New York
  • 6Department of Psychiatry, University of Rochester School of Medicine and Dentistry, Rochester, New York
  • 7Office of Research and Development, Veterans Health Administration, Department of VA, Washington, DC
JAMA Netw Open. 2023;6(12):e2347616. doi:10.1001/jamanetworkopen.2023.47616
Key Points

Question Is the risk of suicide among US veterans who served in the Vietnam War between 1961 and 1975 greater than that of veterans who served in the military at the same time but not in the Vietnam War theater?

Findings In this cohort analysis of approximately 9.6 million Vietnam War–era veterans, there was no increased risk of suicide mortality among veterans who served in the Vietnam War theater and veterans who served during the Vietnam War but were not deployed to Vietnam.

Meaning Findings of this study suggest that although deployment and service during the Vietnam War were not associated with an increased risk of suicide, the high number of suicides among US military personnel between 1979 and 2019 is noteworthy and merits the ongoing attention of health policymakers and mental health professionals.

Abstract

Importance There are persistent questions about suicide deaths among US veterans who served in the Vietnam War. It has been believed that Vietnam War veterans may be at an increased risk for suicide.

Objective To determine whether military service in the Vietnam War was associated with an increased risk of suicide, and to enumerate the number of suicides and analyze patterns in suicides among Vietnam War theater veterans compared with the US population.

Design, Setting, and Participants This cohort study compiled a roster of all Vietnam War–era veterans and Vietnam War theater veterans who served between February 28, 1961, and May 7, 1975. The 2 cohorts included theater veterans, defined as those who were deployed to the Vietnam War, and nontheater veterans, defined as those who served during the Vietnam War era but were not deployed to the Vietnam War. Mortality in these 2 cohorts was monitored from 1979 (beginning of follow-up) through 2019 (end of follow-up). Data analysis was performed between January 2022 and July 2023.

Main Outcomes and Measures The outcome of interest was death by suicide occurring between January 1, 1979, and December 31, 2019. Suicide mortality was ascertained from the National Death Index. Hazard ratios (HRs) that reflected adjusted associations between suicide risk and theater status were estimated with Cox proportional hazards regression models. Standardized mortality rates (SMRs) were calculated to compare the number of suicides among theater and nontheater veterans with the expected number of suicides among the US population.

Results This study identified 2 465 343 theater veterans (2 450 025 males [99.4%]; mean [SD] age at year of entry, 33.8 [6.7] years) and 7 122 976 nontheater veterans (6 874 606 males [96.5%]; mean [SD] age at year of entry, 33.3 [8.2] years). There were 22 736 suicides (24.1%) among theater veterans and 71 761 (75.9%) among nontheater veterans. After adjustments for covariates, Vietnam War deployment was not associated with an increased risk of suicide (HR, 0.94; 95% CI, 0.93-0.96). There was no increased risk of suicide among either theater (SMR, 0.97; 95% CI, 0.96-0.99) or nontheater (SMR, 0.97; 95% CI, 0.97-0.98) veterans compared with the US population.

Conclusions and Relevance This cohort study found no association between Vietnam War–era military service and increased risk of suicide between 1979 and 2019. Nonetheless, the 94 497 suicides among all Vietnam War–era veterans during this period are noteworthy and merit the ongoing attention of health policymakers and mental health professionals.

Introduction

Suicide is a leading cause of death in the US, with close to 50 000 deaths annually.1,2 In 2016, the US Department of Veterans Affairs (VA) began publishing national data on veteran suicides in annual suicide prevention reports.3 In the 2022 National Veteran Suicide Prevention report, the VA reported that although veterans composed only 7.6% of the US population, they accounted for almost 14% of US suicides.4 This report compared sex- and age-adjusted suicide rates in veterans vs nonveterans for each year from 2001 to 2020 and found that suicide rates among veterans remained consistently higher than among nonveterans.4 In 2001, the unadjusted annual population-based estimate of suicide mortality was 23.3 per 100 000 persons for veterans vs 12.6 per 100 000 persons for nonveterans. By 2020, the suicide mortality rate had increased to 31.7 per 100 000 persons for veterans vs 16.1 per 100 000 persons for nonveterans.4

In the 1980s, media accounts cited figures of more than 50 000 Vietnam War theater veteran suicides in 1981 and 58 000 in 1985.5,6 Although these estimates were later found to be incorrect,7 questions regarding the association of Vietnam War service with suicide risk emerged and remain. Most published research assessing postservice mortality of Vietnam War theater veterans relied on samples and specific subsets and did not find an association between serving in the Vietnam War and suicide risk.8-16 Other studies examined theater veteran cohorts with characteristics believed to be associated with greater risk of suicide, including veterans diagnosed with posttraumatic stress disorder17 and veterans who were wounded during military service.18 Both studies reported an increased risk of suicide among their specific groups of theater veterans. Limitations of the previous literature included a focus on population subsets of veterans from specific military branches8,9,11 or veterans who resided in specific states.12-14 Observation periods in these studies were highly variable.

To address these limitations, this Vietnam era veteran mortality study identified all veterans during the Vietnam War era, between 1961 and 1975, and divided them into 2 cohorts: (1) theater veterans, defined as those who were deployed to the Vietnam War, and (2) nontheater veterans, defined as those who served during the Vietnam War era but were not deployed to the Vietnam War. The suicide risk among these cohorts was evaluated over a 41-year period (1979-2019). The study objectives were to examine (1) baseline statistics on suicide by theater service (ie, theater vs nontheater); (2) crude suicide mortality rates for all theater and nontheater veterans combined, for each theater service group, and for the US male population; (3) adjusted associations between suicide mortality and theater service; and (4) suicide risk among Vietnam War–era veterans stratified by theater service compared with the expected suicide risk for the US general population. The latter analysis differed methodologically from the annual suicide prevention reports by the VA, which compare adjusted suicide rates by calendar year between all veterans and nonveterans.3,4,19,20 In contrast, the present study used a cumulative, stratified assessment of suicide risk across calendar years as well as calculated crude rates (CRs) by calendar year among all Vietnam War–era veterans collectively and, when stratified, by theater vs nontheater service. We believe this study makes a unique and comprehensive contribution to existing Vietnam veteran suicide mortality research by examining the universe of Vietnam War–era veterans starting from 1979, with a particular focus on those who served in the Vietnam War theater.

Methods

The VA Central Institutional Review Board approved this cohort study and waived the informed consent requirement because the study involved minimal risk to study participants. We followed the Strengthening the Reporting of Observational Studies in Epidemiology () reporting guideline.

Vietnam War–Era Cohort Selection

The US Veterans Eligibility Trends and Statistics (USVETS) database21,22 was used to identify all veterans (n = 9 826 155) (Figure 1) who served in the US military between February 28, 1961, and May 7, 1975.23 Theater service was determined using the Defense Manpower Data Center (DMDC) Vietnam File (n = 2 828 834),24 which recorded all veterans who served in Southeast Asia during the Vietnam War. The USVETS database and the Vietnam File included Social Security numbers, demographics, and military service characteristics. Race and ethnicity data were obtained from the USVETS database, which compiled the data from the electronic records of various federal sources, including the US Department of Defense (DOD) and VA. It is not known whether race and ethnicity were self-reported by the veterans at the time of initial contact with the VA or DOD or were based on observation.

Records with a date of birth prior to 1900 or after 1958 and/or a date of death prior to January 1, 1979 (beginning of follow-up for this study), in the USVETS database were excluded. Merging of the records from USVETS database and the Vietnam File resulted in the identification of 2 518 221 theater veterans and 7 307 819 nontheater veterans (Figure 1). After exclusions for birth year and date of death, final counts were 2 465 343 for eligible theater veterans and 7 122 976 for eligible nontheater veterans.

Vital Status and Cause of Death Ascertainment

The VA/DOD Mortality Data Repository (MDR)25 and the USVETS database provided vital status and date of death. The MDR also provided cause of death obtained from annual searches of the National Death Index (NDI) of the National Center for Health Statistics.26 At the time of analysis, MDR data were available from 1979 through 2019.

Suicide analyses were limited to records that had been searched in the NDI to enable identification of suicide deaths based on International Classification of Diseases (ICD) codes. Veterans were excluded if their Social Security numbers were not searched in the NDI. Applying these criteria yielded 2 443 238 eligible theater and 7 051 418 eligible nontheater veterans for inclusion in suicide risk models (Figure 1).

Statistical Analysis

Follow-up began on January 1, 1979, the earliest date for which cause of death was available, and ended on date of death or December 31, 2019, whichever occurred first. Suicide deaths included International Classification of Diseases, Ninth Revision (ICD-9) codes E950 to E959 for deaths occurring between 1979 and 1998 and International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10) codes X60 to X84, Y87.0, and U03 for deaths after 1998.27,28 The CRs per 100 000 person-years at risk approximated risks of suicide among theater and nontheater veterans separately and all Vietnam War–era veterans collectively.

The outcome of interest for this study was death by suicide occurring between January 1, 1979, the beginning of follow-up, and December 31, 2019, the end of follow-up. Hazard ratios (HRs) and 95% CIs that reflected adjusted associations between suicide risk and theater status were estimated using Cox proportional hazards regression models generated by the PHREG procedure (SAS Institute Inc).29 Model covariates were age at entry to follow-up (January 1, 1979), sex (male or female), and race (White or other). Race was based on a composite of race and ethnicity (Hispanic: yes, no, or unknown). White referred to White race and non-Hispanic ethnicity. Other races referred to Hispanic, non-Hispanic Black, and other racial groups (American Indian or Alaska Native, Asian, Hawaiian, and Pacific Islander). White race and unknown Hispanic ethnicity were treated as unknown and excluded from Cox proportional hazards regression models; this coding scheme followed the algorithm used by the Centers for Disease Control and Prevention (CDC) in a study assessing suicide rates by race and ethnicity.30

Observations missing 1 or more of the covariates included in the Cox proportional hazards regression models were excluded by SAS from the analysis. Cox proportional hazards regression model analyses were also stratified by branch of service and method-specific suicide risks. The latter included the following suicide methods: (1) intentional poisoning, including intentional drug overdoses (ICD-9 codes E950-E952; ICD-10 codes X60-X69); (2) suffocation, including hanging (ICD-9 code E953; ICD-10 code X70); (3) firearm (ICD-9 codes E955.0-E955.4; ICD-10 codes X72-X74); and (4) all other suicides. The underlying assumption of the Cox proportional hazards regression model (ie, the relative hazard remained constant over time with different variable or covariate levels) was tested both graphically31 and statistically32 (eAppendix in Supplement 1).

Standardized mortality rates (SMRs) were calculated to compare the observed counts of suicide for all Vietnam War–era veterans collectively and for theater and nontheater veterans separately with the expected counts for the US population (eAppendix in Supplement 1). Specifically, the SMR was the ratio of the number of deaths observed in the veteran cohort to the number of deaths that would be expected if the veteran cohort had the same age-, race and ethnicity–, and sex-specific mortality rates as the US population33 (eAppendix in Supplement 1).

The PROC LIFETEST procedure (SAS Institute Inc)29 examined differences in crude suicide hazard rates (CRs) for veterans over time by theater service. Hazard rates for both veterans and similarly aged US males, derived from the CDC WONDER database,34 were calculated per 100 000 persons at risk (eAppendix in Supplement 1). Changes in temporal patterns of suicide for the veteran cohorts and the similarly aged US male population were further examined using Joinpoint Trend Analysis Software, version 4.0 (National Cancer Institute).35 Joinpoint determined the annual percent change (APC) and the average annual percent change (AAPC) in suicide rates since entry into follow-up for veterans and by calendar year for US males based on the best-fit joinpoint regression models selected at an overall α = .05. For veterans, the number of years since entry to follow-up was used as the input for joinpoint and was based on midpoints of calendar years. The initial period of January 1979 to mid-June 1979 (5.5 months) was excluded from the joinpoint analysis because deaths occurring during this period with less than 1 full year of follow-up would not have sufficient follow-up time to support statistical calculations.

Unlike veteran rates, US male suicide rates were enumerated for each calendar year from 1979 to 2019. Because follow-up began on the same date for all veterans (January 1, 1979), calendar year was approximately analogous to the elapsed number of years since entry to follow-up; therefore, the rates for both veterans and male US population were plotted by calendar years.

Two-sided P < .05 indicated statistical significance. Data analysis was performed between January 2022 and July 2023, using SAS, version 9.3 (SAS Institute Inc).

Results
Demographic and Military Service Characteristics

A total of 2 465 343 theater veterans and 7 122 976 nontheater veterans were identified. The theater veterans included 2 450 025 males (99.4%) and 15 277 females (0.6%) with a mean (SD) age at year of entry of 33.8 (6.7) years. The nontheater veterans included 6 874 606 males (96.5%) and 245 684 females (3.5%) with a mean (SD) age at year of entry of 33.3 (8.2) years (Table 1). Most Vietnam War–era veterans were identified as non-Hispanic White individuals (78.4% theater vs 76.0% nontheater). Theater veterans compared with their nontheater counterparts were more likely to have been enlisted (90.4% vs 59.7%) and to have served in the Army (56.3% vs 41.5%).

Because there were 2 data sources for theater veterans (USVETS database and DMDC Vietnam File) and only 1 source for nontheater veterans (USVETS database), the percentages of missing data were higher among nontheater veterans (Table 1). A substantial portion of theater and nontheater veterans were missing race and ethnicity data (Table 1). Both the Cox proportional hazards regression model and SMR analyses excluded veterans with a missing value for a covariate or a stratifier respectively; thus, the Cox proportional hazards regression model of suicide risk associated with deployment that included race and ethnicity as a covariate excluded 29.8% of theater and 34.2% of nontheater veteran suicides that were missing race and ethnicity (eAppendix in Supplement 1). The SMR analyses also lost observations due to missing race and ethnicity, although the numbers were smaller due to the different coding of race and ethnicity for SMR analysis, which followed the coding scheme for US expected rates. Sensitivity analyses conducted with Cox proportional hazards regression models that allowed for a greater inclusion of observations as a result of differential coding of race and ethnicity produced varied results regarding suicide risk associated with Vietnam War service (eAppendix in Supplement 1). Given the large numbers of both theater and nontheater veterans and the small numbers of veterans missing sex and age, missing sex and age data would have no implications for the reported risk of suicide.

Overall and Method-Specific Risks of Suicide

A total of 94 497 suicides were identified, of which 22 736 (24.1%) were among theater veterans (CR, 25.9 per 100 000 person-years at risk) and 71 761 (75.9%) were among nontheater veterans (CR, 28.6 per 100 000 person-years at risk). In adjusted analyses after exclusions for missing covariates, 15 949 theater veterans and 47 176 nontheater veterans were available for Cox proportional hazards regression model analyses (Table 2). For the 41 years between 1979 and 2019, Vietnam War theater status was associated with a decreased risk of suicide mortality among all Vietnam War–era veterans (HR, 0.94; 95% CI, 0.93-0.96) (Table 2). Assessing branch specific suicide risks, only Marines theater veterans did not have a significantly lower risk of suicide (HR, 1.05; 95% CI, 1.00-1.12). Moreover, risk of firearm suicide was not significantly different between theater vs nontheater veterans (HR, 0.99; 95% CI, 0.97-1.01). The leading method of suicide for each of the cohorts based on crude counts (including those excluded from Cox proportional hazards regression models) was by firearm. While graphically, it does not appear that the Cox proportional hazards regression model assumption was violated during the first 35 years of follow-up, statistical analysis did indicate a violation (P < .05). However, in a data set as large as this study’s (over 9 million records) and with small correlations (<0.1), statistically significant violations of proportional hazards should be of little practical significance.36

Suicide Mortality Risk Among Veterans vs US Population

None of the veteran cohorts showed increased risk for suicide mortality compared with the US population after adjusting for sex, race and ethnicity, age, and calendar year. There were fewer than expected suicides among the veteran groups separately (theater: SMR, 0.97 [95% CI, 0.96-0.99]; nontheater: SMR, 0.97 [95% CI, 0.97-0.98]) and all Vietnam War–era veterans collectively (SMR, 0.97; 95% CI, 0.97-0.98).

CRs of Suicide Mortality Over Time

Figure 2 presents the temporal patterns of suicide rates for Vietnam War–era veterans collectively, for each theater status, and for the US male population by calendar years 1979 to 2019. Deaths during the first period of less than a year were excluded; these 5.5 months included 964 suicides among nontheater veterans and 293 suicides among theater veterans. In 1979, after 1 year of follow-up, the CR of suicide among theater veterans was 21.9 per 100 000 persons, increasing to 34.9 per 100 000 persons in 2019 after 41 years. During the first 38 years of follow-up (1979-2016), there was a significant increase in suicide rates (APC, 0.4; 95% CI, 0.2-0.6), and the AAPC for the entire 41 years also showed a significant increase in rates (AAPC, 1.1; 95% CI, 0.3-1.8).

After 1 year of follow-up, the suicide CR among nontheater veterans was 27.9 per 100 000 persons, increasing to 33.0 per 100 000 persons after 41 years. In comparison with theater veterans, over the entire 41 years of follow-up, there was no overall increase in suicide rates (AAPC, 0.4; 95% CI, –0.2 to 1.1).

Among all Vietnam War–era veterans after 1 year of follow-up, the suicide CR was 26.4 per 100 000 persons, increasing to 33.5 after 41 years of follow-up. There were 3 periods of distinct trends representing an increase in suicide rates: 1979 to 1987 (APC, 0.8; 95% CI, 0.4-1.3), 2005 to 2009 (APC, 4.5; 95% CI, 2.1-7.0), and 2012 to 2019 (APC, 2.1; 95% CI, 1.3-2.8). Over the entire 41 years of follow-up, there was also a significant increasing trend in suicide rates (AAPC, 0.6; 95% CI, 0.1-1.0).

Among US males, the suicide rate in 1979, corresponding to approximately 1 year of follow-up for veterans in this study, was 25.2 per 100 000 persons and increased to 26.4 in 2019, corresponding to 41 years of follow-up for veterans in this study, with no calendar year having a significant change in suicide rates. The suicide rates for veterans and US males exhibited similar patterns of increase and decrease over time as they aged (Figure 2).

Discussion

To our knowledge, this cohort study is the most comprehensive examination to date of the association between Vietnam War military service and suicide risk. The results showed no association between either Vietnam War–era military service or theater service and increased risk of suicide.

Compared with the expected US population rates using SMRs, suicide risk across 41 years of observation was not increased for theater, nontheater, and all Vietnam War–era veterans collectively. All 3 groups of veterans had a similar deficit of suicides compared with the US population (eAppendix in Supplement 1). In veteran-to-veteran comparisons using Cox proportional hazards regression models, theater veterans had modestly decreased suicide risk vs nontheater veterans, which was robust to adjustment across all branches of military service. There was also no increased risk in method-specific suicide associated with theater status.

Comparing the AAPCs between theater and nontheater veterans, there was no statistical difference between the overall trends between the 2 groups. However, among theater veterans there was a significant increase in AAPC in the 41 years of post-1979 follow-up, and risk increased from age 33 years (mean age at year of entry) through age 71 years at follow-up. Comparing the veteran cohorts to US males revealed that the CRs of suicides among veterans remained higher than the suicide CRs of the US population for each of the 41 calendar years, but the overall trends were similar over time.

Of interest was the finding of no increased suicide risk after 41 years of follow-up among veterans compared with the US population based on the SMR, although higher CRs of suicide by calendar year were also found among the veterans compared with US males. The higher rates by calendar year in veterans are consistent with findings in annual VA suicide reports that compare adjusted suicide rates by calendar year among all veterans to the US population. A recent study37 that compared veteran and nonveteran suicide rates using both direct and indirect standardization methods similarly found that estimates of suicide risk among veterans were lower based on indirect standardization (ie, SMR). The CRs by calendar year reported in the present study, the direct adjusted rates by calendar year in the VA suicide prevention reports, and the cumulative estimate of suicide risk (as approximated by the SMR) are all valid estimates of suicide risk. The CRs provided snapshots of rates for each year, whereas the SMR provided cumulative risk estimates wherein the expected number of events was standardized to the veteran cohort (eAppendix in Supplement 1). Furthermore, SMR analyses in this study differed from the annual VA data reports in that they examined rates across 41 years rather than over single years, used indirect vs direct adjustment, and adjusted for race and ethnicity, which were not accounted for in the VA suicide prevention reports. A complete explanation of the methodological differences between these 2 estimates of risk is provided in the eAppendix in Supplement 1.

Strengths and Limitations

This study has 2 major strengths that distinguish it from prior studies assessing suicide mortality among Vietnam War–era veterans. First, the entire population of Vietnam War–era veterans, rather than samples or subsets, was included. Second, 41 years of follow-up were available, adding 19 additional years to the Vietnam Experience Study on postservice mortality.9

The most consequential limitation of the present Vietnam Era veteran mortality study was the absence of mortality data before January 1979. Given that the Vietnam Experience Study,9 a 30-year follow-up study of Vietnam War–era veterans, reported the highest suicide rate within the first 5 years after military discharge, the current study may have underestimated the number of suicides among all Vietnam War–era veterans and missed an excess of theater veteran suicides that might have occurred soon after discharge from the military and/or deactivation after the end of deployment. Given the modest decreased risks among theater veterans compared with nontheater veterans (HR, 0.94; 95% CI, 0.93-0.96) and with the US population (SMR, 0.97; 95% CI, 0.96-0.99) as well as the lack of mortality data for years reportedly with the highest suicide rates among Vietnam War theater veterans, any statement of a decreased risk of suicide among theater veterans, although accurate in the context of this study, may not accurately reflect the risk when pre-1979 suicides were included.

Another limitation was the lack of data on potential confounders or exposures (eg, serving in a combat role while deployed). Previous research showed that certain subgroups of theater veterans were at an increased risk for suicide (eg, those diagnosed with posttraumatic stress disorder and those wounded during service).17,18 However, the present study does answer the question of whether Vietnam War service in general is associated with suicide risk. Underreporting of suicides on death certificates has been documented.38,39 However, underreporting should be random regarding theater and nontheater status.

Results of suicide risk estimates varied by how missing data for race and ethnicity were handled (eAppendix in Supplement 1). Excluding rather than including race and ethnicity (Table 2) from models produced similar suicide risk estimates (ie, no increased risk of suicide associated with deployment). Suicide risks differed primarily by how those with White race or unknown Hispanic ethnicity were identified (ie, missing vs White). Given the large number of observations with missing race and ethnicity and the possibility that missing data may not be random regarding theater vs nontheater status, it was unwise to attempt to impute missing values.

Conclusions

To our knowledge, this was the first study to include all Vietnam War–era veterans in an assessment of suicide dichotomized by service in the Vietnam War theater. Evaluation of mortality data from 1979 to 2019 revealed no association between theater service and risk of suicide. While individual calendar year suicide rates among veterans were higher than that for the US male population, after 41 years of follow-up, we found that all Vietnam War–era veterans collectively or by theater status did not have increased risks of suicide compared with the US population.

Despite the absence of an increased risk of suicide among theater veterans or among all Vietnam War–era veterans, it is nonetheless sobering to reflect on the 94 497 suicide deaths that have occurred among all Vietnam War–era veterans since 1979. This loss of life deserves not only to be noted on behalf of these veterans and their survivors but also merits the ongoing attention of health policymakers and mental health professionals, especially given that suicide rates have increased over 41 years among all Vietnam War–era veterans, veterans of other eras of military service, and the wider US population.

Back to top
Article Information

Accepted for Publication: October 31, 2023.

Published: December 28, 2023. doi:10.1001/jamanetworkopen.2023.47616

Correction: This article was corrected on May 17, 2024, to fix errors in the Results section and Tables 1 and 2.

Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2023 Bullman TA et al. vlog Open.

Corresponding Author: Tim A. Bullman, MA, Health Outcomes Military Exposures (133), Veterans Health Administration, Department of Veterans Affairs, 810 Vermont Ave NW, Washington, DC 20420 (tim.bullman@va.gov).

Author Contributions: Mr Bullman had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Concept and design: Bullman, Morley, Weitlauf, Culpepper, Britton, Davey.

Acquisition, analysis, or interpretation of data: Bullman, Akhtar, Morley, Weitlauf, Cypel, Culpepper, Schneiderman, Davey.

Drafting of the manuscript: Bullman, Morley, Weitlauf, Cypel, Davey.

Critical review of the manuscript for important intellectual content: All authors.

Statistical analysis: Bullman, Akhtar, Cypel, Culpepper.

Obtained funding: Schneiderman, Davey.

Administrative, technical, or material support: Weitlauf, Schneiderman, Davey.

Supervision: Culpepper, Schneiderman, Davey.

Conflict of Interest Disclosures: Dr. Weitlauf reported receiving grants from the US Department of Veterans Affairs (VA) Center for Implementation to Innovation and the Veterans Integrated Service Network 21 Sierra Pacific Mental Illness, Research, Education, and Clinical Center during the conduct of the study. No other disclosures were reported.

Funding/Support: This study was funded through in-kind support from the VA Office of Research and Development and the Health Outcomes Military Exposures.

Role of the Funder/Sponsor: The funder had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.

Data Sharing Statement: See Supplement 2.

Additional Contributions: Dennis Fried, PhD, provided initial cohort programming. George Fitzelle, PhD, assisted with initial study design and data acquisition. Rowan Poehler, MHA, Aracely Suarez, MHA, and Kadion Thompson, BS, provided project support and editorial assistance. Pheakday Lim, MHA, and Charles Lin, PhD, provided access to data and understanding of available data. Ying Lu, PhD, provided statistical consultation. These individuals received no additional compensation, outside of their usual salary, for their contributions.

References
1.
Centers for Disease Control and Prevention. Preventing suicide. 2022. Accessed July 7, 2023.
2.
Curtin SC, Hedegaard H, Ahmad FB. Provisional numbers and rates of suicide by month and demographic characteristics: United States, 2020: Vital Statistics Rapid Release, no 16. National Center for Health Statistics. November 2021. Accessed July 12, 2023.
3.
US Department of Veterans Affairs, Office of Mental Health and Suicide Prevention. Suicide among veterans and other Americans 2001-2014. 2016. Accessed July 20, 2023.
4.
US Department of Veterans Affairs. 2022 National Veteran Suicide Prevention annual report. Accessed July 12, 2023.
5.
Anderson R. Vietnam legacy: veterans’ suicide toll may top war casualties. Seattle Times. March 18, 1981.
6.
Langone J. The war that has no ending. پDZ. 1985;(June):44-54.
7.
Pollock DA, Rhodes P, Boyle CA, Decoufle P, McGee DL. Estimating the number of suicides among Vietnam veterans. Am J Psychiatry. 1990;147(6):772-776. doi:
8.
Centers for Disease Control Vietnam Experience Study. Postservice mortality among Vietnam veterans. Ѵ. 1987;257(6):790-795. doi:
9.
Boehmer TK, Flanders WD, McGeehin MA, Boyle C, Barrett DH. Postservice mortality in Vietnam veterans: 30-year follow-up. Arch Intern Med. 2004;164(17):1908-1916. doi:
10.
Watanabe KK, Kang HK, Thomas TL. Mortality among Vietnam veterans: with methodological considerations. J Occup Med. 1991;33(7):780-785. doi:
11.
Watanabe KK, Kang HK. Military service in Vietnam and the risk of death from trauma and selected cancers. Ann Epidemiol. 1995;5(5):407-412. doi:
12.
Anderson H, Hanrahan LP, Jensen M, Laurin D, Yick W, Weigman P. Wisconsin Vietnam Veteran Mortality. Wisconsin Division of Health; 1986.
13.
Kogan MD, Clapp RW. Mortality Among Vietnam Veterans in Massachusetts, 1972-1983. Massachusetts Department of Public Health, Division of Health Statistics; 1985.
14.
Lawrence CE, Reilly AA, Quickenton P, Greenwald P, Page WF, Kuntz AJ. Mortality patterns of New York State Vietnam Veterans. Am J Public Health. 1985;75(3):277-279. doi:
15.
Forsberg CW, Estrada SA, Baraff A, et al. Risk factors for suicide in the Vietnam-era twin registry. Suicide Life Threat Behav. 2022;52(4):631-641. doi:
16.
Kang HK, Cypel Y, Kilbourne AM, et al. HealthViEWS: mortality study of female US Vietnam era veterans, 1965-2010. Am J Epidemiol. 2014;179(6):721-730. doi:
17.
Bullman TA, Kang HK. Posttraumatic stress disorder and the risk of traumatic deaths among Vietnam veterans. J Nerv Ment Dis. 1994;182(11):604-610. doi:
18.
Bullman TA, Kang HK. The risk of suicide among wounded Vietnam veterans. Am J Public Health. 1996;86(5):662-667. doi:
19.
US Department of Veterans Affairs, Office of Mental Health and Suicide Prevention. 2020 National Veteran Suicide Prevention annual report. Accessed July 20, 2023.
20.
US Department of Veterans Affairs, Office of Mental Health and Suicide Prevention. 2021 National Veteran Suicide Prevention annual report. Accessed July 26, 2023.
21.
VHA Data Portal. United States Veterans Eligibility Trends and Statistics (USVETS). Accessed July 19, 2023.
22.
VA Information Resource Center (VIReC). United States Veterans Eligibility Trends and Statistics (USVETS): a new data source with socioeconomic variables. 2019. Accessed July 18, 2023.
23.
Title 10–Armed Forces: this title was enacted by act Aug. 10, 1956, ch. 1041, § 1, 70A Stat. 1. Accessed July 25, 2023.
24.
Defense Manpower Data Center. Vietnam File. July 1999. Accessed July 14, 2023.
25.
US Department of Veterans Affairs. Using data to prevent veteran mortality. 2022. Accessed July 10, 2023.
26.
Hoffmire CABS, Barth SK, Bossarte RM. Reevaluating suicide mortality for veterans with data from the VA-DoD mortality data repository, 2000-2010. Psychiatr Serv. 2020;71(6):612-615. doi:
27.
World Health Organization. Manual of the International Statistical Classification of Diseases: Inquiries, and Causes of Death, Ninth Revision. World Health Organization; 1977.
28.
World Health Organization. International Statistical Classification of Diseases and Related Health Problems, Tenth Revision. World Health Organization; 1992.
29.
SAS Institute SAS/STAT. 9.3 Users Guide. SAS Institute Inc; 2011.
30.
Stone DM, Mack KA, Qualters J. Notes from the field: recent changes in suicide rates, by race and ethnicity and age group—United States, 2021. MMWR Morb Mortal Wkly Rep. 2023;72(6):160-162. doi:
31.
Allison P. Survival Analysis Using SAS: A Practical Guide. 2nd ed. SAS Institute Inc; 2010.
32.
Kleinbaum DG, Klein M. Survival Analysis: A Self-Learning Text. 3rd ed. Springer; 2012. doi:
33.
Breslow NE, Day NE. Statistical methods in cancer research: volume II—the design and analysis of cohort studies. IARC Sci Publ. 1987;(82):1-406.
34.
Centers for Disease Control and Prevention. CDC WONDER: compressed mortality file. May 12, 2023. Accessed July 21, 2023.
35.
Kim HJ, Fay MP, Feuer EJ, Midthune DN. Permutation tests for joinpoint regression with applications to cancer rates. Stat Med. 2000;19(3):335-351. doi:
36.
Therneau TM, Grambsch PM. Testing proportional hazards. In: Modeling Survival Data: Extending the Cox Model. Statistics for Biology and Health. Springer; 2000. doi:
37.
Morral AR, Schell TL, Smart R. Comparison of suicide rates among US veteran and nonveteran populations. Ѵ Netw Open. 2023;6(7):e2324191. doi:
38.
Bakst SS, Braun T, Zucker I, Amitai Z, Shohat T. The accuracy of suicide statistics: are true suicide deaths misclassified? Soc Psychiatry Psychiatr Epidemiol. 2016;51(1):115-123. doi:
39.
Snowdon J, Choi NG. Undercounting of suicides: where suicide data lie hidden. Glob Public Health. 2020;15(12):1894-1901. doi:
×