COVID-19 dynamics in an Ohio prison.

Incarcerated individuals are a highly vulnerable population for infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Understanding the transmission of respiratory infections within prisons and between prisons and surrounding communities is a crucial component of pandemic preparedness and response. Here, we use mathematical and statistical models to analyze publicly available data on the spread of SARS-CoV-2 reported by the Ohio Department of Rehabilitation and Corrections (ODRC). Results from mass testing conducted on April 16, 2020 were analyzed together with time of first reported SARS-CoV-2 infection among Marion Correctional Institution (MCI) inmates. Extremely rapid, widespread infection of MCI inmates was reported, with nearly 80% of inmates infected within 3 weeks of the first reported inmate case. The dynamical survival analysis (DSA) framework that we use allows the derivation of explicit likelihoods based on mathematical models of transmission. We find that these data are consistent with three non-exclusive possibilities: (i) a basic reproduction number >14 with a single initially infected inmate, (ii) an initial superspreading event resulting in several hundred initially infected inmates with a reproduction number of approximately three, or (iii) earlier undetected circulation of virus among inmates prior to April. All three scenarios attest to the vulnerabilities of prisoners to COVID-19, and the inability to distinguish among these possibilities highlights the need for improved infection surveillance and reporting in prisons.

Social Media Sentiment about COVID-19 Vaccination Predicts Vaccine Acceptance among Peruvian Social Media Users the Next Day.

Drawing upon theories of risk and decision making, we present a theoretical framework for how the emotional attributes of social media content influence risk behaviors. We apply our framework to understanding how COVID-19 vaccination Twitter posts influence acceptance of the vaccine in Peru, the country with the highest relative number of COVID-19 excess deaths. By employing computational methods, topic modeling, and vector autoregressive time series analysis, we show that the prominence of expressed emotions about COVID-19 vaccination in social media content is associated with the daily percentage of Peruvian social media survey respondents who are vaccine-accepting over 231 days. Our findings show that net (positive) sentiment and trust emotions expressed in tweets about COVID-19 are positively associated with vaccine acceptance among survey respondents one day after the post occurs. This study demonstrates that the emotional attributes of social media content, besides veracity or informational attributes, may influence vaccine acceptance for better or worse based on its valence.

Projecting COVID-19 cases and hospital burden in Ohio.

As the Coronavirus 2019 disease (COVID-19) started to spread rapidly in the state of Ohio, the Ecology, Epidemiology and Population Health (EEPH) program within the Infectious Diseases Institute (IDI) at The Ohio State University (OSU) took the initiative to offer epidemic modeling and decision analytics support to the Ohio Department of Health (ODH). This paper describes the methodology used by the OSU/IDI response modeling team to predict statewide cases of new infections as well as potential hospital burden in the state. The methodology has two components: (1) A Dynamical Survival Analysis (DSA)-based statistical method to perform parameter inference, statewide prediction and uncertainty quantification. (2) A geographic component that down-projects statewide predicted counts to potential hospital burden across the state. We demonstrate the overall methodology with publicly available data. A Python implementation of the methodology is also made publicly available. This manuscript was submitted as part of a theme issue on "Modelling COVID-19 and Preparedness for Future Pandemics".

Catch the tweet to fight the flu: Using Twitter to promote flu shots on a college campus.

Objective: Over the 2018-2019 flu season we conducted a randomized controlled trial examining the efficacy of a Twitter campaign on vaccination rates. Concurrently we investigated potential interactions between digital social network structure and vaccination status. Participants: Undergratuates at a large midwestern public university were randomly assigned to an intervention (n = 353) or control (n = 349) group. Methods: Vaccination data were collected via monthly surveys. Participant Twitter data were collected through the public-facing Twitter API. Intervention impact was assessed with logistic regression. Standard network science tools examined vaccination coverage over online social networks. Results: The campaign had no effect on vaccination outcome. Receiving a flu shot the prior year had a positive impact on participant vaccination. Evidence of an interaction between digital social network structure and vaccination status was detected. Conclusions: Social media campaigns may not be sufficient for increasing vaccination rates. There may be potential for social media campaigns that leverage network structure.

Host movement, transmission hot spots, and vector-borne disease dynamics on spatial networks.

We examine how spatial heterogeneity combines with mobility network structure to influence vector-borne disease dynamics. Specifically, we consider a Ross-Macdonald-type disease model on n spatial locations that are coupled by host movement on a strongly connected, weighted, directed graph. We derive a closed form approximation to the domain reproduction number using a Laurent series expansion, and use this approximation to compute sensitivities of the basic reproduction number to model parameters. To illustrate how these results can be used to help inform mitigation strategies, as a case study we apply these results to malaria dynamics in Namibia, using published cell phone data and estimates for local disease transmission. Our analytical results are particularly useful for understanding drivers of transmission when mobility sinks and transmission hot spots do not coincide.

Projecting COVID-19 Cases and Subsequent Hospital Burden in Ohio.

As the Coronavirus 2019 (COVID-19) disease started to spread rapidly in the state of Ohio, the Ecology, Epidemiology and Population Health (EEPH) program within the Infectious Diseases Institute (IDI) at the Ohio State University (OSU) took the initiative to offer epidemic modeling and decision analytics support to the Ohio Department of Health (ODH). This paper describes the methodology used by the OSU/IDI response modeling team to predict statewide cases of new infections as well as potential hospital burden in the state. The methodology has two components: 1) A Dynamic Survival Analysis (DSA)-based statistical method to perform parameter inference, statewide prediction and uncertainty quantification. 2) A geographic component that down-projects statewide predicted counts to potential hospital burden across the state. We demonstrate the overall methodology with publicly available data. A Python implementation of the methodology has been made available publicly. Highlights: We present a novel statistical approach called Dynamic Survival Analysis (DSA) to model an epidemic curve with incomplete data. The DSA approach is advantageous over standard statistical methods primarily because it does not require prior knowledge of the size of the susceptible population, the overall prevalence of the disease, and also the shape of the epidemic curve.The principal motivation behind the study was to obtain predictions of case counts of COVID-19 and the resulting hospital burden in the state of Ohio during the early phase of the pandemic.The proposed methodology was applied to the COVID-19 incidence data in the state of Ohio to support the Ohio Department of Health (ODH) and the Ohio Hospital Association (OHA) with predictions of hospital burden in each of the Hospital Catchment Areas (HCAs) of the state.

Relating Eulerian and Lagrangian spatial models for vector-host disease dynamics through a fundamental matrix.

We explore the relationship between Eulerian and Lagrangian approaches for modeling movement in vector-borne diseases for discrete space. In the Eulerian approach we account for the movement of hosts explicitly through movement rates captured by a graph Laplacian matrix L. In the Lagrangian approach we only account for the proportion of time that individuals spend in foreign patches through a mixing matrix P. We establish a relationship between an Eulerian model and a Lagrangian model for the hosts in terms of the matrices L and P. We say that the two modeling frameworks are consistent if for a given matrix P, the matrix L can be chosen so that the residence times of the matrix P and the matrix L match. We find a sufficient condition for consistency, and examine disease quantities such as the final outbreak size and basic reproduction number in both the consistent and inconsistent cases. In the special case of a two-patch model, we observe how similar values for the basic reproduction number and final outbreak size can occur even in the inconsistent case. However, there are scenarios where the final sizes in both approaches can significantly differ by means of the relationship we propose.

Sentinel node approach to monitoring online COVID-19 misinformation.

Understanding how different online communities engage with COVID-19 misinformation is critical for public health response. For example, misinformation confined to a small, isolated community of users poses a different public health risk than misinformation being consumed by a large population spanning many diverse communities. Here we take a longitudinal approach that leverages tools from network science to study COVID-19 misinformation on Twitter. Our approach provides a means to examine the breadth of misinformation engagement using modest data needs and computational resources. We identify a subset of accounts from different Twitter communities discussing COVID-19, and follow these 'sentinel nodes' longitudinally from July 2020 to January 2021. We characterize sentinel nodes in terms of a linked domain preference score, and use a standardized similarity score to examine alignment of tweets within and between communities. We find that media preference is strongly correlated with the amount of misinformation propagated by sentinel nodes. Engagement with sensationalist misinformation topics is largely confined to a cluster of sentinel nodes that includes influential conspiracy theorist accounts. By contrast, misinformation relating to COVID-19 severity generated widespread engagement across multiple communities. Our findings indicate that misinformation downplaying COVID-19 severity is of particular concern for public health response. We conclude that the sentinel node approach can be an effective way to assess breadth and depth of online misinformation penetration.

Health Behaviors and Health Status among Middle-Aged and Older Adults with Chronic Diseases in Taiwan.

Changes in lifestyle behaviors may effectively maintain or improve the health status of individuals with chronic diseases. However, such health behaviors adopted by individuals are unlikely to demonstrate similar patterns. This study analyzed the relationship between the heterogeneous latent classes of health behavior and health statuses among middle-aged and older adults with hypertension, diabetes, or hyperlipidemia in Taiwan. After selecting 2103 individuals from the 2005 and 2009 Taiwan National Health Interview Survey (NHIS), we first identified heterogeneous groups of health behaviors through latent class analysis (LCA). We further explored the relationship between each latent class of health behavior and health status through ordered logit regression. We identified the following five distinct health behavior classes: the all-controlled, exercise and relaxation, healthy diet and reduced smoking or drinking, healthy diet, and least-controlled classes. Regression results indicated that individuals in classes other than the all-controlled class all reported poor health statuses. We also found great magnitude of the coefficient estimates for individuals who reported their health status to be poor or very poor for the least-controlled class. Therefore, health authorities and medical providers may develop targeted policies and interventions that address multiple modifiable health behaviors in each distinct latent class of health behavior.

Modeling the role of carrier and mobile herds on foot-and-mouth disease virus endemicity in the Far North Region of Cameroon.

Foot and mouth disease virus (FMDV) is an RNA virus that infects cloven-hoofed animals, often produces either epidemic or endemic conditions, and negatively affects agricultural economies worldwide. FMDV epidemic dynamics have been extensively studied, but understanding of drivers of disease persistence in areas in which FMDV is endemic, such as most of sub-Saharan Africa, is lacking. We present a spatial stochastic model of disease dynamics that incorporates a spatial transmission kernel in a modified Gillespie algorithm, and use it to evaluate two hypothesized drivers of endemicity: asymptomatic carriers and the movement of mobile herds. The model is parameterized using data from the pastoral systems in the Far North Region of Cameroon. Our computational study provides evidence in support of the hypothesis that asymptomatic carriers, but not mobile herds, are a driver of endemicity.

Patients with minor diseases who access high-tier medical care facilities: New evidence from classification and regression trees.

OBJECTIVES: Patients in Taiwan's National Health Insurance (NHI) program can choose a medical care facility of any tier for outpatient visits, without a referral. However, this system results in high medical expenditures and costs of outpatient visits. In this study, patients who had only minor diseases but who accessed high-tier medical care facilities were investigated using classification and regression trees. METHODS: For this study, data were obtained from the Taiwan NHI Research Database. First, 280 diseases, coded according to the Clinical Classification Software (CCS), were examined to determine whether patients chose the most appropriate facility when seeking medical care. After controlling for the CCS codes, an investigation into the types of patients who visit high-tier medical care facilities was conducted. RESULTS: Chronic disease status and CCS code were critical for constructing the classification trees. Male patients living in urban areas and earning a higher income were more likely to access high-tier medical care facilities. However, changes to the NHI copayment policies have significantly reduced the probability of utilizing high-tier medical care facilities. CONCLUSIONS: Factors relevant to patients' selection of high-tier medical care facilities were identified. Overall, increasing patients' out-of-pocket payments significantly reduced the probability of accessing high-tier medical facilities.

Modeling outbreak data: Analysis of a 2012 Ebola virus disease epidemic in DRC.

We describe two approaches to modeling data from a small to moderate-sized epidemic outbreak. The first approach is based on a branching process approximation and direct analysis of the transmission network, whereas the second one is based on a survival model derived from the classical SIR equations with no explicit transmission information. We compare these approaches using data from a 2012 outbreak of Ebola virus disease caused by Bundibugyo ebolavirus in city of Isiro, Democratic Republic of the Congo. The branching process model allows for a direct comparison of disease transmission across different environments, such as the general community or the Ebola treatment unit. However, the survival model appears to yield parameter estimates with more accuracy and better precision in some circumstances.

Complex contagion leads to complex dynamics in models coupling behaviour and disease.

Models coupling behaviour and disease as two unique but interacting contagions have existed since the mid 2000s. In these coupled contagion models, behaviour is typically treated as a 'simple contagion'. However, the means of behaviour spread may in fact be more complex. We develop a family of disease-behaviour coupled contagion compartmental models in order to examine the effect of behavioural contagion type on disease-behaviour dynamics. Coupled contagion models treating behaviour as a simple contagion and a complex contagion are investigated, showing that behavioural contagion type can have a significant impact on dynamics. We find that a simple contagion behaviour leads to simple dynamics, while a complex contagion behaviour supports complex dynamics with the possibility of bistability and periodic orbits.

The large graph limit of a stochastic epidemic model on a dynamic multilayer network.

We consider a Markovian SIR-type (Susceptible → Infected → Recovered) stochastic epidemic process with multiple modes of transmission on a contact network. The network is given by a random graph following a multilayer configuration model where edges in different layers correspond to potentially infectious contacts of different types. We assume that the graph structure evolves in response to the epidemic via activation or deactivation of edges of infectious nodes. We derive a large graph limit theorem that gives a system of ordinary differential equations (ODEs) describing the evolution of quantities of interest, such as the proportions of infected and susceptible vertices, as the number of nodes tends to infinity. Analysis of the limiting system elucidates how the coupling of edge activation and deactivation to infection status affects disease dynamics, as illustrated by a two-layer network example with edge types corresponding to community and healthcare contacts. Our theorem extends some earlier results describing the deterministic limit of stochastic SIR processes on static, single-layer configuration model graphs. We also describe precisely the conditions for equivalence between our limiting ODEs and the systems obtained via pair approximation, which are widely used in the epidemiological and ecological literature to approximate disease dynamics on networks. The flexible modeling framework and asymptotic results have potential application to many disease settings including Ebola dynamics in West Africa, which was the original motivation for this study.

Model distinguishability and inference robustness in mechanisms of cholera transmission and loss of immunity.

Mathematical models of cholera and waterborne disease vary widely in their structures, in terms of transmission pathways, loss of immunity, and a range of other features. These differences can affect model dynamics, with different models potentially yielding different predictions and parameter estimates from the same data. Given the increasing use of mathematical models to inform public health decision-making, it is important to assess model distinguishability (whether models can be distinguished based on fit to data) and inference robustness (whether inferences from the model are robust to realistic variations in model structure). In this paper, we examined the effects of uncertainty in model structure in the context of epidemic cholera, testing a range of models with differences in transmission and loss of immunity structure, based on known features of cholera epidemiology. We fit these models to simulated epidemic and long-term data, as well as data from the 2006 Angola epidemic. We evaluated model distinguishability based on fit to data, and whether the parameter values, model behavior, and forecasting ability can accurately be inferred from incidence data. In general, all models were able to successfully fit to all data sets, both real and simulated, regardless of whether the model generating the simulated data matched the fitted model. However, in the long-term data, the best model fits were achieved when the loss of immunity structures matched those of the model that simulated the data. Two parameters, one representing person-to-person transmission and the other representing the reporting rate, were accurately estimated across all models, while the remaining parameters showed broad variation across the different models and data sets. The basic reproduction number (R0) was often poorly estimated even using the correct model, due to practical unidentifiability issues in the waterborne transmission pathway which were consistent across all models. Forecasting efforts using noisy data were not successful early in the outbreaks, but once the epidemic peak had been achieved, most models were able to capture the downward incidence trajectory with similar accuracy. Forecasting from noise-free data was generally successful for all outbreak stages using any model. Our results suggest that we are unlikely to be able to infer mechanistic details from epidemic case data alone, underscoring the need for broader data collection, such as immunity/serology status, pathogen dose response curves, and environmental pathogen data. Nonetheless, with sufficient data, conclusions from forecasting and some parameter estimates were robust to variations in the model structure, and comparative modeling can help to determine how realistic variations in model structure may affect the conclusions drawn from models and data.

Network-based analysis of a small Ebola outbreak.

We present a method for estimating epidemic parameters in network-based stochastic epidemic models when the total number of infections is assumed to be small. We illustrate the method by reanalyzing the data from the 2014 Democratic Republic of the Congo (DRC) Ebola outbreak described in Maganga et al. (2014).

Modeling the trade-off between transmissibility and contact in infectious disease dynamics.

Symptom severity affects disease transmission both by impacting contact rates, as well as by influencing the probability of transmission given contact. This involves a trade-off between these two factors, as increased symptom severity will tend to decrease contact rates, but increase the probability of transmission given contact (as pathogen shedding rates increase with symptom severity). This paper explores this trade-off between contact and transmission given contact, using a simple compartmental susceptible-infected-recovered type model. Under mild assumptions on how contact and transmission probability vary with symptom severity, we give sufficient, biologically intuitive criteria for when the basic reproduction number varies non-monotonically with symptom severity. Multiple critical points are possible. We give a complete characterization of the region in parameter space where multiple critical points are located in the special case where contact rate decreases exponentially with symptom severity. We consider a multi-strain version of the model with complete cross-immunity and no super-infection. In this model, we prove that the strain with highest basic reproduction number drives the other strains to extinction. This has both evolutionary and epidemiological implications, including the possibility of an intervention paradoxically resulting in increased infection prevalence.

The impact of spatial arrangements on epidemic disease dynamics and intervention strategies.

The role of spatial arrangements on the spread and management strategies of a cholera epidemic is investigated. We consider the effect of human and pathogen movement on optimal vaccination strategies. A metapopulation model is used, incorporating a susceptible-infected-recovered system of differential equations coupled with an equation modelling the concentration of Vibrio cholerae in an aquatic reservoir. The model compared spatial arrangements and varying scenarios to draw conclusions on how to effectively manage outbreaks. The work is motivated by the 2010 cholera outbreak in Haiti. Results give guidance for vaccination strategies in response to an outbreak.

Measuring inequality in physician distributions using spatially adjusted Gini coefficients.

OBJECTIVE: To measure inequality in physician distributions using Gini coefficient and spatially adjusted Gini coefficients. DESIGN: Measurements were based on the distribution of physician data from the Taiwan National Health Insurance Research Database (NHIRD) and population data from the Ministry of the Interior in Taiwan. SETTINGS: The distribution of population and physicians in Taiwan from 2001 to 2010. PARTICIPANTS: This study considered 35 000 physicians who are registered in Taiwan. MAIN OUTCOME MEASURES: To calculate the Gini coefficient and spatially adjusted Gini coefficients in Taiwan from 2001 to 2010. RESULTS: The Gini coefficient for each year, from 2001 to 2010, ranged from 0.5128 to 0.4692, while the spatially adjusted Gini coefficients based on travel time and travel distance ranged, respectively, from 0.4324 to 0.4066 and from 0.4408 to 0.4178. We found that, in each year, irrespective of the type of spatial adjustment, the spatially adjusted Gini coefficient was smaller than the Gini coefficient itself. Our empirical findings support that the Gini coefficient may overestimate the maldistribution of physicians. CONCLUSIONS: Our simulations demonstrate that increasing the number of physicians in medium-sized cities (such as capitals of counties or provinces), and/or improving the transportation time between medium-sized cities and rural areas, could be feasible solutions to mitigate the problem of geographical maldistribution of physicians.