Scalable Bayesian inference for high-dimensional mixed-type multivariate spatial data
Fuente:
arXiv
Saved in:
| Main Authors: | , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866913052637528064 |
|---|---|
| author | Mukherjee, Arghya Hazra, Arnab Vats, Dootika |
| author_facet | Mukherjee, Arghya Hazra, Arnab Vats, Dootika |
| contents | Spatial generalized linear mixed-effects models are popularly used to analyze 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 for jointly modeling such mixed-type multivariate spatial responses are rare. Using multivariate Gaussian processes (GPs) in the latent layer, we present a class of Bayesian spatial methods applicable to 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 uses elliptical slice sampling within 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. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2510_13233 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Scalable Bayesian inference for high-dimensional mixed-type multivariate spatial data Mukherjee, Arghya Hazra, Arnab Vats, Dootika Methodology Computation 62M30 Spatial generalized linear mixed-effects models are popularly used to analyze 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 for jointly modeling such mixed-type multivariate spatial responses are rare. Using multivariate Gaussian processes (GPs) in the latent layer, we present a class of Bayesian spatial methods applicable to 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 uses elliptical slice sampling within 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. |
| title | Scalable Bayesian inference for high-dimensional mixed-type multivariate spatial data |
| topic | Methodology Computation 62M30 |
| url | https://arxiv.org/abs/2510.13233 |