Open Access
November 2017 Accelerated Gibbs sampling of normal distributions using matrix splittings and polynomials
Colin Fox, Albert Parker
Bernoulli 23(4B): 3711-3743 (November 2017). DOI: 10.3150/16-BEJ863


Standard Gibbs sampling applied to a multivariate normal distribution with a specified precision matrix is equivalent in fundamental ways to the Gauss–Seidel iterative solution of linear equations in the precision matrix. Specifically, the iteration operators, the conditions under which convergence occurs, and geometric convergence factors (and rates) are identical. These results hold for arbitrary matrix splittings from classical iterative methods in numerical linear algebra giving easy access to mature results in that field, including existing convergence results for antithetic-variable Gibbs sampling, REGS sampling, and generalizations. Hence, efficient deterministic stationary relaxation schemes lead to efficient generalizations of Gibbs sampling. The technique of polynomial acceleration that significantly improves the convergence rate of an iterative solver derived from a symmetric matrix splitting may be applied to accelerate the equivalent generalized Gibbs sampler. Identicality of error polynomials guarantees convergence of the inhomogeneous Markov chain, while equality of convergence factors ensures that the optimal solver leads to the optimal sampler. Numerical examples are presented, including a Chebyshev accelerated SSOR Gibbs sampler applied to a stylized demonstration of low-level Bayesian image reconstruction in a large 3-dimensional linear inverse problem.


Download Citation

Colin Fox. Albert Parker. "Accelerated Gibbs sampling of normal distributions using matrix splittings and polynomials." Bernoulli 23 (4B) 3711 - 3743, November 2017.


Received: 1 April 2015; Revised: 1 December 2015; Published: November 2017
First available in Project Euclid: 23 May 2017

zbMATH: 06778301
MathSciNet: MR3654821
Digital Object Identifier: 10.3150/16-BEJ863

Keywords: Bayesian inference , Gaussian Markov random field , Gibbs sampling , matrix splitting , multivariate normal distribution , non-stationary stochastic iteration , polynomial acceleration

Rights: Copyright © 2017 Bernoulli Society for Mathematical Statistics and Probability

Vol.23 • No. 4B • November 2017
Back to Top