I have done many bayesian models using the MCMCglmm package in R, like this one:
model=MCMCglmm(lifespan ~ weight * littersize, random=~idv(DNA1