Simulation of the posterior distribution