posterior_interval.optim_fit function

Extract posterior intervals for parameters of the model