get_mmrm_sample function

Fit MMRM and returns parameter estimates