re_estimate function

Re-estimate the posterior distributions with different burn-in