BBreg function

Fit a beta-binomial logistic regression model