Time Series Analysis of fMRI Data: Spatial Modelling and Bayesian Computation

Ming Teng,Timothy Johnson,Farouk Nathoo
DOI: https://doi.org/10.48550/arXiv.1609.02123
2018-03-04
Abstract:Time series analysis of fMRI data is an important area of medical statistics for neuroimaging data. The neuroimaging community has embraced mean-field variational Bayes (VB) approximations, which are implemented in Statistical Parametric Mapping (SPM) software. While computationally efficient, the quality of VB approximations remains unclear even though they are commonly used in the analysis of neuroimaging data. For reliable statistical inference, it is important that these approximations be accurate and that users understand the scenarios under which they may not be accurate. We consider this issue for a particular model that includes spatially-varying coefficients. To examine the accuracy of the VB approximation we derive Hamiltonian Monte Carlo (HMC) for this model and conduct simulation studies to compare its performance with VB. As expected we find that the computation time required for VB is considerably less than that for HMC. In settings involving a high or moderate signal-to-noise ratio (SNR) we find that the two approaches produce very similar results suggesting that the VB approximation is useful in this setting. On the other hand, when one considers a low SNR, substantial differences are found, suggesting that the approximation may not be accurate in such cases and we demonstrate that VB produces Bayes estimators with larger mean squared error (MSE). A real application related to face perception is also carried out. Overall, our work clarifies the usefulness of VB for the spatiotemporal analysis of fMRI data, while also pointing out the limitation of VB when the SNR is low and the utility of HMC in this case.
Computation
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to evaluate the accuracy of the variational Bayes (VB) approximation in the time - series analysis of functional magnetic resonance imaging (fMRI) data, especially its performance under low signal - to - noise ratio (SNR) conditions. Specifically, the author focuses on models containing spatially varying coefficients and derives the Hamiltonian Monte Carlo (HMC) algorithm for comparison with the VB method. The main objectives of the study are: 1. **Develop the HMC algorithm**: For the first time, an HMC algorithm is derived for Bayesian inference in the time - series analysis of fMRI data. 2. **Compare HMC and VB**: Through detailed simulation studies and the application of actual fMRI data sets, compare the performance of the HMC algorithm and the VB algorithm currently used in the SPM software. The comparison is mainly carried out in terms of estimation accuracy, posterior variability, spatial smoothness of the estimated image, and computation time. Through these comparisons, the author hopes to determine in which cases the VB approximation is accurate and in which cases it may be inaccurate, so as to provide guidance for practitioners to better select and implement these models. In addition, the author also provides the C code for the HMC implementation for researchers to use.