glm_grouped function

Fit a grouped binomial outcome from a data set with case weights