Abstract:Nonstationary and non-Gaussian spatial data are common in various fields, including ecology (e.g., counts of animal species), epidemiology (e.g., disease incidence counts in susceptible regions), and environmental science (e.g., remotely-sensed satellite imagery). Due to modern data collection methods, the size of these datasets have grown considerably. Spatial generalized linear mixed models (SGLMMs) are a flexible class of models used to model nonstationary and non-Gaussian datasets. Despite their utility, SGLMMs can be computationally prohibitive for even moderately large datasets (e.g., 5,000 to 100,000 observed locations). To circumvent this issue, past studies have embedded nested radial basis functions into the SGLMM. However, two crucial specifications (knot placement and bandwidth parameters), which directly affect model performance, are typically fixed prior to model-fitting. We propose a novel approach to model large nonstationary and non-Gaussian spatial datasets using adaptive radial basis functions. Our approach: (1) partitions the spatial domain into subregions; (2) employs reversible-jump Markov chain Monte Carlo (RJMCMC) to infer the number and location of the knots within each partition; and (3) models the latent spatial surface using partition-varying and adaptive basis functions. Through an extensive simulation study, we show that our approach provides more accurate predictions than competing methods while preserving computational efficiency. We demonstrate our approach on two environmental datasets - incidences of plant species and counts of bird species in the United States.

deform: An R Package for Nonstationary Spatial Gaussian Process Models by Deformations and Dimension Expansion

Flexible models for nonstationary dependence: Methodology and examples

Nonrigid Shape Recovery by Gaussian Process Regression

Per-Gaussian Embedding-Based Deformation for Deformable 3D Gaussian Splatting

Modeling Big, Heterogeneous, Non-Gaussian Spatial and Spatio-Temporal Data Using FRK

Modelling of Non-Stationary Spatial Structure Using Parametric Radial Basis Deformations

Estimation of Spatial Deformation for Nonstationary Processes via Variogram Alignment

Modelling Big, Heterogeneous, Non-Gaussian Spatial and Spatio-Temporal Data using FRK

Bayesian Modeling for Nonstationary Spatial Point Process via Spatial Deformations

Flexible Basis Representations for Modeling Large Non-Gaussian Spatial Data

Gaussian Process for Functional Data Analysis: The GPFDA Package for R

Gaussian Markov random field spatial models in GAMLSS

Geostatistics for Large Datasets on Riemannian Manifolds: A Matrix-Free Approach

Mixed-Stationary Gaussian Process for Flexible Non-Stationary Modeling of Spatial Outcomes

Estimating deformations of isotropic Gaussian random fields on the plane

Deformation prediction of foundation pit using Gaussian Process machine learning

Gaussian Process Regression for Transportation System Estimation and Prediction Problems: the Deformation and a Hat Kernel

Estimates of Non-Stationary Spatial Structure Using Gaussian Radial Basis Deformations

spNNGP R package for Nearest Neighbor Gaussian Process models

Manifold Gaussian Processes for Regression

Non-stationary Gaussian Process Surrogates