Hello -    Thank you for this suggestion. I have the NLMIXED model working for a beta-binomial  with random effects (intercept) but for some datasets where the N is large (>150) it is failing to compute I think due to the factorial and beta functions in the beta-binomial distribution.    The non-random effects Beta-binomial models works in PROC FMM but not in PROC NLMIXED (no random effects) when counts (N) are large. I have tried different optimization methods and nothing works.    Any suggestions ?   ... View more