High-dimensional Multivariate Geostatistics: A Bayesian Matrix-Normal Approach

Abstract

Joint modeling of spatially-oriented dependent variables is commonplace in the environmental sciences, where scientists seek to estimate the relationships among a set of environmental outcomes accounting for dependence among these outcomes and the spatial dependence for each outcome. Such modeling is now sought for massive data sets with variables measured at a very large number of locations. Bayesian inference, while attractive for accommodating uncertainties through hierarchical structures, can become computationally onerous for modeling massive spatial data sets because of its reliance on iterative estimation algorithms. This manuscript develops a conjugate Bayesian framework for analyzing multivariate spatial data using analytically tractable posterior distributions that obviate iterative algorithms. We discuss differences between modeling the multivariate response itself as a spatial process and that of modeling a latent process in a hierarchical model. We illustrate the computational and inferential benefits of these models using simulation studies and analysis of a Vegetation Index data set with spatially dependent observations numbering in the millions.

arXiv comment22 pages, 12 figures

References in corpus (6)

• Cross-Covariance Functions for Multivariate Geostatistics
• Permutation and Grouping Methods for Sharpening Gaussian Process Approximations
• Multivariate spatio-temporal models for high-dimensional areal data with application to Longitudinal Employer-Household Dynamics
• Discussion on Computationally Efficient Multivariate Spatio-Temporal Models for High-Dimensional Count-Valued Data by Bradley et al
• Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains
• Practical Bayesian Modeling and Inference for Massive Spatial Datasets On Modest Computing Environments

Cited by in corpus (8)

• Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data
• Modeling massive highly-multivariate nonstationary spatial data with the basis graphical lasso
• Correlation-based sparse inverse Cholesky factorization for fast Gaussian-process inference
• A Framework of Zero-Inflated Bayesian Negative Binomial Regression Models For Spatiotemporal Data
• Bayesian Signal Matching for Transfer Learning in ERP-Based Brain Computer Interface
• Gridding and Parameter Expansion for Scalable Latent Gaussian Models of Spatial Multivariate Data
• Applications of Conjugate Gradient in Bayesian computation
• Inside-out cross-covariance for spatial multivariate data