branching process
Neher RA & Shraiman BI 2011 Genetic draft and quasi-neutrality in large facultatively sexual populations. Genetics 188:975-996.
- simulation of the branching process
- the simulation of the population dynamics is complemented by a simulation of the branching process that can be directly compared to analytic calculations
- we simulate the process described by Equation 1, using an event-driven algorithm
- the simulation keeps track of all individuals that currently carry the allele in question
- for each individual, the time t + Δt of the next event is determined by drawing a Δt from an exponential distribution with parameter B(t) + D + r
- B and D are the birth and death rates
- at time t + Δt, a birth, death, or recombination event is performed with probabilities proportional to B(t + Δt), D, and r, respectively
- the waiting-time distribution for the next event is not exactly exponential, since the birth rate is time dependent
- this, however, amounts only to a correction of order σ2 ≪ 1
