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.

Exploring the Distributional Properties of the Non-Gaussian Random Field Models

Spatial Hyperspheric Models for Compositional Data

A multivariate Gaussian random field prior against spatial confounding

Modeling spatial extremes using normal mean-variance mixtures

Additive Models with Spatio-Temporal Data

Flexible Basis Representations for Modeling Large Non-Gaussian Spatial Data

Nearest-Neighbor Mixture Models for Non-Gaussian Spatial Processes

Gaussian Markov random field spatial models in GAMLSS

Modelling Point Referenced Spatial Count Data: A Poisson Process Approach

Hybrid Parametric Classes of Isotropic Covariance Functions for Spatial Random Fields

Estimating covariance functions of multivariate skew-Gaussian random fields on the sphere

Birnbaum–Saunders spatial modelling and diagnostics applied to agricultural engineering data

Multivariate nearest‐neighbors Gaussian processes with random covariance matrices

Multivariate Matérn Models -- A Spectral Approach

A flexible Bayesian non-confounding spatial model for analysis of dispersed count data in clinical studies

Nonstationary Spatial Process Models with Spatially Varying Covariance Kernels

Bayesian Nonstationary Spatial Modeling for Very Large Datasets

Bayesian nonparametric location-scale-shape mixtures

Highly Multivariate High-dimensionality Spatial Stochastic Processes-A Mixed Conditional Approach

Clustering using skewed multivariate heavy tailed distributions with flexible tail behaviour

A method of generating multivariate non-normal random numbers with desired multivariate skewness and kurtosis