A Time-Series Model for Areal Data Using Spatially Correlated Gaussian Processes

Traditional spatio-temporal models for areal data typically begin with spatial structure imposed at the level of random effects and later extend to include temporal dynamics. We propose an alternative hierarchical modeling framework that captures temporal trends in areal data through Gaussian processes that share spatial information via correlated variance components. This allows the model to better capture shared patterns of variability across regions while preserving local temporal dynamics, offering a more flexible representation of spatio-temporal processes. Specifically, we extend independent Gaussian-process models for time-series data to a spatially correlated framework by placing a conditional autoregressive (CAR) prior on the parameters governing the temporal variability and imposing a conditional dependence of the temporal range on the temporal variance. We apply this model to two case studies: monthly malaria incidence in Niassa, Mozambique, and weekly food insecurity prevalence in Cameroon. Inference is conducted within a Bayesian framework using approximate posterior sampling. Given the hierarchical structure of the model, we employ a combination of Markov chain Monte Carlo (MCMC) techniques, including the Metropolis-adjusted Langevin algorithm (MALA), Metropolis-Hastings, and Gibbs sampling. In both applications, the model demonstrates strong in-sample performance with narrow credible intervals and outperforms established spatio-temporal approaches in many regions when forecasting. These results underscore the model's ability to capture complex spatio-temporal dependencies while maintaining interpretability, key in settings with sparse data and policy relevance. By accounting for spatio-temporal variation through the evolution of temporal dynamics themselves, our approach offers a flexible and principled tool for many applied contexts. Supplementary material included