I'd start by identifying a good example to generate simulated data on, and then use rstanarm::stan_glm as backend