LMMsolve function

Solve Linear Mixed Models