get_MAP function

Find the maximum a posteriori (MAP) estimate of the MCMC chain