Negative binomial maximum likelihood estimate implementation using L-BFGS-B