PT  - JOURNAL ARTICLE
AU  - Kandula, Sasikiran
AU  - Martinez-Alés, Gonzalo
AU  - Rutherford, Caroline
AU  - Gimbrone, Catherine
AU  - Olfson, Mark
AU  - Gould, Madelyn S.
AU  - Keyes, Katherine M.
AU  - Shaman, Jeffrey
TI  - Modeling suicide mortality in US counties using population socioeconomic indicators
AID  - 10.1101/2022.06.06.22275887
DP  - 2022 Jan 01
TA  - medRxiv
PG  - 2022.06.06.22275887
4099  - http://medrxiv.org/content/early/2022/06/06/2022.06.06.22275887.short
4100  - http://medrxiv.org/content/early/2022/06/06/2022.06.06.22275887.full
AB  - Background Suicide is one of the leading causes of death in the United States and population risk prediction models can inform the type, location, and timing of public health interventions. Here, we report the development of a prediction model of suicide risk using population characteristics.Methods All suicide deaths reported to the Nation Vital Statistics System between 2005-2019 were identified, and age, sex, race, and county-of-residence of the decedents were extracted to calculate baseline risk. County-wise annual measures of socioeconomic predictors of suicide risk — unemployment, weekly wage, poverty prevalence, median household income, and population density — along with two state-wise measures of prevalence of major depressive disorder and firearm ownership were compiled from public sources. Conditional autoregressive (CAR) models, which account for spatiotemporal autocorrelation in response and predictors, were used to estimate county-level risk.Results Estimates derived from CAR models were more accurate than from models not adjusted for spatiotemporal autocorrelation. Inclusion of suicide risk/protective covariates further reduced errors. Suicide risk was estimated to increase with each standard deviation increase in firearm ownership (2.8%), prevalence of major depressive episode (1%) and unemployment (2.8%). Conversely, risk was estimated to decrease by 4.3% for each standard deviation increase in both median household income and population density. Increased heterogeneity of risk across counties was also noted.Conclusions Area-level characteristics and the CAR model structure can estimate population-level suicide risk and thus inform decisions on resource allocation and focused interventions during outbreaks.Competing Interest StatementKeyes and Rutherford have been compensated for expert witness work in litigation; Shaman and Columbia University declare partial ownership of SK Analytics; Shaman was a consultant for BNI; Kandula was a consultant for SK Analytics. Martinez-Alés, Gimbrone, Gould and Olfson have no disclosures.Funding StatementThis work is funded by a grant from the National Institute of Mental Health (R01-MH121410)Author DeclarationsI confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.YesI confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.YesI understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).YesI have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.YesCovariate data were gathered from public datasets and references to the primary sources have been provided in the manuscript. Mortality data was acquired under a restricted use agreement that does not allow sharing data publicly.