fit_glm function

Fit GLM on Covariate