Home > Research > Publications & Outputs > Direct simulation for discrete mixture distribu...
View graph of relations

Direct simulation for discrete mixture distributions.

Research output: Contribution to Journal/MagazineJournal articlepeer-review

Published
<mark>Journal publication date</mark>04/2005
<mark>Journal</mark>Statistics and Computing
Issue number2
Volume15
Number of pages9
Pages (from-to)125-133
Publication StatusPublished
<mark>Original language</mark>English

Abstract

We demonstrate how to perform direct simulation for discrete mixture models. The approach is based on directly calculating the posterior distribution using a set of recursions which are similar to those of the Forward-Backward algorithm. Our approach is more practicable than existing perfect simulation methods for mixtures. For example, we analyse 1096 observations from a 2 component Poisson mixture, and 240 observations under a 3 component Poisson mixture (with unknown mixture proportions and Poisson means in each case). Simulating samples of 10,000 perfect realisations took about 17 minutes and an hour respectively on a 900 MHz ultraSPARC computer. Our method can also be used to perform perfect simulation from Markov-dependent mixture models. A byproduct of our approach is that the evidence of our assumed models can be calculated, which enables different models to be compared.