Bayesian Hierarchical Model In R