binomial likelihood function in r