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.

Fast Bayesian Inference for Spatial Mean-Parameterized Conway–Maxwell–Poisson Models

Bayesian statistics in spatial epidemiology]

Fast Bayesian estimation of spatial count data models

A bayesian multivariate mixture model for high throughput spatial transcriptomics

Bayesian estimation for longitudinal data in a joint model with HPCs

Bayesian nonparametric models for spatially indexed data of mixed type

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

Dimension Reduction and Alleviation of Confounding for Spatial Generalized Linear Mixed Models

Efficient Estimation of COM-Poisson Regression and Generalized Additive Model

Efficient Estimation of COM-Poisson Regression and a Generalized Additive Model

Spatial Factor Modeling: A Bayesian Matrix-Normal Approach for Misaligned Data

New biased estimators for the Conway–Maxwell-Poisson Model

Bayesian inference for high-dimensional discrete-time epidemic models: spatial dynamics of the UK COVID-19 outbreak

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

Bayesian modeling of spatial molecular profiling data via Gaussian process

Computationally Scalable Bayesian SPDE Modeling for Censored Spatial Responses

A Tree-Based Semi-Varying Coefficient Model for the COM-Poisson Distribution

Bayesian Bivariate Conway-Maxwell-Poisson Regression Model for Correlated Count Data in Sports

A Bayesian Zero-Inflated Dirichlet-Multinomial Regression Model for Multivariate Compositional Count Data

Bayesian Inference for Finite Populations Under Spatial Process Settings

Flexible Basis Representations for Modeling Large Non-Gaussian Spatial Data