Moment-closure approximations for mass-action models

Although stochastic population models have proved to be a powerful tool in the study of process generating mechanisms across a wide range of disciplines, all too often the associated mathematical development involves nonlinear mathematics, which immediately raises difficult and challenging analytic problems that need to be solved if useful progress is to be made. One approximation that is often employed to estimate the moments of a stochastic process is moment closure. This approximation essentially truncates the moment equations of the stochastic process. A general expression for the marginal- and joint-moment equations for a large class of stochastic population models is presented. The generalisation of the moment equations allows this approximation to be applied easily to a wide range of models. Software is available from http://pysbml.googlecode.com/ to implement the techniques presented here.