Improving the accuracy of brain activation maps in the group-level analysis of fMRI data utilizing spatiotemporal Gaussian process model

Abstract

ObjectiveAccuracy and precision of the statistical analysis methods used for brain activation maps are essential. Adjusting models to consider spatiotemporal correlation embedded in fMRI data may increase their accuracy, but it also introduces a high computational cost. The present study aimed to apply and assess the spatiotemporal Gaussian process (STGP) model to improve accuracy and reduce cost. MethodsWe applied the spatiotemporal Gaussian process (STGP) model for both simulated and experimental memory tfMRI data and compared the findings with fast, fully Bayesian, and General Linear Models (GLM). To assess their accuracy and precision, the models were fitted to the simulated data (1000 voxels,100 times point for 50 people), and an average of accuracy indexes of 100 repetitions was computed. Functional and activation maps for all models were calculated in experimental data analysis. ResultsSTGP model resulted in a higher Z-score in the whole brain, in the 1000 most activated voxels, and in the frontal lobe as the approved memory area. Based on the simulated data, the STGP model showed more accuracy and precision than the other two models. However, its computational time was more than the GLM, as the price of model correction, but much less than that of the fast, fully Bayesian model. ConclusionSpatiotemporal correlation further improved the accuracy of the STGP compared to the GLM and fast, fully Bayesian model. This can result in more accurate activation maps. Moreover, the STGP model’s computational speed appears to be reasonable for model application.