Fitting a GLM model