Fit a linear mixed model