Simulate the Moran process for evolutionary dynamics. Calculate fixation probabilities, expected fixation times, and fitness-dependent selection in finite populations.
The Moran Process Calculator models evolutionary dynamics in finite populations using the classic Moran process — a fundamental model in mathematical biology. It computes fixation probability (the chance a mutant type takes over the entire population) and expected fixation time under frequency-dependent or constant selection. That makes it a useful bridge between qualitative evolutionary intuition and the quantitative behavior of finite populations.
In the Moran process, at each time step one individual is chosen to reproduce (proportional to fitness) and one random individual dies. This birth-death process models genetic drift and selection in populations of constant size N. A single mutant with fitness r in a population of N-1 wild-type individuals fixes with probability (1 - 1/r) / (1 - 1/rᴺ) for r ≠ 1.
Enter the population size, initial mutant count, relative fitness, and explore how selection strength and population size affect evolutionary outcomes. Compare neutral drift (r=1) to advantageous, deleterious, and strongly selected mutations.
Use this calculator when you want fixation probability and takeover time from a finite-population model instead of relying on qualitative intuition alone. It is useful for population genetics, evolutionary game theory, and mathematical biology work where drift and selection need to be compared explicitly. That makes the long-run behavior easier to discuss in a quantitative way.
Fixation probability (r ≠ 1): ρᵢ = (1 - 1/rⁱ) / (1 - 1/rᴺ). Neutral (r = 1): ρᵢ = i/N. Expected fixation time (conditional, r = 1): t = -N × Σ_{j=1}^{N-1} [(N-j)/j] × ln(1 - j/N) (approximation).
Result: Fixation probability ≈ 9.5%
A single mutant (i=1) with 5% fitness advantage (r=1.05) in a population of 100: ρ₁ = (1 - 1/1.05) / (1 - 1/1.05^100) = 0.0476 / 0.5017 ≈ 0.095 or 9.5%. Under neutral drift (r=1), the probability would be only 1/100 = 1%.
The Moran process is one of two canonical models for finite-population evolution (the other being the Wright-Fisher model). Unlike Wright-Fisher, which has non-overlapping generations, the Moran process uses overlapping generations with continuous birth-death events. This makes it more realistic for many biological populations.
The key insight is that even beneficial mutations don't always fix — they can be lost by random drift, especially when rare. Conversely, deleterious mutations can fix, especially in small populations. This tension between selection (deterministic) and drift (stochastic) is central to evolutionary biology.
The Moran process on graphs extends the model to structured populations where individuals interact only with neighbors. The population structure can dramatically affect fixation probabilities. "Amplifiers of selection" (like the star graph) increase the fixation probability of beneficial mutants. "Suppressors of selection" make evolution more neutral-like.
This has implications for understanding evolution in spatially structured environments such as tissues (cancer), biofilms (bacteria), and social networks (cultural evolution).
Tumors evolve by a process similar to the Moran process: cells divide, acquire mutations, and compete for space. Because tumor populations can be relatively small (especially early), genetic drift plays a significant role alongside selection. Understanding fixation probabilities helps predict resistance evolution and design treatment strategies.
A stochastic model of evolution in a finite population of constant size N. Each step: one individual reproduces (fitness-weighted) and one dies (random). Over time, one type eventually fixates (reaches 100%). It's used to model genetic drift and selection.
The probability that a mutant type starting at frequency i/N eventually replaces all wild-type individuals. For a neutral mutation (r=1), this equals i/N (e.g., 1/N for a single mutant). Beneficial mutations fix with higher probability.
A mutant with relative fitness r > 1 has enhanced fixation probability: ρ₁ ≈ 1 - 1/r for large N. Even a 1% advantage (r=1.01) dramatically increases fixation probability in large populations. Deleterious mutants (r < 1) can still fix by drift, especially in small populations.
When r = 1, all individuals have equal fitness, and evolution is purely random (genetic drift). A single mutant in a population of N fixes with probability 1/N. In a population of 1000, there's a 0.1% chance — small but nonzero.
Larger populations make selection more efficient and drift less important. In large N: beneficial mutations fix more reliably, deleterious ones are eliminated more effectively, and neutral drift takes longer. In small N: drift dominates and even harmful mutations can fix.
Modeling evolution on graphs (structured populations), cancer tumor dynamics, antibiotic resistance evolution, game theory (evolutionary games), immune system dynamics, and cultural evolution. It's the simplest birth-death process with selection.