posterior.mode function

Estimates the marginal parameter modes using kernel density estimation