Fully Bayesian inference for structural MRI: application to segmentation and statistical analysis of T2-hypointensities.

PLoS ONE. 2013;8(7):e68196 DOI 10.1371/journal.pone.0068196


Journal Homepage

Journal Title: PLoS ONE

ISSN: 1932-6203 (Online)

Publisher: Public Library of Science (PLoS)

LCC Subject Category: Medicine | Science

Country of publisher: United States

Language of fulltext: English

Full-text formats available: PDF, HTML, XML



Paul Schmidt
Volker J Schmid
Christian Gaser
Dorothea Buck
Susanne Bührlen
Annette Förschler
Mark Mühlau


Peer review

Editorial Board

Instructions for authors

Time From Submission to Publication: 24 weeks


Abstract | Full Text

Aiming at iron-related T2-hypointensity, which is related to normal aging and neurodegenerative processes, we here present two practicable approaches, based on Bayesian inference, for preprocessing and statistical analysis of a complex set of structural MRI data. In particular, Markov Chain Monte Carlo methods were used to simulate posterior distributions. First, we rendered a segmentation algorithm that uses outlier detection based on model checking techniques within a Bayesian mixture model. Second, we rendered an analytical tool comprising a Bayesian regression model with smoothness priors (in the form of Gaussian Markov random fields) mitigating the necessity to smooth data prior to statistical analysis. For validation, we used simulated data and MRI data of 27 healthy controls (age: [Formula: see text]; range, [Formula: see text]). We first observed robust segmentation of both simulated T2-hypointensities and gray-matter regions known to be T2-hypointense. Second, simulated data and images of segmented T2-hypointensity were analyzed. We found not only robust identification of simulated effects but also a biologically plausible age-related increase of T2-hypointensity primarily within the dentate nucleus but also within the globus pallidus, substantia nigra, and red nucleus. Our results indicate that fully Bayesian inference can successfully be applied for preprocessing and statistical analysis of structural MRI data.