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

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.
Type de document :
Article dans une revue
IEEE Journal of Selected Topics in Signal Processing, IEEE, 2016, 10 (2), pp.416-427. <10.1109/JSTSP.2015.2502553>
Liste complète des métadonnées


https://hal.inria.fr/hal-01221115
Contributeur : Pablo Mesejo Santiago <>
Soumis le : mardi 8 mars 2016 - 17:00:59
Dernière modification le : vendredi 11 mars 2016 - 00:06:04
Document(s) archivé(s) le : dimanche 13 novembre 2016 - 11:54:04

Fichiers

IEEE_JSTSP.pdf
Fichiers produits par l'(les) auteur(s)

Identifiants

Collections

Citation

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, IEEE, 2016, 10 (2), pp.416-427. <10.1109/JSTSP.2015.2502553>. <hal-01221115v2>

Partager

Métriques

Consultations de
la notice

460

Téléchargements du document

456