The validity of pharmacokinetic parameters estimated by the maximum a posteriori probability (MAP) Bayesian method was investigated by simulation studies. A 1-compartment model with bolus intravenous administration was used as a pharmacokinetic model, and the coefficients of variation for the parameters and residual error were set at 30% and 10%, respectively. The accuracy of the posterior modes of pharmacokinetic parameters estimated by the MAP Bayesian method was assessed by the difference between the true value and the estimated value. The results showed that the accuracy of the Bayesian estimation depended on sampling times and on the differences between the prior means and individual true parameter values. For assessing the reliability and accuracy of the Bayesian estimation, the authors suggest using the whole posterior distribution of the pharmacokinetic parameters to describe the 95th percentile range for predicted blood concentration profiles. The authors believe that the proposed procedures provide helpful information for evaluating the Bayesian estimation of pharmacokinetic profiles.