fit_nbinom
