A wavelet-based Bayesian approach to regression models with long memory errors and its application to FMRI data

This article considers linear regression models with long memory errors. These models have been proven useful for application in many areas, such as medical imaging, signal processing, and econometrics. Wavelets, being self-similar, have a strong connection to long memory data. Here we employ discrete wavelet transforms as whitening filters to simplify the dense variance-covariance matrix of the data. We then adopt a Bayesian approach for the estimation of the model parameters. Our inferential procedure uses exact wavelet coefficients variances and leads to accurate estimates of the model parameters. We explore performances on simulated data and present an application to an fMRI data set. In the application we produce posterior probability maps (PPMs) that aid interpretation by identifying voxels that are likely activated with a given confidence.