log_likelihood function

Calculate the negative log-likelihood for a beta binomial regression model