Adaptive sampling for Bayesian geospatial models

Journal Article (Journal Article)

Bayesian hierarchical modeling with Gaussian process random effects provides a popular approach for analyzing point-referenced spatial data. For large spatial data sets, however, generic posterior sampling is infeasible due to the extremely high computational burden in decomposing the spatial correlation matrix. In this paper, we propose an efficient algorithm—the adaptive griddy Gibbs (AGG) algorithm—to address the computational issues with large spatial data sets. The proposed algorithm dramatically reduces the computational complexity. We show theoretically that the proposed method can approximate the real posterior distribution accurately. The sufficient number of grid points for a required accuracy has also been derived. We compare the performance of AGG with that of the state-of-the-art methods in simulation studies. Finally, we apply AGG to spatially indexed data concerning building energy consumption.

Full Text

Duke Authors

Cited Authors

  • Yang, H; Liu, F; Ji, C; Dunson, D

Published Date

  • November 1, 2014

Published In

Volume / Issue

  • 24 / 6

Start / End Page

  • 1101 - 1110

Electronic International Standard Serial Number (EISSN)

  • 1573-1375

International Standard Serial Number (ISSN)

  • 0960-3174

Digital Object Identifier (DOI)

  • 10.1007/s11222-013-9422-4

Citation Source

  • Scopus