Figure 1

Pseudocode for sampling algorithm. A biochemical signal x(t) obtained from Gillespie SSA simulation can be characterized using two vectors x = [x₁, x₂, ..., x _M ] and t = [t₁, t₂, ..., t _M ]. Time points t _i , i = 1, 2, ..., M, define when the concentration has changed and x _i , i = 1, 2, ..., M, describe the concentration within the interval [t _i , t _i +1), i = 1, 2, ..., M -1. These vectors, as well as the starting time for sampling (T), the number of samples (N), and the time step (Δt), are given as parameters for the algorithm.