fitNL function

Fit Negative Binomial Mixture Model