Abstract:Stochastic reaction network models arise in intracellular chemical reactions, epidemiological models and other population process models, and are a class of continuous time Markov chains which have the nonnegative integer lattice as state space. We consider the problem of estimating the conditional probability distribution of a stochastic reaction network given exact partial state observations in time snapshots. We propose a particle filtering method called the targeting method. Our approach takes into account that the reaction counts in between two observation snapshots satisfy linear constraints and also uses inhomogeneous Poisson processes as proposals for the reaction counts to facilitate exact interpolation. We provide rigorous analysis as well as numerical examples to illustrate our method and compare it with other alternatives.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: given discrete - time snapshot observations of some states, estimate the conditional probability distribution of a stochastic reaction network. Specifically, the author focuses on how to estimate the conditional probability distribution of a stochastic reaction network when the states of some species are accurately observed. ### Problem Background Stochastic reaction network models are widely used in intracellular chemical reactions, epidemiological models, and other population process models. These models are a class of continuous - time Markov chains (CTMC), and their state space is the non - negative integer lattice \( \mathbb{Z}^n_+ \). Each reaction event will change the state of the system by a fixed amount \( \nu_j\in\mathbb{Z}^n \). Due to the wide application of these models in many fields, the problem of inferring system states and parameters from partial observation data has become crucial. ### Specific Problems In this paper, the author pays special attention to the following problems: - **Partial State Observation**: Assume that the state of the system \( Z(t)=(X(t), Y(t)) \), where \( Y(t) \) is the partial state accurately observed in discrete - time snapshots \( t_1 < t_2 <\cdots< t_{T_s} \). - **Conditional Probability Distribution Estimation**: The goal is to estimate the conditional probability distribution \( \pi(t, z) = P\{Z(t) = z|Y(t_l)=y_l, l = 1,\ldots, T_s\} \) of the system state \( Z(t) \) given the partial state observations \( Y(t_l)=y_l \). ### Solutions To meet this challenge, the author proposes a particle filtering algorithm called the "targeting method". The main features of this method include: 1. **Linear Constraints**: Considering that the number of reactions between two observation snapshots satisfies linear constraints. 2. **Inhomogeneous Poisson Process**: Use the inhomogeneous Poisson process as the proposed distribution of the number of reactions for accurate interpolation. ### Method Advantages Compared with traditional particle filtering methods, the targeting method improves the accuracy and efficiency of estimation by forcing the filtering process to satisfy the next observed value. In addition, the author also provides strict mathematical analysis and numerical experiments to prove the superior performance of this method. ### Summary The core problem of this paper is to solve the problem of estimating the conditional probability distribution of a stochastic reaction network under partial - state discrete - time snapshot observations, and proposes the targeting method based on linear constraints and inhomogeneous Poisson processes to achieve this goal.

Stochastic Filtering of Reaction Networks Partially Observed in Time Snapshots

State and parameter estimation from exact partial state observation in stochastic reaction networks

Stochastic modeling of non-linear adsorption with Gaussian kernel density estimators

Stochastic filters based on hybrid approximations of multiscale stochastic reaction networks

Stochastic filtering for multiscale stochastic reaction networks based on hybrid approximations

Fast reactions with non-interacting species in stochastic reaction networks

Conditional Monte Carlo for Reaction Networks

Filtered finite state projection method for the analysis and estimation of stochastic biochemical reaction networks

Stochastic reaction networks with input processes: Analysis and applications to reporter gene systems

Convergence of Regularized Particle Filters for Stochastic Reaction Networks

Bayesian inference of Stochastic reaction networks using Multifidelity Sequential Tempered Markov Chain Monte Carlo

Bayesian Inference for Jump-Diffusion Approximations of Biochemical Reaction Networks

Block particle filters for state estimation of stochastic reaction-diffusion systems

Synthetic likelihood method for reaction network inference

Neural-network solutions to stochastic reaction networks

Inference for reaction networks using the Linear Noise Approximation

Optimal and $H_\infty$ Control of Stochastic Reaction Networks

Stochastic reaction networks within interacting compartments

Intermittent State Observer Design for Neural Networks with Reaction–Diffusion Terms Using Partial Measurements

Reduction for stochastic biochemical reaction networks with multiscale conservations

A Scalable Computational Framework for Establishing Long-Term Behavior of Stochastic Reaction Networks