fit_nb function

Fit a negative binomial to a vector of observations and expectations