Robust and efficient estimation of the parameters of the negative binomial distribution. Two-step estimation procedure, with computation of an initial estimate which is used to identify and downweight outliers. Then a weighted maximum likelihood estimate is computed. The initial estimator is a minimum disparity estimator.