ISM Research Memorandum
No.
906
Title:
fMRI activation maps based on the NN-ARx model
Author(s):
Riera, J.(NICHe, Tohoku Univ.), Bosch, J.(Cuban Neuroscience Center), Yamashita, O.(The Institute of Statis tical Mathematics), Kawashima, R.(NICHe, Tohoku Univ.), Sadato, N.(National Institute for Physiological Science), Okada, T.(Institute of Biomedical Research and Innovation), Ozaki, T.(The Institute of Statistical Mathematics)
Key words:
Hemodynamics Approach; Autoregressive Models; fMRI; Neuronal Activation.
Abstract:
The most significant progresses in the understanding of human brain functions have been possible due to the use of fMRI, which, when used in combination with other standard neuroimaging techniques (i.e. EEG), provides researchers with a potential tool to elucidate many biophysical principles, established previously by means of animal comparative studies. However, to date most of the methods proposed in the literature seeking fMRI signs have been limited to the use of a top-down data analysis approach, thus ignoring a pool of physiological facts. In spite of the important contributions achieved by applying these methods to actual data, there is a disproportionate gap between theoretical models and data-analysis strategies while trying to focus on several new prospects, like for example fMRI/EEG data fusion, non-linear BOLD signal dynamics and causality/connectivity patterns. In this paper, we propose a new approach which will allow many of the above mentioned hot topics to be addressed in the near future with an underlying interpretability based on bottom-up modeling. In particular, the a-MAP presented in the paper to test brain activation corresponds very well with the standardized T-test of the SPM99 toolbox. Additionally, a new Impulse Response Function has been formulated, directly related to the sell-established concept of the Hemodynamics Response Function (HRF). The model uses not only the information contained in the signal but also that in the structure of the background noise to simultaneously estimate the HRF and the autocorrelation function by using an autoregressive model with a filtered Poisson process driving the dynamics. The short-range contributions of near-neighbor voxels are also included, and the potential drift was characterized by a polynomial series. Since our model originated from an immediate extension of the hemodynamics approach, a natural interpretability of the results is feasible.