Spatial Multivariate Trees for Big Data Bayesian Regression
arXiv:2012.00943
Abstract
High resolution geospatial data are challenging because standard geostatistical models based on Gaussian processes are known to not scale to large data sizes. While progress has been made towards methods that can be computed more efficiently, considerably less attention has been devoted to big data methods that allow the description of complex relationships between several outcomes recorded at high resolutions by different sensors. Our Bayesian multivariate regression models based on spatial multivariate trees (SpamTrees) achieve scalability via conditional independence assumptions on latent random effects following a treed directed acyclic graph. Information-theoretic arguments and considerations on computational efficiency guide the construction of the tree and the related efficient sampling algorithms in imbalanced multivariate settings. In addition to simulated data examples, we illustrate SpamTrees using a large climate data set which combines satellite data with land-based station data. Source code is available at https://github.com/mkln/spamtree
References in corpus (7)
- Cross-Covariance Functions for Multivariate Geostatistics
- Permutation and Grouping Methods for Sharpening Gaussian Process Approximations
- Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains
- Multiresolution Gaussian Processes
- Modeling Massive Spatial Datasets Using a Conjugate Bayesian Linear Regression Framework
- Fast increased fidelity approximate Gibbs samplers for Bayesian Gaussian process regression
- Conjugate Nearest Neighbor Gaussian Process Models for Efficient Statistical Interpolation of Large Spatial Data