We develop a stochastic, agent-based model to study how genetic traits and experiential changes in the state of agents and available resources influence individuals’ foraging and movement behaviors. These behaviors are manifest as decisions on when to stay and exploit a current resource patch or move to a particular neighboring patch, based on information of the resource qualities of the patches and the anticipated level of intraspecific competition within patches. We use a genetic algorithm approach and an individual’s biomass as a fitness surrogate to explore the foraging strategy diversity of evolving guilds under clonal versus hermaphroditic sexual reproduction. We first present the resource exploitation processes, movement on cellular arrays, and genetic algorithm components of the model. We then discuss their implementation on the Nova software platform. This platform seamlessly combines the dynamical systems modeling of consumer-resource interactions with agent-based modeling of individuals moving over a landscapes, using an architecture that lays transparent the following four hierarchical simulation levels: 1.) within-patch consumer-resource dynamics, 2.) within-generation movement and competition mitigation processes, 3.) across-generation evolutionary processes, and 4.) multiple runs to generate the statistics needed for comparative analyses. The focus of our analysis is on the question of how the biomass production efficiency and the diversity of guilds of foraging strategy types, exploiting resources over a patchy landscape, evolve under clonal versus random hermaphroditic sexual reproduction. Our results indicate greater biomass production efficiency under clonal reproduction only at higher population densities, and demonstrate that polymorphisms evolve and are maintained under random mating systems. The latter result questions the notion that some type of associative mating structure is needed to maintain genetic polymorphisms among individuals exploiting a common patchy resource on an otherwise spatially homogeneous landscape.

Epidemiological models are dominated by compartmental models, of which SIR formulations are the most commonly used. These formulations can be continuous or discrete (in either the state-variable values or time), deterministic or stochastic, or spatially homogeneous or heterogeneous, the latter often embracing a network formulation. Here we review the continuous and discrete deterministic and discrete stochastic formulations of the SIR dynamical systems models, and we outline how they can be easily and rapidly constructed using Numerus Model Builder, a graphically-driven coding platform. We also demonstrate how to extend these models to a metapopulation setting using NMB network and mapping tools.

One contribution of 15 to a theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes'.Subject Areas: computational biology, health and disease and epidemiology Dynamic SEIR (Susceptible, Exposed, Infectious, Removed) compartmental models provide a tool for predicting the size and duration of both unfettered and managed outbreaks-the latter in the context of interventions such as case detection, patient isolation, vaccination and treatment. The reliability of this tool depends on the validity of key assumptions that include homogeneity of individuals and spatio-temporal homogeneity. Although the SEIR compartmental framework can easily be extended to include demographic (e.g. age) and additional disease (e.g. healthcare workers) classes, dependence of transmission rates on time, and metapopulation structure, fitting such extended models is hampered by both a proliferation of free parameters and insufficient or inappropriate data. This raises the question of how effective a tool the basic SEIR framework may actually be. We go some way here to answering this question in the context of the 2014-2015 outbreak of Ebola in West Africa by comparing fits of an SEIR time-dependent transmission model to both country-and district-level weekly incidence data. Our novel approach in estimating the effective-size-of-the-populations-at-risk (N eff ) and initial number of exposed individuals (E 0 ) at both district and country levels, as well as the transmission function parameters, including a time-to-halving-the-force-ofinfection (t f/2 ) parameter, provides new insights into this Ebola outbreak. It reveals that the estimate R 0 1.7 from country-level data appears to seriously underestimate R 0 3.3 2 4.3 obtained from more spatially homogeneous district-level data. Country-level data also overestimate t f/2 22 weeks, compared with 8-10 weeks from district-level data. Additionally, estimates for the duration of individual infectiousness is around two weeks from spatially inhomogeneous country-level data compared with 2.4-4.5 weeks from spatially more homogeneous district-level data, which estimates are rather high compared with most values reported in the literature.This article is part of the theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: approaches and important themes'. This issue is linked with the subsequent theme issue 'Modelling infectious disease outbreaks in humans, animals and plants: epidemic forecasting and control'.

We present a stochastic transmission chain simulation model for Ebola viral disease (EVD) in West Africa, with the salutary result that the virus may be more controllable than previously suspected. The ongoing tactics to detect cases as rapidly as possible and isolate individuals as safely as practicable is essential to saving lives in the current outbreaks in Guinea, Liberia, and Sierra Leone. Equally important are educational campaigns that reduce contact rates between susceptible and infectious individuals in the community once an outbreak occurs. However, due to the relatively low R 0 of Ebola (around 1.5 to 2.5 next generation cases are produced per current generation case in naïve populations), rapid isolation of infectious individuals proves to be highly efficacious in containing outbreaks in new areas, while vaccination programs, even with low efficacy vaccines, can be decisive in curbing future outbreaks in areas where the Ebola virus is maintained in reservoir populations.

Background No versatile web app exists that allows epidemiologists and managers around the world to comprehensively analyze the impacts of COVID-19 mitigation. The http://covid-webapp.numerusinc.com/ web app presented here fills this gap. Methods Our web app uses a model that explicitly identifies susceptible, contact, latent, asymptomatic, symptomatic and recovered classes of individuals, and a parallel set of response classes, subject to lower pathogen-contact rates. The user inputs a CSV file of incidence and, if of interest, mortality rate data. A default set of parameters is available that can be overwritten through input or online entry, and a user-selected subset of these can be fitted to the model using maximum-likelihood estimation (MLE). Model fitting and forecasting intervals are specifiable and changes to parameters allow counterfactual and forecasting scenarios. Confidence or credible intervals can be generated using stochastic simulations, based on MLE values, or on an inputted CSV file containing Markov chain Monte Carlo (MCMC) estimates of one or more parameters. Results We illustrate the use of our web app in extracting social distancing, social relaxation, surveillance or virulence switching functions (i.e., time varying drivers) from the incidence and mortality rates of COVID-19 epidemics in Israel, South Africa, and England. The Israeli outbreak exhibits four distinct phases: initial outbreak, social distancing, social relaxation, and a second wave mitigation phase. An MCMC projection of this latter phase suggests the Israeli epidemic will continue to produce into late November an average of around 1500 new case per day, unless the population practices social-relaxation measures at least 5-fold below the level in August, which itself is 4-fold below the level at the start of July. Our analysis of the relatively late South African outbreak that became the world’s fifth largest COVID-19 epidemic in July revealed that the decline through late July and early August was characterised by a social distancing driver operating at more than twice the per-capita applicable-disease-class (pc-adc) rate of the social relaxation driver. Our analysis of the relatively early English outbreak, identified a more than 2-fold improvement in surveillance over the course of the epidemic. It also identified a pc-adc social distancing rate in early August that, though nearly four times the pc-adc social relaxation rate, appeared to barely contain a second wave that would break out if social distancing was further relaxed. Conclusion Our web app provides policy makers and health officers who have no epidemiological modelling or computer coding expertise with an invaluable tool for assessing the impacts of different outbreak mitigation policies and measures. This includes an ability to generate an epidemic-suppression or curve-flattening index that measures the intensity with which behavioural responses suppress or flatten the epidemic curve in the region under consideration.