Abstract:Spatial modelling often uses Gaussian random fields to capture the stochastic nature of studied phenomena. However, this approach incurs significant computational burdens (O(n3)), primarily due to covariance matrix computations. In this study, we propose to use a low-rank approximation of a thin plate spline as a spatial random effect in Bayesian spatial models. We compare its statistical performance and computational efficiency with the approximated Gaussian random field (by the SPDE method). In this case, the dense matrix of the thin plate spline is approximated using a truncated spectral decomposition, resulting in computational complexity of O(kn2) operations, where k is the number of knots. Bayesian inference is conducted via the Hamiltonian Monte Carlo algorithm of the probabilistic software Stan, which allows us to evaluate performance and diagnostics for the proposed models. A simulation study reveals that both models accurately recover the parameters used to simulate data. However, models using a thin plate spline demonstrate superior execution time to achieve the convergence of chains compared to the models utilizing an approximated Gaussian random field. Furthermore, thin plate spline models exhibited better computational efficiency for simulated data coming from different spatial locations. In a real application, models using a thin plate spline as spatial random effect produced similar results in estimating a relative index of abundance for a benthic marine species when compared to models incorporating an approximated Gaussian random field. Although they were not the more computational efficient models, their simplicity in parametrization, execution time and predictive performance make them a valid alternative for spatial modelling under Bayesian inference.

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

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

Flexible Basis Representations for Modeling Large Non-Gaussian Spatial Data

Estimation of Spatial Deformation for Nonstationary Processes via Variogram Alignment

Flexible models for nonstationary dependence: Methodology and examples

Online Spatiotemporal Modeling for High Spatial-Dimensional DPSs under Nonstationary Sensor Layout

Spatial Functional Data analysis: Irregular spacing and Bernstein polynomials

Nonstationary Spatial Gaussian Markov Random Fields

Bayesian Nonstationary Spatial Modeling for Very Large Datasets

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

Improving Piecewise Linear Snow Density Models through Hierarchical Spatial and Orthogonal Functional Smoothing

Bayesian Modeling with Spatial Curvature Processes

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

Bayesian Modeling for Nonstationary Spatial Point Process via Spatial Deformations

Bayesian Inference of Spatially Varying Parameters in Soil Constitutive Models by Using Deformation Observation Data

Bayesian Function-on-Function Regression for Spatial Functional Data

Modeling random isotropic vector fields on the sphere: theory and application to the noise in GNSS station position time series

Nonparametric Spherical Regression Using Diffeomorphic Mappings

Why not a thin plate spline for spatial models? A comparative study using Bayesian inference

Spatial Hyperspheric Models for Compositional Data

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