Self-correlated spatial random variables: From an auto- to a sui- model respecification

Abstract

This paper marks the 50-year publication anniversary of Besag's seminal spatial auto- models paper. His classic article synthesizes generic autoregressive specifications (i.e., a response variable appears on both sides of a regression equation and/or probability function equal sign) for the following six popular random variables: normal, logistic (i.e., Bernoulli), binomial, Poisson, exponential, and gamma. Besag dismisses these last two while recognizing failures of both as well as the more scientifically critical counts-oriented auto-Poisson. His initially unsuccessful subsequent work first attempted to repair them (e.g., pseudo-likelihood estimation), and then successfully revise them within the context of mixed models, formulating a spatially structured random effects term that effectively and efficiently absorbs and accounts for spatial autocorrelation in geospatial data. One remaining weakness of all but the auto-normal is a need to resort to Markov chain Monte Carlo (MCMC) techniques for legitimate estimation purposes. Recently, Griffith succeeded in devising an innovative uniform distribution genre—sui-uniform random variables—that accommodates spatial autocorrelation, too. Its most appealing feature is that, by applying two powerful mathematical statistical theorems (i.e., the probability integral transform, and the quantile function), it redeems Besag's auto- model failures. This paper details conversion of Besag's initial six modified variates, exemplifying them with both simulation experiments and publicly accessible real-world georeferenced data. The principal outcome is valuable spatial statistical advancements, with special reference to Moran eigenvector spatial filtering.