GillespieSSA2: 'GillespieSSA2': Gillespie's Stochastic Simulation Algorithm...

Description Details The stochastic simulation algorithm Contents of this package See Also


GillespieSSA2 is a fast, scalable, and versatile framework for simulating large systems with Gillespie's Stochastic Simulation Algorithm (SSA). This package is the spiritual successor to the GillespieSSA package originally written by Mario Pineda-Krch.


GillespieSSA2 has the following added benefits:

The SSA methods currently implemented are: Exact (ssa_exact()), Explicit tau-leaping (ssa_etl()), and the Binomial tau-leaping (ssa_btl()).

The stochastic simulation algorithm

The stochastic simulation algorithm (SSA) is a procedure for constructing simulated trajectories of finite populations in continuous time. If X_i(t) is the number of individuals in population i (i = 1,...,N) at time t, the SSA estimates the state vector X(t) = (X_1(t),...,X_N(t)), given that the system initially (at time t_0) was in state X(t_0) = x_0.

Reactions are single instantaneous events changing at least one of the populations (e.g. birth, death, movement, collision, predation, infection, etc). These cause the state of the system to change over time.

The SSA procedure samples the time tau to the next reaction R_j (j = 1,...,M) and updates the system state X(t) accordingly.

Each reaction R_j is characterized mathematically by two quantities; its state-change vector \bm{ν_j} and its propensity function a_j(\mathbf{x}). The state-change vector is defined as nu_j = (nu_1j,...,nu_Nj), where nu_ij is the change in the number of individuals in population i caused by one reaction of type j. The propensity function is defined as a_j(\mathbf{x}), where a_j(x)dt is the probability that a particular reaction j will occur in the next infinitesimal time interval [t,t+dt].

Contents of this package

See Also

ssa() for more explanation on how to use GillespieSSA2

GillespieSSA2 documentation built on May 18, 2021, 5:06 p.m.