1
|
Shaman J, Kandula S, Pei S, Galanti M, Olfson M, Gould M, Keyes K. Quantifying suicide contagion at population scale. SCIENCE ADVANCES 2024; 10:eadq4074. [PMID: 39083618 PMCID: PMC11290520 DOI: 10.1126/sciadv.adq4074] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2024] [Accepted: 06/27/2024] [Indexed: 08/02/2024]
Abstract
The spread of suicidal behavior among individuals is often described as a contagion; however, rigorous modeling of suicide as a dynamic, contagious process is minimal. Here, we develop and validate a model-inference system depicting suicide ideation and death and use it to quantify the contagion processes in the US associated with two prominent celebrity suicide events: Robin Williams during 2014 and Kate Spade and Anthony Bourdain, which occurred 3 days apart during 2018. We show that both events produced large transient increases of suicide contagion contact rates, i.e., the spread of suicidal thought and behavior, and a period of elevated suicidal ideation in the general population. Our modeling approach provides a framework for quantifying suicidal contagion and better understanding, preventing, and containing its spread.
Collapse
Affiliation(s)
- Jeffrey Shaman
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
- Columbia Climate School, Columbia University, New York, NY 10025, USA
| | - Sasikiran Kandula
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
| | - Sen Pei
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
| | - Marta Galanti
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
| | - Mark Olfson
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
- Department of Psychiatry, Columbia University, New York, NY 10032, USA
| | - Madelyn Gould
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
- Department of Psychiatry, Columbia University, New York, NY 10032, USA
| | - Katherine Keyes
- Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, NY 10032, USA
| |
Collapse
|
2
|
Fonseca LL, Böttcher L, Mehrad B, Laubenbacher RC. Surrogate modeling and control of medical digital twins. ARXIV 2024:arXiv:2402.05750v2. [PMID: 38827450 PMCID: PMC11142319] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 06/04/2024]
Abstract
The vision of personalized medicine is to identify interventions that maintain or restore a person's health based on their individual biology. Medical digital twins, computational models that integrate a wide range of health-related data about a person and can be dynamically updated, are a key technology that can help guide medical decisions. Such medical digital twin models can be high-dimensional, multi-scale, and stochastic. To be practical for healthcare applications, they often need to be simplified into low-dimensional surrogate models that can be used for optimal design of interventions. This paper introduces surrogate modeling algorithms for the purpose of optimal control applications. As a use case, we focus on agent-based models (ABMs), a common model type in biomedicine for which there are no readily available optimal control algorithms. By deriving surrogate models that are based on systems of ordinary differential equations, we show how optimal control methods can be employed to compute effective interventions, which can then be lifted back to a given ABM. The relevance of the methods introduced here extends beyond medical digital twins to other complex dynamical systems.
Collapse
Affiliation(s)
- Luis L. Fonseca
- Laboratory for Systems Medicine, Department of Medicine, University of Florida, Gainesville, FL, USA
| | - Lucas Böttcher
- Laboratory for Systems Medicine, Department of Medicine, University of Florida, Gainesville, FL, USA
- Department of Computational Science and Philosophy, Frankfurt School of Finance and Management, 60322 Frankfurt am Main, Germany
| | - Borna Mehrad
- Laboratory for Systems Medicine, Department of Medicine, University of Florida, Gainesville, FL, USA
| | - Reinhard C. Laubenbacher
- Laboratory for Systems Medicine, Department of Medicine, University of Florida, Gainesville, FL, USA
| |
Collapse
|
3
|
He K, Foerster S, Vora NM, Blaney K, Keeley C, Hendricks L, Varma JK, Long T, Shaman J, Pei S. Evaluating completion rates of COVID-19 contact tracing surveys in New York City. BMC Public Health 2024; 24:414. [PMID: 38331772 PMCID: PMC10854191 DOI: 10.1186/s12889-024-17920-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Accepted: 01/29/2024] [Indexed: 02/10/2024] Open
Abstract
IMPORTANCE Contact tracing is the process of identifying people who have recently been in contact with someone diagnosed with an infectious disease. During an outbreak, data collected from contact tracing can inform interventions to reduce the spread of infectious diseases. Understanding factors associated with completion rates of contact tracing surveys can help design improved interview protocols for ongoing and future programs. OBJECTIVE To identify factors associated with completion rates of COVID-19 contact tracing surveys in New York City (NYC) and evaluate the utility of a predictive model to improve completion rates, we analyze laboratory-confirmed and probable COVID-19 cases and their self-reported contacts in NYC from October 1st 2020 to May 10th 2021. METHODS We analyzed 742,807 case investigation calls made during the study period. Using a log-binomial regression model, we examined the impact of age, time of day of phone call, and zip code-level demographic and socioeconomic factors on interview completion rates. We further developed a random forest model to predict the best phone call time and performed a counterfactual analysis to evaluate the change of completion rates if the predicative model were used. RESULTS The percentage of contact tracing surveys that were completed was 79.4%, with substantial variations across ZIP code areas. Using a log-binomial regression model, we found that the age of index case (an individual who has tested positive through PCR or antigen testing and is thus subjected to a case investigation) had a significant effect on the completion of case investigation - compared with young adults (the reference group,24 years old < age < = 65 years old), the completion rate for seniors (age > 65 years old) were lower by 12.1% (95%CI: 11.1% - 13.3%), and the completion rate for youth group (age < = 24 years old) were lower by 1.6% (95%CI: 0.6% -2.6%). In addition, phone calls made from 6 to 9 pm had a 4.1% (95% CI: 1.8% - 6.3%) higher completion rate compared with the reference group of phone calls attempted from 12 and 3 pm. We further used a random forest algorithm to assess its potential utility for selecting the time of day of phone call. In counterfactual simulations, the overall completion rate in NYC was marginally improved by 1.2%; however, certain ZIP code areas had improvements up to 7.8%. CONCLUSION These findings suggest that age and time of day of phone call were associated with completion rates of case investigations. It is possible to develop predictive models to estimate better phone call time for improving completion rates in certain communities.
Collapse
Affiliation(s)
- Kaiyu He
- Department of Biostatistics, Mailman School of Public Health, Columbia University, New York, NY, 10032, USA
| | - Steffen Foerster
- New York City Department of Health and Mental Hygiene (DOHMH), Long Island City, NY, 11001, USA
| | - Neil M Vora
- New York City Department of Health and Mental Hygiene (DOHMH), Long Island City, NY, 11001, USA
| | - Kathleen Blaney
- New York City Department of Health and Mental Hygiene (DOHMH), Long Island City, NY, 11001, USA
| | | | | | - Jay K Varma
- Department of Population Health Sciences, Weill Cornell Medical College, New York, NY, 10065, USA
| | - Theodore Long
- NYC Health + Hospitals, New York, NY, USA
- Department of Population Health, New York University, New York, NY, 10016, USA
| | - Jeffrey Shaman
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY, 10032, USA
- Columbia Climate School, Columbia University, New York, NY, 10025, USA
| | - Sen Pei
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, NY, 10032, USA.
| |
Collapse
|
4
|
Böttcher L, Chou T, D’Orsogna MR. Forecasting drug-overdose mortality by age in the United States at the national and county levels. PNAS NEXUS 2024; 3:pgae050. [PMID: 38725534 PMCID: PMC11079616 DOI: 10.1093/pnasnexus/pgae050] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Accepted: 01/25/2024] [Indexed: 05/12/2024]
Abstract
The drug-overdose crisis in the United States continues to intensify. Fatalities have increased 5-fold since 1999 reaching a record high of 108,000 deaths in 2021. The epidemic has unfolded through distinct waves of different drug types, uniquely impacting various age, gender, race, and ethnic groups in specific geographical areas. One major challenge in designing interventions and efficiently delivering treatment is forecasting age-specific overdose patterns at the local level. To address this need, we develop a forecasting method that assimilates observational data obtained from the CDC WONDER database with an age-structured model of addiction and overdose mortality. We apply our method nationwide and to three select areas: Los Angeles County, Cook County, and the five boroughs of New York City, providing forecasts of drug-overdose mortality and estimates of relevant epidemiological quantities, such as mortality and age-specific addiction rates.
Collapse
Affiliation(s)
- Lucas Böttcher
- Department of Computational Science and Philosophy, Frankfurt School of Finance and Management, 60322 Frankfurt am Main, Germany
| | - Tom Chou
- Department of Computational Medicine, University of California, Los Angeles, Los Angeles, CA 90095-1766, USA
| | - Maria R D’Orsogna
- Department of Computational Medicine, University of California, Los Angeles, Los Angeles, CA 90095-1766, USA
- Department of Mathematics, California State University at Northridge, Los Angeles, CA 91330-8313, USA
| |
Collapse
|
5
|
Zhang R, Tai J, Pei S. Ensemble inference of unobserved infections in networks using partial observations. PLoS Comput Biol 2023; 19:e1011355. [PMID: 37549190 PMCID: PMC10434926 DOI: 10.1371/journal.pcbi.1011355] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2023] [Revised: 08/17/2023] [Accepted: 07/12/2023] [Indexed: 08/09/2023] Open
Abstract
Undetected infections fuel the dissemination of many infectious agents. However, identification of unobserved infectious individuals remains challenging due to limited observations of infections and imperfect knowledge of key transmission parameters. Here, we use an ensemble Bayesian inference method to infer unobserved infections using partial observations. The ensemble inference method can represent uncertainty in model parameters and update model states using all ensemble members collectively. We perform extensive experiments in both model-generated and real-world networks in which individuals have differential but unknown transmission rates. The ensemble method outperforms several alternative approaches for a variety of network structures and observation rates, despite that the model is mis-specified. Additionally, the computational complexity of this algorithm scales almost linearly with the number of nodes in the network and the number of observations, respectively, exhibiting the potential to apply to large-scale networks. The inference method may support decision-making under uncertainty and be adapted for use for other dynamical models in networks.
Collapse
Affiliation(s)
- Renquan Zhang
- School of Mathematical Sciences, Dalian University of Technology, Dalian, China
| | - Jilei Tai
- School of Mathematical Sciences, Dalian University of Technology, Dalian, China
| | - Sen Pei
- Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, New York, United States of America
| |
Collapse
|
6
|
Huang S, Sun J, Feng L, Xie J, Wang D, Hu Y. Identify hidden spreaders of pandemic over contact tracing networks. Sci Rep 2023; 13:11621. [PMID: 37468540 DOI: 10.1038/s41598-023-32542-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2022] [Accepted: 03/29/2023] [Indexed: 07/21/2023] Open
Abstract
The COVID-19 infection cases have surged globally, causing devastations to both the society and economy. A key factor contributing to the sustained spreading is the presence of a large number of asymptomatic or hidden spreaders, who mix among the susceptible population without being detected or quarantined. Due to the continuous emergence of new virus variants, even if vaccines have been widely used, the detection of asymptomatic infected persons is still important in the epidemic control. Based on the unique characteristics of COVID-19 spreading dynamics, here we propose a theoretical framework capturing the transition probabilities among different infectious states in a network, and extend it to an efficient algorithm to identify asymptotic individuals. We find that using pure physical spreading equations, the hidden spreaders of COVID-19 can be identified with remarkable accuracy, even with incomplete information of the contract-tracing networks. Furthermore, our framework can be useful for other epidemic diseases that also feature asymptomatic spreading.
Collapse
Affiliation(s)
- Shuhong Huang
- School of Data and Computer Science, Sun Yat-sen University, Guangzhou, 510006, China
- Institute of Neuroscience, Technical University of Munich, Munich, 80802, Germany
| | | | - Ling Feng
- Institute of High Performance Computing, Agency for Science, Technology and Research (A*STAR), Singapore, 138632, Singapore
- Department of Physics, National University of Singapore, Singapore, 117551, Singapore
| | - Jiarong Xie
- School of Data and Computer Science, Sun Yat-sen University, Guangzhou, 510006, China
| | - Dashun Wang
- Kellogg School of Management, Northwestern University, Evanston, IL, USA
| | - Yanqing Hu
- Department of Statistics and Data Science, College of Science, Southern University of Science and Technology, 518055, Shenzhen, China.
| |
Collapse
|
7
|
Böttcher L, Chou T, D’Orsogna MR. Modeling and forecasting age-specific drug overdose mortality in the United States. THE EUROPEAN PHYSICAL JOURNAL. SPECIAL TOPICS 2023; 232:1-10. [PMID: 37359186 PMCID: PMC10132445 DOI: 10.1140/epjs/s11734-023-00801-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/30/2022] [Accepted: 02/24/2023] [Indexed: 06/28/2023]
Abstract
Drug overdose deaths continue to increase in the United States for all major drug categories. Over the past two decades the total number of overdose fatalities has increased more than fivefold; since 2013 the surge in overdose rates is primarily driven by fentanyl and methamphetamines. Different drug categories and factors such as age, gender, and ethnicity are associated with different overdose mortality characteristics that may also change in time. For example, the average age at death from a drug overdose has decreased from 1940 to 1990 while the overall mortality rate has steadily increased. To provide insight into the population-level dynamics of drug overdose mortality, we develop an age-structured model for drug addiction. Using an augmented ensemble Kalman filter (EnKF), we show through a simple example how our model can be combined with synthetic observation data to estimate mortality rate and an age-distribution parameter. Finally, we use an EnKF to combine our model with observation data on overdose fatalities in the United States from 1999 to 2020 to forecast the evolution of overdose trends and estimate model parameters.
Collapse
Affiliation(s)
- Lucas Böttcher
- Department of Computational Science and Philosophy, Frankfurt School of Finance and Management, 60322 Frankfurt am Main, Germany
| | - Tom Chou
- Department of Computational Medicine, University of California, Los Angeles, Los Angeles, 90095 CA USA
- Department of Mathematics, University of California, Los Angeles, Los Angeles, 90095 CA USA
| | - Maria R. D’Orsogna
- Department of Computational Medicine, University of California, Los Angeles, Los Angeles, 90095 CA USA
- Department of Mathematics, California State University at Northridge, Los Angeles, 91330 CA USA
| |
Collapse
|
8
|
Masoomy H, Chou T, Böttcher L. Impact of random and targeted disruptions on information diffusion during outbreaks. CHAOS (WOODBURY, N.Y.) 2023; 33:033145. [PMID: 37003816 DOI: 10.1063/5.0139844] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/23/2022] [Accepted: 03/02/2023] [Indexed: 06/19/2023]
Abstract
Outbreaks are complex multi-scale processes that are impacted not only by cellular dynamics and the ability of pathogens to effectively reproduce and spread, but also by population-level dynamics and the effectiveness of mitigation measures. A timely exchange of information related to the spread of novel pathogens, stay-at-home orders, and other measures can be effective at containing an infectious disease, particularly during the early stages when testing infrastructure, vaccines, and other medical interventions may not be available at scale. Using a multiplex epidemic model that consists of an information layer (modeling information exchange between individuals) and a spatially embedded epidemic layer (representing a human contact network), we study how random and targeted disruptions in the information layer (e.g., errors and intentional attacks on communication infrastructure) impact the total proportion of infections, peak prevalence (i.e., the maximum proportion of infections), and the time to reach peak prevalence. We calibrate our model to the early outbreak stages of the SARS-CoV-2 pandemic in 2020. Mitigation campaigns can still be effective under random disruptions, such as failure of information channels between a few individuals. However, targeted disruptions or sabotage of hub nodes that exchange information with a large number of individuals can abruptly change outbreak characteristics, such as the time to reach the peak of infection. Our results emphasize the importance of the availability of a robust communication infrastructure during an outbreak that can withstand both random and targeted disruptions.
Collapse
Affiliation(s)
- Hosein Masoomy
- Department of Physics, Shahid Beheshti University, 1983969411 Tehran, Iran
| | - Tom Chou
- Department of Computational Medicine and Department of Mathematics, UCLA, Los Angeles, California 90095, USA
| | - Lucas Böttcher
- Department of Computational Science and Philosophy, Frankfurt School of Finance and Management, 60322 Frankfurt am Main, Germany
| |
Collapse
|