Transformed Gaussian Markov Random Fields and Spatial Modeling

TL;DR

This paper introduces transformed Gaussian Markov random fields (TGMRFs) that allow for asymmetric and heavy-tailed spatial dependence modeling, improving upon traditional Gaussian-based models in ecological applications.

Contribution

It proposes a novel class of spatial models using transformations of GMRFs with Gaussian copulas, enabling flexible marginal distributions and better ecological data modeling.

Findings

TGMRFs outperform traditional models in ecological data analysis

Bayesian inference is effectively applied to the new models

Simulation studies validate the models' accuracy and robustness

Abstract

The Gaussian random field (GRF) and the Gaussian Markov random field (GMRF) have been widely used to accommodate spatial dependence under the generalized linear mixed model framework. These models have limitations rooted in the symmetry and thin tail of the Gaussian distribution. We introduce a new class of random fields, termed transformed GRF (TGRF), and a new class of Markov random fields, termed transformed GMRF (TGMRF). They are constructed by transforming the margins of GRFs and GMRFs, respectively, to desired marginal distributions to accommodate asymmetry and heavy tail as needed in practice. The Gaussian copula that characterizes the dependence structure facilitates inferences and applications in modeling spatial dependence. This construction leads to new models such as gamma or beta Markov fields with Gaussian copulas, which can be used to model Poisson intensity or Bernoulli rate in a spatial generalized linear mixed model. The method is naturally implemented in a Bayesian framework. We illustrate the utility of the methodology in an ecological application with spatial count data and spatial presence/absence data of some snail species, where the new models are shown to outperform the traditional spatial models. The validity of Bayesian inferences and model selection are assessed through simulation studies for both spatial Poisson regression and spatial Bernoulli regression.