Specifying a Gaussian Markov Random Field by a Sparse Cholesky Triangle

ABSTRACT This note discusses the approach of specifying a Gaussian Markov random field (GMRF) by the Cholesky triangle of the precision matrix. A such representation can be made extremely sparse using numerical techniques for incomplete sparse Cholesky factorization, and provide very computational efficient representation for simulating from the GMRF. However, we provide theoretical and empirical justification showing that the sparse Cholesky triangle representation is fragile when conditioning a GMRF on a subset of the variables or observed data, meaning that the computational cost increases.

[1]  H. Rue,et al.  Norges Teknisk-naturvitenskapelige Universitet Approximating Hidden Gaussian Markov Random Fields Approximating Hidden Gaussian Markov Random Fields , 2003 .

[2]  Peter J. Diggle,et al.  Bayesian analysis of agricultural field experiments - Discussion , 1999 .

[3]  Timothy C. Coburn,et al.  Hierarchical Modeling and Analysis for Spatial Data , 2007 .

[4]  I. Gustafsson.  A class of first order factorization methods , 1978 .

[5]  J. Besag.  Statistical Analysis of Non-Lattice Data , 1975 .

[6]  J. Besag,et al.  Bayesian analysis of agricultural field experiments , 1999 .

[7]  N. Munksgaard,et al.  Solving Sparse Symmetric Sets of Linear Equations by Preconditioned Conjugate Gradients , 1980, TOMS.

[8]  Peter Congdon,et al.  Gaussian Markov Random Fields: Theory and Applications , 2007 .

[9]  L. Held,et al.  Gaussian Markov Random Fields: Theory And Applications (Monographs on Statistics and Applied Probability) , 2005 .

[10]  H. Rue,et al.  On Block Updating in Markov Random Field Models for Disease Mapping , 2002 .

[11]  J. Besag.  Spatial Interaction and the Statistical Analysis of Lattice Systems , 1974 .

[12]  H. Rue,et al.  Fitting Gaussian Markov Random Fields to Gaussian Fields , 2002 .

[13]  Susan A. Murphy,et al.  Monographs on statistics and applied probability , 1990 .

[14]  H. Rue.  Fast sampling of Gaussian Markov random fields , 2000 .

[15]  Håvard Rue,et al.  On block updating in Markov random field models for disease mapping. (REVISED, May 2001) , 2000 .

[16]  H. Künsch.  Gaussian Markov random fields , 1979 .

[17]  J. Meijerink,et al.  An iterative solution method for linear systems of which the coefficient matrix is a symmetric -matrix , 1977 .

[18]  L. Fahrmeir,et al.  Bayesian inference for generalized additive mixed models based on Markov random field priors , 2001 .

[19]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .