BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Christian A. Yates

Christian A. Yates appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
Quantitative MethodsPopulations and Evolutioncond-mat.stat-mechphysics.soc-phRoboticsSystems and Control

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2026arXiv

An agent-based modelling approach to investigate the impact of gender on tuberculosis transmission in Uganda

Tuberculosis (TB) is an airborne disease caused by the pathogen Mycobacterium tuberculosis. In 2023, it returned to being the leading cause of death from an infectious agent globally, replacing COVID-19; in the nineteenth century, one in seven of all humans died of tuberculosis. More than 10 million people are diagnosed with TB every year. The majority of cases in adults occur in males (62.5% of all global adult cases in 2023, compared to 37.5% in females). The main reasons for males suffering from a higher burden of global TB cases, compared to females, may be in large part due to population-scale factors, such as employment type, the quantity and type of social contacts they make, and their health-seeking behaviours (e.g. differences in diagnostic and treatment delays between genders). To investigate which population-scale factors are most important in determining this higher TB burden in males, we have developed an age- and gender-stratified, spatially heterogeneous epidemiological agent-based model. We have focused specifically on Kampala, the capital of Uganda, which is a high-burden TB country. We considered counterfactual scenarios to elucidate the impact of gender on the epidemiology of TB. Setting disease progression parameters equal between the genders leads to a reduction in both male-to-female case ratio and total case numbers.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Stochastic drift in discrete waves of non-locally interacting-particles

In this paper, we investigate a generalised model of $N$ particles undergoing second-order non-local interactions on a lattice. Our results have applications across many research areas, including the modelling of migration, information dynamics and Muller's ratchet -- the irreversible accumulation of deleterious mutations in an evolving population. Strikingly, numerical simulations of the model are observed to deviate significantly from its mean-field approximation even for large population sizes. We show that the disagreement between deterministic and stochastic solutions stems from finite-size effects that change the propagation speed and cause the position of the wave to fluctuate. These effects are shown to decay anomalously as $(\log N)^{-2}$ and $(\log N)^{-3}$, respectively -- much slower than the usual $N^{-1/2}$ factor. Our results suggest that the accumulation of deleterious mutations in a Muller's ratchet and the loss of awareness in a population may occur much faster than predicted by the corresponding deterministic models. The general applicability of our model suggests that this unexpected scaling could be important in a wide range of real-world applications.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An equivalence framework for an age-structured multi-stage representation of the cell cycle

We develop theoretical equivalences between stochastic and deterministic models for populations of individual cells stratified by age. Specifically, we develop a hierarchical system of equations describing the full dynamics of an age-structured multi-stage Markov process for approximating cell cycle time distributions. We further demonstrate that the resulting mean behaviour is equivalent, over large timescales, to the classical McKendrick-von Foerster integro-partial differential equation. We conclude by extending this framework to a spatial context, facilitating the modelling of travelling wave phenomena and cell-mediated pattern formation. More generally, this methodology may be extended to myriad reaction-diffusion processes for which the age of individuals is relevant to the dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Critical weaknesses in shielding strategies for COVID-19

The COVID-19 pandemic, caused by the coronavirus SARS-CoV-2, has led to a wide range of non-pharmaceutical interventions being implemented around the world to curb transmission. However, the economic and social costs of some of these measures, especially lockdowns, has been high. An alternative and widely discussed public health strategy for the COVID-19 pandemic would have been to 'shield' those most vulnerable to COVID-19 (minimising their contacts with others), while allowing infection to spread among lower risk individuals with the aim of reaching herd immunity. Here we retrospectively explore the effectiveness of this strategy using a stochastic SEIR framework, showing that even under the unrealistic assumption of perfect shielding, hospitals would have been rapidly overwhelmed with many avoidable deaths among lower risk individuals. Crucially, even a small (20%) reduction in the effectiveness of shielding would have likely led to a large increase (>150%) in the number of deaths compared to perfect shielding. Our findings demonstrate that shielding the vulnerable while allowing infections to spread among the wider population would not have been a viable public health strategy for COVID-19 and is unlikely to be effective for future pandemics.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Efficiently simulating discrete-state models with binary decision trees

Stochastic simulation algorithms (SSAs) are widely used to numerically investigate the properties of stochastic, discrete-state models. The Gillespie Direct Method is the pre-eminent SSA, and is widely used to generate sample paths of so-called agent-based or individual-based models. However, the simplicity of the Gillespie Direct Method often renders it impractical where large-scale models are to be analysed in detail. In this work, we carefully modify the Gillespie Direct Method so that it uses a customised binary decision tree to trace out sample paths of the model of interest. We show that a decision tree can be constructed to exploit the specific features of the chosen model. Specifically, the events that underpin the model are placed in carefully-chosen leaves of the decision tree in order to minimise the work required to keep the tree up-to-date. The computational efficencies that we realise can provide the apparatus necessary for the investigation of large-scale, discrete-state models that would otherwise be intractable. Two case studies are presented to demonstrate the efficiency of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Incorporating domain growth into hybrid methods for reaction-diffusion systems

Reaction--diffusion mechanism are a robust paradigm that can be used to represent many biological and physical phenomena over multiple spatial scales. Applications include intracellular dynamics, the migration of cells and the patterns formed by vegetation in semi-arid landscapes. Moreover, domain growth is an important process for embryonic growth and wound healing. There are many numerical modelling frameworks capable of simulating such systems on growing domains, however each of these may be well suited to different spatial scales and particle numbers. Recently, spatially extended hybrid methods on static domains have been produced in order to bridge the gap between these different modelling paradigms in order to represent multiscale phenomena. However, such methods have not been developed with domain growth in mind. In this paper, we develop three hybrid methods on growing domains, extending three of the prominent static domain hybrid methods. We also provide detailed algorithms to allow others to employ them. We demonstrate that the methods are able to accurately model three representative reaction-diffusion systems accurately and without bias.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A hybrid algorithm for coupling PDE and compartment-based dynamics

Stochastic simulation methods can be applied successfully to model exact spatio-temporally resolved reaction-diffusion systems. However, in many cases, these methods can quickly become extremely computationally intensive with increasing particle numbers. An alternative description of many of these systems can be derived in the diffusive limit as a deterministic, continuum system of partial differential equations. Although the numerical solution of such partial differential equations is, in general, much more efficient than the full stochastic simulation, the deterministic continuum description is generally not valid when copy numbers are low and stochastic effects dominate. Therefore, to take advantage of the benefits of both of these types of models, each of which may be appropriate in different parts of a spatial domain, we have developed an algorithm that can be used to couple these two types of model together. This hybrid coupling algorithm uses an overlap region between the two modelling regimes. By coupling fluxes at one end of the interface and using a concentration-matching condition at the other end, we ensure that mass is appropriately transferred between PDE- and compartment-based regimes. Our methodology gives notable reductions in simulation time in comparison with using a fully stochastic model, whilst maintaining the important stochastic features of the system and providing detail in appropriate areas of the domain. We test our hybrid methodology robustly by applying it to several biologically motivated problems including diffusion and morphogen gradient formation. Our analysis shows that the resulting error is small, unbiased and does not grow over time.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Coupling volume-excluding compartment-based models of diffusion at different scales: Voronoi and pseudo-compartment approaches

Numerous processes across both the physical and biological sciences are driven by diffusion. Partial differential equations (PDEs) are a popular tool for modelling such phenomena deterministically, but it is often necessary to use stochastic models to accurately capture the behaviour of a system, especially when the number of diffusing particles is low. The stochastic models we consider in this paper are `compartment-based': the domain is discretized into compartments, and particles can jump between these compartments. Volume-excluding effects (crowding) can be incorporated by blocking movement with some probability. Recent work has established the connection between fine-grained models and coarse-grained models incorporating volume exclusion, but only for uniform lattices. In this paper we consider non-uniform, hybrid lattices that incorporate both fine- and coarse-grained regions, and present two different approaches to describing the interface of the regions. We test both techniques in a range of scenarios to establish their accuracy, benchmarking against fine-grained models, and show that the hybrid models developed in this paper can be significantly faster to simulate than the fine-grained models in certain situations, and are at least as fast otherwise.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Extending the multi-level method for the simulation of stochastic biological systems

The multi-level method for discrete state systems, first introduced by Anderson and Higham [Multiscale Model. Simul. 10:146--179, 2012], is a highly efficient simulation technique that can be used to elucidate statistical characteristics of biochemical reaction networks. A single point estimator is produced in a cost-effective manner by combining a number of estimators of differing accuracy in a telescoping sum, and, as such, the method has the potential to revolutionise the field of stochastic simulation. The first term in the sum is calculated using an approximate simulation algorithm, and can be calculated quickly but is of significant bias. Subsequent terms successively correct this bias by combining estimators from approximate stochastic simulations algorithms of increasing accuracy, until a desired level of accuracy is reached. In this paper we present several refinements of the multi-level method which render it easier to understand and implement, and also more efficient. Given the substantial and complex nature of the multi-level method, the first part of this work (Sections 2 - 5) is written as a tutorial, with the aim of providing a practical guide to its use. The second part (Sections 6 - 8) takes on a form akin to a research article, thereby providing the means for a deft implementation of the technique, and concludes with a discussion of a number of open problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Reconciling transport models across scales: the role of volume exclusion

Diffusive transport is a universal phenomenon, throughout both biological and physical sciences, and models of diffusion are routinely used to interrogate diffusion-driven processes. However, most models neglect to take into account the role of volume exclusion, which can significantly alter diffusive transport, particularly within biological systems where the diffusing particles might occupy a significant fraction of the available space. In this work we use a random walk approach to provide a means to reconcile models that incorporate crowding effects on different spatial scales. Our work demonstrates that coarse-grained models incorporating simplified descriptions of excluded volume can be used in many circumstances, but that care must be taken in pushing the coarse-graining process too far.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The pseudo-compartment method for coupling PDE and compartment-based models of diffusion

Spatial reaction-diffusion models have been employed to describe many emergent phenomena in biological systems. The modelling technique most commonly adopted in the literature implements systems of partial differential equations (PDEs), which assumes there are sufficient densities of particles that a continuum approximation is valid. However, due to recent advances in computational power, the simulation, and therefore postulation, of computationally intensive individual-based models has become a popular way to investigate the effects of noise in reaction-diffusion systems in which regions of low copy numbers exist. The stochastic models with which we shall be concerned in this manuscript are referred to as `compartment-based'. These models are characterised by a discretisation of the computational domain into a grid/lattice of `compartments'. Within each compartment particles are assumed to be well-mixed and are permitted to react with other particles within their compartment or to transfer between neighbouring compartments. We develop two hybrid algorithms in which a PDE is coupled to a compartment-based model. Rather than attempting to balance average fluxes, our algorithms answer a more fundamental question: `how are individual particles transported between the vastly different model descriptions?' First, we present an algorithm derived by carefully re-defining the continuous PDE concentration as a probability distribution. Whilst this first algorithm shows strong convergence to analytic solutions of test problems, it can be cumbersome to simulate. Our second algorithm is a simplified and more efficient implementation of the first, it is derived in the continuum limit over the PDE region alone. We test our hybrid methods for functionality and accuracy in a variety of different scenarios by comparing the averaged simulations to analytic solutions of PDEs for mean concentrations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An adaptive multi-level simulation algorithm for stochastic biological systems

Discrete-state, continuous-time Markov models are widely used in the modeling of biochemical reaction networks. Their complexity often precludes analytic solution, and we rely on stochastic simulation algorithms to estimate system statistics. The Gillespie algorithm is exact, but computationally costly as it simulates every single reaction. As such, approximate stochastic simulation algorithms such as the tau-leap algorithm are often used. Potentially computationally more efficient, the system statistics generated suffer from significant bias unless tau is relatively small, in which case the computational time can be comparable to that of the Gillespie algorithm. The multi-level method (Anderson and Higham, Multiscale Model. Simul. 2012) tackles this problem. A base estimator is computed using many (cheap) sample paths at low accuracy. The bias inherent in this estimator is then reduced using a number of corrections. Each correction term is estimated using a collection of paired sample paths where one path of each pair is generated at a higher accuracy compared to the other (and so more expensive). By sharing random variables between these paired paths the variance of each correction estimator can be reduced. This renders the multi-level method very efficient as only a relatively small number of paired paths are required to calculate each correction term. In the original multi-level method, each sample path is simulated using the tau-leap algorithm with a fixed value of $τ$. This approach can result in poor performance when the reaction activity of a system changes substantially over the timescale of interest. By introducing a novel, adaptive time-stepping approach where $τ$ is chosen according to the stochastic behaviour of each sample path we extend the applicability of the multi-level method to such cases. We demonstrate the efficiency of our method using a number of examples.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mathematical Modelling of Turning Delays in Swarm Robotics

We investigate the effect of turning delays on the behaviour of groups of differential wheeled robots and show that the group-level behaviour can be described by a transport equation with a suitably incorporated delay. The results of our mathematical analysis are supported by numerical simulations and experiments with e-puck robots. The experimental quantity we compare to our revised model is the mean time for robots to find the target area in an unknown environment. The transport equation with delay better predicts the mean time to find the target than the standard transport equation without delay.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Ergodic directional switching in mobile insect groups

We obtain a Fokker-Planck equation describing experimental data on the collective motion of locusts. The noise is of internal origin and due to the discrete character and finite number of constituents of the swarm. The stationary probability distribution shows a rich phenomenology including non-monotonic behavior of several order/disorder transition indicators in noise intensity. This complex behavior arises naturally as a result of the randomness in the system. Its counterintuitive character challenges standard interpretations of noise induced transitions and calls for an extension of this theory in order to capture the behavior of certain classes of biologically motivated models. Our results suggest that the collective switches of the group's direction of motion might be due to a random ergodic effect and, as such, they are inherent to group formation.

0 citations0 reviewsNo code linkedOpen access