9867

Stochastic Rock-Paper-Scissors Population Dynamics

When three strategies (A, B, and C) compete for the same resources, a property called cyclic dominance can often be observed: the presence of As hinders the reproduction of Bs, Bs have an adverse effect on the reproduction of Cs, and Cs have a similar advantage over As. Real-life examples have been studied in colonies of bacteria (Kerr et al., 2002) and in lizard populations (Sinervo and Lively, 1996). You can examine the oscillating population structure that results from this rule and find out why random effects have a critical influence, as derived in (Reichenbach et al., 2006; 2008), by running this real-time simulation.

SNAPSHOTS

  • [Snapshot]
  • [Snapshot]
  • [Snapshot]

DETAILS

The core assumption of the rock-paper-scissors model is that members of the A species survive longer in populations with higher B counts ("A beats B") and lower C counts ("C beats A"), and that cyclic equivalents of this rule hold for the B and C survival rates.
In an ecosystem dominated by As, for example, the B population will dwindle. Since this makes it easier for Cs to multiply, however, the A count will suffer in the long term, helping the B species to recover, etc. These oscillations would continue ad infinitum in a simple deterministic model, although with ever-increasing amplitudes.
To mirror the random influences that always affect real ecosystems, evolutionary selection (death) and reproduction can both be modeled as separate "chemical reactions" that occur randomly in a four-type population (A, B, C, and "empty" for any unused resources in the ecosystem):
1. In a reproduction process, two random individuals are chosen. If either one is "empty", it is replaced with the species of the other individual.
2. A selection process only occurs if two randomly chosen individuals belong to different species. Then, the weaker one (in the sense of the rock-paper-scissors principle described above) is replaced with "empty".
The resulting deviations from the deterministic case described above may be small, especially for large populations, but eventually a species with low population count will become extinct. As one of the remaining species has a clear advantage over the other one due to the rock-paper-scissors rules, the game now has an obvious winner that will be able to take over all of the system.
In this Demonstration, the recent population history is visualized as a trajectory on a triangular simplex: close proximity to one of its vertices signals that a large share of the total population belongs to the species indicated next to the vertex. How close the system is to its maximum carrying capacity is shown by the slider below the simplex.
Initially the population consists of an equal number of As, Bs, and Cs. As no species is favored in this situation, the trajectory remains at the fixed point in the simplex center until random deviations push it toward one of the corners. Having orbited around the simplex center—corresponding to the oscillations of the different species—the trajectory will reach an edge when one species becomes extinct. It will then approach the vertex of the winning species and stay there forever.
You can set the initial conditions by clicking into the simplex or dragging the population slider, or use the sliders to the left of the simplex to change the overall rates at which selection and reproduction occur. By choosing the "instant reproduction" option, you can couple both processes so that As directly replace dead Bs, Bs replace dead Cs, and Cs replace dead As. Instant reproduction keeps the total population count at the initial value and slows down the outward movement of the trajectory because species with a large population share are now less likely to profit from deaths in other subpopulations than in the original case.
The spiral patterns that emerge when a spatial version of this scenario is played out on a grid lattice are shown in another Demonstrations project, "Biodiversity in Spatial Rock-Paper-Scissors Games".
References:
B. Kerr, M. A. Riley, M. W. Feldman, and B. J. M. Bohannan, "Local Dispersal Promotes Biodiversity in a Real-Life Game of Rock-Paper-Scissors," Nature, 418, 2002 pp. 171–174.
B. Sinervo and C. M. Lively, "The Rock-Scissors-Paper Game and the Evolution of Alternative Male Strategies," Nature, 380, 1996 pp. 240-243.
T. Reichenbach, M. Mobilia, and E. Frey, "Coexistence versus Extinction in the Stochastic Cyclic Lotka-Volterra Model," Phys. Rev. E, 74(5), 2006.
T. Reichenbach, M. Mobilia, and E. Frey, "Self-Organization of Mobile Populations in Cyclic Competition," J. Theor. Biol., 254(2), 2008 pp. 368–383.
    • Share:

Embed Interactive Demonstration New!

Just copy and paste this snippet of JavaScript code into your website or blog to put the live Demonstration on your site. More details »

Files require Wolfram CDF Player or Mathematica.









 
RELATED RESOURCES
Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+