Moran Process Calculator

Simulate the Moran process for evolutionary dynamics. Calculate fixation probabilities, expected fixation times, and fitness-dependent selection in finite populations.

Moran Process Calculator

Total number of individuals
Number of mutant individuals at start
1 = neutral, >1 = advantageous, <1 = deleterious
Fixation Probability
4.7984%
ρ = 0.047984
Neutral Fixation
1.0000%
i/N = 1/100
Selection Advantage
4.80×
Fixation boost vs neutral drift
Selection Coefficient
+0.0500
Advantageous
Loss Probability
95.20%
Probability mutant goes extinct
Initial Frequency
1.0%
1 of 100 individuals
Fixation: 4.80% | Extinction: 95.20%

Fitness vs Fixation Probability (i=1, N=100)

Fitness (r)Fixation %vs NeutralProbability
0.500.000%0.00×
0.800.000%0.00×
0.900.000%0.00×
0.950.031%0.03×
0.990.583%0.58×
1.001.000%1.00×
1.011.571%1.57×
1.054.798%4.80×
1.109.092%9.09×
1.2016.667%16.67×
1.5033.333%33.33×
2.0050.000%50.00×

Population Size Effect (r=1.05, i=1)

NNeutralSelectedRatio
520.00%21.998%1.10×
1010.00%12.334%1.23×
205.00%7.642%1.53×
502.00%5.217%2.61×
1001.00%4.798%4.80×
2000.50%4.762%9.52×
5000.20%4.762%23.81×
10000.10%4.762%47.62×

Initial Count vs Fixation (r=1.05, N=100)

iFreqFixation %Probability
11.0%4.80%
22.0%9.37%
33.0%13.72%
44.0%17.87%
55.0%21.81%
66.0%25.57%
77.0%29.15%
88.0%32.56%
99.0%35.81%
1010.0%38.90%
1111.0%41.85%
1212.0%44.66%
1313.0%47.33%
1414.0%49.87%
1515.0%52.30%
1616.0%54.60%
1717.0%56.80%
1818.0%58.90%
1919.0%60.89%
2020.0%62.79%
Planning notes, formulas, and examples

About the Moran Process Calculator

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.

When This Page Helps

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.

How to Use the Inputs

  1. Enter the total population size (N).
  2. Enter the initial number of mutant individuals (i).
  3. Enter the relative fitness of the mutant type (r).
  4. View fixation probability and expected fixation time.
  5. Compare different fitness values in the sweep table.
  6. Use presets for common evolutionary scenarios.
  7. Examine how population size affects fixation probability.
Formula used
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).

Example Calculation

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%.

Tips & Best Practices

  • For a single mutant (i=1) with large N: ρ₁ ≈ 1 - 1/r when r > 1 (strongly selected).
  • The "1/3 rule" in evolutionary game theory: cooperation is favored if benefit/cost > N/3.
  • Population structure (graphs) can amplify or suppress selection compared to well-mixed populations.
  • Fixation time scales as N² for neutral drift and N·ln(N) for strongly advantageous mutations.
  • In cancer evolution, tumor populations are small enough that drift significantly affects which mutations fix.
  • Compare neutral (r=1) to selected (r≠1) to understand the relative roles of drift and selection.

Evolutionary Dynamics Fundamentals

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.

Structured Populations and Evolutionary Graph Theory

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).

Applications in Cancer Biology

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.

Sources & Methodology

Last updated:

Frequently Asked Questions

  • 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.