Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2601.12712 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866917209646825472 |
|---|---|
| author | Kimpson, Tom Germano, Domenic P. J. Flegg, Jennifer A. Flegg, Mark B. |
| author_facet | Kimpson, Tom Germano, Domenic P. J. Flegg, Jennifer A. Flegg, Mark B. |
| contents | Many biological systems exhibit multiscale dynamics, where some species occur in high copy numbers while others remain rare. This heterogeneity necessitates hybrid modelling approaches: deterministic models are computationally efficient but inaccurate for low-count species, while fully stochastic simulations are accurate but prohibitively expensive. Hybrid methods like the Jump-Switch-Flow (JSF) algorithm address this by simulating low-count species stochastically and high-count species deterministically. However, selecting regime-switching thresholds to control errors for specific observables remains an open challenge. We develop a principled framework for threshold selection targeting extinction probability. We formalise JSF as a piecewise-deterministic Markov process and derive backward equations for extinction under exact and hybrid dynamics. Near extinction boundaries, complex nonlinear dynamics reduce to tractable time-inhomogeneous linear birth-death processes. This structure yields a rigorous error decomposition based on early and late excursions. Isolating the dominant error term motivates a fast, actionable heuristic. We demonstrate via Monte Carlo studies on a stochastic Lotka-Volterra model that our heuristic reliably upper-bounds empirical errors in extinction probability. This enables users to select the smallest threshold that satisfies a target error tolerance. This work paves the way for principled, efficient multiscale modelling and simulation in stochastic biological systems. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2601_12712 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Multiscale Modelling of Birth-Death Processes Kimpson, Tom Germano, Domenic P. J. Flegg, Jennifer A. Flegg, Mark B. Populations and Evolution Many biological systems exhibit multiscale dynamics, where some species occur in high copy numbers while others remain rare. This heterogeneity necessitates hybrid modelling approaches: deterministic models are computationally efficient but inaccurate for low-count species, while fully stochastic simulations are accurate but prohibitively expensive. Hybrid methods like the Jump-Switch-Flow (JSF) algorithm address this by simulating low-count species stochastically and high-count species deterministically. However, selecting regime-switching thresholds to control errors for specific observables remains an open challenge. We develop a principled framework for threshold selection targeting extinction probability. We formalise JSF as a piecewise-deterministic Markov process and derive backward equations for extinction under exact and hybrid dynamics. Near extinction boundaries, complex nonlinear dynamics reduce to tractable time-inhomogeneous linear birth-death processes. This structure yields a rigorous error decomposition based on early and late excursions. Isolating the dominant error term motivates a fast, actionable heuristic. We demonstrate via Monte Carlo studies on a stochastic Lotka-Volterra model that our heuristic reliably upper-bounds empirical errors in extinction probability. This enables users to select the smallest threshold that satisfies a target error tolerance. This work paves the way for principled, efficient multiscale modelling and simulation in stochastic biological systems. |
| title | Multiscale Modelling of Birth-Death Processes |
| topic | Populations and Evolution |
| url | https://arxiv.org/abs/2601.12712 |