Home > Research > Publications & Outputs > Computational Methods for Complex Stochastic Sy...

Electronic data

Links

Text available via DOI:

View graph of relations

Computational Methods for Complex Stochastic Systems: A Review of Some Alternatives to MCMC.

Research output: Contribution to Journal/MagazineJournal articlepeer-review

Published
<mark>Journal publication date</mark>06/2008
<mark>Journal</mark>Statistics and Computing
Issue number2
Volume18
Number of pages21
Pages (from-to)151-171
Publication StatusPublished
<mark>Original language</mark>English

Abstract

We consider analysis of complex stochastic models based upon partial information. MCMC and reversible jump MCMC are often the methods of choice for such problems, but in some situations they can be difficult to implement; and suffer from problems such as poor mixing, and the difficulty of diagnosing convergence. Here we review three alternatives to MCMC methods: importance sampling, the forward-backward algorithm, and sequential Monte Carlo (SMC). We discuss how to design good proposal densities for importance sampling, show some of the range of models for which the forward-backward algorithm can be applied, and show how resampling ideas from SMC can be used to improve the efficiency of the other two methods. We demonstrate these methods on a range of examples, including estimating the transition density of a diffusion and of a discrete-state continuous-time Markov chain; inferring structure in population genetics; and segmenting genetic divergence data.