Bayesian methods for pharmacokinetic models in dynamic contrast-enhanced magnetic resonance imaging

This paper proposes a new method for estimating kinetic parameters of dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) based on adaptive Gaussian Markov random fields. Kinetic parameter estimates using neighboring voxels reduce the observed variability in local tumor regions while preserving sharp transitions between heterogeneous tissue boundaries. Asymptotic results for standard errors from likelihood-based nonlinear regression are compared with those derived from the posterior distribution using Bayesian estimation with and without neighborhood information. Application of the method to the analysis of breast tumors based on kinetic parameters has shown that the use of Bayesian analysis combined with adaptive Gaussian Markov random fields provides improved convergence behavior and more consistent morphological and functional statistics.