TY - JOUR

T1 - An efficient forward–reverse expectation-maximization algorithm for statistical inference in stochastic reaction networks

AU - Bayer, Christian

AU - Moraes, Alvaro

AU - Tempone, Raul

AU - Vilanova, Pedro

N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The work described in this paper was supported by King Abdullah University of Science and Technology (KAUST). A. Moraes, R. Tempone and P. Vilanova are members of the KAUST SRI Center for Uncertainty Quantification at the Computer, Electrical and Mathematical Science and Engineering Division, KAUST.

PY - 2016/2/18

Y1 - 2016/2/18

N2 - © 2016 Taylor & Francis Group, LLC. ABSTRACT: In this work, we present an extension of the forward–reverse representation introduced by Bayer and Schoenmakers (Annals of Applied Probability, 24(5):1994–2032, 2014) to the context of stochastic reaction networks (SRNs). We apply this stochastic representation to the computation of efficient approximations of expected values of functionals of SRN bridges, that is, SRNs conditional on their values in the extremes of given time intervals. We then employ this SRN bridge-generation technique to the statistical inference problem of approximating reaction propensities based on discretely observed data. To this end, we introduce a two-phase iterative inference method in which, during phase I, we solve a set of deterministic optimization problems where the SRNs are replaced by their reaction-rate ordinary differential equations approximation; then, during phase II, we apply the Monte Carlo version of the expectation-maximization algorithm to the phase I output. By selecting a set of overdispersed seeds as initial points in phase I, the output of parallel runs from our two-phase method is a cluster of approximate maximum likelihood estimates. Our results are supported by numerical examples.

AB - © 2016 Taylor & Francis Group, LLC. ABSTRACT: In this work, we present an extension of the forward–reverse representation introduced by Bayer and Schoenmakers (Annals of Applied Probability, 24(5):1994–2032, 2014) to the context of stochastic reaction networks (SRNs). We apply this stochastic representation to the computation of efficient approximations of expected values of functionals of SRN bridges, that is, SRNs conditional on their values in the extremes of given time intervals. We then employ this SRN bridge-generation technique to the statistical inference problem of approximating reaction propensities based on discretely observed data. To this end, we introduce a two-phase iterative inference method in which, during phase I, we solve a set of deterministic optimization problems where the SRNs are replaced by their reaction-rate ordinary differential equations approximation; then, during phase II, we apply the Monte Carlo version of the expectation-maximization algorithm to the phase I output. By selecting a set of overdispersed seeds as initial points in phase I, the output of parallel runs from our two-phase method is a cluster of approximate maximum likelihood estimates. Our results are supported by numerical examples.

UR - http://hdl.handle.net/10754/621486

UR - http://www.tandfonline.com/doi/full/10.1080/07362994.2015.1116396

UR - http://www.scopus.com/inward/record.url?scp=84959118707&partnerID=8YFLogxK

U2 - 10.1080/07362994.2015.1116396

DO - 10.1080/07362994.2015.1116396

M3 - Article

VL - 34

SP - 193

EP - 231

JO - Stochastic Analysis and Applications

JF - Stochastic Analysis and Applications

SN - 0736-2994

IS - 2

ER -