ml_glm function

A function to fit generalized linear models using maximum likelihood.