“…The convergence property of the sampling chains will be checked by the Gelman-Rubin index, which is the degree of approximating 1 [47]. Applications can be found in [41] and [48]. When it is converged, a fixed number of samples can be generated from the posterior functions of the parameter vector after a burn-in period (e.g., the first 1000 samples), i.e., θ Mentioned that parameter of the priors in (26) are calibrated using MLE estimates based on actual data, and then updated on the basis of the same data used to formulate the likelihood in (27), which is not fully Bayesian approach.…”