scholarly journals Learning differential equation models from stochastic agent-based model simulations

2021 ◽  
Vol 18 (176) ◽  
Author(s):  
John T. Nardini ◽  
Ruth E. Baker ◽  
Matthew J. Simpson ◽  
Kevin B. Flores

Agent-based models provide a flexible framework that is frequently used for modelling many biological systems, including cell migration, molecular dynamics, ecology and epidemiology. Analysis of the model dynamics can be challenging due to their inherent stochasticity and heavy computational requirements. Common approaches to the analysis of agent-based models include extensive Monte Carlo simulation of the model or the derivation of coarse-grained differential equation models to predict the expected or averaged output from the agent-based model. Both of these approaches have limitations, however, as extensive computation of complex agent-based models may be infeasible, and coarse-grained differential equation models can fail to accurately describe model dynamics in certain parameter regimes. We propose that methods from the equation learning field provide a promising, novel and unifying approach for agent-based model analysis. Equation learning is a recent field of research from data science that aims to infer differential equation models directly from data. We use this tutorial to review how methods from equation learning can be used to learn differential equation models from agent-based model simulations. We demonstrate that this framework is easy to use, requires few model simulations, and accurately predicts model dynamics in parameter regions where coarse-grained differential equation models fail to do so. We highlight these advantages through several case studies involving two agent-based models that are broadly applicable to biological phenomena: a birth–death–migration model commonly used to explore cell biology experiments and a susceptible–infected–recovered model of infectious disease spread.

2019 ◽  
Vol 16 (159) ◽  
pp. 20190421 ◽  
Author(s):  
Nabil T. Fadai ◽  
Ruth E. Baker ◽  
Matthew J. Simpson

Understanding how cells proliferate, migrate and die in various environments is essential in determining how organisms develop and repair themselves. Continuum mathematical models, such as the logistic equation and the Fisher–Kolmogorov equation, can describe the global characteristics observed in commonly used cell biology assays, such as proliferation and scratch assays. However, these continuum models do not account for single-cell-level mechanics observed in high-throughput experiments. Mathematical modelling frameworks that represent individual cells, often called agent-based models, can successfully describe key single-cell-level features of these assays but are computationally infeasible when dealing with large populations. In this work, we propose an agent-based model with crowding effects that is computationally efficient and matches the logistic and Fisher–Kolmogorov equations in parameter regimes relevant to proliferation and scratch assays, respectively. This stochastic agent-based model allows multiple agents to be contained within compartments on an underlying lattice, thereby reducing the computational storage compared to existing agent-based models that allow one agent per site only. We propose a systematic method to determine a suitable compartment size. Implementing this compartment-based model with this compartment size provides a balance between computational storage, local resolution of agent behaviour and agreement with classical continuum descriptions.


2019 ◽  
Author(s):  
Nabil T. Fadai ◽  
Ruth E. Baker ◽  
Matthew J. Simpson

AbstractUnderstanding how cells proliferate, migrate, and die in various environments is essential in determining how organisms develop and repair themselves. Continuum mathematical models, such as the logistic equation and the Fisher-Kolmogorov equation, can describe the global characteristics observed in commonly-used cell biology assays, such as proliferation and scratch assays. However, these continuum models do not account for single-cell-level mechanics observed in high-throughput experiments. Mathematical modelling frameworks that represent individual cells, often called agent-based models, can successfully describe key single-cell-level features of these assays, but are computationally infeasible when dealing with large populations. In this work, we propose an agent-based model with crowding effects that is computationally efficient and matches the logistic and Fisher-Kolmogorov equations in parameter regimes relevant to proliferation and scratch assays, respectively. This stochastic agent-based model allows multiple agents to be contained within compartments on an underlying lattice, thereby reducing the computational storage compared to existing agent-based models that allow one agent per site only. We propose a systematic method to determine a suitable compartment size. Implementing this compartment-based model with this compartment size provides a balance between computational storage, local resolution of agent behaviour, and agreement with classical continuum descriptions.


2015 ◽  
Author(s):  
Heiko Enderling

For quantitative cancer models to be meaningful and interpretable the number of unknown parameters must be kept minimal. Experimental data can be utilized to calibrate model dynamics rates or rate constants. Proper integration of experimental data, however, depends on the chosen theoretical framework. Using live imaging of cell proliferation as an example, we show how to derive cell cycle distributions in agent-based models and averaged proliferation rates in differential equation models. We focus on a tumor hierarchy of cancer stem and progenitor non-stem cancer cells.


2018 ◽  
Vol 10 (1) ◽  
Author(s):  
Sultanah Alshammari ◽  
Armin Mikler

ObjectiveTo develop a computational model to assess the risk of epidemics in global mass gatherings and evaluate the impact of various measures of prevention and control of infectious diseases.IntroductionGlobal Mass gatherings (MGs) such as Olympic Games, FIFA World Cup, and Hajj (Muslim pilgrimage to Makkah), attract millions of people from different countries. The gathering of a large population in a proximity facilitates transmission of infectious diseases [1]. Attendees arrive from different geographical areas with diverse disease history and immune responses. The associated travel patterns with global events can contribute to a further disease spread affecting a large number of people within a short period and lead to a potential pandemic. Global MGs pose serious health threats and challenges to the hosting countries and home countries of the participants [2]. Advanced planning and disease surveillance systems are required to control health risks in these events. The success of computational models in different areas of public health and epidemiology motivates using these models in MGs to study transmission of infectious diseases and assess the risk of epidemics. Computational models enable simulation and analysis of different disease transmission scenarios in global MGs. Epidemic models can be used to evaluate the impact of various measures of prevention and control of infectious diseases.MethodsThe annual event of the Hajj is selected to illustrate the main aspects of the proposed model and to address the associated challenges. Every year, more than two million pilgrims from over 186 countries arrive in Makkah to perform Hajj with the majority arriving by air. Foreign pilgrims can stay at one of the holy cities of Makkah and Madinah up to 30-35 days prior the starting date of the Hajj. The long duration of the arrival phase of the Hajj allows a potential epidemic to proceed in the population of international pilgrims. Stochastic SEIR (Susceptible−Exposed−Infected−Recovered) agent-based model is developed to simulate the disease transmission among pilgrims. The agent-based model is used to simulate pilgrims and their interactions during the various phases of the Hajj. Each agent represents a pilgrim and maintains a record of demographic data (gender, country of origin, age), health data (infectivity, susceptibility, number of days being exposed or infected), event related data (location, arrival date and time), and precautionary or health-related behaviors.Each pilgrim can be either healthy but susceptible to a disease, exposed who are infected but cannot transmit the infection, or infectious (asymptomatic or symptomatic) who are infected and can transmit the disease to other susceptibles. Exposed individuals transfer to the infectious compartment after 1/α days, and infectious individuals will recover and gain immunity to that disease after 1/γ days. Where α is the latent period and γ is the infectious period. Moving susceptible individuals to exposed compartment depends on a successful disease transmission given a contact with an infectious individual. The disease transmission rate is determined by the contact rate and thetransmission probability per contact. Contact rate and mixing patterns are defined by probabilistic weights based on the features of infectious pilgrims and the duration and setting of the stage where contacts are taking place. The initial infections are seeded in the population using two scenarios (Figure 1) to measure the effects of changing, the timing for introducing a disease into the population and the likelihood that a particular flight will arrive with one or more infected individuals.ResultsThe results showed that the number of initial infections is influenced by increasing the value of λ and selecting starting date within peak arrival days. When starting from the first day, the average size of the initial infectious ranges from 0.05% to 1% of the total arriving pilgrims. Using the SEIR agent-based model, a simulation of the H1N1 Influenza epidemic was completed for the 35-days arrival stage of the Hajj. The epidemic is initiated with one infectious pilgrim per flight resulting in infected 0.5% of the total arriving pilgrims. As pilgrims spend few hours at the airport, the results obtained from running the epidemic model showed only new cases of susceptible individuals entering the exposed state in a range of 0.20% to 0.35% of total susceptibles. The number of new cases is reduced by almost the same rate of the number of infectious individuals following precautionary behaviors.ConclusionsA data-driven stochastic SEIR agent-based model is developed to simulate disease spread at global mass gatherings. The proposed model can provide initial indicators of infectious disease epidemic at these events and evaluate the possible effects of intervention measures and health-related behaviors. The proposed model can be generalized to model the spread of various diseases in different mass gatherings, as it allows different factors to vary and entered as parameters.References1. Memish ZA, Stephens GM, Steffen R, Ahmed QA. Emergence of medicine for mass gatherings: lessons from the Hajj. The Lancet infectious diseases. 2012 Jan 31;12(1):56-65.2. Chowell G, Nishiura H, Viboud C. Modeling rapidly disseminating infectious disease during mass gatherings. BMC medicine. 2012 Dec 7;10(1):159.


Author(s):  
Herbert Dawid ◽  
Simon Gemkow ◽  
Philipp Harting ◽  
Sander van der Hoog ◽  
Michael Neugart

This chapter introduces the Eurace@Unibi model, one of the agent-based simulation models that are relatively new additions to the toolbox of macroeconomists, and the research that has been done within this framework. It shows how an agent-based model can be used to identify economic mechanisms and how it can be applied to spatial policy analysis. The assessment is that agent-based models in economics have passed the proof-of-concept phase and it is now time to move beyond that stage. It has been shown that new kinds of insights can be obtained that complement established modeling approaches. The chapter concludes by pointing toward some potentially fruitful areas of agent-based macroeconomic research.


2018 ◽  
Vol 151 ◽  
pp. 21-28 ◽  
Author(s):  
Kelsey L. Spence ◽  
Terri L. O’Sullivan ◽  
Zvonimir Poljak ◽  
Amy L. Greer

Sign in / Sign up

Export Citation Format

Share Document