Male and Female Physician Suicidality: A Systematic Review and Meta-analysis

Key Points

Question  Are age-standardized suicide mortality ratios (SMRs) for male and female physicians high worldwide after 1980, and did they decrease from before to after 1980?

Findings  In a meta-analysis of 9 studies and databases, female physicians’ SMRs were high, while male physicians’ SMRs were significantly lower since 1980, and both SMRs decreased over time. A systematic review showed increased risk associated with male sex, relationship difficulties, youthful or elderly age ranges, and especially career difficulties; suicide methods were by poisoning, firearms, and asphyxiation.

Meaning  More research is critical to address modifiable risks of physician suicide and understand what causes vulnerability in physician subpopulations.

Importance  Population-based findings on physician suicide are of great relevance because this is an important and understudied topic.

Objective  To evaluate male and female physician suicide risks compared with the general population from 1980 to date and test whether there is a reduction of SMR in cohorts after 1980 compared with before 1980 via a meta-analysis, modeling studies, and a systematic review emphasizing physician suicide risk factors.

Data Sources  This study uses studies retrieved from PubMed, Scielo, PsycINFO, and Lilacs for human studies published by October 3, 2019, using the search term “(((suicide) OR (self-harm) OR (suicidality)) AND ((physicians) OR (doctors))).” Databases were also searched from countries listed in articles selected for review. Data were also extracted from an existing article by other authors to facilitate comparisons of the pre-1980 suicide rate with the post-1980 changes.

Study Selection  Original articles assessing male and/or female physician suicide were included; for the meta-analysis, only cohorts from 1980 to the present were included.

Data Extraction and Synthesis  The preregistered systematic review and meta-analysis followed Cochrane, PRISMA, and MOOSE guidelines. Data were extracted into standardized tables per a prespecified structured checklist, and quality scores were added. Heterogeneity was tested via Q test, I², and τ². For pooled effect estimates, we used random-effects models. The Begg and Egger tests, sensitivity analyses, and meta-regression were performed. Proportional mortality ratios were presented when SMR data could not be extracted.

Main Outcomes and Measures  Suicide SMRs for male and female physicians from 1980 to the present and changes over time (before and after 1980).

Results  Of 7877 search results, 32 articles were included in the systematic review and 9 articles and data sets in the meta-analysis. Meta-analysis showed a significantly higher suicide SMR in female physicians compared with women in general (1.46 [95% CI, 1.02-1.91]) and a significantly lower suicide SMR in male physicians compared with men in general (0.67 [95% CI, 0.55-0.79]). Male and female physician SMRs significantly decreased after 1980 vs before 1980 (male physicians: SMR, −0.84 [95% CI, −1.26 to −0.42]; P < .001; female physicians: SMR, −1.96 [95% CI, −3.09 to −0.84]; P = .002). No evidence of publication bias was found.

Conclusions and Relevance  In this study, suicide SMR was found to be high in female physicians and low in male physicians since 1980 but also to have decreased over time in both groups. Physician suicides are multifactorial, and further research into these factors is critical.

Suicide leads to more annual deaths globally than natural disasters, violence inflicted by others, war, and conflict combined.^1,2 Suicide caused almost 800 000 deaths worldwide in 2016, down from nearly 850 000 suicide deaths in 2004.^1,2 In 2016, the age-standardized suicide rates were 10.53 per 100 000 persons globally.^1,2 Crude suicide rates for both men and women decreased between 2000 and 2016 in Europe and all World Health Organization regions other than the Americas.^1,2 In the United States, age-standardized suicide rates rose for both men and women over the same period, reaching 21.1 and 6.4 deaths per 100 000 persons, respectively.^1,2

Suicide risk factors are male sex (for completed suicides),^3,4 younger age, fewer years of formal education, unmarried status, and the presence of mental disorders.^5-7 Occupational hazards are also associated with suicide risk^8-10; physicians typically have higher suicide rates than the general population, including the military.⁹ Overall, age-standardized suicide mortality ratios (SMRs) for physicians were significantly higher (ie, higher suicide rates in physicians compared with the general population) according to a meta-analysis, including much higher findings for female physicians and moderately higher findings for male physicians.¹¹ However, heterogeneity was considerable for the meta-analysis overall, and data quality was higher for the suicides of male physicians, possibly because the data ranged from 1910 to 1998 and multiple articles had no data on female physicians.

Meanwhile, from 1980 to 2015, the proportion of female US medical school graduates increased from 23.3% to 47.6%,¹² but female participants in the US civilian labor force rose only from 42.5% to 46.8%.¹³ This implies that gender dynamics in medicine changed much more rapidly than in the general population during this period, which may have affected female and male physicians’ SMRs. However, many additional cultural, socioeconomic, and political changes influenced health care during this period,^14,15 reinforcing the need for action regarding physician suicides.

As to the general population, the rates of college and university graduation (a protective factor against suicide) are increasing among women but decreasing among men.¹⁶ Additionally, unemployment, which affects suicide risk among men more than women,¹⁷ has been increasing worldwide.¹⁸ Even among individuals who are employed, suicide rates are highest among working-class men.¹⁹ These factors are consistent with increasing suicide rates among men in general, which are less likely to affect physicians in recent times.

In light of the literature, we hypothesized a higher risk of male and female physician suicide compared with the general population from the 1980s to date; we also hypothesized that the suicide rate ratio of both female and male physicians seen before 1980 would be reduced after 1980.¹¹ In this systematic review and meta-analysis, we aimed to answer the following research questions: are age-standardized suicide SMRs for male and female physicians high worldwide after 1980, and did they decrease after 1980 compared with before 1980?

This systematic review and meta-analysis were conducted in accordance with recommendations of the Cochrane group, the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) Statement,²⁰ and the Meta-analyses of Observational Studies in Epidemiology Checklist (MOOSE)²¹ guidelines. The study is registered at the International Prospective Register of Ongoing Systematic Reviews (CRD42019124985).

We used Boolean operators on PubMed, Scielo, PsycINFO, and Lilacs for human research studies published from database inception to October 3, 2019, using the following keywords: “(((suicide) OR (self-harm) OR (suicidality)) AND ((physicians) OR (doctors)).” In PubMed, the Medical Subject Heading (MeSH) term suicide was combined with the MeSH term physicians. Additionally, we largely followed the methods used by Schernhammer and Colditz¹¹ to make our meta-analysis comparable with that study (eMethods 1 in the Supplement). We also searched for databases from the countries listed in the articles selected for review, including the US National Institute for Occupational Safety and Health, the US National Outcomes Measurement System, and the UK Office of National Statistics (ONS) (eMethods 1 in the Supplement).

Inclusion criteria were original articles in English, Portuguese, Spanish, or French that assessed male and/or female physician suicide as an outcome. For the meta-analysis, we included only cohort periods from 1980 to the present. We also allowed for up to 1 year of data prior to 1980.

Exclusion criteria included studies lacking abstracts; nonoriginal articles; studies missing completed suicide data; studies without physicians’ suicide data or with only physicians of some medical specialties included; and articles covering similar periods plus or minus 2 years, if the source material was exactly the same, but not if additional source material was used (eg, an additional database or different US states).

For the meta-analysis, we excluded articles that did not compare male and female physician age-standardized suicide rates with the general population (observed/expected suicide rates), if this information could not be calculated from relevant sources (ie, the databases listed). We also excluded studies where cohort data could not be extracted from 1980 and after and studies with a substantial overlap (of 2 or more years) that used the same source material.

After excluding ineligible records, 4 authors (D.D., M.M.E., W.G., and T.C.C.) extracted data into standardized tables per a prespecified structured checklist. Three authors (D.D., M.M.E., and W.G.) further reviewed the extracted data, entered the relevant numbers into the final extraction form, and added the quality scores.

We replicated the methods of Schernhammer and Colditz,¹¹ including their quality scores. Additionally, we performed analyses based on approaches from recent studies.^22,23

Following the meta-analysis approach if the primary outcome (the SMR) was not provided by the articles, we calculated SMR values from the observed (O) and expected (E) suicide deaths (SMR = O/E, where O equaled total physician suicides divided by the total physician population and E equaled total population suicides divided by the total population, using the same age groups, locations, and year ranges; eMethods 2 in the Supplement). That is, SMRs on suicides in physicians within each country were based on the comparison with the general population of that country of the same age group.

The calculation for SMRs of the period being studied was the male or female physician suicide mortality rate per 100 000 person-years, divided by the suicide mortality rate of male or female individuals in the general population per 100 000 persons per year. If 95% CIs were missing, we estimated them using the logarithm of the proportions, assuming an approximate gaussian distribution. We also log-transformed the main SMR data using upper and lower CIs, accounting for the effect of potential extreme outliers, to check if the transformation changes the results.

We calculated all analyses using the statistical packages for meta-analysis of Stata version 15.1 (StataCorp). We tested for heterogeneity by giving commands for the Q test, I², and also τ² to estimate heterogeneity variance.²⁴ Values of P < .10 for the Q test and greater than 35% for I² were deemed indicative of study heterogeneity, according to the Cochrane Handbook.²⁵ We used random-effects models for the analysis of male and female physicians, for comparability and because we assumed that the cohort populations across time and countries have differences that could affect the true results; this model addresses effect-size distributions and is better for potentially heterogeneous populations.^26,27

The Begg and Egger tests,^28,29 sensitivity analyses, meta-regression (including before and after 1980), and cumulative meta-analysis were performed (eMethods 3 in the Supplement). We also present proportionate mortality ratios (PMRs) when SMR data could not be extracted, such as National Institute for Occupational Safety and Health–National Outcomes Measurement System (eMethods 3 in the Supplement).

Furthermore, to evaluate the effect of quality scores and duration of cohort on SMR, we performed a meta-regression with each variable in a separate model and then ran a model with all variables (with the SMR as the dependent variable and the above-mentioned variables as the independent variables) by using a modified model.³⁰ We also performed a meta-regression testing the effect of time (before and after 1980), sex, and the interaction term (time × sex) on SMR and on observed and expected suicide rates. For this meta-regression, pre-1980 data was extracted from Schernhammer and Colditz¹¹ and post-1980 data was mainly from our meta-analysis.

The initial search led to 7877 results. On applying our criteria, we ended up with 32 articles^8-10,31-61 (eFigure 1 and eTable 2 in the Supplement). We excluded 13 studies from the systematic review (eTable 1 in the Supplement) and meta-analysis.

Twenty-one articles^{8,9,31,32,34-48,54,61,62} showed significant differences between suicide rates in physicians and those of the general population; 11 articles^{10,49-53,55-57,61,62} showed none. Most articles separated physicians by sex, with 11 studies^8,31-40 and 9 articles^{8,33,35,36,39,41-44} showing significantly higher suicide risk in female and male physicians, respectively, than the general population. One article⁴⁵ and 3 articles^9,33,45 showed significantly lower rates among female and male physicians compared with the general population. Of articles that combined male and female individuals, 2 of 4 reported lower suicide rates among physicians,^45,46 and the other half reported higher rates^47,48 (eTable 2 in the Supplement).

Psychiatric illnesses were reported as risk factors in only 4 articles,^44,49-51 mainly depression, substance abuse, or both. Prior suicide attempts and suicidal ideation were not discussed.

Eight studies^{34,36,37,44,47,50-52} reported the most prominent demographic risk factors for suicide. These were sex (although results were mixed), followed by increasing age.^34,37,44,51 Additionally, there were mixed findings for race/ethnicity^36,44,47 and consistent findings of increased risk for interpersonal relationship factors,^47,50-52 especially being divorced,^47,51,52 but also being widowed,⁵¹ single,⁵² or having relationship problems.⁵⁰

Specialties were addressed in 7 articles.^{37-39,41,47,53,54} Physician specialties associated with higher risk of suicide included psychiatry^37,41,53,54 and anesthesiology,^39,41,54 followed by radiology,⁴¹ rehabilitation medicine,⁴¹ community health,⁵⁴ and general practice.⁵⁴ As to surgery, results were mixed.^38,53 One study⁴⁷ reported no differences between specialties.

Eight studies^{8,10,31,33,36,38,43,47} compared physician suicides with those of nonphysician occupations, with mixed results. Some found that physicians had higher suicide rates than other professions (including teachers, academics, and veterinarians),^{8,31,33,38,43,47} while other studies found increased risk in other professions (including dentists, pharmacists, and nurses).^10,31,33,43 Some studies reported on a specific race or sex^8,36,38,43 or found different occupation associations by sex.³¹

Work or training demands were indirectly hinted at by 6 studies.^{42,44,45,49,50,55} In the United States, male physicians mostly died by suicide when professionally productive,⁵⁵ during early training,⁴⁵ or following competitive schooling, when they had work problems, had work incapacity, or were working in isolation.⁴² Heavy workloads or isolation also contributed to suicidality.^44,49,50

Sixteen articles^{9,31,33,37,39-41,44,45,47,48,50-53,55} included at least partial descriptions of physician suicide methods (eTable 2 in the Supplement), typically by poison, firearms, or suffocation. In the United States, Brazil, and South Africa, the common suicide methods were firearms and poisoning and were modulated by geographic location and sex (eTable 2 in the Supplement).^{41,44,45,51,52} Poisoning was also the most common choice in Europe,^{33,37,39,48,53} Australia,⁵⁰ and New Zealand.⁹ Asphyxiation by hanging or other methods was often mentioned^{9,31,33,41,45,50,52,53} but was never reported as the most common method.

In our meta-analysis of 9 studies and data sets, there was a total of 547 male physician suicides and 162 female physician suicides (Table). Similar to Schernhammer and Colditz,¹¹ consistency in quality scoring between reviewers was 97%.

Male physicians’ SMR was significantly lower than the SMR of the general male population (ie, there were lower suicide rates in male physicians than men in general; 0.67 [95% CI, 0.55-0.79]; Figure 1A). Heterogeneity testing results included a Cochran Q test of 9.08 (P < .001), an I² of 11.9% (indicating low heterogeneity), and a τ²of 0.4%. Additionally, all studies’ 95% CIs overlapped.

To better evaluate heterogeneity, we then performed sensitivity analyses on the male physician cohort by removing individual studies, which did not change the direction or significance of the pooled results (Figure 2B). When analyzing the subset of studies with higher quality scores (>5), the forest plot and sensitivity analysis were similar (eFigure 2 and eFigure 3 in the Supplement). Further evaluating heterogeneity, 2 univariate meta-regressions each showed no significant association between quality score (β = −0.067; P = .41) and duration of cohort (β = −0.004; P = .85) with SMR. A multiple meta-regression model including quality score and duration of cohort as variables and SMR as an outcome was not significant. Finally, log-transformed data showed a significant pooled effect of −0.29 (95% CI, −0.43 to −0.14) in the same direction as the main SMR analysis.

Female physicians’ SMR was significantly higher than the SMR of the general female population (ie, higher suicide rates in female physicians than women in general; SMR, 1.46 [95% CI, 1.02-1.91]; Figure 1B). Heterogeneity testing results included a Cochran Q test of 16.02 (P = .04), an I² of 50.1% (moderate heterogeneity), and a τ² of 22.6%.

To better evaluate heterogeneity, we then performed a sensitivity analysis by removing individual studies (Figure 2C). We found that on excluding several articles,^{33,34,37,50,52} the 95% CI crossed the null; that is, the risk of suicide among female physicians became not significantly greater than that of the general female population. Additionally, higher-quality studies^{31,34,37,50,54,57-59,63} had a nonsignificant SMR for female physicians, and the sensitivity analysis has a pattern similar to that of main analysis (eFigure 4 and eFigure 5 in the Supplement). Further evaluating heterogeneity, 2 univariate meta-regressions each showed no significant association between quality score (β = −0.081; P = .76) and duration of cohort (β = −0.042; P = .47) with SMR. A multiple meta-regression model including quality score and duration of cohort as variables and SMR as an outcome was not significant. Finally, log-transformed data showed a pooled effect size of 0.43 (95% CI, 0.13-0.72), which was significant and in the same direction as the SMR.

There was no evidence of publication bias in studies on male and female physician suicides, as shown by respectively regressing suicide SMR with the inverse of study variance; in other words, the Egger test (intercept estimates: male physicians, −0.29; P = .78; female physicians, 0.85; P = .42), the Begg test (male physicians: continuity-corrected z score, 0.94; P = .35; female physicians: continuity-corrected z score, 0.63; P = .53), and the visual assessment of the funnel plot showed relative symmetry (Figure 2B).

We performed a cumulative meta-analysis (a statistical procedure to retrospectively calculate the sequence of meta-analyses by starting with a single study—the one with the earliest cohort—and adding the later studies one at a time); it shows how the overall estimate changes over time as each study is added to the pool. Calculating the cumulative evidence over time revealed significantly decreasing male physician SMR and increasing precision (effect sizes from 0.58 [95% CI, −0.20 to 1.36] to 0.67 [95% CI, 0.56-0.78]; eFigure 6 in the Supplement). This goes in the opposite direction of the previous cumulative meta-analysis.¹¹ The cumulative meta-analysis among female physicians showed an early increase (effect sizes from 0.62 [95% CI, -0.23 to 1.47] to 1.76 [95% CI, 1.34-2.18]), followed by a trend toward decreasing suicide risk (effect size, 1.43 [95% CI, 1.12-1.73]) that is still high compared with women in the general population (eFigure 7 in the Supplement).

In male physicians, the meta-regression of SMR over time (dichotomized by before and after 1980) led to a significant result of −0.84 (95% CI, −1.26 to −0.42; P < .001) and explained 36.9% (by adjusted R²) of the variance in the model. This represents an SMR decrease of 0.84 from the cohorts before 1980 to those after 1980.

In female physicians, the meta-regression of SMR over time (dichotomized to before and after 1980) led to a significant result of −1.96 (95% CI, −3.09 to −0.84; P = .002) and explained 43.9% (by adjusted R²) of the variance in the model. This represents an SMR decrease of 1.96.

Using a stepwise approach, we then tested a multivariate meta-regression model including time (before and after 1980) and sex (male and female physicians) and another model adding the interaction term time × sex. The meta-regression model better reflected the data and was significant (P < .001) with several terms (time, −0.85 [95% CI, −1.40 to –0.29]; P = .003; sex, 1.94 [95% CI, 1.25-2.62]; P < .001; interaction term time × sex, −1.06 [95% CI, −2.03 to −0.10]; P = .03). This is consistent with the greatest change being seen in female physicians after 1980 (Figure 3).

We performed separate meta-regressions of observed suicide rates in male and female physicians each, as well as separate meta-regressions of expected suicide rates in men and women of the general population each. The meta-regressions help us better understand whether the SMR changes over time (before and after 1980) were driven by observed suicide rates in physicians (male and female) or alternatively by expected suicide rates in men and women of the general population.

We found that on meta-regression assessing the association of time with expected suicide rates in men in the general population and the association of time with observed suicide rates in male physicians, both meta-regression results were in the direction of decreased suicides after 1980, but this was only significant in male physicians (β = −35.37 [95% CI, −62.33 to −8.41]; P = .01), and time explained 20.21% of the variance in the data. Meanwhile, the results were nonsignificant in each of the meta-regressions for expected suicide rates among women in the general population) and observed suicide rates in female physicians (eMethods 5 in the Supplement).

Physician SMRs compare observed suicide rates (suicides in physicians divided by physician populations) to expected suicide rates (suicides in the general population divided by the general population). Conversely, physician-suicide PMRs indicate the ratio of physician-suicide proportionate mortality (ie, suicides divided by all-cause mortality within physician populations) to general population–suicide proportionate mortality (ie, suicides divided by all-cause mortality within general populations). In the United States (per National Outcomes Measurement System–National Institute for Occupational Safety and Health data), physician PMRs respectively increased in white male physicians and especially black male physicians from the first recorded time span (1985-1998) to the second (roughly 1999-2014) (white male physicians, 1.77 [95% CI, 1.60-1.96] vs 2.03 [95% CI, 1.81-2.27]; black male physicians, 2.50 [95% CI, 1.08-4.92] vs 4.24 [95% CI, 2.12-7.59]) but decreased for white female physicians during those periods (2.66 [95% CI, 1.92-3.60] vs 2.42 [95% CI, 1.80-3.18]) (Table). (There were insufficient suicide data on black female physicians.) The US PMR data indicate that suicide makes up a higher proportion of all-cause mortality within all these physician groups than it does for the comparable general populations of the same ages, sexes, and races—that is, suicide causes more of the deaths within the physician populations than within the general population.

As to the English data set (2011-2015; Table), the suicide PMR was also high for male physicians (1.85 [95% CI, 1.41-2.38]) and higher for female physicians (1.97 [95% CI, 1.23-2.98]), which was significant in both cases despite the much lower and nonsignificant risks suggested by their suicide SMRs (0.63 [95% CI, 0.43-2.28] and 1.01 [95% CI, 0.44-2.31] for male and female physicians, respectively). This apparent discrepancy within the same sample reflects the differences in how SMRs and PMRs are calculated.

Suicide rates in male physicians (observed rate, 12.4 per 100 000 population) were lower than rates in men in the general population (expected rate, 18.9 per 100 000 population), which was consistent with the low SMR for male physicians and implied that the elevated PMR was driven by relatively lower all-cause mortality rather than more suicides in male physicians compared with men in general. Likewise, female physicians and women in general had similar suicide rates (5.5 per 100 000 population and 5.4 per 100 000 population, respectively); therefore, the elevated PMR in female physicians may be driven by lower all-cause mortality relative to women in general.

The latest meta-analysis of physician suicide was published in 2004¹¹ and included cohort data mainly before 1980. That study found suicide SMRs in male physicians (1.41 [95% CI, 1.21-1.65]) and especially female physicians (2.27 [95% CI, 1.90-2.73]) significantly higher than the expected population suicide rates. We had therefore hypothesized a higher risk of physician suicides compared with the general population from the 1980s to date, and our data confirmed a high SMR for female physicians (1.46) but a low SMR for male physicians (0.67). We had also hypothesized that the suicide SMRs of both female and male physicians seen before 1980 would decrease after 1980, and this was confirmed by the models; additionally, female physicians had the greatest decrease in SMR after 1980. These population-based findings are of great relevance, because physician suicide is an important and understudied topic.

Our meta-regression suggested that the decrease in male physician suicide SMR over time was driven by the rate of physician suicides rather than population suicides. This may be because male physicians are relatively protected from workforce or unemployment factors affecting men of lower socioeconomic status, which may have obscured any burnout-associated outcomes. Conversely, we could not make the same claim for female physicians vs women in the general population, potentially because the pre-1980 data was underpowered. Yet the SMR for female physicians, while decreasing, was still high, implying that greater female representation in the physician workforce may not have overcome the magnitude of their increased risk compared with women in general. It also raises the question of a lag in improved workforce conditions compared with workforce numbers.

Meanwhile, the PMR and SMR data within a single sample provide different views of physician suicides. For example, SMR data suggest that male physicians in England are at lower risk of suicide than the male population while female physicians had risk comparable with that of women in general, yet the PMR data suggest that suicide makes up a relatively large proportion of mortality in both male and female physicians, perhaps because physicians die less frequently of other (eg, cardiovascular, pulmonary) causes. This is consistent with studies in Denmark⁴⁸ and Maryland⁴¹ suggesting lower overall physician mortality yet greater suicide mortality compared with what would be expected. As to the United States, we find that PMRs are high in black and white male physicians and white female physicians. It is important to highlight that SMRs and PMRs are different health indicators with potentially different applications—for example, a physician might care more about risk of suicide compared with dying of other causes, while a policy maker might prioritize how suicide risk in physicians compares with that of the general population. Overall, suicide is among the top 10 mortality causes in physicians and general population per the US Centers for Disease Control and Prevention and National Institute for Mental Health.⁶⁴ We further discuss important risk factors.

Generally speaking, major risk factors for suicide in the population are previous suicide attempts; a family history of suicide; mental disorders; homosexuality; male sex; divorced, widowed, or single marital or relationship status; young adult or elderly age ranges; acute stressors; and childhood trauma.^7,65-67 However, we only found reports associated with sex, relationship status, and age factors. Suicide methods were by poisoning or drugs, firearms, and suffocation, but the frequency varied by geographic location (and accessibility), sex, and other factors. In this review, career-associated factors seemed to be the most prominent risk, while psychiatric illnesses seemed to be underreported.^52,68,69

The most common disorders in this review were depression and/or substance abuse. Unfortunately, physicians are less likely to seek mental health services⁷⁰ out of career concerns, culture, and/or a predisposition toward self-reliance. Additionally, retrospective toxicology screening of suicide data finds that physicians are more likely than nonphysicians to have positive results for antipsychotics, benzodiazepines, and barbiturates but not antidepressants.⁶⁸ This links to physicians being undertreated for depression despite current effective treatments,⁷¹ whereas the positive toxicology findings may indicate substance abuse, self-medication, or intentional toxicity.

Psychiatry and anesthesiology were the top 2 specialties at risk in this review. However, we excluded studies analyzing only specialist data, some of which reported higher suicide rates for anesthesiologists, typically by drugs.^72-75 It is well established that anesthesiologists tend to have much higher rates of substance abuse than other physicians.⁷⁶ Psychiatrists are also known to have more mental distress, mental illness, and burnout compared with other physician groups, and 1 systematic review⁷⁷ found concerning rates of depression, psychotropic drug use, and “almost-clinical-level emotional disturbance caused by patient suicide.”^77(p212) Our review hints at career or workload burden contributions to suicidality.^39,43,51 Some associated factors could be the judicialization of medicine, organization oversight, and burnout, topics that are beyond the scope of this article.

Some of the most consistent physician suicide risk factors were increasing age, and divorce or relationships, although marriage is reportedly protective in men and not clearly so in women or female physicians.⁷⁰ Additionally, there is some suggestion of an association between sex, race, and suicidality (eg, an elevated risk in black male physicians).

Physician suicide is suboptimally reported and addressed; our meta-analysis showed decreasing suicide rates, but the data sets showed paradoxical trends that may be affected by gender, race, location, and time. What is driving the data? Studies need to be systematized to clarify underlying patterns, identify vulnerable subsets of physicians at risk, and explore suicide prevention by addressing the modifiable risk factors suggested by this research.

A probable bias in physician suicidality studies is underreporting, because colleagues may be unwilling to classify deaths as suicides.^38,43,53 For example, most unintentional poisonings involved prescription drugs, yet physicians probably would not unintentionally overdose more than the population.⁵³ Another limitation is broad study locations, time span, and methods (eg, different International Classification of Diseases codes and data sources). However, quality scores, analysis of heterogeneity, publication bias, and meta-regression were used to increase validity.

Similar to the previous meta-analysis,¹¹ our study primarily reflects northern European and North American populations (eMethods 6 in the Supplement). The National Outcomes Measurement System data sets do not include races other than black and white. However, other racial/ethnic groups, particularly Asian populations, are more highly represented in medicine than in the general population.¹²

Our population-based meta-analysis showed suicide SMRs that were high in female physicians and low in male physicians after 1980, and both SMRs decreased over time (from before to after 1980). Nevertheless, suicide PMRs highlight the importance of suicide relative to other causes of mortality in physicians. This qualitative and quantitative work illustrates physician suicide patterns and characteristics for future studies and policies.

Accepted for Publication: December 12, 2019.

Corresponding Author: Dante Duarte, MD, PhD, Neuromodulation Center, Spaulding Rehabilitation Hospital, Department of Physical Medicine and Rehabilitation, Harvard Medical School, 96 13th St, Boston, MA 02129 (dantegalileu@gmail.com).

Published Online: March 4, 2020. doi:10.1001/jamapsychiatry.2020.0011

Author Contributions: Drs Duarte and El-Hagrassy had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Drs Duarte and El-Hagrassy are equally contributing authors.

Concept and design: Duarte, El-Hagrassy, Couto, Correa.

Acquisition, analysis, or interpretation of data: Duarte, El-Hagrassy, Couto, Gurgel, Fregni.

Drafting of the manuscript: Duarte, El-Hagrassy, Couto, Gurgel.

Critical revision of the manuscript for important intellectual content: Duarte, El-Hagrassy, Gurgel, Fregni, Correa.

Statistical analysis: Duarte, El-Hagrassy, Gurgel, Fregni.

Obtained funding: Correa.

Administrative, technical, or material support: Duarte, Correa.

Supervision: Duarte, Correa.

Conflict of Interest Disclosures: None reported.

Funding/Support: This work was conducted with statistical consultation support from Harvard Catalyst/The Harvard Clinical and Translational Science Center (National Center for Advancing Translational Sciences/National Institutes of Health award UL 1TR002541) and financial contributions from Harvard University and its affiliated academic health care centers.

Role of the Funder/Sponsor: The sponsor 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.

Disclaimer: The content is solely the responsibility of the authors and does not necessarily represent the official views of Harvard Catalyst, Harvard University and its affiliated academic health care centers, the National Institutes of Health, or Federal University of Minas Gerais.

Additional Contributions: We acknowledge Kevin Pacheco-Barrios, MD, Spaulding Neuromodulation Center, for his help on the graphs and figures. He was not compensated for his contributions.