Inside-out cross-covariance for spatial multivariate data

Abstract

As the spatial features of multivariate data are increasingly central in researchers' applied problems, there is a growing demand for novel spatially-aware methods that are flexible, easily interpretable, and scalable to large data. We develop inside-out cross-covariance (IOX) models for multivariate spatial likelihood-based inference. IOX leads to valid cross-covariance matrix functions which we interpret as inducing spatial dependence on independent replicates of a correlated random vector. The resulting sample cross-covariance matrices are "inside-out" relative to the ubiquitous linear model of coregionalization (LMC). However, unlike LMCs, our methods offer direct marginal inference, easy prior elicitation of covariance parameters, the ability to model outcomes with unequal smoothness, and flexible dimension reduction. As a covariance model for a q-variate Gaussian process, IOX leads to scalable models for noisy vector data as well as flexible latent models. For large n cases, IOX complements Vecchia approximations and related process-based methods based on sparse graphical models. We demonstrate superior performance of IOX on synthetic datasets as well as on colorectal cancer proteomics data. An R package implementing the proposed methods is available at github.com/mkln/spiox.

References in corpus (9)

• Cross-Covariance Functions for Multivariate Geostatistics
• A general framework for Vecchia approximations of Gaussian processes
• A multi-resolution approximation for massive spatial datasets
• 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
• Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data
• High-dimensional Multivariate Geostatistics: A Bayesian Matrix-Normal Approach
• Modeling massive highly-multivariate nonstationary spatial data with the basis graphical lasso