Dans cet article, nous synthétisons la méthodologie usuelle et des variantes plus élaborées pour analyser des données individuelles d’Imagerie par Résonance Magnétique fonctionnelle (IRMf). De telles données sont acquises au cours d’expériences sensorielles ou cognitives dont le but est de stimuler le sujet dans le scanner afin d’évoquer une activité cérébrale typique. À partir de données quadri-dimensionnelles (trois dans l’espace et le temps), le but est double : d’abord, détecter les régions impliquées dans les processus sensoriels ou cognitifs que le protocole expérimental manipule ; ensuite, estimer la dynamique cérébrale sous-jacente. La première question est généralement abordée dans le contexte du Modèle Linéaire Général (MLG) qui suppose a priori connue la réponse impulsionnelle du filtre hémodynamique. La seconde question traite de l’estimation de la forme de cette réponse qui fait sens uniquement dans les régions activées. Durant ces cinq dernières années, un nouveau cadre de détection-estimation conjointe traitant ces deux questions simultanément a été proposé successivement dans [ 59 , 60 , 102 ]. Nous montrons ici jusqu’à quel point cette méthodologie étend les approches à base de MLG, les améliore en termes de sensibilité et spécificité statistique, quelles questions supplémentaires ce cadre permet d’investiguer théoriquement et comment il fournit un moyen aux utilisateurs de traiter de façon adaptée leur données non-lissées spatialement, aussi bien dans le répère 3D de l’acquisition que sur la surface corticale.
In this paper, we review classical and advanced methodologies for analysing within-subject functional Magnetic Resonance Imaging (fMRI) data. Such data are usually acquired during sensory or cognitive experiments that aims at stimulating the subject in the scanner and eliciting evoked brain activity. From such four-dimensional datasets (three in space, one in time), the goal is twofold: first, detecting brain regions involved in the sensory or cognitives processes that the experimental design manipulates; second, estimating the underlying activation dynamics. The first issue is usually addressed in the context of the General Linear Model (GLM), which a priori assumes a X form for the impulse response of the hemodynamic filter. The second question aims at estimating this shape which makes sense in activating regions only. In the last five years, a novel Joint Detection-Estimation (JDE) framework addressing these two questions simultaneously has been proposed in [ 59 , 60 , 102 ]. We show to which extent this methodology outperforms the GLM approach in terms of statistical sensitivity and specificity, which additional questions it allows us to investigate theoretically and how it provides us with a well-adapted framework to treat spatially unsmoothed real fMRI data both in the 3D acquisition volume and on the cortical surface.
Mot clés : IRMf, neuro-imagerie, hémodynamique non-paramétrique, inférence bayésienne, MCMC, sélection de modèles, champs de Markov, fonction de partition, champs de Pots
@article{JSFS_2010__151_1_58_0, author = {Ciuciu, Philippe and Vincent, Thomas and Risser, Laurent and Donnet, Sophie}, title = {A joint detection-estimation framework for analysing within-subject {fMRI} data}, journal = {Journal de la soci\'et\'e fran\c{c}aise de statistique}, pages = {58--89}, publisher = {Soci\'et\'e fran\c{c}aise de statistique}, volume = {151}, number = {1}, year = {2010}, mrnumber = {2679310}, zbl = {1316.62039}, language = {en}, url = {http://www.numdam.org/item/JSFS_2010__151_1_58_0/} }
TY - JOUR AU - Ciuciu, Philippe AU - Vincent, Thomas AU - Risser, Laurent AU - Donnet, Sophie TI - A joint detection-estimation framework for analysing within-subject fMRI data JO - Journal de la société française de statistique PY - 2010 SP - 58 EP - 89 VL - 151 IS - 1 PB - Société française de statistique UR - http://www.numdam.org/item/JSFS_2010__151_1_58_0/ LA - en ID - JSFS_2010__151_1_58_0 ER -
%0 Journal Article %A Ciuciu, Philippe %A Vincent, Thomas %A Risser, Laurent %A Donnet, Sophie %T A joint detection-estimation framework for analysing within-subject fMRI data %J Journal de la société française de statistique %D 2010 %P 58-89 %V 151 %N 1 %I Société française de statistique %U http://www.numdam.org/item/JSFS_2010__151_1_58_0/ %G en %F JSFS_2010__151_1_58_0
Ciuciu, Philippe; Vincent, Thomas; Risser, Laurent; Donnet, Sophie. A joint detection-estimation framework for analysing within-subject fMRI data. Journal de la société française de statistique, Tome 151 (2010) no. 1, pp. 58-89. http://www.numdam.org/item/JSFS_2010__151_1_58_0/
[1] The variability of human BOLD hemodynamic responses, Neuroimage, Volume 7 (1998), 574 pages
[2] Variational algorithms for approximate Bayesian inference, University College of London, London, United Kingdom, May (2003) (PhD thesis)
[3] Estimation of the Hemodynamic Response Function in Event-Related Functional MRI: Directed Acyclic Graphs for a General Bayesian Inference Framework, Proc. 18th IPMI (Taylor, Chris; Noble, J. Alison, eds.) (LNCS-2732), Springer Verlag, Ambleside, United Kingdom (2003), pp. 635-646
[4] Unsupervised robust non-parametric estimation of the hemodynamic response function for any fMRI experiment, IEEE Trans. Med. Imag., Volume 22 (2003) no. 10, pp. 1235-1251
[5] Estimation of the Hemodynamic Response Function in Event-Related Functional MRI: Bayesian networks as a framework for efficient Bayesian modeling and inference, IEEE Trans. Med. Imag., Volume 23 (2004) no. 8, pp. 959-967
[6] An ARX model-based approach to trial by trial identification of fMRI-BOLD responses, Neuroimage, Volume 37 (2007), pp. 189-201
[7] Modeling the hemodynamic response to brain activation, Neuroimage, Volume 23, Supplement 1 (2004), p. S220-S233
[8] Fast approximate energy minimization via graph cuts, IEEE Trans. Pattern Anal. Mach. Intell., Volume 23 (2001) no. 11, pp. 1222-1239
[9] Dynamics of blood flow and oxygenation changes during Brain activation: the balloon model, Magn. Reson. Med., Volume 39 (1998), pp. 855-864
[10] Generalizing Swendsen-Wang to sampling arbitrary posterior probabilities, IEEE Trans. Pattern Anal. Mach. Intell., Volume 27 (2005) no. 8, pp. 1239-1253
[11] Deterministic edge-preserving regularization in computed imaging, IEEE Trans. Image Processing, Volume 6 (1997) no. 2, pp. 298-311
[12] A new cortical surface parcellation model and its automatic implementation, Proc. 9th MICCAI (LNCS 4191), Springer Verlag, Copenhaegen, Denmark (2006), pp. 193-200
[13] Marginal likelihood from the Gibbs output, J. Amer. Statist. Assoc., Volume 90 (1995), pp. 1313-1321 | MR | Zbl
[14] Bilinear extrapolation scheme For Fast Estimation of 3D Ising Field Partition Function. Application to fMRI Time Course Analysis., 16th Proc. IEEE ICIP, Cairo, Egypt (2009), pp. 833-836
[15] Schéma d’extrapolation de fonctions de partition de champs de Potts. Application à l’analyse d’images en IRMf, Actes du 22 colloque GRETSI, Dijon, France (2009)
[16] Spatially adaptive mixture modeling for analysis of within-subject fMRI time series, IEEE Trans. Med. Imag., Volume 29 (2010) no. 4, pp. 1059-1074
[17] Marginal likelihood from the Metropolis-Hastings output, J. Amer. Statist. Assoc., Volume 96 (2001) no. 453, pp. 270-281 | MR | Zbl
[18] Parametric analysis of XMRI data using linear systems methods, Neuroimage, Volume 6 (1997), pp. 93-103
[19] The impact of temporal regularization on estimates of the BOLD hemodynamic response function: a comparative analysis., Neuroimage, Volume 40 (2008) no. 4, pp. 1606-1618 | DOI
[20] Single-trial variability in event-related BOLD signals, Neuroimage, Volume 15 (2002) no. 4, p. 823-35
[21] Spatio-Temporal fMRI Analysis Using Markov Random Fields, IEEE Trans. Med. Imag., Volume 17 (1998) no. 6, pp. 1028-1039
[22] Bayesian model comparison and selection in neuroimaging from JDE results (2010) (Research Rep.)
[23] Mixture model mapping of brain activation in functional magnetic resonance images, Hum. Brain Mapp., Volume 7 (1999), pp. 1-14
[24] Modelling spatially correlated data via mixtures: a Bayesian approach, J. R. Statist. Soc. B, Volume 64 (2002) no. 4, pp. 805-826 | MR | Zbl
[25] Dynamic causal modelling., Neuroimage, Volume 19 (2003) no. 4, p. 1273-302
[26] Variational free energy and the Laplace approximation, Neuroimage, Volume 33 (2007), pp. 220-234
[27] Nonlinear responses in fMRI: the balloon model, Volterra kernels, and other hemodynamics, Neuroimage, Volume 12 (2000), pp. 466-477
[28] Hidden Markov Random Field Model Selection Criteria Based on Mean Field-Like Approximations, IEEE Trans. Pattern Anal. Mach. Intell., Volume 25 (2003) no. 9, pp. 1089-1101
[29] Statistical Parametric Mapping, Functional Neuroimaging : Technical Foundations (Thatcher, R.W.; Hallet, M.; Zeffiro, T.; John, E.R.; Huerta, M., eds.) (1994), pp. 79-93
[30] Bayesian spatio-temporal modeling of the hemodynamic response function in BOLD fMRI, Biometrics, Volume 57 (2001), pp. 554-562 | MR | Zbl
[31] Bayesian model choice: asymptotic and exact calculations, J. R. Statist. Soc. B, Volume 56 (1994), pp. 501-514 | MR | Zbl
[32] A Bayesian time-course model for functional magnetic resonance imaging data (with discussion), J. Amer. Statist. Assoc., Volume 95 (2000), pp. 691-719
[33] Bayesian fMRI data analysis with sparse spatial basis function priors, Neuroimage, Volume 34 (2007) no. 3, pp. 1108-1125
[34] Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images, IEEE Trans. Pattern Anal. Mach. Intell., Volume PAMI-6 (1984) no. 6, pp. 721-741 | Zbl
[35] Dissociating Neural Correlates of Cognitive Components in Mental Calculation, Cereb. Cortex, Volume 11 (2001) no. 4, pp. 350-359 | DOI
[36] Deconvolution of impulse response in event-related BOLD fMRI, Neuroimage, Volume 9 (1999), pp. 416-429
[37] Simulating normalizing constants: from importance sampling to bridge sampling to path sampling, Statistical Science, Volume 13 (1998), pp. 163-185 | MR | Zbl
[38] Modeling the haemodynamic response in fXMRI using smooth filters, IEEE Trans. Med. Imag., Volume 19 (2000) no. 12, pp. 1188-1201
[39] Hidden Markov models and desease mapping, J. Amer. Statist. Assoc., Volume 97 (2002) no. 460, pp. 1-16 | MR | Zbl
[40] Constrained Restoration and the Recovery of Discontinuities, IEEE Trans. Pattern Anal. Mach. Intell., Volume 14 (1992) no. 3, pp. 367-383
[41] Reversible jump MCMC computation and Bayesian model determination, Biometrika, Volume 82 (1995), pp. 711-732 | MR | Zbl
[42] Markov Chain Monte Carlo in practice, Chapman and Hall, London, United Kingdom, 1996 | MR | Zbl
[43] Nonlinear Image Recovery with Half-Quadratic Regularization, IEEE Trans. Image Processing, Volume 4 (1995) no. 7, pp. 932-946
[44] A Bayesian approach to significance testing for statistical images from PET, Quantification of Brain Function, Tracer Kinematics in Image Analysis in Brain PET (Excerpta Medica International Congress Series (1030)), Elsevier Science, New York, USA (1993), pp. 521-531
[45] Auxiliary variable methods for Markov chain Monte Carlo with applications, J. Amer. Statist. Assoc., Volume 93 (1998) no. 442, pp. 585-595 | Zbl
[46] Variation of BOLD hemodynamic responses across subjects and brain regions and their effects on statistical analyses, Neuroimage, Volume 21 (2004), pp. 1639-1651
[47] Diffusion-based spatial priors for functional magnetic resonance images., Neuroimage, Volume 41 (2008) no. 2, pp. 408-423 | DOI
[48] Detecting latency differences in event-related BOLD responses: application to words versus nonwords and initial versus repeated face presentations, Neuroimage, Volume 15 (2002) no. 1, pp. 83-97
[49] The choice of basis function in event-related fMRI, Volume 13 (2001), 149 pages
[50] A Half-Quadratic block-coordinate descent method for spectral estimation, Signal Processing, Volume 82 (2002) no. 7, pp. 941-959 | Zbl
[51] Bayesian approach to Inverse problems, ISTE Ltd and John Wiley & Sons Inc, 2008 | MR | Zbl
[52] Robust extrapolation scheme for fast estimation of 3D Ising field partition functions. Application to within-subject fMRI data analysis., 12thProc. MICCAI’09 (Yang, G.-Z., ed.) (LNCS 5761), Springer Verlag Berlin Heidelberg, London, UK (2009), pp. 975-983
[53] Spatially unsupervised analysis of within-subject fMRI data using multiple extrapolations of 3D Ising field partition functions, 2009 IEEE international workshop on Machine Learning for Signal Processing, Grenoble, France (2009)
[54] Nonlinear estimation of the BOLD signal, Magn. Reson. Med., Volume 40 (2008) no. 2, pp. 504-514
[55] Compound Gauss-Markov random fields for image estimation, IEEE Trans. Signal Processing, Volume 39 (1991) no. 3, pp. 683-697
[56] Variational Bayesian inversion of the equivalent current dipole model in EEG/MEG, Neuroimage, Volume 39 (2008), pp. 728-741
[57] Bayes factors, J. Amer. Statist. Assoc., Volume 90 (1995), pp. 773-795 | MR | Zbl
[58] Modeling the hemodynamic response in single-trial functional MRI experiments, Magn. Reson. Med., Volume 42 (1999), pp. 787-797
[59] Empirical and Substantive Models, the Bayesian Paradigm, and Meta-Analysis in Functional Brain Imaging, Hum. Brain Mapp., Volume 5 (1997), pp. 259-263
[60] Monte Carlo strategies in scientific computing, Springer series in Statistics, Springer-Verlag, New-York, 2001 | MR
[61] Min-max extrapolation scheme for fast estimation of 3D Potts field partition functions. Application to the joint detection-estimation of brain activity in fMRI., Journal of Signal Processing Systems, Volume 65 (2010), pp. 325-338 | DOI
[62] Bayesian deconvolution fMRI data using bilinear dynamical systems, Neuroimage, Volume 42 (2008) no. 4, pp. 1381-1396 | DOI
[63] A dynamic causal model study of neuronal population dynamics, Neuroimage, Volume 51 (2010) no. 1, pp. 91-101
[64] Characterizing the hemodynamic response: effects of presentation rate, sampling procedure, and the possibility of ordering brain activity based on relative timing, Neuroimage, Volume 11 (2000), pp. 735-759
[65] Simulating ratios of normalizing constants via a simple identity: a theoretical exploration, Statistica Sinica, Volume 6 (1996), pp. 831-860 | Zbl
[66] Comparison of two convolution models for fMRI time series, Neuroimage, Volume 5 (1997)
[67] Approximate Bayesian inference by the weighted likelihood bootstrap (with discussion), J. R. Statist. Soc. B, Volume 56 (1994), pp. 3-48 | Zbl
[68] Anatomically informed convolution kernels for the projection of fMRI data on the cortical surface, Proc. 9th MICCAI (LNCS 4191), Springer Verlag, Copenhaegen, Denmark (2006), pp. 300-307
[69] From Spatial Regularization to Anatomical Priors in fMRI Analysis, IPMI, Glenwood Springs, Colorado (2005)
[70] Combining spatial priors and anatomical information for fMRI detection, Medical Image Analysis, Volume 14 (2010) no. 3, pp. 318-331
[71] Brain Magnetic Resonance Imaging with contrast dependent on blood oxygenation, Proc. Natl. Acad. Sci. USA, Volume 87 (1990) no. 24, pp. 9868-9872
[72] A two-dimensional model with an order-disorder transition, Phys. Rev., Volume 65 (1944) no. 3& 4, pp. 117-149 | Zbl
[73] Modelling the neurovascular habituation effect on fMRI time series, 34th Proc. IEEE ICASSP, Taipei, Taiwan (2009), pp. 433-436
[74] Improved fMRI group studies based on spatially varying non-parametric BOLD signal modeling, 5th Proc. IEEE ISBI, Paris, France (2008), pp. 1263-1266
[75] Bayesian Comparison of Spatially Regularised General Linear Models, Human Brain Mapping, Volume 28 (2007) no. 4, pp. 275-293 | DOI
[76] Bilinear dynamical systems, Philos Trans R Soc Lond B Biol Sci, Volume 360 (2005) no. 1457, pp. 983-993 | DOI
[77] Improved detection sensitivity of functional MRI data using a brain parcellation technique, Proc. 5th MICCAI (LNCS 2488 (Part I)), Springer Verlag, Tokyo, Japan (2002), pp. 467-474 | Zbl
[78] Robust Bayesian estimation of the hemodynamic response function in event-related BOLD MRI using basic physiological information, Hum. Brain Mapp., Volume 19 (2003) no. 1, pp. 1-17
[79] Joint detection-estimation of brain activity in functional MRI: a multichannel deconvolution solution, IEEE Trans. Signal Processing, Volume 53 (2005) no. 9, pp. 3488-3502 | Zbl
[80] Dealing with the shortcomings of spatial normalization: Multi-subject parcellation of fMRI datasets, Hum. Brain Mapp., Volume 27 (2006) no. 8, pp. 678-693
[81] Are fMRI event-related response constant in time? A model selection answer., Neuroimage (2006), pp. 1169-1176
[82] Probabilistic anatomo-functional parcellation of the cortex: how many regions?, 11thProc. MICCAI, LNCS Springer Verlag, New-York, USA (2008)
[83] Comparing families of dynamic causal models, PLoS Comput Biol, Volume 6 (2010) no. 3
[84] Bayesian fMRI time series analysis with spatial priors, Neuroimage, Volume 23 (2005) no. 2, pp. 350-362
[85] fMRI activation maps based on the NN-ARx model, Neuroimage, Volume 23 (2004), pp. 680-697
[86] Multivariate spatial Gaussian mixture modeling for statistical clustering of hemodynamic parameters in functional MRI, 34th Proc. IEEE ICASSP, Taipei, Taiwan (2009), pp. 445-448
[87] Modeling Hemodynamic Response for Analysis of Functional XMRI Time-Series, Hum. Brain Mapp., Volume 6 (1998), pp. 283-300
[88] Estimating the integrated likelihood via posterior simulation using the harmonic mean identity, Bayesian statistics 8 (Bernardo, J.M.; Bayarri, M.J.; Berger, O.; David, A.P.; Heckermann, D.; Smith, A.F.M.; West, M., eds.), Oxford University Press (2007), pp. 1-45 | Zbl
[89] Bayesian approach to segmentation of statistical parametric maps, IEEE Trans. Biomed. Eng., Volume 48 (2001), pp. 1186-1194
[90] Exploratory analysis of brain connectivity with ICA, IEEE Eng Med Biol Mag, Volume 25 (2006) no. 2, pp. 102-111
[91] A state-space model of the hemodynamic approach: nonlinear filtering of BOLD signal, Neuroimage, Volume 21 (2004), pp. 547-567
[92] Negative functional MRI response correlates with decreases in neuronal activity in monkey visual area V1., Nat Neurosci, Volume 9 (2006) no. 4, pp. 569-577
[93] Dynamic causal models of neural system dynamics:current state and future extensions, J Biosci, Volume 32 (2007) no. 1, pp. 129-144
[94] Biophysical models of fMRI responses, Curr Opin Neurobiol, Volume 14 (2004) no. 5, pp. 629-635 | DOI
[95] Bayesian blind separation of generalized hyperbolic processes in noisy and underdeterminate mixtures, IEEE Trans. Signal Processing, Volume 54 (2006) no. 9, pp. 3257-3269 | Zbl
[96] Probabilistic modeling of single-trial fMRI data, IEEE Trans. Med. Imag., Volume 19 (2000), pp. 19-35
[97] A fully Bayesian approach to the parcel-based detection-estimation of brain activity in fMRI, Neuroimage, Volume 41 (2008) no. 3, pp. 941-969 | DOI
[98] Assessing brain activity through spatial Bayesian variable selection, Neuroimage, Volume 20 (2003), pp. 802-815
[99] Ten simple rules for dynamic causal modeling, Neuroimage, Volume 49 (2010) no. 4, pp. 3099-3109 | DOI
[100] Estimation of Binary Markov Random Fields Using Markov Chain Monte Carlo, J. Comput. and Graph. Stats., Volume 15 (2006) no. 1, pp. 207-227
[101] Modelling the role of excitatory and inhibitory neuronal activity in the generation of the BOLD signal, Neuroimage, Volume 35 (2007), pp. 149-165
[102] Identification and comparison of stochastic metabolic/hemodynamic models (sMHM) for the generation of the BOLD signal, J. Comput. NeuroSci., Volume 26 (2009), pp. 251-269
[103] Spatial mixture modeling of fMRI data, Hum. Brain Mapp., Volume 11 (2000) no. 4, pp. 233-248
[104] Variational Bayes inference of spatial mixture models for segmentation, IEEE Trans. Med. Imag., Volume 25 (2006) no. 10, pp. 1380-1391
[105] Mixture models with adaptive spatial regularization for segmentation with an application to fMRI data, IEEE Trans. Med. Imag., Volume 24 (2005) no. 1, pp. 1-11
[106] Fully Bayesian spatio-temporal modelling of fMRI data, IEEE Trans. Med. Imag., Volume 23 (2004) no. 2, pp. 213-231
[107] Constrained linear basis set for HRF modelling using variational Bayes, Neuroimage, Volume 21 (2004) no. 4, pp. 1748-1761
[108] A general statistical analysis for fMRI data, Neuroimage, Volume 15 (2002) no. 1, pp. 1-15