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

Abstract

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.

Published in PRE. 14 pages, 6 figures

References in corpus (4)

• Evolutionary games on graphs
• Fixation of strategies for an evolutionary game in finite populations
• Fixation times in evolutionary games under weak selection
• Mixing times in evolutionary game dynamics

Cited by in corpus (7)

• Fitness dependence of the fixation-time distribution for evolutionary dynamics on graphs
• Takeover times for a simple model of network infection
• Reaction rates and the noisy saddle-node bifurcation: Renormalization group for barrier crossing
• Decentralized decision making by an ant colony: drift-diffusion model of individual choice, quorum and collective decision
• Optimal control of stochastic reaction networks with entropic control cost and emergence of mode-switching strategies
• First-passage on disordered intervals
• Intermediate time scale in the first product formation time distribution of Michaelis-Menten kinetics with inhibitors