pith. sign in

arxiv: 2504.15166 · v5 · submitted 2025-04-21 · 🧬 q-bio.QM · math.PR· physics.chem-ph· q-bio.MN

Simulating stochastic population dynamics: The Linear Noise Approximation can capture non-linear phenomena

Frederick Truman-Williams , Giorgos Minas This is my paper

Pith reviewed 2026-05-22 18:52 UTC · model grok-4.3

classification 🧬 q-bio.QM math.PRphysics.chem-phq-bio.MN
keywords Linear Noise Approximationstochastic population dynamicscentre manifold theorynon-linear dynamicsoscillatory systemsbistabilitygene regulationcomputational efficiency
0
0 comments X

The pith

The Linear Noise Approximation can be modified using centre manifold theory to accurately simulate long-term non-linear stochastic population dynamics.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper shows how to adapt the Linear Noise Approximation, a fast method for modeling random changes in population sizes, so it works for systems with non-linear features like cycles and multiple steady states. The authors use centre manifold theory to spot the slow parts of the dynamics and create simple, system-specific fixes for whole classes of similar behaviors. These fixes produce reliable long-term forecasts while keeping the speed and tractability that make the original approximation useful. Readers working on gene networks, epidemics, or ecological populations would care because the result gives a practical middle path between quick but limited models and slow but detailed ones.

Core claim

The paper claims that specific modifications to the Linear Noise Approximation, identified using centre manifold theory from non-linear dynamical systems, allow it to capture non-linear phenomena such as oscillations and bistability in stochastic population processes, resulting in accurate long-term simulations that preserve the method's computational efficiency and analytical advantages.

What carries the argument

The modified Linear Noise Approximation using centre manifold theory, which identifies simple system-specific corrections for classes of oscillatory and bi-stable systems.

If this is right

  • The modified LNA produces accurate long-term trajectories for classes of oscillatory systems.
  • The same approach works for bi-stable systems that switch between multiple stable states.
  • Computational cost stays low enough for repeated simulations, sensitivity analysis, and parameter estimation.
  • Concrete examples from molecular population dynamics confirm the gains in accuracy and speed.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The method could be tried on larger gene regulatory networks to test whether the speed advantage scales.
  • Similar centre-manifold corrections might apply to stochastic models in ecology or epidemiology.
  • Direct runtime and error comparisons against full stochastic simulators would quantify the practical benefit.

Load-bearing premise

That centre manifold reductions developed for deterministic systems transfer directly to the stochastic LNA setting and remain accurate for entire classes of oscillatory and bi-stable systems.

What would settle it

Run both the modified LNA and an exact stochastic simulation method on a known oscillatory gene circuit for many thousands of time units and check whether the long-term amplitude and period match; large divergence would show the claim is false.

Figures

Figures reproduced from arXiv: 2504.15166 by Frederick Truman-Williams, Giorgos Minas.

Figure 1
Figure 1. Figure 1: A. Phase portraits of solutions of the RRE in (6) (black) and a SSA simulation (red) [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: LNA method summary. A. LNA is a stochastic approximation of the stochastic [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: General pcLNA algorithm over steps i − 1, i, i + 1. For simplicity we illustrate the case where j − 1 = j = j + 1, so the same RRE solution is used over these steps. pcLNA proceeds with standard LNA steps interrupted by phase corrections with x (ji) (ti) replaced by G(X(ti)) = x (ji) (si), and the perturbations ξ(ti) replaced by κ(si). All that is preventing one from readily applying the steps in the pcLNA… view at source ↗
Figure 4
Figure 4. Figure 4: The dynamics of a generic normal form of the supercritical Hopf bifurcation. The two [PITH_FULL_IMAGE:figures/full_fig_p014_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Phase correction. The phase-correcting map [PITH_FULL_IMAGE:figures/full_fig_p016_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Phase correction (PC) for Hopf bifurcation systems. (Left) The stochastic trajectory [PITH_FULL_IMAGE:figures/full_fig_p017_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Phase Correction (PC) for bistable systems. (Left) The stochastic trajectory (red) [PITH_FULL_IMAGE:figures/full_fig_p019_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Comparison between SSA (red) and pcLNA (blue) for (A) the NF- [PITH_FULL_IMAGE:figures/full_fig_p022_8.png] view at source ↗
read the original abstract

Population dynamics in fields such as molecular biology, epidemiology, and ecology exhibit highly stochastic and non-linear behaviour. In gene regulatory systems in particular, oscillations and multi-stability are especially common. Despite this, none of the currently available stochastic models for population dynamics are both accurate and computationally efficient for long-term predictions. A prominent model in this field, the Linear Noise Approximation (LNA), is computationally efficient for tasks such as simulation, sensitivity analysis, and parameter estimation; however, it is only accurate for linear systems and short-time predictions. Other models may achieve greater accuracy across a broader range of systems, but they sacrifice computational efficiency and analytical tractability. This paper demonstrates that, with specific modifications, the LNA can accurately capture non-linear dynamics in population processes. We introduce a new framework based on centre manifold theory, a classical concept from non-linear dynamical systems. This approach enables the identification of simple, system-specific modifications to the LNA, tailored to classes of qualitatively similar non-linear dynamical systems. With these modifications, the LNA can achieve accurate long-term simulations without compromising computational efficiency. We apply our methodology to classes of oscillatory and bi-stable systems, and present multiple examples from molecular population dynamics that demonstrate accurate long-term simulations alongside significant improvements in computational efficiency.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The paper claims that specific modifications to the Linear Noise Approximation (LNA), derived from centre manifold theory applied to the deterministic drift, enable accurate long-term stochastic simulations of non-linear population dynamics (oscillatory and bi-stable systems) in molecular biology without sacrificing the computational efficiency of the standard LNA.

Significance. If the central construction holds, the work would offer a practical advance for long-horizon stochastic modeling in systems biology and epidemiology by retaining the analytical tractability and speed of the LNA while extending its validity to qualitatively non-linear regimes; the explicit use of centre manifold reduction and system-specific tailoring for entire classes of models is a notable strength.

major comments (2)
  1. [Framework section (centre manifold construction)] The manuscript does not supply a stochastic error bound or Lyapunov function controlling distance to the identified centre manifold under the full chemical master equation when noise intensity is O(1/√Ω). This is load-bearing for the long-term accuracy claim in bi-stable and oscillatory examples, because noise-induced excursions can occur on timescales comparable to the simulation horizon.
  2. [Results / examples section] Quantitative comparison tables, error metrics, and full derivations comparing the modified LNA against exact methods (e.g., Gillespie) over long times are absent from the abstract and would be required to substantiate the claim that post-hoc system-specific corrections preserve marginal statistics without hidden fitting.
minor comments (1)
  1. [Methods] Notation for the modified drift and diffusion terms should be introduced with explicit equations rather than descriptive text to improve reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments. We address each major comment below and outline the revisions we intend to incorporate.

read point-by-point responses

  1. Referee: [Framework section (centre manifold construction)] The manuscript does not supply a stochastic error bound or Lyapunov function controlling distance to the identified centre manifold under the full chemical master equation when noise intensity is O(1/√Ω). This is load-bearing for the long-term accuracy claim in bi-stable and oscillatory examples, because noise-induced excursions can occur on timescales comparable to the simulation horizon.

    Authors: We appreciate the referee pointing out the absence of a rigorous stochastic error bound. The present framework first reduces the deterministic drift via centre manifold theory and then applies the LNA to the reduced system; the stochastic justification rests on the standard LNA scaling (fluctuations of order 1/√Ω) together with numerical evidence that trajectories remain close to the manifold over long horizons in the chosen examples. A full Lyapunov analysis for the chemical master equation lies outside the scope of the current work. In the revision we will add an explicit discussion of this limitation, including references to existing results on stochastic centre manifolds, and we will qualify the long-term claims accordingly. revision: partial

  2. Referee: [Results / examples section] Quantitative comparison tables, error metrics, and full derivations comparing the modified LNA against exact methods (e.g., Gillespie) over long times are absent from the abstract and would be required to substantiate the claim that post-hoc system-specific corrections preserve marginal statistics without hidden fitting.

    Authors: The abstract is a concise summary and therefore does not contain tables or numerical metrics. The body of the manuscript already presents direct comparisons with Gillespie trajectories for both oscillatory and bistable systems over extended times. To strengthen the presentation we will insert additional quantitative error measures (time-averaged L2 errors on means and variances, and distributional distances) into the results section and will revise the abstract to reference these comparisons. The corrections themselves are obtained deterministically from the centre-manifold reduction and are not the result of statistical fitting to stochastic trajectories. revision: yes

Circularity Check

0 steps flagged

No circularity: centre manifold reduction is imported from classical deterministic theory

full rationale

The paper's core step applies standard centre manifold theory (a classical result from deterministic nonlinear dynamics) to derive system-specific corrections to the LNA drift/diffusion terms. This is an external mathematical import rather than a self-definition, fitted parameter, or self-citation chain. The abstract and description give no equations that equate a prediction to its own input by construction, nor any load-bearing uniqueness theorem authored by the same team. The derivation therefore remains self-contained against external benchmarks and does not reduce to tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the applicability of deterministic centre manifold theory to the stochastic LNA and on the existence of simple modifications that work uniformly for classes of oscillatory and bi-stable systems. No free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption Centre manifold theory from non-linear dynamical systems can be used to identify simple modifications to the LNA that capture non-linear phenomena in stochastic population processes.
    Invoked in the abstract as the basis for the new framework enabling accurate long-term simulations.

pith-pipeline@v0.9.0 · 5766 in / 1301 out tokens · 34044 ms · 2026-05-22T18:52:06.422540+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

79 extracted references · 79 canonical work pages

  1. [1]

    Anderson and T.G

    D.F. Anderson and T.G. Kurtz.Continuous Time Markov Chain Models for Chemical Reaction Networks, pages 3–42. Springer New York, 2011

    work page 2011

  2. [2]

    Ferrell, and Eduardo D

    David Angeli, James E. Ferrell, and Eduardo D. Sontag. Detection of multistability, bifur- cations, and hysteresis in a large class of biological positive-feedback systems.Proceedings of the National Academy of Sciences, 101:1822–1827, 2 2004

    work page 2004

  3. [3]

    Fernández Arancibia, Hernán E

    Sol M. Fernández Arancibia, Hernán E. Grecco, and Luis G. Morelli. Effective description of bistability and irreversibility in apoptosis.Phys. Rev. E, 104:064410, Dec 2021

    work page 2021

  4. [4]

    Ashall, C.A

    L. Ashall, C.A. Horton, D.E. Nelson, P. Paszek, C.V. Harper, K. Sillitoe, S. Ryan, D.G. Spiller, J.F. Unitt, D.S. Broomhead, D.B. Kell, D.A. Rand, V. Sée, and M.R.H. White. Pulsatile stimulation determines timing and specificity of nf-κb-dependent transcription. Science (American Association for the Advancement of Science), 324(5924):242–246, 2009

    work page 2009

  5. [5]

    Boland, T

    R.P. Boland, T. Galla, and A.J. McKane. How limit cycles and quasi-cycles are related in systems with intrinsic noise.Journal of Statistical Mechanics: Theory and Experiment, 2008, 9 2008

    work page 2008

  6. [6]

    Alexander Cao, Benjamin Lindner, and Peter J. Thomas. A partial differential equation for the mean–return-time phase of planar stochastic oscillators.SIAM journal on applied mathematics, 80(1):422–447, 2020

    work page 2020

  7. [7]

    Linear mapping approximation of gene regulatory networks with stochastic dynamics.Nature Communications, 9:3305, 8 2018

    Zhixing Cao and Ramon Grima. Linear mapping approximation of gene regulatory networks with stochastic dynamics.Nature Communications, 9:3305, 8 2018

    work page 2018

  8. [8]

    Rationalised experiment design for parameter estimation with sensitivity clustering.Scientific Reports, 14:25864, 10 2024

    Harsh Chhajer and Rahul Roy. Rationalised experiment design for parameter estimation with sensitivity clustering.Scientific Reports, 14:25864, 10 2024

    work page 2024

  9. [9]

    Geometry, epistasis, and developmental patterning

    Francis Corson and Eric Dean Siggia. Geometry, epistasis, and developmental patterning. Proceedings of the National Academy of Sciences, 109:5568–5575, 4 2012

    work page 2012

  10. [10]

    Brown, Boris V

    Mirela Domijan, Paul E. Brown, Boris V. Shulgin, and David A. Rand. Pettsy: a compu- tational tool for perturbation analysis of complex systems biology models.BMC Bioinfor- matics, 17:124, 3 2016

    work page 2016

  11. [11]

    A random dynamical systems perspective on isochronicity for stochastic oscillations.Communications in Mathematical Physics, 386(3):1603–1641, 2021

    Maximilian Engel and Christian Kuehn. A random dynamical systems perspective on isochronicity for stochastic oscillations.Communications in Mathematical Physics, 386(3):1603–1641, 2021

    work page 2021

  12. [12]

    Woodcock, Michal Komorowski, Claire V

    B¨ arbel Finkenst¨ adt, Dan J. Woodcock, Michal Komorowski, Claire V. Harper, Julian R. E. Davis, Mike R. H. White, David a. Rand, Barbel Finkenstadt, Dan J. Woodcock, Michal Komorowski, Claire V. Harper, Julian R. E. Davis, Mike R. H. White, and David a. Rand. Quantifying intrinsic and extrinsic noise in gene transcription using the linear noise approx- ...

    work page 1960

  13. [13]

    Prediction of predator–prey populations modelled by perturbed odes.Journal of Mathematical Biology, 54(3):407–451, 2007

    Sorana Froda and Sévérien Nkurunziza. Prediction of predator–prey populations modelled by perturbed odes.Journal of Mathematical Biology, 54(3):407–451, 2007

    work page 2007

  14. [14]

    Gardner, Charles R

    Timothy S. Gardner, Charles R. Cantor, and James J. Collins. Construction of a genetic toggle switch in escherichia coli.Nature, 403:339–342, 1 2000. 24

    work page 2000

  15. [15]

    Gillespie

    D.T. Gillespie. A general method for numerically simulating the stochastic time evolution of coupled chemical reactions.Journal of Computational Physics, 22(4):403–434, 1976

    work page 1976

  16. [16]

    Gillespie

    D.T. Gillespie. Exact stochastic simulation of coupled chemical reactions.The Journal of Physical Chemistry, 81(25):2340–2361, 1977

    work page 1977

  17. [17]

    Gillespie

    D.T. Gillespie. A rigorous derivation of the chemical master equation.Physica A, 188(1):404–425, 1992

    work page 1992

  18. [18]

    Gillespie

    D.T. Gillespie. The chemical Langevin equation.The Journal of Chemical Physics, 113(1):297–306, 07 2000

    work page 2000

  19. [19]

    Gillespie

    D.T. Gillespie. Approximate accelerated stochastic simulation of chemically reacting sys- tems.The Journal of Chemical Physics, 115(4):1716–1733, 07 2001

    work page 2001

  20. [20]

    Gillespie and L.R

    D.T. Gillespie and L.R. Petzold. Improved leap-size selection for accelerated stochastic simulation.The Journal of Chemical Physics, 119(16):8229–8234, 10 2003

    work page 2003

  21. [21]

    Sharp developmental thresholds defined through bistability by antagonistic gradients of retinoic acid and fgf signaling.De- velopmental Dynamics, 236:1495–1508, 6 2007

    Albert Goldbeter, Didier Gonze, and Olivier Pourquié. Sharp developmental thresholds defined through bistability by antagonistic gradients of retinoic acid and fgf signaling.De- velopmental Dynamics, 236:1495–1508, 6 2007

    work page 2007

  22. [22]

    Wadkin, Sam A

    Andrew Golightly, Laura E. Wadkin, Sam A. Whitaker, Andrew W. Baggaley, Nick G. Parker, and Theodore Kypraios. Accelerating bayesian inference for stochastic epidemic models using incidence data.Statistics and Computing, 33:134, 12 2023

    work page 2023

  23. [23]

    Stochastic mod- els for circadian rhythms: effect of molecular noise on periodic and chaotic behaviour

    Didier Gonze, José Halloy, Jean-Christophe Leloup, and Albert Goldbeter. Stochastic mod- els for circadian rhythms: effect of molecular noise on periodic and chaotic behaviour. Comptes Rendus Biologies, 326(2):189 – 203, 2003

    work page 2003

  24. [24]

    Grunberg and Domitilla Del Vecchio

    Theodore W. Grunberg and Domitilla Del Vecchio. A stein’s method approach to the linear noise approximation for stationary distributions of chemical reaction networks, 2024

    work page 2024

  25. [25]

    An efficient and unbiased method for sensitivity anal- ysis of stochastic reaction networks.Journal of the Royal Society Interface, 11:20140979– 20140979, 2014

    Ankit Gupta and Mustafa Khammash. An efficient and unbiased method for sensitivity anal- ysis of stochastic reaction networks.Journal of the Royal Society Interface, 11:20140979– 20140979, 2014

    work page 2014

  26. [26]

    A.V. Hill. The possible effects of the aggregation of the molecules of hæmoglobin on its dissociation curves.The Journal of Physiology, 40:i–vii, January 1910

    work page 1910

  27. [27]

    Matlab version: 9.13.0 (r2022b), 2022

    The MathWorks Inc. Matlab version: 9.13.0 (r2022b), 2022

    work page 2022

  28. [28]

    Ishii, K

    N. Ishii, K. Nakahigashi, T. Baba, M. Robert, T. Soga, A. Kanai, T. Hirasawa, M. Naba, K. Hirai, A. Hoque, P.Y. Ho, Y. Kakazu, K. Sugawara, S. Igarashi, S. Harada, T. Masuda, N. Sugiyama, T. Togashi, M. Hasegawa, Y. Takai, K. Yugi, K. Arakawa, N. Iwata, Y. Toya, Y. Nakayama, T. Nishioka, K. Shimizu, H. Mori, and M. Tomita. Multiple high-throughput analyse...

    work page 2007

  29. [29]

    Ito and K

    Y. Ito and K. Uchida. Formulas for intrinsic noise evaluation in oscillatory genetic networks. Journal of Theoretical Biology, 267:223–234, 11 2010

    work page 2010

  30. [30]

    Holimap: an accurate and efficient method for solving stochas- tic gene network dynamics.Nature Communications, 15:6557, 8 2024

    Chen Jia and Ramon Grima. Holimap: an accurate and efficient method for solving stochas- tic gene network dynamics.Nature Communications, 15:6557, 8 2024. 25

    work page 2024

  31. [31]

    Johnson and R.S

    K.A. Johnson and R.S. Goody. The original michaelis constant: Translation of the 1913 michaelis–menten paper.Biochemistry (Easton), 50(39):8264–8269, 2011

    work page 1913

  32. [32]

    Michal Komorowski, Maria J Costa, David a Rand, and Michael P H Stumpf. Sensitivity, robustness, and identifiability in stochastic chemical kinetics models.Proceedings of the National Academy of Sciences of the United States of America, 108:8645–8650, 5 2011

    work page 2011

  33. [33]

    Bayesian inference of biochemical kinetic parameters using the linear noise approximation.BMC Bioinformatics, 10:343, 10 2009

    Michal Komorowski, Barbel Finkenstadt, Claire V Harper, David a Rand, B¨ arbel Finkenst¨ adt, Claire V Harper, David a Rand, Barbel Finkenstadt, Claire V Harper, and David a Rand. Bayesian inference of biochemical kinetic parameters using the linear noise approximation.BMC Bioinformatics, 10:343, 10 2009

    work page 2009

  34. [34]

    T.G. Kurtz. Limit theorems for sequences of jump markov processes approximating ordinary differential processes.Journal of Applied Probability, 8(2):344–356, 1971

    work page 1971

  35. [35]

    Kurtz, Society for Industrial, and Applied Mathematics.Approximation of Popula- tion Processes

    T.G. Kurtz, Society for Industrial, and Applied Mathematics.Approximation of Popula- tion Processes. Number nos. 36-40 in Approximation of Population Processes. Society for Industrial and Applied Mathematics, 1981

    work page 1981

  36. [36]

    Thomas G. Kurtz. Strong approximation theorems for density dependent markov chains. Stochastic Processes and their Applications, 6(3):223–240, 1978

    work page 1978

  37. [37]

    Kuznetsov.Elements of Applied Bifurcation Theory, 3rd Edic

    Y. Kuznetsov.Elements of Applied Bifurcation Theory, 3rd Edic. Society for Industrial and Applied Mathematics, 2011

    work page 2011

  38. [38]

    J. F. Le Gall.Brownian motion, martingales, and stochastic calculus. Graduate texts in mathematics, 274. Springer, Switzerland, 2016

    work page 2016

  39. [39]

    Lefever and G

    R. Lefever and G. Nicolis. Chemical instabilities and sustained oscillations.Journal of Theoretical Biology, 30(2):267–284, 1971

    work page 1971

  40. [40]

    Ross, and Glenn Vinnicombe

    Ioannis Lestas, Johan Paulsson, Nicholas E. Ross, and Glenn Vinnicombe. Noise in gene regulatory networks.IEEE Transactions on Automatic Control, 53:189–200, 1 2008

    work page 2008

  41. [41]

    A framework for parameter estimation and model selection from experimental data in systems biology using approximate bayesian computation.Nature Protocols, 9:439–456, 2 2014

    Juliane Liepe, Paul Kirk, Sarah Filippi, Tina Toni, Chris P Barnes, and Michael P H Stumpf. A framework for parameter estimation and model selection from experimental data in systems biology using approximate bayesian computation.Nature Protocols, 9:439–456, 2 2014

    work page 2014

  42. [42]

    Mangan and U

    S. Mangan and U. Alon. Structure and function of the feed-forward loop network motif. Proceedings of the National Academy of Sciences, 100:11980–11985, 10 2003

    work page 2003

  43. [43]

    Adamson, and Nancy Papalopulu

    Elli Marinopoulou, Veronica Biga, Nitin Sabherwal, Anzy Miller, Jayni Desai, Antony D. Adamson, and Nancy Papalopulu. Hes1 protein oscillations are necessary for neural stem cells to exit from quiescence.iScience, 24:103198, 10 2021

    work page 2021

  44. [44]

    Markov chain monte carlo without likelihoods.Proceedings of the National Academy of Sciences, 100:15324– 15328, 12 2003

    Paul Marjoram, John Molitor, Vincent Plagnol, and Simon Tavaré. Markov chain monte carlo without likelihoods.Proceedings of the National Academy of Sciences, 100:15324– 15328, 12 2003

    work page 2003

  45. [45]

    Vestergaard.Gillespie Algorithms for Stochastic Multiagent Dynamics in Populations and Networks

    Naoki Masuda and Christian L. Vestergaard.Gillespie Algorithms for Stochastic Multiagent Dynamics in Populations and Networks. Cambridge University Press, 1 2023

    work page 2023

  46. [46]

    Minas and D.A

    G. Minas and D.A. Rand. Long-time analytic approximation of large stochastic oscillators: Simulation, analysis and inference.PLoS Computational Biology, 13, 7 2017. 26

    work page 2017

  47. [47]

    Minas and D.A

    G. Minas and D.A. Rand. Parameter sensitivity analysis for biochemical reaction networks. Mathematical biosciences and engineering : MBE, 16(5):3965–3987, 2019

    work page 2019

  48. [48]

    Minas, D.J

    G. Minas, D.J. Woodcock, L. Ashall, C. Harper, M.R.H. White, and D.A. Rand. Multi- plexing information flow through dynamic signalling systems.PLoS computational biology, 16(8):e1008076–e1008076, 2020

    work page 2020

  49. [49]

    Parameter sensitivity analysis for biochemical reaction networks.Mathematical Biosciences and Engineering, 16:3965–3987, 4 2019

    Giorgos Minas and David A Rand. Parameter sensitivity analysis for biochemical reaction networks.Mathematical Biosciences and Engineering, 16:3965–3987, 4 2019

    work page 2019

  50. [50]

    Bela Novak and John J. Tyson. Numerical analysis of a comprehensive model of m-phase control in xenopus oocyte extracts and intact embryos.Journal of Cell Science, 106:1153– 1168, 12 1993

    work page 1993

  51. [51]

    Palis and F

    J. Palis and F. Takens. Stability of parametrized families of gradient vector fields.Annals of Mathematics, 118(3):383–421, 1983

    work page 1983

  52. [52]

    Pomerening, Sun Young Kim, and James E

    Joseph R. Pomerening, Sun Young Kim, and James E. Ferrell. Systems-level dissection of the cell-cycle oscillator: Bypassing positive feedback produces damped oscillations.Cell, 122:565–578, 8 2005

    work page 2005

  53. [53]

    Ivanova, and Freddie Bickford Smith

    Tom Rainforth, Adam Foster, Desi R. Ivanova, and Freddie Bickford Smith. Modern bayesian experimental design.Statistical Science, 39, 2 2024

    work page 2024

  54. [54]

    Negative autoregulation speeds the response times of transcription networks.Journal of Molecular Biology, 323:785–793, 11 2002

    Nitzan Rosenfeld, Michael B Elowitz, and Uri Alon. Negative autoregulation speeds the response times of transcription networks.Journal of Molecular Biology, 323:785–793, 11 2002

    work page 2002

  55. [55]

    Approximation and inference meth- ods for stochastic biochemical kinetics - a tutorial review.Journal of Physics A: Mathemat- ical and Theoretical, 50, 8 2016

    David Schnoerr, Guido Sanguinetti, and Ramon Grima. Approximation and inference meth- ods for stochastic biochemical kinetics - a tutorial review.Journal of Physics A: Mathemat- ical and Theoretical, 50, 8 2016

    work page 2016

  56. [56]

    J. T. C Schwabedal and A Pikovsky. Effective phase description of noise-perturbed and noise-induced oscillations: The dynamics of nonlinear stochastic systems.The European physical journal. ST, Special topics, 187:63–76, 2010

    work page 2010

  57. [57]

    Scott, B

    M. Scott, B. Ingalls, and M. Kærn. Estimations of intrinsic and extrinsic noise in models of nonlinear genetic networks.Chaos, 16, 2006

    work page 2006

  58. [58]

    Ali Shahmohammadi and Kimberley B. McAuley. Using prior parameter knowledge in model-based design of experiments for pharmaceutical production.AIChE Journal, 66, 11 2020

    work page 2020

  59. [59]

    A. N. Shoshitaishvili. Bifurcations of topological type at singular points of parametrized vector fields.Functional analysis and its applications, 6(2):169–170, 1972

    work page 1972

  60. [60]

    S. A. Sisson, Y. Fan, and Mark M. Tanaka. Sequential monte carlo without likelihoods. Proceedings of the National Academy of Sciences, 104:1760–1765, 2 2007

    work page 2007

  61. [61]

    Swallow, D.A

    B. Swallow, D.A. Rand, and G. Minas. Bayesian inference for stochastic oscillatory systems using the phase-corrected linear noise approximation.Bayesian Analysis, pages 1 – 30, 2024

    work page 2024

  62. [62]

    Partially hyperbolic fixed points.Topology (Oxford), 10(2):133–147, 1971

    Floris Takens. Partially hyperbolic fixed points.Topology (Oxford), 10(2):133–147, 1971

    work page 1971

  63. [63]

    Inferring coalescence times from dna sequence data.Genetics, 145:505–518, 2 1997

    Simon Tavaré, David J Balding, R C Griffiths, and Peter Donnelly. Inferring coalescence times from dna sequence data.Genetics, 145:505–518, 2 1997. 27

    work page 1997

  64. [64]

    Thomas, H

    P. Thomas, H. Matuschek, and R. Grima. How reliable is the linear noise approximation of gene regulatory networks?BMC Genomics, 14, 2013

    work page 2013

  65. [65]

    Ap- proximate bayesian computation scheme for parameter inference and model selection in dynamical systems.Journal of The Royal Society Interface, 6:187–202, 2 2009

    Tina Toni, David Welch, Natalja Strelkowa, Andreas Ipsen, and Michael P.H Stumpf. Ap- proximate bayesian computation scheme for parameter inference and model selection in dynamical systems.Journal of The Royal Society Interface, 6:187–202, 2 2009

    work page 2009

  66. [66]

    van Kampen.Stochastic Processes in Physics and Chemistry

    N.G. van Kampen.Stochastic Processes in Physics and Chemistry. Elsevier Science Pub- lishers, Amsterdam, 1992

    work page 1992

  67. [67]

    Optimal experiment design for model selection in biochemical networks.BMC Systems Biology, 8:20, 2014

    Joep Vanlier, Christian A Tiemann, Peter AJ Hilbers, and Natal A W van Riel. Optimal experiment design for model selection in biochemical networks.BMC Systems Biology, 8:20, 2014

    work page 2014

  68. [68]

    Waage and C.M

    P. Waage and C.M. Gulberg. Studies concerning affinity.Journal of chemical education, 63(12):1044–, 1986

    work page 1986

  69. [69]

    Wallace, D.T

    E.W.J. Wallace, D.T. Gillespie, K.R. Sanft, and L.R. Petzold. Linear noise approximation is valid over limited times for any chemical system that is sufficiently large.IET Systems Biology, 6:102–115, 8 2012

    work page 2012

  70. [70]

    Weitz, Sang Woo Park, Ceyhun Eksin, and Jonathan Dushoff

    Joshua S. Weitz, Sang Woo Park, Ceyhun Eksin, and Jonathan Dushoff. Awareness-driven behavior changes can shift the shape of epidemics away from peaks and toward plateaus, shoulders, and oscillations.Proceedings of the National Academy of Sciences, 117:32764– 32771, 12 2020

    work page 2020

  71. [71]

    Wilhelm and R

    T. Wilhelm and R. Heinrich. Smallest Chemical-Reaction System With Hopf-Bifurcation. J Math Chem, 17:1–14, 1995

    work page 1995

  72. [72]

    Stochastic modelling for systems biology, 2012

    Darren James Wilkinson. Stochastic modelling for systems biology, 2012. 28 Supporting Information (SI) The SI contains further mathematical details and numerical investigations referenced in the main paper. We used PeTTSy [10] implemented in the MATLAB [27] environment and available athttps://wrap.warwick.ac.uk/id/eprint/77654/to produce all results. S1 T...

    work page 2012

  73. [76]

    Outputs:(t j,y j),j= 1,

    updatei←i+ 1and storetint i andyiny i. Outputs:(t j,y j),j= 1, . . . , i Algorithm 5SSA with thinning Inputs:Initial conditions:t 0 = 0,Y(0) =y 0,i= 0,m= 1. ParameterT >0,δt >0,M positive integer. Steps:Whilet < Tandm≤M

    work page

  74. [77]

    generate two random numbers,r 1 andr 2, from the uniform distribution on the unit interval

    work page

  75. [78]

    sample ˆτaccording to (34) andkaccording to (35)

    work page

  76. [79]

    updatet←t+ ˆτandy←y+ν k

    work page

  77. [80]

    increaseδtand run again

    updatei←i+ 1and ift∈(mδt,(m+ 1)δt), then storetint m andyiny m, and updatem←m+ 1, ift >(m+ 1)δt, end procedure with message “increaseδtand run again.”. Outputs:(t j,y j),j= 1, . . . , m When performing the SSA simulations, we used the so-calledthinningmethod to reduce the computational memory used. That is, we used the Algorithm 5. This records the SSA ge...

    work page

  78. [81]

    Initial values Fig. S10

    In particular, we provide the details of the systems used and additional comparisons between SSA and pcLNA. S4.1 The NF-κB reaction network The NF-κB network consists of11species detailed in Table S2. 36 Table S2: Details of the11species characterising the NF-κB network and their initial concen- trations in three simulation settings with different value o...

    work page

  79. [82]

    Panels (A)–(D) correspond to timest= 0.5000,9.9999,49.9996,100.0000, respectively. Each panel shows (left) the states recorded at timetof the simulations with solutions of the corresponding RRE, and (right) the corresponding empirical probability density functions of each system variable. The initial conditions for the two ODE solutions are provided in Ta...

    work page 2000