scholarly journals Continuum Approximation of Invasion Probabilities

2017 ◽  
Author(s):  
Rebecca K. Borchering ◽  
Scott A. McKinley

AbstractIn the last decade there has been growing criticism of the use of Stochastic Differential Equations (SDEs) to approximate discrete state-space, continuous-time Markov chain population models. In particular, several authors have demonstrated the failure of Diffusion Approximation, as it is often called, to approximate expected extinction times for populations that start in a quasi-stationary state.In this work we investigate a related, but distinct, population dynamics property for which Diffusion Approximation fails: invasion probabilities. We consider the situation in which a few individual are introduced into a population and ask whether their collective lineage can successfully invade. Because the population count is so small during the critical period of success or failure, the process is intrinsically stochastic and discrete. In addition to demonstrating how and why the Diffusion Approximation fails in the large population limit, we contrast this analysis with that of a sometimes more successful alternative WKB-like approach. Through numerical investigations, we also study how these approximations perform in an important intermediate regime. In a surprise, we find that there are times when the Diffusion Approximation performs well: particularly when parameters are near-critical and the population size is small to intermediate.

Author(s):  
Tabea Waizmann ◽  
Luca Bortolussi ◽  
Andrea Vandin ◽  
Mirco Tribastone

Stochastic reaction networks are a fundamental model to describe interactions between species where random fluctuations are relevant. The master equation provides the evolution of the probability distribution across the discrete state space consisting of vectors of population counts for each species. However, since its exact solution is often elusive, several analytical approximations have been proposed. The deterministic rate equation (DRE) gives a macroscopic approximation as a compact system of differential equations that estimate the average populations for each species, but it may be inaccurate in the case of nonlinear interaction dynamics. Here we propose finite-state expansion (FSE), an analytical method mediating between the microscopic and the macroscopic interpretations of a stochastic reaction network by coupling the master equation dynamics of a chosen subset of the discrete state space with the mean population dynamics of the DRE. An algorithm translates a network into an expanded one where each discrete state is represented as a further distinct species. This translation exactly preserves the stochastic dynamics, but the DRE of the expanded network can be interpreted as a correction to the original one. The effectiveness of FSE is demonstrated in models that challenge state-of-the-art techniques due to intrinsic noise, multi-scale populations and multi-stability.


1989 ◽  
Vol 26 (04) ◽  
pp. 880-885 ◽  
Author(s):  
K. Kämmerle

In this paper a bisexual Moran model is introduced. The population consists of N pairs of individuals. At times t = 1, 2, ·· ·two individuals are born, who ‘choose their parents randomly' and independently of each other. Then one of the pairs is removed and replaced by the two individuals born at that instant. The extinction probability of the descendants of a single pair and the number of ancestors of a whole generation are studied. A limit result for large population sizes has been derived by diffusion approximation methods.


2019 ◽  
Vol 340 ◽  
pp. 222-232 ◽  
Author(s):  
Liang-Qun Li ◽  
Xiao-Li Wang ◽  
Wei-Xin Xie ◽  
Zong-Xiang Liu

Sign in / Sign up

Export Citation Format

Share Document