fit_glm_method function

Perform GLM (Generalised Linear Model) fitting