Convergence of regularized particle filters for stochastic reaction networks

TL;DR

This paper investigates the convergence of regularized particle filters for stochastic reaction networks, establishing conditions under which they reliably estimate hidden states despite artificial noise bias.

Contribution

The paper proves convergence conditions for RPFs applied to SRNs, using various dynamical models and sensitivity analysis, ensuring reliable filtering in biological systems.

Findings

RPFs converge under linear propensity growth conditions

Multiple dynamical models support the convergence results

Numerical examples validate theoretical findings

Abstract

Filtering for stochastic reaction networks (SRNs) is an important problem in systems/synthetic biology aiming to estimate the state of unobserved chemical species. A good solution to it can provide scientists valuable information about the hidden dynamic state and enable optimal feedback control. Usually, the model parameters need to be inferred simultaneously with state variables, and a conventional particle filter can fail to solve this problem accurately due to sample degeneracy. In this case, the regularized particle filter (RPF) is preferred to the conventional ones, as the RPF can mitigate sample degeneracy by perturbing particles with artificial noise. However, the artificial noise introduces an additional bias to the estimate, and, thus, it is questionable whether the RPF can provide reliable results for SRNs. In this paper, we aim to identify conditions under which the RPF converges to the exact filter in the filtering problem determined by a bimolecular network. First, we establish computationally efficient RPFs for SRNs on different scales using different dynamical models, including the continuous-time Markov process, tau-leaping model, and piecewise deterministic process. Then, by parameter sensitivity analyses, we show that the established RPFs converge to the exact filters if all reactions leading to an increase of the molecular population have linearly growing propensities and some other mild conditions are satisfied simultaneously. This ensures the performance of the RPF for a large class of SRNs, and several numerical examples are presented to illustrate our results.