University of Cambridge > > Isaac Newton Institute Seminar Series > Efficiency of Stochastic Simulations

Efficiency of Stochastic Simulations

Add to your list(s) Download to your calendar using vCal

If you have a question about this talk, please contact INI IT.

SDBW01 - Opening workshop

Co-authors: David F. Anderson (University of Wisconsin-Madison), Yu Sun (University of Wisconsin-Madison)

I will analyze and compare the computational complexity of different simulation strategies for continuous time Markov chains. I consider the task of approximating the expected value of some functional of the state of the system over a compact time interval. This task is a bottleneck in many large-scale computations arising in biochemical kinetics and cell biology. In this context, the terms 'Gillespie's method', 'The Stochastic Simulation Algorithm' and 'The Next Reaction Method' are widely used to describe exact simulation methods. For example, Google Scholar records more than 6,000 citations to Gillespie's seminal 1977 paper. I will look at the use of standard Monte Carlo when samples are produced by exact simulation and by approximation with tau-leaping or an Euler-Maruyama discretization of a diffusion approximation. In particular, I will point out some possible pitfalls when computational complexity is analysed. Appropriate modifications o f recently proposed multilevel Monte Carlo algorithms will then be studied for the tau-leaping and Euler-Maruyama approaches. I will pay particular attention to a parameterization of the problem that, in the mass action chemical kinetics setting, corresponds to the classical system size scaling.

Related Links

This talk is part of the Isaac Newton Institute Seminar Series series.

Tell a friend about this talk:

This talk is included in these lists:

Note that ex-directory lists are not shown.


© 2006-2023, University of Cambridge. Contact Us | Help and Documentation | Privacy and Publicity