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.
Original language | English (US) |
---|---|
Pages (from-to) | 161-176 |
Number of pages | 16 |
Journal | Communications in Statistics: Simulation and Computation |
Volume | 35 |
Issue number | 1 |
DOIs | |
State | Published - Jan 1 2006 |
Externally published | Yes |
Keywords
- Gaussian Markov random field
- Incomplete Cholesky factorization
- Parameterization
- Sparse matrices
ASJC Scopus subject areas
- Statistics and Probability
- Modeling and Simulation