A marginal sampler for -Stable Poisson-Kingman mixture models
arXiv:1407.4211 · doi:10.1080/10618600.2015.1110526
Abstract
We investigate the class of -stable Poisson-Kingman random probability measures (RPMs) in the context of Bayesian nonparametric mixture modeling. This is a large class of discrete RPMs which encompasses most of the the popular discrete RPMs used in Bayesian nonparametrics, such as the Dirichlet process, Pitman-Yor process, the normalized inverse Gaussian process and the normalized generalized Gamma process. We show how certain sampling properties and marginal characterizations of -stable Poisson-Kingman RPMs can be usefully exploited for devising a Markov chain Monte Carlo (MCMC) algorithm for making inference in Bayesian nonparametric mixture modeling. Specifically, we introduce a novel and efficient MCMC sampling scheme in an augmented space that has a fixed number of auxiliary variables per iteration. We apply our sampling scheme for a density estimation and clustering tasks with unidimensional and multidimensional datasets, and we compare it against competing sampling schemes.
New algorithmic performance comparisons were added
References in corpus (6)
- Are Gibbs-type priors the most natural generalization of the Dirichlet process?
- MCMC for Normalized Random Measure Mixture Models
- Poisson Process Partition Calculus with applications to Exchangeable models and Bayesian Nonparametrics
- A simple example of Dirichlet process mixture inconsistency for the number of components
- Modeling with Normalized Random Measure Mixture Models
- Gibbs Partitions (EPPF's) Derived From a Stable Subordinator are Fox H and Meijer G Transforms
Cited by in corpus (4)
- Completely random measures for modelling block-structured networks
- A hybrid sampler for Poisson-Kingman mixture models
- A priori truncation method for posterior sampling from homogeneous normalized completely random measure mixture models
- Gibbs Partitions, Riemann-Liouville Fractional Operators, Mittag-Leffler Functions, and Fragmentations Derived From Stable Subordinators