Experimental fMRI Data Research Articles

Spatiotemporal clustering of fMRI time series in the spectral domain

We propose a new method for the analysis of functional magnetic resonance images (fMRI). The decision that a voxel v 0 is activated is based not solely on the value of the fMRI signal at v 0, but rather on the comparison of all time series s v ( t) in a small neighborhood N ( v 0 ) around v 0. Our approach explicitly takes into account the intrinsic spatiotemporal correlations that exist in the data. We focus on experimental designs with periodic stimuli, and therefore we can capture most of the features of the BOLD signal with a low dimensional subspace in the frequency domain. The presence of activated time series can be detected by partitioning the time series in this low dimensional space. Experiments with simulated data, and experimental fMRI data, demonstrate that our approach can outperform standard methods of analysis, such as the t-test.

A signal subspace approach for modeling the hemodynamic response function in fMRI.

Many fMRI analysis methods use a model for the hemodynamic response function (HRF). Common models of the HRF, such as the Gaussian or Gamma functions, have parameters that are usually selected a priori by the data analyst. A new method is presented that characterizes the HRF over a wide range of parameters via three basis signals derived using principal component analysis (PCA). Covering the HRF variability, these three basis signals together with the stimulation pattern define signal subspaces which are applicable to both linear and nonlinear modeling and identification of the HRF and for various activation detection strategies. Analysis of simulated fMRI data using the proposed signal subspace showed increased detection sensitivity compared to the case of using a previously proposed trigonometric subspace. The methodology was also applied to activation detection in both event-related and block design experimental fMRI data using both linear and nonlinear modeling of the HRF. The activated regions were consistent with previous studies, indicating the ability of the proposed approach in detecting brain activation without a priori assumptions about the shape parameters of the HRF. The utility of the proposed basis functions in identifying the HRF is demonstrated by estimating the HRF in different activated regions.

Wavelet-based estimation of a semiparametric generalized linear model of fMRI time-series

This paper addresses the problem of detecting significant changes in fMRI time series that are correlated to a stimulus time course. This paper provides a new approach to estimate the parameters of a semiparametric generalized linear model of fMRI time series. The fMRI signal is described as the sum of two effects: a smooth trend and the response to the stimulus. The trend belongs to a subspace spanned by large scale wavelets. The wavelet transform provides an approximation to the Karhunen-Loève transform for the long memory noise and we have developed a scale space regression that permits to carry out the regression in the wavelet domain while omitting the scales that are contaminated by the trend. In order to demonstrate that our approach outperforms the state-of-the art detrending technique, we evaluated our method against a smoothing spline approach. Experiments with simulated data and experimental fMRI data, demonstrate that our approach can infer and remove drifts that cannot be adequately represented with splines.

Multivariate Model Specification for fMRI Data

We present a general method—denoted MoDef—to help specify (or define) the model used to analyze brain imaging data. This method is based on the use of the multivariate linear model on a training data set. We show that when the a priori knowledge about the expected brain response is not too precise, the method allows for the specification of a model that yields a better sensitivity in the statistical results. This obviously relies on the validity of the a priori information, in our case the representativity of the training set, an issue addressed using a cross-validation technique. We propose a fast implementation that allows the use of the method on large data sets as found with functional Magnetic Resonance Images. An example of application is given on an experimental fMRI data set that includes nine subjects who performed a mental computation task. Results show that the method increases the statistical sensitivity of fMRI analyses.

Colored noise and computational inference in neurophysiological (fMRI) time series analysis: resampling methods in time and wavelet domains.

Even in the absence of an experimental effect, functional magnetic resonance imaging (fMRI) time series generally demonstrate serial dependence. This colored noise or endogenous autocorrelation typically has disproportionate spectral power at low frequencies, i.e., its spectrum is (1/f)-like. Various pre-whitening and pre-coloring strategies have been proposed to make valid inference on standardised test statistics estimated by time series regression in this context of residually autocorrelated errors. Here we introduce a new method based on random permutation after orthogonal transformation of the observed time series to the wavelet domain. This scheme exploits the general whitening or decorrelating property of the discrete wavelet transform and is implemented using a Daubechies wavelet with four vanishing moments to ensure exchangeability of wavelet coefficients within each scale of decomposition. For (1/f)-like or fractal noises, e.g., realisations of fractional Brownian motion (fBm) parameterised by Hurst exponent 0 < H < 1, this resampling algorithm exactly preserves wavelet-based estimates of the second order stochastic properties of the (possibly nonstationary) time series. Performance of the method is assessed empirically using (1/f)-like noise simulated by multiple physical relaxation processes, and experimental fMRI data. Nominal type 1 error control in brain activation mapping is demonstrated by analysis of 13 images acquired under null or resting conditions. Compared to autoregressive pre-whitening methods for computational inference, a key advantage of wavelet resampling seems to be its robustness in activation mapping of experimental fMRI data acquired at 3 Tesla field strength. We conclude that wavelet resampling may be a generally useful method for inference on naturally complex time series.