Modeling the haemodynamic response in fMRI using smooth FIR filters

Area of Science:

• Neuroimaging
• Statistical Modeling
• Biophysics

Background:

• Accurate modeling of the haemodynamic response is crucial for analyzing functional magnetic resonance imaging (fMRI) data.
• Previous methods relied on limited parametric response functions.
• A more flexible approach is needed to capture the complex haemodynamic response.

Purpose of the Study:

• To develop a semi-parametric model for the haemodynamic response function (HRF) in fMRI.
• To incorporate Gaussian process priors to manage the increased degrees of freedom.
• To compare the proposed model with existing methods.

Main Methods:

• Utilized a semi-parametric approach based on finite impulse response (FIR) filters.
• Introduced Gaussian process priors on FIR filter parameters.
• Employed evidence framework for hyper-parameter optimization and Markov Chain Monte Carlo (MCMC) for sampling.
• Validated the model on simulated and real visual stimulation fMRI data.

Main Results:

• The proposed Gaussian process-FIR model offers a flexible alternative to standard haemodynamic response kernels.
• The model effectively handles the increased complexity of the haemodynamic response.
• Demonstrated successful application in analyzing visual stimulation fMRI data.

Conclusions:

• The Gaussian process-FIR model provides a robust and flexible framework for haemodynamic response modeling in fMRI.
• This approach enhances the analysis of functional neuroimages by offering improved modeling capabilities.
• The study validates the utility of advanced statistical methods in neuroimaging research.