Estimation of the hemodynamic response in event-related functional MRI: Bayesian networks as a framework for efficient Bayesian modeling and inference.
A convenient way to analyze blood-oxygen-level-dependent functional magnetic resonance imaging data consists of modeling the whole brain as a stationary, linear system characterized by its transfer function: the hemodynamic response function (HRF). HRF estimation, though of the greatest interest, is still under investigation, for the problem is ill-conditioned. In this paper, we recall the most general Bayesian model for HRF estimation and show how it can beneficially be translated in terms of Bayesian graphical models, leading to 1) a clear and efficient representation of all structural and functional relationships entailed by the model, and 2) a straightforward numerical scheme to approximate the joint posterior distribution, allowing for estimation of the HRF, as well as all other model parameters. We finally apply this novel technique on both simulations and real data.
Marrelec
Guillaume
G
INSERM U494, CHU Pitié-Salpêtrière, 91 boulevard de l'Hôpital, 75634 Paris Cedex 13, France.
Ciuciu
Philippe
P
Pélégrini-Issac
Mélanie
M
Benali
Habib
H
D000293:Adolescent; D000328:Adult; D000465:Algorithms; D001499:Bayes Theorem; D001921:Brain / Q000098:blood supply* / Q000502:physiology*; D001931:Brain Mapping / Q000379:methods; D002560:Cerebrovascular Circulation / Q000502:physiology; D003071:Cognition / Q000502:physiology; D003198:Computer Simulation; D005071:Evoked Potentials / Q000502:physiology*; D006801:Humans; D007090:Image Interpretation, Computer-Assisted / Q000379:methods*; D008279:Magnetic Resonance Imaging / Q000379:methods*; D008955:Models, Cardiovascular; D008959:Models, Neurological; D015233:Models, Statistical; D015203:Reproducibility of Results; D012680:Sensitivity and Specificity