Abstract
This paper demonstrates how Gaussian Markov random fields (conditional autoregressions) can be sampled quickly by using numerical techniques for sparse matrices. The algorithm is general and efficient, and expands easily to various forms for conditional simulation and evaluation of normalization constants. We demonstrate its use by constructing efficient block updates in Markov chain Monte Carlo algorithms for disease mapping.
Original language | English (US) |
---|---|
Pages (from-to) | 325-338 |
Number of pages | 14 |
Journal | Journal of the Royal Statistical Society. Series B: Statistical Methodology |
Volume | 63 |
Issue number | 2 |
DOIs | |
State | Published - 2001 |
Externally published | Yes |
Keywords
- Block sampling
- Conditional autoregressive model
- Divide-and-conquer strategy
- Gaussian Markov random field
- Markov chain Monte Carlo methods
ASJC Scopus subject areas
- Statistics and Probability
- Statistics, Probability and Uncertainty