Abstract
In the recent COVID-19 pandemic, mathematical modeling constitutes an important tool to evaluate the prospective effectiveness of non-pharmaceutical interventions (NPIs) and to guide policy-making. Most research is, however, centered around characterizing the epidemic based on point estimates like the average infectiousness or the average number of contacts.
In this work, we use stochastic simulations to investigate the consequences of a population’s heterogeneity regarding connectivity and individual viral load levels.
Therefore, we translate a COVID-19 ODE model to a stochastic multi-agent system. We use contact networks to model complex interaction structures and a probabilistic infection rate to model individual viral load variation.
We observe a large dependency of the dispersion and dynamical evolution on the population’s heterogeneity that is not adequately captured by point estimates, for instance, used in ODE models. In particular, models that assume the same clinical and transmission parameters may lead to different conclusions, depending on different types of heterogeneity in the population. For instance, the existence of hubs in the contact network leads to an initial increase of dispersion and the effective reproduction number, but to a lower herd immunity threshold (HIT) compared to homogeneous populations or a population where the heterogeneity stems solely from individual infectivity variations.
Author summary Computational modeling can support decision-making in the face of pandemics like COVID-19. Models help to understand transmission data and predict important epidemiological properties (e.g., When will herd immunity be reached?). They can also examine the effectiveness of certain measures, and—to a limited extent—extrapolate the dynamics under specific assumptions. In all these cases, the heterogeneity of the population plays an important role. For instance, it is known that connectivity differences in (and among) age groups influence the dynamics of epidemic propagation. Here we focus on two types of differences among individuals: their social interactions and on how infectious they are. We show that only considering population averages (e.g., What is the average number of contacts of an individual?) may lead to misleading conclusions, because the individual differences (such as those related to the epidemic (over-)dispersion) play an important role in shaping the epidemic dynamics. Many commonly used model classes, such as SEIR-type ODE compartmental models, ignore differences within a population to a large extent. This omission bears the potential of misleading conclusions.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This work was partially supported by the DFG project MULTIMODE.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
There were no ethics committee approvals needed.
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I 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).
Yes
I 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.
Yes
Data Availability
This is a simulation study. Code is available online.