Abstract:Spatial models for areal data are often constructed such that all pairs of adjacent regions are assumed to have near‐identical spatial autocorrelation. In practice, data can exhibit dependence structures more complicated than can be represented under this assumption. In this article, we develop a new model for spatially correlated data observed on graphs, which can flexibly represented many types of spatial dependence patterns while retaining aspects of the original graph geometry. Our method implies an embedding of the graph into Euclidean space wherein covariance can be modeled using traditional covariance functions, such as those from the Matérn family. We parameterize our model using a class of graph metrics compatible with such covariance functions, and which characterize distance in terms of network flow, a property useful for understanding proximity in many ecological settings. By estimating the parameters underlying these metrics, we recover the "intrinsic distances" between graph nodes, which assist in the interpretation of the estimated covariance and allow us to better understand the relationship between the observed process and spatial domain. We compare our model to existing methods for spatially dependent graph data, primarily conditional autoregressive models and their variants, and illustrate advantages of our method over traditional approaches. We fit our model to bird abundance data for several species in North Carolina, and show how it provides insight into the interactions between species‐specific spatial distributions and geography.

A full scale approximation of covariance functions for large spatial data sets

Iterative Methods for Full-Scale Gaussian Process Approximations for Large Spatial Data

Covariance approximation for large multivariate spatial data sets with an application to multiple climate model errors

A Vecchia Approximation for High-Dimensional Gaussian Cumulative Distribution Functions Arising from Spatial Data

Rapid Numerical Approximation Method for Integrated Covariance Functions Over Irregular Data Regions

Spatial Process Approximations: Assessing Their Necessity

Visibility graph-based covariance functions for scalable spatial analysis in non-convex partially Euclidean domains

Hybrid Parametric Classes of Isotropic Covariance Functions for Spatial Random Fields

Likelihood Approximation With Hierarchical Matrices For Large Spatial Datasets

Modelling local and global phenomena with sparse Gaussian processes

A flexible and interpretable spatial covariance model for data on graphs

A Variational Approach for Modeling High-dimensional Spatial Generalized Linear Mixed Models

Normalizing Basis Functions: Approximate Stationary Models for Large Spatial Data

Semiparametric Estimation of Cross-covariance Functions for Multivariate Random Fields

Hierarchical low rank approximation of likelihoods for large spatial datasets

Kryging: Geostatistical analysis of large-scale datasets using Krylov subspace methods

Inference for Gaussian Processes with Matérn Covariogram on Compact Riemannian Manifolds

Exploring the Efficacy of Statistical and Deep Learning Methods for Large Spatial Datasets: A Case Study

Posterior Covariance Structures in Gaussian Processes

High-dimensional copula variational approximation through transformation