Fast Bayesian Functional Regression for Non-Gaussian Spatial Data

Abstract

Functional generalized linear models (FGLM) have been widely used to study the relations between non-Gaussian response and functional covariates. However, most existing works for FGLM assume independence among observations and therefore they are of limited applicability for correlated data. A particularly important example is spatial functional data, where we observe functions over spatial domains, such as the age population curve or temperature curve at each areal unit. In this paper, we extend FGLM by incorporating spatial random effects. Especially, we study the relationship between the non-Gaussian response variable and functional covariates that are spatially observed. However, such model has computational and inferential challenges. The high-dimensional spatial random effects cause the slow mixing of Markov chain Monte Carlo (MCMC) algorithms. Furthermore, spatial confounding can lead to bias in parameter estimates and inflate their variances. To address these issues, we propose an efficient Bayesian method using a sparse reparameterization of high-dimensional random effects. We also study the average coverage probabilities of the credible intervals of functional parameters. We apply our methods to simulated and real data examples, including malaria incidence data and US COVID-19 data. The proposed method is fast while providing accurate functional estimates.