Graphical Gaussian Process Models for Highly Multivariate Spatial Data
For multivariate spatial (Gaussian) process models, common cross-covariance functions do not exploit graphical models to ensure process-level conditional independence among the variables. This is undesirable, especially for highly multivariate settings, where popular cross-covariance functions such as the multivariate Matérn suffer from a "curse of dimensionality" as the number of parameters and floating point operations scale up in quadratic and cubic order, respectively, in the number of variables. We propose a class of multivariate "graphical Gaussian Processes" using a general construction called "stitching" that crafts cross-covariance functions from graphs and ensure process-level conditional independence among variables. For the Matérn family of functions, stitching yields a multivariate GP whose univariate components are exactly Matérn GPs, and conforms to process-level conditional independence as specified by the graphical model. For highly multivariate settings and decomposable graphical models, stitching offers massive computational gains and parameter dimension reduction. We demonstrate the utility of the graphical Matérn GP to jointly model highly multivariate spatial data using simulation examples and an application to air-pollution modelling.
READ FULL TEXT