Stochastic dynamics in finite populations: Difference between revisions

Stochastic differential equations (SDE) provide a general framework to describe the evolutionary dynamics of an arbitrary number of types in finite populations, which results in demographic noise, and to incorporate mutations. For large, but finite populations this allows to include demographic noise without requiring explicit simulations. Instead, the population size only rescales the amplitude of the noise. Moreover, this framework admits the inclusion of mutations between different types, provided that mutation rates, \(\mu\), are not too small compared to the inverse population size \(1/N\). This ensures that all types are almost always represented in the population and that the occasional extinction of one type does not result in an extended absence of that type. For \(\mu N\ll1\) this limits the use of SDE’s, but in this case well established alternative approximations are available based on time scale separation. We illustrate our approach by a [[Rock-Scissors-Paper game]] with mutations, where we demonstrate excellent agreement with simulation based results for sufficiently large populations. In the absence of mutations the excellent agreement extends to small population sizes.

[[Image:Stochastic dynamics - noise term Cxx, no mutations.png|300px|thumb|Value of the element \(\mathcal C_{xx}(x,y,z)\) of the noise matrix \(\mathcal C(\mathbf x)\) for \(d = 3\) strategies and \(\mu = 0\). \(\mathcal C_{xx}(x,y,z)\) determines how the noise in the \(x\)-direction affects the \(x\)-coordinate. In the case of \(\mu = 0\), this noise vanishes for \(x\to0\). For \(y\to0\) and \(z\to0\) we recover the usual multiplicative noise from one-dimensional evolutionary processes.]]

Comparisons between the deterministic dynamics in infinite populations, the stochastic dynamics in finite populations and individual based simulations focus on the [[Rock-Paper-Scissors game]] with a generic payoff

\begin{pmatrix}0 & \frac{s}{2} & -1 \\

\frac{1+s}{3} & -1 & 0\end{pmatrix}\end{matrix}.

In the limit \(N\to\infty\) with \(s=1.4\) and without mutations, \(\mu=0\), \(\hat x\) is an attractor of the replicator dynamics. The figure shows a sample trajectory that spirals towards the interior fixed point \(\hat x\).

# Traulsen, A., Claussen, J. C. & Hauert, C. (2012) Stochastic differential equations for evolutionary dynamics with demographic noise and mutations. ''Phys. Rev. E'' in print.