Are there any packages such as brms or blme?
They allow to use bayesian computation (Stan) through a simple formula.
blme: Maximum a posteriori estimation for linear and generalized linear mixed-effects models in a Bayesian setting.
brms: Provides an interface to fit Bayesian generalized (non-)linear multilevel models using Stan.
A wide range of distributions and link functions are supported, allowing users to fit linear, robust linear, count data, survival, response times, ordinal, zero-inflated, hurdle, and even self-defined mixture models all in a multilevel context. Auto-correlation structures, censored data, meta-analytic standard errors…