log_likelihood2 function

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