Advanced level

The binomial and Poisson distributions

The number of successes in a Bernoulli scheme is a random variable, and its distribution is called binomial. It has two parameters and two short formulas for its mean and its spread. When the trials are very many and the success very rare, that distribution turns into the Poisson one — the model of breakdowns, calls and typos.

Before you start

This topic builds on earlier ideas. Before you start, it's worth working through the lessons below — they'll make everything click:

Where this is used

Real situations where you count exactly the way this lesson teaches:

  • Overbooking a flight
    An airline sells 210 tickets for 200 seats because a passenger shows up with probability 0.95. The number who show up is binomial with E(X) = 210 · 0.95 = 199.5 and variance 210 · 0.95 · 0.05 = 9.975, that is a standard deviation of 3.16 passengers. On average half a seat stays free, but a spread of three passengers means overbooking happens on a sizeable share of flights — and that is the number the airline is knowingly buying.
  • Server failures per month
    If a server fails on average 3 times a month, the count of failures is Poisson with λ = 3. A month with no failure at all occurs with probability e to the minus three, that is 0.0498 — under 5 per cent of months. At most two failures has probability 0.0498 · (1 + 3 + 4.5) = 0.423, so more than half of all months bring at least three.
  • Typos in a book
    A proofreader estimates that 0.5 typos per page survive their work. The count on a random page is Poisson with λ = 0.5, so a clean page occurs with probability e to the minus one half, that is 0.6065. Across a 300-page book that gives about 300 · 0.3935 = 118 pages carrying at least one typo — a realistic forecast rather than pessimism.
  • Prize caps on bottles
    A bottle cap hides a prize with probability 1/50. Buying 20 bottles wins nothing with probability 0.98 to the twentieth power, that is 0.668, so the chance of at least one prize is 33 per cent. To push past one half you have to buy 35 — because 0.98 to the thirty-fifth power is already 0.493.

All formulas

  • Binomial distribution

    XB(n,p),P(X=k)=(nk)pk(1p)nkX \sim B(n, p), \qquad P(X = k) = \binom{n}{k} p^k (1-p)^{n-k}

    the formula derived in the Bernoulli lesson

  • Expected value

    E(X)=npE(X) = np

    trials times the chance of success

  • Variance

    σ2(X)=np(1p)\sigma^2(X) = np(1-p)

    largest at p = 1/2

  • At least one success

    P(X1)=1(1p)nP(X \geq 1) = 1 - (1-p)^n

    through the complement

  • Poisson distribution

    XPoiss(λ),P(X=k)=λkeλk!X \sim \text{Poiss}(\lambda), \qquad P(X = k) = \frac{\lambda^k e^{-\lambda}}{k!}

    rare events in a fixed window

  • Poisson parameters

    E(X)=σ2(X)=λE(X) = \sigma^2(X) = \lambda

    one number describes the whole distribution

  • The limiting passage

    λ=np\lambda = np

    n large, p small, the product fixed

In the Bernoulli scheme we computed the probability of exactly kk successes in nn trials. The previous lesson gave us the notion of a random variable. Put the two together: the number of successes is a random variable, and the set of its probabilities is its distribution.

The binomial distribution

The variable XX counting successes in a Bernoulli scheme is said to have a binomial distribution, written

XB(n,p),P(X=k)=(nk)pk(1p)nkX \sim B(n, p), \qquad P(X = k) = \binom{n}{k} p^k (1-p)^{n-k}

The formula is not derived a second time — it is exactly the one from before, binomial coefficient counting paths and all. What is new is the viewpoint: instead of asking about one kk, we look at the whole table at once.

−10123456789101100.10.20.3kP(X = k)p = 0.5p = 0.2
Ten trials, two different success probabilities. The wide bars are p = 0.5 — symmetric about 5. The narrow ones are p = 0.2, leaning left with a peak at 2. Both have exactly eleven bars, from k = 0 to k = 10.

The drawing shows two things at once. The peak sits near npnp: at 55 for p=0.5p = 0.5, at 22 for p=0.2p = 0.2. And symmetry is the exception rather than the rule — it belongs to p=12p = \tfrac{1}{2} alone; for smaller pp the distribution has a longer right tail.

The two parameters

Rather than compute an expected value from a table of eleven terms, we have short formulas:

E(X)=np,σ2(X)=np(1p)E(X) = np, \qquad \sigma^2(X) = np(1-p)

The first has a one-line justification with no binomial coefficient in it. A single trial is a variable taking 11 with probability pp and 00 with probability qq, so its expected value is

1p+0q=p1 \cdot p + 0 \cdot q = p

The number of successes is the sum of nn such variables, and the expected value of a sum is the sum of the expected values — hence npnp.

The variance carries an observation of its own: the factor p(1p)p(1-p) is largest at p=12p = \tfrac{1}{2}, where it equals 14\tfrac{1}{4}. An outcome is hardest to predict when success and failure are equally likely; for pp near zero or one almost everything is settled and the spread disappears.

We inspect 200 units from a production run in which 3% are faulty. How many faults do we expect to find, and with what spread?

At least once

The question asked most often about a binomial distribution is not “how many exactly” but “will it happen at all”. The answer goes through the complement and fits on one line:

P(X1)=1P(X=0)=1(1p)nP(X \geq 1) = 1 - P(X = 0) = 1 - (1-p)^n
A bottle cap hides a prize with probability 1/50. How many bottles must you buy for the chance of any prize to exceed 50%?

The Poisson distribution

Now take a situation where the trials are very many and the success very rare: millions of ad impressions and a few dozen purchases; 2.6 million seconds in a month and three failures. There nn and pp separately stop being known or meaningful — what stays meaningful is their product, the average number of events.

A variable with a Poisson distribution has one parameter λ\lambda and the distribution

P(X=k)=λkeλk!,k=0,1,2,P(X = k) = \frac{\lambda^k e^{-\lambda}}{k!}, \qquad k = 0, 1, 2, \dots
−1012345678900.10.20.3kP(X = k)
The Poisson distribution with λ = 2. The two tallest bars, at k = 1 and k = 2, are exactly equal — a property of this distribution at whole λ. To the right the tail never ends, it only shrinks fast.

The parameters here are remarkably short:

E(X)=σ2(X)=λE(X) = \sigma^2(X) = \lambda

The equality of mean and variance is a practical test of the model: if the spread in the data is clearly larger than the mean, the events are not independent and a Poisson distribution does not fit them.

The limiting passage

Where does the formula with ee come from? From the binomial distribution, in the limit nn \to \infty, p0p \to 0, with the product held fixed:

λ=np\lambda = np

You need not take that on faith — the numbers are enough.

A thousand tickets are sold and each wins with probability 0.002. Compare the exact probability of three wins with the Poisson approximation.

In practice the approximation is used for n50n \geq 50 and np10np \leq 10. This is exactly why the Poisson distribution is called the distribution of rare events: not because the events are few, but because each single opportunity for one is negligible while the opportunities are very many.

Exercises

The set asks three things and does not repeat the question about P(X=k)P(X = k) for the binomial — that is the same formula you practised in the Bernoulli lesson.

The P(X1)P(X \geq 1) prompt gives nn and pp and asks about at least one success; the answer is a fraction in lowest terms. The E(X)E(X) or σ2(X)\sigma^2(X) prompt gives the same two parameters and asks for the expected value or the variance — the answer is a plain number, because that is a quantity on the variable’s own scale rather than a probability. The P(X=k)P(X = k) prompt with a parameter λ\lambda is about the Poisson distribution: there eλe^{-\lambda} is irrational, so the answer is given rounded to four decimal places.

Practice

Work through a set of exercises — they get harder as you go. At the end you'll see your score and the mistakes worth reviewing.

Exercise 1 of 8Score: 0
P(X ≥ 1): n = 3, p = 1/2

The generator does not ask you to judge whether a Poisson approximation is allowed here — that is a verdict, not a number. The criterion is stated above: many trials, small pp, a fixed product.

Common mistakes

  • Reading npnp as the outcome rather than the averageE(X)=6E(X) = 6 does not mean there will be six faults; it means six on average, with a deviation of 2.412.41.
  • Applying the binomial distribution to dependent trials — drawing without replacement from a small set is not a Bernoulli scheme, so none of these formulas holds.
  • Dropping the factor eλe^{-\lambda} — without it the probabilities do not add up to one and can exceed one as early as k=2k = 2.
  • Substituting nn and pp into the Poisson formula instead of λ\lambda — the parameter is their product npnp, and it goes into both the power and the exponent.
  • Approximating with Poisson at large pp — the limit needs small pp; at p=0.3p = 0.3 the approximation already parts company in the second decimal.
  • Expecting the Poisson variance to be smaller than the mean — both equal λ\lambda, and a spread larger than the mean in real data is a sign that the model is wrong.

Formula card

Topic: Binomial and Poisson

  • Binomial distribution

    XB(n,p),P(X=k)=(nk)pk(1p)nkX \sim B(n, p), \qquad P(X = k) = \binom{n}{k} p^k (1-p)^{n-k}

    the formula derived in the Bernoulli lesson

  • Expected value

    E(X)=npE(X) = np

    trials times the chance of success

  • Variance

    σ2(X)=np(1p)\sigma^2(X) = np(1-p)

    largest at p = 1/2

  • At least one success

    P(X1)=1(1p)nP(X \geq 1) = 1 - (1-p)^n

    through the complement

  • Poisson distribution

    XPoiss(λ),P(X=k)=λkeλk!X \sim \text{Poiss}(\lambda), \qquad P(X = k) = \frac{\lambda^k e^{-\lambda}}{k!}

    rare events in a fixed window

  • Poisson parameters

    E(X)=σ2(X)=λE(X) = \sigma^2(X) = \lambda

    one number describes the whole distribution

  • The limiting passage

    λ=np\lambda = np

    n large, p small, the product fixed

−10123456789101100.10.20.3kP(X = k)p = 0.5p = 0.2
Two binomial distributions with the same number of trials and different p. At p = 0.5 the distribution is symmetric about 5; at p = 0.2 it leans left with its peak at 2 — because E(X) = np moves with p.
−1012345678900.10.20.3kP(X = k)
The Poisson distribution with λ = 2. The peak sits at 1 and 2, and the tail runs right without end — a count of rare events has no upper bound, only ever smaller probabilities.

Frequently asked questions

Related articles