cnorm.betabinomial1 function

Fit a beta binomial regression model