stats_glm_logit function

Fit outcome models using stats::glm