glm_fit function

Fit design matrices to one or multiple outcomes