brm_formula function

Model formula