Maximum likelihood estimation of log-affine models using detailed-balanced reaction networks

Abstract

A fundamental question in the field of molecular computation is what computational tasks a biochemical system can carry out. In this work, we focus on the problem of finding the maximum likelihood estimate (MLE) for log-affine models. We revisit a construction due to Gopalkrishnan of a mass-action system with the MLE as its unique positive steady state, which is based on choosing a basis for the kernel of the design matrix of the model. We extend this construction to allow for any finite spanning set of the kernel, and explore how the choice of spanning set influences the dynamics of the resulting network, including the existence of boundary steady states, the deficiency of the network, and the rate of convergence. In particular, we prove that using a Markov basis as the spanning set guarantees global stability of the MLE steady state.

25 pages, 5 figures. Improved exposition. To appear in J. Math. Biol

References in corpus (21)

• SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
• Array Programming with NumPy
• On the toric algebra of graphical models
• A geometric approach to the Global Attractor Conjecture
• Generalized mass action systems: Complex balancing equilibria and sign vectors of the stoichiometric and kinetic-order subspaces
• Product-form Poisson-like distributions and complex balanced reaction systems
• Nonpolynomial vector fields under the Lotka-Volterra normal form
• New Approach to the Stability of Chemical Reaction Networks: Piecewise Linear in Rates Lyapunov Functions
• The Maximum Likelihood Degree of Toric Varieties
• Catalysis in Reaction Networks
• Autocatalysis in Reaction Networks
• Stochastic Chemical Reaction Networks for Robustly Approximating Arbitrary Probability Distributions
• Multistationarity in the space of total concentrations for systems that admit a monomial parametrization
• Maximum Likelihood Estimation of Toric Fano Varieties
• Hierarchical subspace models for contingency tables
• Rate-Independent Computation in Continuous Chemical Reaction Networks
• Toric invariant theory for maximum likelihood estimation in log-linear models
• A reaction network scheme which implements inference and learning for Hidden Markov Models
• Families of polytopes with rational linear precision in higher dimensions
• The structure of the moduli space of toric dynamical systems of a reaction network
• Reaction extent or advancement of the reaction: A new general definition