non-exact calculations in mem_exact

[lbau7]

Hi there,

I currently use the mem_exact function a lot for simulations with a relatively small number of baskets (3-5). I noticed that the function only calculates the posterior probabilities exactly, but the posterior mean, median and HPD-intervals are still estimated from MCMC samples. This slows down computation a lot, since 100.000 samples are drawn and there is no way to change this. Drawing the posterior samples takes much longer than calculating the exact posterior probabilities for a small number of baskets.

Is there a specific reason that the rest of the output isn't also calculated exactly? Since the posterior weights are already calculated in the function, it seems that it would be much faster to also calculate the posterior mean and median exactly. It would be awesome if this could be changed!