Extract point estimates of compositions from fitted model.
rmspe_calc(par, n_bins, n_draws, target)
The parameter (standard deviation) to be searched over to find a Dirichlet equivalent
Bins for the Dirichlet distribution
Numbers of samples to use for doing calculation
The goal of the specified prior, e.g. 1 or 1/n_bins