Hierarchical Bayesian Approach In R