The main function MMEst() performs (Restricted) Maximum Likelihood in a variance component mixed models using a Min-Max (MM) algorithm (Hunter, D. R., & Lange, K. (2004) <doi:10.1198/0003130042836>).
