Gradient Scan Gibbs Sampler: an efficient algorithm for high-dimensional Gaussian distributions

Abstract : This paper deals with Gibbs samplers that include high dimensional conditional Gaussian distributions. It proposes an efficient algorithm that avoids the high dimensional Gaussian sampling and relies on a random excursion along a small set of directions. The algorithm is proved to converge, i.e. the drawn samples are asymptotically distributed according to the target distribution. Our main motivation is in inverse problems related to general linear observation models and their solution in a hierarchical Bayesian framework implemented through sampling algorithms. It finds direct applications in semi-blind / unsupervised methods as well as in some non-Gaussian methods. The paper provides an illustration focused on the unsupervised estimation for super-resolution methods.
Complete list of metadatas

Cited literature [44 references]  Display  Hide  Download

https://hal-supelec.archives-ouvertes.fr/hal-01252598
Contributor : François Orieux <>
Submitted on : Thursday, January 7, 2016 - 5:24:34 PM
Last modification on : Monday, October 14, 2019 - 5:02:03 PM
Long-term archiving on : Friday, April 8, 2016 - 1:34:56 PM

File

gsgs-paper.pdf
Files produced by the author(s)

Identifiers

Citation

Olivier Féron, François Orieux, Jean-François Giovannelli. Gradient Scan Gibbs Sampler: an efficient algorithm for high-dimensional Gaussian distributions. IEEE Journal of Selected Topics in Signal Processing, IEEE, 2016, 10 (2), pp.343-352. ⟨10.1109/JSTSP.2015.2510961⟩. ⟨hal-01252598⟩

Share

Metrics

Record views

403

Files downloads

331