hanova1 function

Posterior sampling for a 1-way hierarchical ANOVA