Abstract
AbstractCellular processes are open systems, situated in a heterogeneous context, rather than operating in isolation. Chemical reaction networks (CRNs) whose reaction rates are modelled as external stochastic processes account for the heterogeneous environment when describing the embedded process. A marginal description of the embedded process is of interest for (i) fast simulations that bypass the co-simulation of the environment, (ii) obtaining new process equations from which moment equations can be derived, (iii) the computation of information-theoretic quantities, and (iv) state estimation. It is known since Snyder’s and related works that marginalization over a stochastic intensity turns point processes into self-exciting ones. While the Snyder filter specifies the exact history-dependent propensities in the framework of CRNs in Markov environment, it was recently suggested to use approximate filters for the marginal description. By regarding the chemical reactions as events, we establish a link between CRNs in a linear random environment and Hawkes processes, a class of self-exciting counting processes widely used in event analysis. The Hawkes approximation can be obtained via moment closure scheme or as the optimal linear approximation under the quadratic criterion. We show the equivalence of both approaches. Furthermore, we use martingale techniques to provide results on the agreement of the Hawkes process and the exact marginal process in their second order statistics, i.e., covariance, auto/cross-correlation. We introduce an approximate marginal simulation algorithm and illustrate it in case studies.AMS subject classifications37M05, 60G35, 60G55, 60J28, 60K37, 62M15
Publisher
Cold Spring Harbor Laboratory
Reference59 articles.
1. P. K. Andersen , O. Borgan , R. D. Gill , and N. Keiding , Statistical models based on counting processes, New York: Springer Verlag, 1993.
2. A modified next reaction method for simulating chemical systems with time dependent propensities and delays
3. Mutual information, relative entropy, and estimation in the Poisson channel;IEEE Transactions on Information Theory,2012
4. The spectral analysis of point processes;Journal of the Royal Statistical Society: Series B (Methodological),1963
5. J. S. Bendat and A. G. Piersol , Random data: analysis and measurement procedures, Wiley series in probability and statistics, Hoboken, N.J. : Wiley, 4th ed., 2010.