Gibbs Sampler for Matrix Generalized Inverse Gaussian Distributions

02/20/2023
by   Yasuyuki Hamura, et al.
0

Sampling from matrix generalized inverse Gaussian (MGIG) distributions is required in Markov Chain Monte Carlo (MCMC) algorithms for a variety of statistical models. However, an efficient sampling scheme for the MGIG distributions has not been fully developed. We here propose a novel blocked Gibbs sampler for the MGIG distributions, based on the Choleski decomposition. We show that the full conditionals of the diagonal and unit lower-triangular entries are univariate generalized inverse Gaussian and multivariate normal distributions, respectively. Several variants of the Metropolis-Hastings algorithm can also be considered for this problem, but we mathematically prove that the average acceptance rates become extremely low in particular scenarios. We demonstrate the computational efficiency of the proposed Gibbs sampler through simulation studies and data analysis.

READ FULL TEXT
research
09/10/2018

The Coordinate Sampler: A Non-Reversible Gibbs-like MCMC Sampler

In this article, we derive a novel non-reversible, continuous-time Marko...
research
05/06/2020

A Bayesian approach for clustering skewed data using mixtures of multivariate normal-inverse Gaussian distributions

Non-Gaussian mixture models are gaining increasing attention for mixture...
research
07/04/2018

Discrete Sampling using Semigradient-based Product Mixtures

We consider the problem of inference in discrete probabilistic models, t...
research
09/19/2010

Robust graphical modeling of gene networks using classical and alternative T-distributions

Graphical Gaussian models have proven to be useful tools for exploring n...
research
08/27/2021

Bayesian Sparse Blind Deconvolution Using MCMC Methods Based on Normal-Inverse-Gamma Prior

Bayesian estimation methods for sparse blind deconvolution problems conv...
research
06/02/2011

Restricted Collapsed Draw: Accurate Sampling for Hierarchical Chinese Restaurant Process Hidden Markov Models

We propose a restricted collapsed draw (RCD) sampler, a general Markov c...
research
07/21/2020

Backfitting for large scale crossed random effects regressions

Regression models with crossed random effect error models can be very ex...

Please sign up or login with your details

Forgot password? Click here to reset