I am fitting a 5 parameters model to some data using Maximum Likelihood and Non Informative Bayesian Inference using Metropolis-Hastings algorithm.
From the maximum likelihood fit, i calculated the expected value of the parameters and their 75% Wald confidence interval and the 75% likelihood ratio based confidence interval.
In the Bayesian framework, I used the Metropolis Hastings algorithm and I calculated the expected value for each chain (which is almost coincident with maximum likelihood estimate. To be honest, I would like to compare the uncertainties in this estimation with the uncertainties in Maximum Likelihood. So what I did is to calculate the 0.125 and 0.875 quantile of the stabilized part of the Markov Chain. The obtained values sometimes differ from the likelihood ratio confidence intervals, especially when the Wald CI and the likelihood ratio CI are not coincident. .
Is the procedure of evaluating quantiles of the converged Markov Chain correct in order to estimate Credible Intervals?