Fitting distributions the Bayesian way