Study of dimerization system, initial values x = [301, 0], parameters k = [1.66 × 10−3, 0.2]. (a) Single SSA realisation. (b) Average of 100 000 SSA simulations. (c) and (d) Histogram of SSA runs (grey bars) and probability density of normal distribution (blue line) calculated from mean and variance of SSA runs corresponding to points c and d in (b). (e) Mean for both variables, calculated using SSA, the moment approximation using 1 moment (deterministic) and two central moments (2m). (f) Variance of x 1 calculated using SSA, two central moments (2m), three central moments with Gaussian closure (3m), and four central moments (4m). (g) Third central moment calculated using SSA and moment approximation method, fourth central moment calculated using SSA, and moment approximation method.
Study of Michaelis-Menten kinetics, with parameters , and initial conditions, S(0) = 301, P(0) = 0, and E(0) = 120. (a) Single SSA realisation. Trajectories calculated using (b) moment approximation including only the mean (deterministic). (c) and (d) Variance of S and covariance between S and P calculated using SSA and approximation using 2 moments. (e) Skewness of S calculated using SSA and approximation with 3 moments. (f) Kurtosis calculated using SSA and approximation up to 4 and 6 moments.
Study of p53 model, parameter set q 1 = [90, 0.002, 1.7, 1.1, 0.93, 0.96, 0.01]. (a) Single SSA realisation. (b) Average of 100 000 SSA simulations. Trajectories calculated using (c) moment approximation including only the mean (deterministic), (d) linear noise approximation, (e) mean and the variance, (f) up to three central moments, (g) five central moments, and (h) six central moments. (i) and (j) Cumulative difference between mean trajectory calculated with SSA and trajectories calculated using (i) 2 moments and (j) 6 moments.
Analysis of distribution for p53 model. (a) Variance calculated based on SSA runs. (b) Variance calculated with SSA, LNA, and moment approximation method. (c) Skewness calculated based on SSA runs (blue line) and skewness for normal distribution (cyan dashed line). (d) and (f) Histograms calculated based on SSA for points d, e, and f in Figure 4(c) , and probability density function of normal distribution calculated using mean and variance based on SSA (cyan line).
Contour plots for the p53 system of the distance d between SSA trajectories and trajectories calculated with the moment expansion method. (a) and (b) Contour for varying k 1 and k 3, calculated using expansion up to 2 and 6 moments, respectively. (c) and (d) Contour for varying k 5 and k 6, calculated using expansion up to 2 and 6 moments, respectively. Red dots indicate parameter values used for SSA simulations.
Assessing the sensitivity to parameter k 1 of the dimerisation system of Sec. III A . Initial conditions and parameters are the same as previously, except that we additionally consider a perturbed value of k 1 = 1.1 × 1.66 × 10−3. (a) Average of 100 000 SSA simulations for the original and perturbed k 1 values. (b) As (a), except the mean is estimated using the proposed moment expansion approximation with 6 moments. (c) The sensitivity coefficients ∂μ1(t, [k 1, k 2])/∂k 1 and ∂μ2(t, [k 1, k 2])/∂k 2 (where μ1 and μ2 represent the means of x 1 and x 2, respectively) estimated from the SSA and 6m output via finite different approximations. (d), (e), and (f) As (a), (b), and (c), but for the second central moment. (g), (h), and (i) As (a), (b), and (c), but for the third central moment.
Study of the deviation from the mean (x − μ). (a) Deterministic mean and single SSA trajectory for dimerisation system. (b) Deterministic mean and 9 points taken from different SSA trajectories for dimerisation system. (c) Deterministic mean and single SSA trajectory for p53 system. (d) Deterministic mean and 9 points taken from different SSA trajectories for p53 system.
Error between mean, second, and third central moment calculated with SSA and approximation methods for the dimerization system.
Number of potentially nonzero central moment terms to include in moment approximation for different numbers of variables (columns) and number of central moments (rows).
Article metrics loading...
Full text loading...