Premium
Fitting Gaussian Markov Random Fields to Gaussian Fields
Author(s) -
RUE HÅÅVARD,
TJELMELAND HÅÅKON
Publication year - 2002
Publication title -
scandinavian journal of statistics
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.359
H-Index - 65
eISSN - 1467-9469
pISSN - 0303-6898
DOI - 10.1111/1467-9469.00058
Subject(s) - markov chain , gaussian , random field , mathematics , markov chain monte carlo , gaussian random field , statistical physics , gaussian process , variable order markov model , markov property , inference , markov model , bayesian probability , algorithm , artificial intelligence , computer science , statistics , physics , quantum mechanics
This paper discusses the following task often encountered in building Bayesian spatial models: construct a homogeneous Gaussian Markov random field (GMRF) on a lattice with correlation properties either as present in some observed data, or consistent with prior knowledge. The Markov property is essential in designing computationally efficient Markov chain Monte Carlo algorithms to analyse such models. We argue that we can restate both tasks as that of fitting a GMRF to a prescribed stationary Gaussian field on a lattice when both local and global properties are important. We demonstrate that using the Kullback–Leibler discrepancy often fails for this task, giving severely undesirable behaviour of the correlation function for lags outside the neighbourhood. We propose a new criterion that resolves this difficulty, and demonstrate that GMRFs with small neighbourhoods can approximate Gaussian fields surprisingly well even with long correlation lengths. Finally, we discuss implications of our findings for likelihood based inference for general Markov random fields when global properties are also important.