Data-driven HRF estimation for encoding and decoding models.

Area of Science:

• Neuroimaging
• Computational Neuroscience
• Biostatistics

Background:

• The canonical hemodynamic response function (HRF) is widely used in BOLD fMRI but shows variability across brain regions and individuals.
• Data-driven HRF estimation can increase statistical power but often leads to unstable results due to high dimensionality.
• Existing methods struggle with accurate and efficient HRF modeling in fMRI data analysis.

Purpose of the Study:

• To develop a robust and efficient method for joint estimation of activation and HRF in BOLD fMRI data.
• To improve the statistical power and accuracy of fMRI data modeling by incorporating data-driven HRF estimation.
• To introduce the Rank-1 GLM (R1-GLM) as a superior alternative to conventional HRF modeling techniques.

Main Methods:

• Developed a joint estimation framework for activation and HRF using a rank constraint, ensuring HRF consistency across conditions while allowing voxel-wise differences.
• Implemented an efficient quasi-Newton optimization method for model estimation, leveraging fast gradient computations.
• Extended the R1-GLM to GLM with separate designs, enhancing brain activity decoding accuracy.

Main Results:

• The R1-GLM model demonstrated superior performance compared to 10 other HRF modeling methods on two independent datasets.
• Achieved higher scores in both encoding and decoding settings, indicating improved accuracy in modeling BOLD fMRI data.
• The proposed method offers a balance of accuracy and computational efficiency for fMRI analysis.

Conclusions:

• The R1-GLM provides a more accurate and statistically powerful approach to HRF modeling in BOLD fMRI.
• This method addresses the limitations of canonical HRF assumptions and unstable unconstrained estimation.
• R1-GLM represents a significant advancement for both encoding and decoding applications in neuroimaging research.