A differential evolution-based approach for fitting a nonlinear biophysical model to fMRI BOLD data - Archive ouverte HAL Access content directly
Journal Articles IEEE Journal of Selected Topics in Signal Processing Year : 2016

A differential evolution-based approach for fitting a nonlinear biophysical model to fMRI BOLD data

(1) , (2) , (2) , (3) , (2) , (1)
1
2
3

Abstract

Physiological and biophysical models have been proposed to link neuronal activity to the Blood Oxygen Level-Dependent (BOLD) signal in functional MRI (fMRI). Those models rely on a set of parameter values that cannot always be extracted from the literature. In some applications, interesting insight into the brain physiology or physiopathology can be gained from an estimation of the model parameters from measured BOLD signals. This estimation is challenging because there are more than 10 potentially interesting parameters involved in nonlinear equations and whose interactions may result in identifiability issues. However, the availability of statistical prior knowledge about these parameters can greatly simplify the estimation task. In this work we focus on the extended Balloon model and propose the estimation of 15 parameters using two stochastic approaches: an Evolutionary Computation global search method called Differential Evolution (DE) and a Markov Chain Monte Carlo version of DE. To combine both the ability to escape local optima and to incorporate prior knowledge, we derive the target function from Bayesian modeling. The general behavior of these algorithms is analyzed and compared with the de facto standard Expectation Maximization Gauss-Newton (EM/GN) approach, providing very promising results on challenging real and synthetic fMRI data sets involving rats with epileptic activity. These stochastic optimizers provided a better performance than EM/GN in terms of distance to the ground truth in 3 out of 6 synthetic data sets and a better signal fitting in 11 out of 12 real data sets. Non-parametric statistical tests showed the existence of statistically significant differences between the real data results obtained by DE and EM/GN. Finally, the estimates obtained from DE for these parameters seem both more realistic and more stable or at least as stable across sessions as the estimates from EM/GN.
Fichier principal
Vignette du fichier
IEEE_JSTSP.pdf (5.07 Mo) Télécharger le fichier
Vignette du fichier
fMRImodelFIT_v1.zip (4.35 Mo) Télécharger le fichier
Origin : Files produced by the author(s)
Format : Other
Origin : Files produced by the author(s)
Comment : Codes used in the paper to perform the estimation of biophysical parameters. This package also includes the 6 synthetic datasets used in the experiments.
Loading...

Dates and versions

hal-01221115 , version 1 (27-10-2015)
hal-01221115 , version 2 (08-03-2016)

Identifiers

Cite

Pablo Mesejo, Sandrine Saillet, Olivier David, Christian Bénar, Jan M. Warnking, et al.. A differential evolution-based approach for fitting a nonlinear biophysical model to fMRI BOLD data. IEEE Journal of Selected Topics in Signal Processing, 2016, 10 (2), pp.416-427. ⟨10.1109/JSTSP.2015.2502553⟩. ⟨hal-01221115v2⟩
603 View
793 Download

Altmetric

Share

Gmail Facebook Twitter LinkedIn More