Score: 0

Scalable Bayesian inference for high-dimensional mixed-type multivariate spatial data

Published: October 15, 2025 | arXiv ID: 2510.13233v1

By: Arghya Mukherjee, Arnab Hazra, Dootika Vats

Potential Business Impact:

Models different kinds of data together in places.

Business Areas:
Geospatial Data and Analytics, Navigation and Mapping

Spatial generalized linear mixed-effects methods are popularly used to model spatially indexed univariate responses. However, with modern technology, it is common to observe vector-valued mixed-type responses, e.g., a combination of binary, count, or continuous types, at each location. Methods that allow joint modeling of such mixed-type multivariate spatial responses are rare. Using latent multivariate Gaussian processes (GPs), we present a class of Bayesian spatial methods that can be employed for any combination of exponential family responses. Since multivariate GP-based methods can suffer from computational bottlenecks when the number of spatial locations is high, we further employ a computationally efficient Vecchia approximation for fast posterior inference and prediction. Key theoretical properties of the proposed model, such as identifiability and the structure of the induced covariance, are established. Our approach employs a Markov chain Monte Carlo-based inference method that utilizes elliptical slice sampling in a blocked Metropolis-within-Gibbs sampling framework. We illustrate the efficacy of the proposed method through simulation studies and a real-data application on joint modeling of wildfire counts and burnt areas across the United States.

Page Count
38 pages

Category
Statistics:
Methodology