deterministic, synchronous simulation methodology  (Rocha labs)

Journal: The American naturalist

Article Title: Evolution at the edge of expanding populations

doi: 10.1086/704594

Figure Lengend Snippet: Faster-dispersing species take over expanding populations despite having a slower growth as long as v2 > v1. A: Phase diagram of evolutionary outcome in r2/r1,D2/D1 space, with 2 subdomains of interest in the success domain: (1) is the domain of success with higher growth rate and slower dispersal. (2) is the domain of success with faster dispersal but lower growth rate. The (+) symbol represents the hyperswarmer phenotype with respect to wild-type phenotype, as measured in the experimental conditions (Table 1). The (*) symbol represents the wild-type phenotype with respect to the hyperswarmer phenotype. B: Fixation probability obtained from the stochastic model with death rate (stochastic simulations performed with S = 1, K = 100, L = 2λ1). C: Takeover time obtained from the deterministic model (simulations performed with L = 2λ1 and S = 0.2). The empty orange circles represent the conditions where species 2 fails to takeover within the duration of the simulation (r1t = 180). In A-C, the red dot depicts the reference point (D2 = D1 and r2 = r1).

Article Snippet: Numerical simulations The deterministic numerical simulations (used for , , and ) were performed in MATLAB (The MathWorks) following Euler’s method, with dx = 0.1, dt = 0.001, D 1 = 1, r 1 = 1, and total spatial range of 400.

Techniques: