corrHLfit function

Fits a mixed model, typically a spatial GLMM.