Efficient Simulations of Individual Based Models for Adaptive Dynamics and the Canonical Equation

Efficient Simulations of Individual Based Models for Adaptive Dynamics and the Canonical Equation
Vaibhav Madhok

We propose a faster algorithm for individual based simulations for adaptive dynamics based on a simple modification to the standard Gillespie Algorithm for simulating stochastic birth-death processes. We provide an analytical explanation that shows that simulations based on the modified algorithm, in the deterministic limit, lead to the same equations of adaptive dynamics as well as same conditions for evolutionary branching as those obtained from the standard Gillespie algorithm. Based on this algorithm, we provide an intuitive and simple interpretation of the canonical equation of adaptive dynamics. With the help of examples we compare the performance of this algorithm to the standard Gillespie algorithm and demonstrate its efficiency. We also study an example using this algorithm to study evolutionary dynamics in a multi-dimensional phenotypic space and study the question of predictability of evolution.

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s