brm_model function

Fit an MMRM.