An enhanced Gibbs sampler algorithm for non-conditional simulation of Gaussian random vectors

Authors:
Daisy Arroyo;Xavier Emery;MaríA PeláEz
Affiliations:
Department of Mathematics, Northern Catholic University, Chile;Department of Mining Engineering, University of Chile, Chile and ALGES Laboratory, Advanced Mining Technology Center, University of Chile, Chile;Department of Mathematics, Northern Catholic University, Chile
Venue:
Computers & Geosciences
Year:
2012

Citing 4
Cited 0

Using the Gibbs sampler for conditional simulation of Gaussian-based random fields

Computers & Geosciences
TBSIM: A computer program for conditional simulation of three-dimensional Gaussian random fields via the turning bands method

Computers & Geosciences
Truncated Gaussian simulation of discrete-valued, ordinal coregionalized variables

Computers & Geosciences
Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images

IEEE Transactions on Pattern Analysis and Machine Intelligence

Quantified Score

Hi-index	0.00

Visualization

Abstract

This paper addresses the problem of simulating a Gaussian random vector with zero mean and given variance-covariance matrix, without conditioning constraints. Variants of the Gibbs sampler algorithm are presented, based on the proposal by Galli and Gao, which do not require inverting the variance-covariance matrix and therefore allow considerable time savings. Numerical experiments are performed to check the accuracy of the algorithm and to determine implementation parameters (in particular, the updating and blocking strategies) that increase the rates of convergence and mixing.