Is ot possible to do Bayesian model fitting with raw R code?