The purpose of this study is to investigate a method, using simulations, to improve contrast agent quantification in Dynamic Contrast Enhanced MRI. Bayesian hierarchical models (BHMs) are applied to smaller images ($10\times10\times10$) such that spatial information can be incorporated. Then exploratory analysis is done for larger images ($64\times64\times64$) by using maximum a posteriori (MAP). For smaller images: the estimators of proposed BHMs show improvements in terms of the root mean squared error compared to the estimators in existing method for a noise level equivalent of a 12-channel head coil at 3T. Moreover, Leroux model outperforms Besag models. For larger images: MAP estimators also show improvements by assigning Leroux prior.