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

Tobias Galla

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

ResearcherUnclaimed source record
cond-mat.stat-mechPopulations and Evolutionphysics.soc-phnlin.AOcond-mat.dis-nnQuantitative MethodsBiological PhysicsSocial and Information NetworksComputer Science and Game Theoryhep-latMolecular Networksnlin.PSquant-phTissues and OrgansApplicationscond-mat.softMethodologynlin.CDnlin.CGphysics.med-phq-fin.TR

Catalog footprint

What is connected

51works
21topics
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

51 published item(s)

preprint2021arXiv

Network meta-analysis and random walks

Network meta-analysis (NMA) is a central tool for evidence synthesis in clinical research. The results of an NMA depend critically on the quality of evidence being pooled. In assessing the validity of an NMA, it is therefore important to know the proportion contributions of each direct treatment comparison to each network treatment effect. The construction of proportion contributions is based on the observation that each row of the hat matrix represents a so-called 'evidence flow network' for each treatment comparison. However, the existing algorithm used to calculate these values is associated with ambiguity according to the selection of paths. In this work we present a novel analogy between NMA and random walks. We use this analogy to derive closed-form expressions for the proportion contributions. A random walk on a graph is a stochastic process that describes a succession of random 'hops' between vertices which are connected by an edge. The weight of an edge relates to the probability that the walker moves along that edge. We use the graph representation of NMA to construct the transition matrix for a random walk on the network of evidence. We show that the net number of times a walker crosses each edge of the network is related to the evidence flow network. By then defining a random walk on the directed evidence flow network, we derive analytically the matrix of proportion contributions. The random-walk approach, in addition to being computationally more efficient, has none of the associated ambiguity of the existing algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Sampling rare trajectories using stochastic bridges

The numerical quantification of the statistics of rare events in stochastic processes is a challenging computational problem. We present a sampling method that constructs an ensemble of stochastic trajectories that are constrained to have fixed start and end points (so-called stochastic bridges). We then show that by carefully choosing a set of such bridges and assigning an appropriate statistical weight to each bridge, one can focus more processing power on the rare events of a target stochastic process while faithfully preserving the statistics of these rate trajectories. Further, we also compare the stochastic bridges produced using our method to the Wentzel-Kramers-Brillouin (WKB) optimal paths of the target process, derived in the limit of low noise. We see that the paths produced using our method, encoding the full statistics of the process, collapse onto the WKB optimal path as the level of noise is reduced. We propose that our method can be used to judge the accuracy of the WKB approximation at finite levels of noise.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Degree irregularity and rank probability bias in network meta-analysis

Network meta-analysis (NMA) is a statistical technique for the comparison of treatment options. The nodes of the network are the competing treatments and edges represent comparisons of treatments in trials. Outcomes of Bayesian NMA include estimates of treatment effects, and the probabilities that each treatment is ranked best, second best and so on. How exactly network geometry affects the accuracy and precision of these outcomes is not fully understood. Here we carry out a simulation study and find that disparity in the number of trials involving different treatments leads to a systematic bias in estimated rank probabilities. This bias is associated with an increased variation in the precision of treatment effect estimates. Using ideas from the theory of complex networks, we define a measure of `degree irregularity' to quantify asymmetry in the number of studies involving each treatment. Our simulations indicate that more regular networks have more precise treatment effect estimates and smaller bias of rank probabilities. We also find that degree regularity is a better indicator of NMA quality than both the total number of studies in a network and the disparity in the number of trials per comparison. These results have implications for planning future trials. We demonstrate that choosing trials which reduce the network's irregularity can improve the precision and accuracy of NMA outcomes.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ecological communities from random generalised Lotka-Volterra dynamics with non-linear feedback

We investigate the outcome of generalised Lotka-Volterra dynamics of ecological communities with random interaction coefficients and non-linear feedback. We show in simulations that the saturation of non-linear feedback stabilises the dynamics. This is confirmed in an analytical generating-functional approach to generalised Lotka-Volterra equations with piecewise linear saturating response. For such systems we are able to derive self-consistent relations governing the stable fixed-point phase, and to carry out a linear stability analysis to predict the onset of unstable behaviour. We investigate in detail the combined effects of the mean, variance and co-variance of the random interaction coefficients, and the saturation value of the non-linear response. We find that stability and diversity increases with the introduction of non-linear feedback, where decreasing the saturation value has a similar effect to decreasing the co-variance. We also find co-operation to no longer have a detrimental effect on stability with non-linear feedback, and the order parameters mean abundance and diversity to be less dependent on the symmetry of interactions with stronger saturation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective diffusion coefficients in reaction-diffusion systems with anomalous transport

We show that the Turing patterns in reaction systems with subdiffusion can be replicated in an effective system with Markovian cross-diffusion. The effective system has the same Turing instability as the original system, and the same patterns. If particles are short-lived, the transient dynamics are captured as well. We use the cross-diffusive system to define effective diffusion coefficients for the system with anomalous transport, and we show how they can be used to efficiently describe the Turing instability. We also demonstrate that the mean squared displacement of a suitably defined ensemble of subdiffusing particles grows linearly with time, with a diffusion coefficient which agrees with our earlier calculations. We verify these deductions by numerically integrating both the fractional reaction-diffusion equation and its normally diffusing counterpart. Our findings suggest that cross-diffusive behaviour can come about as a result of anomalous transport.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Zealots in multi-state noisy voter models

The noisy voter model is a stylised representation of opinion dynamics. Individuals copy opinions from other individuals, and are subject to spontaneous state changes. In the case of two opinion states this model is known to have a noise-driven transition between a unimodal phase, in which both opinions are present, and a bimodal phase in which one of the opinions dominates. The presence of zealots can remove the unimodal and bimodal phases in the model with two opinion states. Here, we study the effects of zealots in noisy voter models with M>2 opinion states on complete interaction graphs. We find that the phase behaviour diversifies, with up to six possible qualitatively different types of stationary states. The presence of zealots removes some of these phases, but not all. We analyse situations in which zealots affect the entire population, or only a fraction of agents, and show that this situation corresponds to a single-community model with a fractional number of zealots, further enriching the phase diagram. Our study is conducted analytically based on effective birth-death dynamics for the number of individuals holding a given opinion. Results are confirmed in numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Motion, fixation probability and the choice of an evolutionary process

Different evolutionary models are known to make disparate predictions for the success of an invading mutant in some situations. For example, some evolutionary mechanics lead to amplification of selection in structured populations, while others suppress it. Here, we use computer simulations to study evolutionary populations moved by flows, and show how the speed of this motion impacts the fixation probability of an invading mutant. Flows of different speeds interpolate between evolutionary dynamics on fixed heterogeneous graphs and in well-stirred populations. We find that the motion has an active role in amplifying or suppressing selection, accomplished by fragmenting and reconnecting the interaction graph. While increasing flow speeds suppress selection for most evolutionary models, we identify characteristic responses to flow for the different update rules we test. We suggest these responses as a potential aid for choosing the most suitable update rule for a given biological system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Birth-death processes with quenched uncertainty and intrinsic noise

The dynamics of populations is frequently subject to intrinsic noise. At the same time unknown interaction networks or rate constants can present quenched uncertainty. Existing approaches often involve repeated sampling of the quenched disorder and then running the stochastic birth-death dynamics on these samples. In this paper we take a different view, and formulate an effective jump process, representative of the ensemble of quenched interactions as a whole. Using evolutionary games with random payoff matrices as an example, we develop an algorithm to simulate this process, and we discuss diffusion approximations in the limit of weak intrinsic noise.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Intrinsic noise in systems with switching environments

We study individual-based dynamics in finite populations, subject to randomly switching environmental conditions. These are inspired by models in which genes transition between on and off states, regulating underlying protein dynamics. Similarly switches between environmental states are relevant in bacterial populations and in models of epidemic spread. Existing piecewise-deterministic Markov process (PDMP) approaches focus on the deterministic limit of the population dynamics while retaining the randomness of the switching. Here we go beyond this approximation and explicitly include effects of intrinsic stochasticity at the level of the linear-noise approximation. Specifically we derive the stationary distributions of a number of model systems, in good agreement with simulations. This improves existing approaches which are limited to the regimes of fast and slow switching.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mechanisms of stochastic onset and termination of atrial fibrillation episodes: Insights using a cellular automaton model

Mathematical models of cardiac electrical excitation are increasingly complex, with multiscale models seeking to represent and bridge physiological behaviours across temporal and spatial scales. The increasing complexity of these models makes it computationally expensive to both evaluate long term (>60 seconds) behaviour and determine sensitivity of model outputs to inputs. This is particularly relevant in models of atrial fibrillation (AF), where individual episodes last from seconds to days, and inter-episode waiting times can be minutes to months. Potential mechanisms of transition between sinus rhythm and AF have been identified but are not well understood, and it is difficult to simulate AF for long periods of time using state-of-the-art models. In this study, we implemented a Moe-type cellular automaton on a novel, topologically correct surface geometry of the left atrium. We used the model to simulate stochastic initiation and spontaneous termination of AF, arising from bursts of spontaneous activation near pulmonary veins. The simplified representation of atrial electrical activity reduced computational cost, and so permitted us to investigate AF mechanisms in a probabilistic setting. We computed large numbers (~10^5) of sample paths of the model, to infer stochastic initiation and termination rates of AF episodes using different model parameters. By generating statistical distributions of model outputs, we demonstrated how to propagate uncertainties of inputs within our microscopic level model up to a macroscopic level. Lastly, we investigated spontaneous termination in the model and found a complex dependence on its past AF trajectory, the mechanism of which merits future investigation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sojourn times and fixation dynamics in multi-player games with fluctuating environments

We study evolutionary multi-player games in finite populations, subject to fluctuating environments. The population undergoes a birth-death process with absorbing states, and the environment follows a Markovian process, resulting in a fluctuating payoff matrix for the evolutionary game. Our focus is on the fixation or extinction of a single mutant in a population of wildtypes. We show that the nonlinear nature of fitnesses in multi-player games gives rise to an intricate interplay of selection, genetic drift and environmental fluctuations. This generates effects not seen in simpler two-player games. To analyse trajectories towards fixation we analytically calculate sojourn times for general birth-death processes in populations of two types of individuals and in fluctuating environments.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic transport in the presence of spatial disorder: fluctuation-induced corrections to homogenization

Motivated by uncertainty quantification in natural transport systems, we investigate an individual-based transport process involving particles undergoing a random walk along a line of point sinks whose strengths are themselves independent random variables. We assume particles are removed from the system via first-order kinetics. We analyse the system using a hierarchy of approaches when the sinks are sparsely distributed, including a stochastic homogenization approximation that yields explicit predictions for the extrinsic disorder in the stationary state due to sink strength fluctuations. The extrinsic noise induces long-range spatial correlations in the particle concentration, unlike fluctuations due to the intrinsic noise alone. Additionally, the mean concentration profile, averaged over both intrinsic and extrinsic noise, is elevated compared with the corresponding profile from a uniform sink distribution, showing that the classical homogenization approximation can be a biased estimator of the true mean.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Assessing Measures of Atrial Fibrillation Clustering via Stochastic Models of Episode Recurrence and Disease Progression

Atrial fibrillation (AF) is a leading cause of morbidity and mortality. AF prevalence increases with age, which is attributed to pathophysiological changes that aid AF initiation and perpetuation. Current state-of-the-art models are only capable of simulating short periods of atrial activity at high spatial resolution, whilst the majority of clinical recordings are based on infrequent temporal datasets of limited spatial resolution. Being able to estimate disease progression informed by both modelling and clinical data would be of significant interest. In addition an analysis of the temporal distribution of recorded fibrillation episodes AF density can provide insights into recurrence patterns. We present an initial analysis of the AF density measure using a simplified idealised stochastic model of a binary time series representing AF episodes. The future aim of this work is to develop robust clinical measures of progression which will be tested on models that generate long-term synthetic data. These measures would then be of clinical interest in deciding treatment strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bursting noise in gene expression dynamics: Linking microscopic and mesoscopic models

The dynamics of short-lived mRNA results in bursts of protein production in gene regulatory networks. We investigate the propagation of bursting noise between different levels of mathematical modelling, and demonstrate that conventional approaches based on diffusion approximations can fail to capture bursting noise. An alternative coarse-grained model, the so-called piecewise deterministic Markov process (PDMP), is seen to outperform the diffusion approximation in biologically relevant parameter regimes. We provide a systematic embedding of the PDMP model into the landscape of existing approaches, and we present analytical methods to calculate its stationary distribution and switching frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generating functionals and Gaussian approximations for interruptible delay reactions

We develop a generating functional description of the dynamics of non-Markovian individual-based systems, in which delay reactions can be terminated before completion. This generalises previous work in which a path-integral approach was applied to dynamics in which delay reactions complete with certainty. We construct a more widely applicable theory, and from it we derive Gaussian approximations of the dynamics, valid in the limit of large, but finite population sizes. As an application of our theory we study predator-prey models with delay dynamics due to gestation or lag periods to reach the reproductive age. In particular we focus on the effects of delay on noise-induced cycles.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Modelling the progression of atrial fibrillation: A stochastic individual-based approach

We propose a stochastic individual-based model of the progression of atrial fibrillation (AF). The model operates at patient level over a lifetime and is based on elements of the physiology and biophysics of AF, making contact with existing mechanistic models. The outputs of the model are times when the patient is in normal rhythm and AF, and we carry out a population-level analysis of the statistics of disease progression. While the model is stylised at present and not directly predictive, future improvements are proposed to tighten the gap between existing mechanistic models of AF, and epidemiological data, with a view towards model-based personalised medicine.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stochastic tunneling and metastable states during the somatic evolution of cancer

Tumors initiate when a population of proliferating cells accumulates a certain number and type of genetic and/or epigenetic alterations. The population dynamics of such sequential acquisition of (epi)genetic alterations has been the topic of much investigation. The phenomenon of stochastic tunneling, where an intermediate mutant in a sequence does not reach fixation in a population before generating a double mutant, has been studied using a variety of computational and mathematical methods. However, the field still lacks a comprehensive analytical description since theoretical predictions of fixation times are only available for cases in which the second mutant is advantageous. Here, we study stochastic tunneling in a Moran model. Analyzing the deterministic dynamics of large populations we systematically identify the parameter regimes captured by existing approaches. Our analysis also reveals fitness landscapes and mutation rates for which finite populations are found in long-lived metastable states. These are landscapes in which the final mutant is not the most advantageous in the sequence, and resulting metastable states are a consequence of a mutation-selection balance. The escape from these states is driven by intrinsic noise, and their location affects the probability of tunneling. Existing methods no longer apply. In these regimes it is the escape from the metastable states that is the key bottleneck; fixation is no longer limited by the emergence of a successful mutant lineage. We used the so-called Wentzel-Kramers-Brillouin method to compute fixation times in these parameter regimes, successfully validated by stochastic simulations. Our work fills a gap left by previous approaches and provides a more comprehensive description of the acquisition of multiple mutations in populations of somatic cells.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

When the mean is not enough: Calculating fixation time distributions in birth-death processes

Studies of fixation dynamics in Markov processes predominantly focus on the mean time to absorption. This may be inadequate if the distribution is broad and skewed. We compute the distribution of fixation times in one-step birth-death processes with two absorbing states. These are expressed in terms of the spectrum of the process, and we provide different representations as forward-only processes in eigenspace. These allow efficient sampling of fixation time distributions. As an application we study evolutionary game dynamics, where invading mutants can reach fixation or go extinct. We also highlight the median fixation time as a possible analog of mixing times in systems with small mutation rates and no absorbing states, whereas the mean fixation time has no such interpretation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fixation in finite populations evolving in fluctuating environments

The environment in which a population evolves can have a crucial impact on selection. We study evolutionary dynamics in finite populations of fixed size in a changing environment. The population dynamics are driven by birth and death events. The rates of these events may vary in time depending on the state of the environment, which follows an independent Markov process. We develop a general theory for the fixation probability of a mutant in a population of wild-types, and for mean unconditional and conditional fixation times. We apply our theory to evolutionary games for which the payoff structure varies in time. The mutant can exploit the environmental noise; a dynamic environment that switches between two states can lead to a probability of fixation that is higher than in any of the individual environmental states. We provide an intuitive interpretation of this surprising effect. We also investigate stationary distributions when mutations are present in the dynamics. In this regime, we find two approximations of the stationary measure. One works well for rapid switching, the other for slowly fluctuating environments.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gaussian approximations for stochastic systems with delay: chemical Langevin equation and application to a Brusselator system

We present a heuristic derivation of Gaussian approximations for stochastic chemical reaction systems with distributed delay. In particular we derive the corresponding chemical Langevin equation. Due to the non-Markovian character of the underlying dynamics these equations are integro-differential equations, and the noise in the Gaussian approximation is coloured. Following on from the chemical Langevin equation a further reduction leads to the linear-noise approximation. We apply the formalism to a delay variant of the celebrated Brusselator model, and show how it can be used to characterise noise-driven quasi-cycles, as well as noise-triggered spiking. We find surprisingly intricate dependence of the typical frequency of quasi-cycles on the delay period.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

When do microscopic assumptions determine the outcome in evolutionary game dynamics?

The modelling of evolutionary game dynamics in finite populations requires microscopic processes that determine how strategies spread. The exact details of these processes are often chosen without much further consideration. Different types of microscopic models, including in particular fitness-based selection rules and imitation-based dynamics, are often used as if they were interchangeable. We challenge this view and investigate how robust these choices on the micro-level really are. Focusing on a key macroscopic observable, the probability for a single mutant to take over a population of wild-type individuals, we show that there is a unique pair of a fitness-based process and an imitation process leading to identical outcomes for arbitrary games and for all intensities of selection. This highlights the perils of making arbitrary choices at the micro-level without regard of the consequences at the macro-level.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Demographic noise and resilience in a semi-arid ecosystem model

The scarcity of water characterising drylands forces vegetation to adopt appropriate survival strategies. Some of these generate water-vegetation feedback mechanisms that can lead to spatial self-organisation of vegetation, as it has been shown with models representing plants by a density of biomass, varying continuously in time and space. However, although plants are usually quite plastic they also display discrete qualities and stochastic behaviour. These features may give rise to demographic noise, which in certain cases can influence the qualitative dynamics of ecosystem models. In the present work we explore the effects of demographic noise on the resilience of a model semi-arid ecosystem. We introduce a spatial stochastic eco-hydrological hybrid model in which plants are modelled as discrete entities subject to stochastic dynamical rules, while the dynamics of surface and soil water are described by continuous variables. The model has a deterministic approximation very similar to previous continuous models of arid and semi-arid ecosystems. By means of numerical simulations we show that demographic noise can have important effects on the extinction and recovery dynamics of the system. In particular we find that the stochastic model escapes extinction under a wide range of conditions for which the corresponding deterministic approximation predicts absorption into desert states.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Pattern formation in individual-based systems with time-varying parameters

We study the patterns generated in finite-time sweeps across symmetry-breaking bifurcations in individual-based models. Similar to the well-known Kibble-Zurek scenario of defect formation, large-scale patterns are generated when model parameters are varied slowly, whereas fast sweeps produce a large number of small domains. The symmetry breaking is triggered by intrinsic noise, originating from the discrete dynamics at the micro-level. Based on a linear-noise approximation, we calculate the characteristic length scale of these patterns. We demonstrate the applicability of this approach in a simple model of opinion dynamics, a model in evolutionary game theory with a time-dependent fitness structure, and a model of cell differentiation. Our theoretical estimates are confirmed in simulations. In further numerical work, we observe a similar phenomenon when the symmetry-breaking bifurcation is triggered by population growth.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stochastic processes with distributed delays: chemical Langevin equation and linear-noise approximation

We develop a systematic approach to the linear-noise approximation for stochastic reaction systems with distributed delays. Unlike most existing work our formalism does not rely on a master equation, instead it is based upon a dynamical generating functional describing the probability measure over all possible paths of the dynamics. We derive general expressions for the chemical Langevin equation for a broad class of non-Markovian systems with distributed delay. Exemplars of a model of gene regulation with delayed auto-inhibition and a model of epidemic spread with delayed recovery provide evidence of the applicability of our results.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Complexity measures, emergence, and multiparticle correlations

We study correlation measures for complex systems. First, we investigate some recently proposed measures based on information geometry. We show that these measures can increase under local transformations as well as under discarding particles, thereby questioning their interpretation as a quantifier for complexity or correlations. We then propose a refined definition of these measures, investigate its properties and discuss its numerical evaluation. As an example, we study coupled logistic maps and study the behavior of the different measures for that case. Finally, we investigate other local effects during the coarse graining of the complex system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Computing complexity measures for quantum states based on exponential families

Given a multiparticle quantum state, one may ask whether it can be represented as a thermal state of some Hamiltonian with k-particle interactions only. The distance from the exponential family defined by these thermal states can be considered as a measure of complexity of a given state. We investigate the resulting optimization problem and show how symmetries can be exploited to simplify the task of finding the nearest thermal state in a given exponential family. We also present an algorithm for the computation of the complexity measure and consider specific examples to demonstrate its applicability.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Demographic noise and piecewise deterministic Markov processes

We explore a class of hybrid (piecewise deterministic) systems characterized by a large number of individuals inhabiting an environment whose state is described by a set of continuous variables. We use analytical and numerical methods from nonequilibrium statistical mechanics to study the influence that intrinsic noise has on the qualitative behavior of the system. We discuss the application of these concepts to the case of semiarid ecosystems. Using a system-size expansion we calculate the power spectrum of the fluctuations in the system. This predicts the existence of noise-induced oscillations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects of communication and utility-based decision making in a simple model of evacuation

We present a simple cellular automaton based model of decision making during evacuation. Evacuees have to choose between two different exit routes, resulting in a strategic decision making problem. Agents take their decisions based on utility functions, these can be revised as the evacuation proceeds, leading to complex interaction between individuals and to jamming transitions. The model also includes the possibility to communicate and exchange information with distant agents, information received may affect the decision of agents. We show that under a wider range of evacuation scenarios performance of the model system as a whole is optimal at an intermediate fraction of evacuees with access to communication.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Effects of noise on convergent game learning dynamics

We study stochastic effects on the lagging anchor dynamics, a reinforcement learning algorithm used to learn successful strategies in iterated games, which is known to converge to Nash points in the absence of noise. The dynamics is stochastic when players only have limited information about their opponents' strategic propensities. The effects of this noise are studied analytically in the case where it is small but finite, and we show that the statistics and correlation properties of fluctuations can be computed to a high accuracy. We find that the system can exhibit quasicycles, driven by intrinsic noise. If players are asymmetric and use different parameters for their learning, a net payoff advantage can be achieved due to these stochastic oscillations around the deterministic equilibrium.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Fluctuations in meta-population exclusion processes

We introduce a meta-population version of models of asymmetric exclusion models, consisting of a spatial arrangement of patches. Patches are of a specific size, indicating the maximal number of particles they can hold. We use an expansion in the inverse patch size to calculate the spectral properties of fluctuations in such systems. This provides a systematic derivation from first principles of effective Langevin descriptions discussed in the literature. We apply our approach to the totally asymmetric simple exclusion process, to variants with an overall constraint on the total particle number and to a two-species exclusion model. The theory provides semi-analytical results, these are confirmed in numerical simulations, and give good approximations to conventional exclusion models. These are recovered when the patch size is set to unity.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Mixing times in evolutionary game dynamics

Without mutation and migration, evolutionary dynamics ultimately leads to the extinction of all but one species. Such fixation processes are well understood and can be characterized analytically with methods from statistical physics. However, many biological arguments focus on stationary distributions in a mutation-selection equilibrium. Here, we address the equilibration time required to reach stationarity in the presence of mutation, this is known as the mixing time in the theory of Markov processes. We show that mixing times in evolutionary games have the opposite behaviour from fixation times when the intensity of selection increases: In coordination games with bistabilities, the fixation time decreases, but the mixing time increases. In coexistence games with metastable states, the fixation time increases, but the mixing time decreases. Our results are based on simulations and the WKB approximation of the master equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Stochastic oscillations of adaptive networks: application to epidemic modelling

Adaptive-network models are typically studied using deterministic differential equations which approximately describe their dynamics. In simulations, however, the discrete nature of the network gives rise to intrinsic noise which can radically alter the system's behaviour. In this article we develop a method to predict the effects of stochasticity in adaptive networks by making use of a pair-based proxy model. The technique is developed in the context of an epidemiological model of a disease spreading over an adaptive network of infectious contact. Our analysis reveals that in this model the structure of the network exhibits stochastic oscillations in response to fluctuations in the disease dynamic.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The mechanics of stochastic slowdown in evolutionary games

We study the stochastic dynamics of evolutionary games, and focus on the so-called `stochastic slowdown' effect, previously observed in (Altrock et. al, 2010) for simple evolutionary dynamics. Slowdown here refers to the fact that a beneficial mutation may take longer to fixate than a neutral one. More precisely, the fixation time conditioned on the mutant taking over can show a maximum at intermediate selection strength. We show that this phenomenon is present in the prisoner's dilemma, and also discuss counterintuitive slowdown and speedup in coexistence games. In order to establish the microscopic origins of these phenomena, we calculate the average sojourn times. This allows us to identify the transient states which contribute most to the slowdown effect, and enables us to provide an understanding of slowdown in the takeover of a small group of cooperators by defectors: Defection spreads quickly initially, but the final steps to takeover can be delayed substantially. The analysis of coexistence games reveals even more intricate behavior. In small populations, the conditional average fixation time can show multiple extrema as a function of the selection strength, e.g., slowdown, speedup, and slowdown again. We classify two-player games with respect to the possibility to observe non-monotonic behavior of the conditional average fixation time as a function of selection strength.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Voter models with conserved dynamics

We propose a modified voter model with locally conserved magnetization and investigate its phase ordering dynamics in two dimensions in numerical simulations. Imposing a local constraint on the dynamics has the surprising effect of speeding up the phase ordering process. The system is shown to exhibit a scaling regime characterized by algebraic domain growth, at odds with the logarithmic coarsening of the standard voter model. A phenomenological approach based on cluster diffusion and similar to Smoluchowski ripening correctly predicts the observed scaling regime. Our analysis exposes unexpected complexity in the phase ordering dynamics without thermodynamic potential.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Complex dynamics in learning complicated games

Game theory is the standard tool used to model strategic interactions in evolutionary biology and social science. Traditional game theory studies the equilibria of simple games. But is traditional game theory applicable if the game is complicated, and if not, what is? We investigate this question here, defining a complicated game as one with many possible moves, and therefore many possible payoffs conditional on those moves. We investigate two-person games in which the players learn based on experience. By generating games at random we show that under some circumstances the strategies of the two players converge to fixed points, but under others they follow limit cycles or chaotic attractors. The dimension of the chaotic attractors can be very high, implying that the dynamics of the strategies are effectively random. In the chaotic regime the payoffs fluctuate intermittently, showing bursts of rapid change punctuated by periods of quiescence, similar to what is observed in fluid turbulence and financial markets. Our results suggest that such intermittency is a highly generic phenomenon, and that there is a large parameter regime for which complicated strategic interactions generate inherently unpredictable behavior that is best described in the language of dynamical systems theory

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cycles of cooperation and defection in imperfect learning

When people play a repeated game they usually try to anticipate their opponents' moves based on past observations, and then decide what action to take next. Behavioural economics studies the mechanisms by which strategic decisions are taken in these adaptive learning processes. We here investigate a model of learning the iterated prisoner's dilemma game. Players have the choice between three strategies, always defect (ALLD), always cooperate (ALLC) and tit-for-tat (TFT). The only strict Nash equilibrium in this situation is ALLD. When players learn to play this game convergence to the equilibrium is not guaranteed, for example we find cooperative behaviour if players discount observations in the distant past. When agents use small samples of observed moves to estimate their opponent's strategy the learning process is stochastic, and sustained oscillations between cooperation and defection can emerge. These cycles are similar to those found in stochastic evolutionary processes, but the origin of the noise sustaining the oscillations is different and lies in the imperfect sampling of the opponent's strategy. Based on a systematic expansion technique, we are able to predict the properties of these learning cycles, providing an analytical tool with which the outcome of more general stochastic adaptation processes can be characterised.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Fixation and escape times in stochastic game learning

Evolutionary dynamics in finite populations is known to fixate eventually in the absence of mutation. We here show that a similar phenomenon can be found in stochastic game dynamical batch learning, and investigate fixation in learning processes in a simple 2x2 game, for two-player games with cyclic interaction, and in the context of the best-shot network game. The analogues of finite populations in evolution are here finite batches of observations between strategy updates. We study when and how such fixation can occur, and present results on the average time-to-fixation from numerical simulations. Simple cases are also amenable to analytical approaches and we provide estimates of the behaviour of so-called escape times as a function of the batch size. The differences and similarities with escape and fixation in evolutionary dynamics are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Individual impact of agent actions in financial markets

We present an analysis of the price impact associated with trades effected by different financial firms. Using data from the Spanish Stock Market, we find a high degree of heterogeneity across different market members, both in the instantaneous impact functions and in the time-dependent market response to trades by individual members. This heterogeneity is statistically incompatible with the existence of market-wide universal impact dynamics which apply uniformly to all trades and suggests that rather, market dynamics emerge from the complex interaction of different behaviors of market participants. Several possible reasons for this are discussed, along with potential extensions one may consider to increase the range of applicability of existing models of market impact.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Learning to play public good games

We extend recent analyses of stochastic effects in game dynamical learning to cases of multi-player games, and to games defined on networked structures. By means of an expansion in the noise strength we consider the weak-noise limit, and present an analytical computation of spectral properties of fluctuations in multi-player public good games. This extends existing work on two-player games. In particular we show that coherent cycles may emerge driven by noise in the adaptation dynamics. These phenomena are not too dissimilar from cyclic strategy switching observed in experiments of behavioural game theory.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Optimizing evacuation flow in a two-channel exclusion process

We use a basic setup of two coupled exclusion processes to model a stylised situation in evacuation dynamics, in which evacuees have to choose between two escape routes. The coupling between the two processes occurs through one common point at which particles are injected, the process can be controlled by directing incoming individuals into either of the two escape routes. Based on a mean-field approach we determine the phase behaviour of the model, and analytically compute optimal control strategies, maximising the total current through the system. Results are confirmed by numerical simulations. We also show that dynamic intervention, exploiting fluctuations about the mean-field stationary state, can lead to a further increase in total current.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stochastic waves in a Brusselator model with nonlocal interaction

We show that intrinsic noise can induce spatio-temporal phenomena such as Turing patterns and travelling waves in a Brusselator model with nonlocal interaction terms. In order to predict and to characterize these quasi-waves we analyze the nonlocal model using a system-size expansion. The resulting theory is used to calculate the power spectra of the quasi-waves analytically, and the outcome is tested successfully against simulations. We discuss the possibility that nonlocal models in other areas, such as epidemic spread or social dynamics, may contain similar stochastically-induced patterns.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Evolutionary dynamics, intrinsic noise and cycles of co-operation

We use analytical techniques based on an expansion in the inverse system size to study the stochastic evolutionary dynamics of finite populations of players interacting in a repeated prisoner's dilemma game. We show that a mechanism of amplification of demographic noise can give rise to coherent oscillations in parameter regimes where deterministic descriptions converge to fixed points with complex eigenvalues. These quasi-cycles between co-operation and defection have previously been observed in computer simulations; here we provide a systematic and comprehensive analytical characterization of their properties. We are able to predict their power spectra as a function of the mutation rate and other model parameters, and to compare the relative magnitude of the cycles induced by different types of underlying microscopic dynamics. We also extend our analysis to the iterated prisoner's dilemma game with a win-stay lose-shift strategy, appropriate in situations where players are subject to errors of the trembling-hand type.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Imitation, internal absorption and the reversal of local drift in stochastic evolutionary games

Evolutionary game dynamics in finite populations is typically subject to noise, inducing effects which are not present in deterministic systems, including fixation and extinction. In the first part of this paper we investigate the phenomenon of drift reversal in finite populations, taking into account that drift is a local quantity in strategy space. Secondly, we study a simple imitation dynamics, and show that it can lead to fixation at internal mixed-strategy fixed points even in finite populations. Imitation in infinite populations is adequately described by conventional replicator dynamics, and these equations are known to have internal fixed points. Internal absorption in finite populations on the other hand is a novel dynamic phenomenon. Due to an outward drift in finite populations this type of dynamic arrest is not found in other commonly studied microscopic dynamics, not even in those with the same deterministic replicator limit as imitation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Intrinsic fluctuations in stochastic delay systems: theoretical description and application to a simple model of gene regulation

The effects of intrinsic noise on stochastic delay systems is studied within an expansion in the inverse system size. We show that the stochastic nature of the underlying dynamics may induce oscillatory behaviour in parameter ranges where the deterministic system does not sustain cycles, and compute the power spectra of these stochastic oscillations analytically, in good agreement with simulations. The theory is developed in the context of a simple one-dimensional toy model, but is applicable more generally. Gene regulatory systems in particular often contain only a small number of molecules, leading to significant fluctuations in mRNA and protein concentrations. As an application we therefore study a minimalistic model of the expression levels of hes1 mRNA and Hes1 protein, representing the simple motif of an auto-inhibitory feedback loop and motivated by its relevance to somite segmentation.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Intrinsic noise in game dynamical learning

Demographic noise has profound effects on evolutionary and population dynamics, as well as on chemical reaction systems and models of epidemiology. Such noise is intrinsic and due to the discreteness of the dynamics in finite populations. We here show that similar noise-sustained trajectories arise in game dynamical learning, where the stochasticity has a different origin: agents sample a finite number of moves of their opponents in-between adaptation events. The limit of infinite batches results in deterministic modified replicator equations, whereas finite sampling leads to a stochastic dynamics. The characteristics of these fluctuations can be computed analytically using methods from statistical physics, and such noise can affect the attractors significantly, leading to noise-sustained cycling or removing periodic orbits of the standard replicator dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Limit cycles, complex Floquet multipliers and intrinsic noise

We study the effects of intrinsic noise on chemical reaction systems, which in the deterministic limit approach a limit cycle in an oscillatory manner. Previous studies of systems with an oscillatory approach to a fixed point have shown that the noise can transform the oscillatory decay into sustained coherent oscillations with a large amplitude. We show that a similar effect occurs when the stable attractors are limit cycles. We compute the correlation functions and spectral properties of the fluctuations in suitably co-moving Frenet frames for several model systems including driven and coupled Brusselators, and the Willamowski-Roessler system. Analytical results are confirmed convincingly in numerical simulations. The effect is quite general, and occurs whenever the Floquet multipliers governing the stability of the limit cycle are complex, with the amplitude of the oscillations increasing as the instability boundary is approached.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Minority games, evolving capitals and replicator dynamics

We discuss a simple version of the Minority Game (MG) in which agents hold only one strategy each, but in which their capitals evolve dynamically according to their success and in which the total trading volume varies in time accordingly. This feature is known to be crucial for MGs to reproduce stylised facts of real market data. The stationary states and phase diagram of the model can be computed, and we show that the ergodicity breaking phase transition common for MGs, and marked by a divergence of the integrated response is present also in this simplified model. An analogous majority game turns out to be relatively void of interesting features, and the total capital is found to diverge in time. Introducing a restraining force leads to a model akin to replicator dynamics of evolutionary game theory, and we demonstrate that here a different type of phase transition is observed. Finally we briefly discuss the relation of this model with one strategy per player to more sophisticated Minority Games with dynamical capitals and several trading strategies per agent.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Effects of noise and confidence thresholds in nominal and metric Axelrod dynamics of social influence

We study the effects of bounded confidence thresholds and of interaction and external noise on Axelrod's model of social influence. Our study is based on a combination of numerical simulations and an integration of the mean-field Master equation describing the system in the thermodynamic limit. We find that interaction thresholds affect the system only quantitatively, but that they do not alter the basic phase structure. The known crossover between an ordered and a disordered state in finite systems subject to external noise persists in models with general confidence threshold. Interaction noise here facilitates the dynamics and reduces relaxation times. We also study Axelrod systems with metric features, and point out similarities and differences compared to models with nominal features. Metric features are used to demonstrate that a small group of extremists can have a significant impact on the opinion dynamics of a population of Axelrod agents.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

Two-population replicator dynamics and number of Nash equilibria in random matrix games

We study the connection between the evolutionary replicator dynamics and the number of Nash equilibria in large random bi-matrix games. Using techniques of disordered systems theory we compute the statistical properties of both, the fixed points of the dynamics and the Nash equilibria. Except for the special case of zero-sum games one finds a transition as a function of the so-called co-operation pressure between a phase in which there is a unique stable fixed point of the dynamics coinciding with a unique Nash equilibrium, and an unstable phase in which there are exponentially many Nash equilibria with statistical properties different from the stationary state of the replicator equations. Our analytical results are confirmed by numerical simulations of the replicator dynamics, and by explicit enumeration of Nash equilibria.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

SUSY Ward identities in 1-loop perturbation theory

We present preliminary results of a study of the supersymmetric (SUSY) Ward identities (WIs) for the N=1 SU(2) SUSY Yang-Mills theory in the context of one-loop lattice perturbation theory. The supersymmetry on the lattice is explicitly broken by the gluino mass and the lattice artifacts. However, the renormalization of the supercurrent can be carried out in a scheme that restores the nominal continuum WIs. The perturbative calculation of the renormalization constants and mixing coefficients for the local supercurrent is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

On the 1-loop lattice perturbation theory of the supersymmetric Ward identities

The one loop corrections to the supersymmetric Ward identities (WIs) in the discretized N=1 SU(2) supersymmetric Yang-Mills theory can be investigated by means of lattice perturbation theory. The supersymmetry (SUSY) is explicitly broken by the lattice discretization as well as by the introduction of Wilson fermions. However, the renormalization of the supercurrent can be carried out in a scheme that restores the nominal continuum WIs. We present our work in progress which is concerned with the 1-loop renormalization of the local supercurrent, i.e. with the perturbative computation of the corresponding renormalization constants and mixing coefficients.

0 citations0 reviewsNo code linkedOpen access