Journal
SCANDINAVIAN JOURNAL OF STATISTICS
Volume 29, Issue 1, Pages 31-49Publisher
BLACKWELL PUBL LTD
DOI: 10.1111/1467-9469.00058
Keywords
conditional autoregressive models; empirical Bayes; Gaussian fields; Gaussian Markov random field; Kullback-Leibler discrepancy; Markov chain Monte Carlo; maximum likelihood
Categories
Ask authors/readers for more resources
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.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available