Bayesian inference in R