Learning Bayesian Models With R