Dear Authors,
Recently, a new package for MMRM models (used widely in clinical trials) was born. It employs the GLS estimation to obtain marginal model with lots of residual covariance structures and both KR and S. degrees of freedom.
I thought you might find it useful, just to explore a new calculation engine. It's always good to know alternatives! It's quite fast and powerful. In case you were interested, please find: https://github.com/openpharma/mmrm