Reversible jump algorithm for analysis of gamma mixtures

Author(s):  
Petar M. Djuric
2019 ◽  
Vol 161 ◽  
pp. 32-51 ◽  
Author(s):  
Philippe Gagnon ◽  
Mylène Bédard ◽  
Alain Desgagné

2009 ◽  
Vol 178 (3) ◽  
pp. 1411-1436 ◽  
Author(s):  
Thomas Bodin ◽  
Malcolm Sambridge

2014 ◽  
Vol 72 ◽  
pp. 298-314 ◽  
Author(s):  
Silvia Pandolfi ◽  
Francesco Bartolucci ◽  
Nial Friel

2020 ◽  
pp. 1-33
Author(s):  
Andriy Norets

This article develops a Markov chain Monte Carlo (MCMC) method for a class of models that encompasses finite and countable mixtures of densities and mixtures of experts with a variable number of mixture components. The method is shown to maximize the expected probability of acceptance for cross-dimensional moves and to minimize the asymptotic variance of sample average estimators under certain restrictions. The method can be represented as a retrospective sampling algorithm with an optimal choice of auxiliary priors and as a reversible jump algorithm with optimal proposal distributions. The method is primarily motivated by and applied to a Bayesian nonparametric model for conditional densities based on mixtures of a variable number of experts. The mixture of experts model outperforms standard parametric and nonparametric alternatives in out of sample performance comparisons in an application to Engel curve estimation. The proposed MCMC algorithm makes estimation of this model practical.


2020 ◽  
Author(s):  
Felipe J Medina-Aguayo ◽  
Xavier Didelot ◽  
Richard G Everitt

AbstractBacteria reproduce clonally but most species recombine frequently, so that the ancestral process is best captured using an ancestral recombination graph. This graph model is often too complex to be used in an inferential setup, but it can be approximated for example by the ClonalOrigin model. Inference in the ClonalOrigin model is performed via a Reversible-Jump Markov Chain Monte Carlo algorithm, which attempts to jointly explore: the recombination rate, the number of recombination events, the departure and arrival points on the clonal genealogy for each recombination event, and the range of genomic sites affected by each recombination event. However, the Reversible-Jump algorithm usually performs poorly due to the complexity of the target distribution since it needs to explore spaces of different dimensions. Recent developments in Bayesian computation methodology have provided ways to improve existing methods and code, but are not well-known outside the statistics community. We show how exploiting one of these new computational methods can lead to faster inference under the ClonalOrigin model.


Biometrika ◽  
2011 ◽  
Vol 98 (1) ◽  
pp. 231-236 ◽  
Author(s):  
M. Papathomas ◽  
P. Dellaportas ◽  
V. G. S. Vasdekis

Sign in / Sign up

Export Citation Format

Share Document