When the *r* parameter is an integer, the negative
binomial pdf is

$$y=f(x|r,p)=\left(\begin{array}{c}r+x-1\\ x\end{array}\right){p}^{r}{q}^{x}{I}_{(0,1,\mathrm{...})}(x)$$

where *q* = 1 – *p*.
When *r* is not an integer, the binomial coefficient
in the definition of the pdf is replaced by the equivalent expression

$$\frac{\Gamma (r+x)}{\Gamma (r)\Gamma (x+1)}$$

In its simplest form (when *r* is an integer),
the negative binomial distribution models the number of failures *x* before
a specified number of successes is reached in a series of independent,
identical trials. Its parameters are the probability of success in
a single trial, *p*, and the number of successes, *r*.
A special case of the negative binomial distribution, when *r* =
1, is the geometric
distribution, which models the number of failures before the
first success.

More generally, *r* can take on non-integer
values. This form of the negative binomial distribution has no interpretation
in terms of repeated trials, but, like the Poisson distribution, it is useful
in modeling count data. The negative binomial distribution is more
general than the Poisson distribution because it has a variance that
is greater than its mean, making it suitable for count data that do
not meet the assumptions of the Poisson distribution. In the limit,
as *r* increases to infinity, the negative binomial
distribution approaches the Poisson distribution.

Suppose you are collecting data on the number of auto accidents on a busy highway, and would like to be able to model the number of accidents per day. Because these are count data, and because there are a very large number of cars and a small probability of an accident for any specific car, you might think to use the Poisson distribution. However, the probability of having an accident is likely to vary from day to day as the weather and amount of traffic change, and so the assumptions needed for the Poisson distribution are not met. In particular, the variance of this type of count data sometimes exceeds the mean by a large amount. The data below exhibit this effect: most days have few or no accidents, and a few days have a large number.

accident = [2 3 4 2 3 1 12 8 14 31 23 1 10 7 0]; m = mean(accident)

m = 8.0667

v = var(accident)

v = 79.3524

The negative binomial distribution is more general than the
Poisson, and is often suitable for count data when the Poisson is
not. The function `nbinfit`

returns
the maximum likelihood estimates (MLEs) and confidence intervals for
the parameters of the negative binomial distribution. Here are the
results from fitting the `accident`

data:

[phat,pci] = nbinfit(accident)

`phat = `*1×2*
1.0060 0.1109

`pci = `*2×2*
0.2152 0.0171
1.7968 0.2046

It is difficult to give a physical interpretation in this case to the individual parameters. However, the estimated parameters can be used in a model for the number of daily accidents. For example, a plot of the estimated cumulative probability function shows that while there is an estimated 10% chance of no accidents on a given day, there is also about a 10% chance that there will be 20 or more accidents.

plot(0:50,nbincdf(0:50,phat(1),phat(2)),'.-'); xlabel('Accidents per Day') ylabel('Cumulative Probability')

Compute and plot the pdf using four different values for the parameter `r`

, the desired number of successes: `.1`

, `1`

, `3`

, and `6`

. In each case, the probability of success `p`

is `.5`

.

x = 0:10; plot(x,nbinpdf(x,.1,.5),'s-', ... x,nbinpdf(x,1,.5),'o-', ... x,nbinpdf(x,3,.5),'d-', ... x,nbinpdf(x,6,.5),'^-'); legend({'r = .1' 'r = 1' 'r = 3' 'r = 6'}) xlabel('x') ylabel('f(x|r,p)')

The plot shows that the negative binomial distribution can take on a variety of shapes, ranging from very skewed to nearly symmetric, depending on the value of `r`

.