bm_glm function

Fit Generalized Linear Models