A differential evolution-based approach for fitting a nonlinear biophysical model to fMRI BOLD data - INRIA - Institut National de Recherche en Informatique et en Automatique Accéder directement au contenu
Article Dans Une Revue IEEE Journal of Selected Topics in Signal Processing Année : 2016

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

Résumé

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
fMRImodelFIT_v1.zip (4.35 Mo) Télécharger le fichier
Origine : Fichiers produits par l'(les) auteur(s)
Format : Autre
Origine : Fichiers produits par l'(les) auteur(s)
Commentaire : 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 et versions

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

Identifiants

Citer

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⟩
636 Consultations
841 Téléchargements

Altmetric

Partager

Gmail Facebook X LinkedIn More