HLfit function

Fit mixed models with given correlation matrix