Binomial to Poisson Approximation

A2 · S2 · 11 min

When the number of trials is large and the probability of success is small, binomial probabilities are tedious to calculate but almost identical to Poisson probabilities with the same mean. The syllabus expects you to recognise when this approximation is appropriate, to justify it, and to use it. On Paper 6 it usually appears as "use a suitable approximation", and choosing the right approximation (Poisson or normal) is part of the mark.

Why the approximation works

The Poisson distribution was built as the limit of binomial distributions in which nn grows and pp shrinks while the mean npnp stays fixed. So a binomial with a large nn and a small pp is already "most of the way" to its Poisson limit.

Compare the two sets of moments:

B(n,p)B(n, p)Po(np)\text{Po}(np)
Meannpnpnpnp
Variancenp(1−p)np(1 - p)npnp

The means agree exactly. The variances differ by the factor 1−p1 - p, which is close to 11 when pp is small. Here is how close the probabilities are for B(100,0.03)B(100, 0.03) and Po(3)\text{Po}(3):

rr0011223344556677
B(100,0.03)B(100, 0.03)0.04760.04760.14710.14710.22520.22520.22750.22750.17060.17060.10130.10130.04960.04960.02060.0206
Po(3)\text{Po}(3)0.04980.04980.14940.14940.22400.22400.22400.22400.16800.16800.10080.10080.05040.05040.02160.0216

Agreement to about two decimal places, with a far simpler formula.

The conditions

Key result

If X∼B(n,p)X \sim B(n, p) with nn large and pp small, then approximately

X∼Po(np).X \sim \text{Po}(np).

The syllabus conditions are

n>50andnp<5(approximately).n > 50 \quad\text{and}\quad np < 5 \qquad \text{(approximately)}.

Both conditions matter. A large nn alone is not enough: B(200,0.4)B(200, 0.4) has np=80np = 80 and is much better approximated by a normal distribution. A small pp alone is not enough either: B(10,0.05)B(10, 0.05) is easy to compute exactly, and nn is too small for the approximation to be close.

Choosing between the Poisson and normal approximations

SituationApproximationConditions
nn large, pp smallPo(np)\text{Po}(np)n>50n > 50, np<5np < 5
nn large, pp not near 00 or 11N(np,npq)N(np, npq) with continuity correctionnp>5np > 5 and nq>5nq > 5

The normal approximation to the binomial is from S1. The deciding number is npnp: below 55, use the Poisson; above 55 (with nq>5nq > 5 too), use the normal.

When pp is close to 1

If almost every trial is a success, the number of failures is the rare event. Define Y=n−XY = n - X, the number of failures. Then Y∼B(n,q)Y \sim B(n, q) with q=1−pq = 1 - p small, and YY can be approximated by Po(nq)\text{Po}(nq). Translate the question about successes into a question about failures before you approximate.

Using a Poisson approximation
  1. Define XX and state the exact distribution: X∼B(n,p)X \sim B(n, p).
  2. Check the conditions: n>50n > 50 and np<5np < 5. If pp is near 11, switch to counting failures.
  3. State the approximating distribution: X∼Po(np)X \sim \text{Po}(np) approximately, giving the numerical value of npnp.
  4. Calculate the required probability with the Poisson formula. No continuity correction: both distributions are discrete, on the same whole numbers.
  5. Give the answer to 3 significant figures.
A routine approximation

On average 2%2\% of the light bulbs produced by a machine are faulty. A random sample of 150150 bulbs is taken.

(a) Use a suitable approximation to find the probability that the sample contains at most 22 faulty bulbs.

(b) Justify your approximation.

Solution

(a) Let XX be the number of faulty bulbs. Then X∼B(150,0.02)X \sim B(150, 0.02).

np=150×0.02=3np = 150 \times 0.02 = 3, so approximately X∼Po(3)X \sim \text{Po}(3).

P(X≤2)≈e−3(1+3+322)=8.5e−3=0.423P(X \le 2) \approx e^{-3}\left(1 + 3 + \frac{3^2}{2}\right) = 8.5e^{-3} = 0.423

(b) n=150>50n = 150 > 50 and np=3<5np = 3 < 5.

For comparison, the exact binomial probability is 0.4210.421, so the approximation is good.

Counting failures

98%98\% of the seeds in a large batch germinate. A gardener plants 200200 seeds. Use a suitable approximation to find the probability that more than 195195 seeds germinate.

Solution

The number germinating is B(200,0.98)B(200, 0.98), but p=0.98p = 0.98 is not small. Count the failures instead.

Let YY be the number of seeds that do not germinate: Y∼B(200,0.02)Y \sim B(200, 0.02).

n=200>50n = 200 > 50 and nq=4<5nq = 4 < 5, so approximately Y∼Po(4)Y \sim \text{Po}(4).

More than 195195 germinate means 196196 to 200200 germinate, so 00 to 44 fail:

P(Y≤4)≈e−4(1+4+422+436+4424)=e−4(1+4+8+10.667+10.667)=0.629P(Y \le 4) \approx e^{-4}\left(1 + 4 + \frac{4^2}{2} + \frac{4^3}{6} + \frac{4^4}{24}\right) = e^{-4}(1 + 4 + 8 + 10.667 + 10.667) = 0.629
Finding the sample size

The probability that a randomly chosen person has a particular rare blood group is 0.0040.004. Use a suitable approximation to find the least number of people who must be tested so that the probability of finding at least one person with this blood group is greater than 0.950.95.

Solution

Let XX be the number with the blood group among nn people: X∼B(n,0.004)X \sim B(n, 0.004). With nn large and pp small, approximately X∼Po(0.004n)X \sim \text{Po}(0.004n).

P(X≥1)=1−e−0.004n>0.95P(X \ge 1) = 1 - e^{-0.004n} > 0.95e−0.004n<0.05  ⇒  −0.004n<ln⁡0.05  ⇒  n>ln⁡200.004=748.9e^{-0.004n} < 0.05 \;\Rightarrow\; -0.004n < \ln 0.05 \;\Rightarrow\; n > \frac{\ln 20}{0.004} = 748.9

The least number is 749749.

Check the approximation is reasonable: n=749>50n = 749 > 50 and np=3.0<5np = 3.0 < 5. (The exact binomial calculation, 1−0.996n>0.951 - 0.996^n > 0.95, gives 748748. The approximation is close but not identical, which is why the question specifies the method.)

An exam-style multi-part question

0.8%0.8\% of the eggs from a farm have a double yolk. The eggs are sold in boxes of 250250.

(a) Use a suitable approximation to find the probability that a box contains more than 33 double-yolked eggs.

(b) Explain why a normal approximation would not be appropriate here.

(c) A shop buys 66 boxes. Find the probability that at most one of these boxes contains more than 33 double-yolked eggs.

Solution

(a) Let XX be the number of double-yolked eggs in a box. X∼B(250,0.008)X \sim B(250, 0.008).

n=250>50n = 250 > 50 and np=2<5np = 2 < 5, so approximately X∼Po(2)X \sim \text{Po}(2).

P(X>3)=1−P(X≤3)≈1−e−2(1+2+2+43)=1−0.85712=0.14288P(X > 3) = 1 - P(X \le 3) \approx 1 - e^{-2}\left(1 + 2 + 2 + \frac{4}{3}\right) = 1 - 0.85712 = 0.14288

So P(X>3)=0.143P(X > 3) = 0.143 (3 s.f.).

(b) For a normal approximation we need np>5np > 5, but np=2np = 2. The distribution is very skewed, not bell-shaped.

(c) Let NN be the number of the 6 boxes with more than 3 double-yolked eggs. N∼B(6,0.14288)N \sim B(6, 0.14288).

P(N≤1)=0.857126+6(0.14288)(0.85712)5=0.39651+0.39658=0.793P(N \le 1) = 0.85712^6 + 6(0.14288)(0.85712)^5 = 0.39651 + 0.39658 = 0.793

Note that NN is not approximated: n=6n = 6 is small, so the exact binomial is used.

Common mistakes
  • Applying a continuity correction. The continuity correction is only for approximating a discrete distribution by a continuous one. Binomial to Poisson is discrete to discrete: no correction.
  • Using the variance npqnpq as λ\lambda. The Poisson parameter is the mean, npnp.
  • Approximating when pp is near 1 without switching. B(200,0.98)B(200, 0.98) has np=196np = 196, nowhere near 55. Count failures.
  • Forgetting to rewrite the event after switching. "More than 195 successes" becomes "at most 4 failures". Check with an extreme case: 200 successes is 0 failures.
  • Using the Poisson when npnp is large. If np>5np > 5 and nq>5nq > 5, the normal is the suitable approximation.
  • Justifying with only one condition. "Because nn is large" is incomplete. Give both n>50n > 50 and np<5np < 5, with numbers.
Exam tip
  • "Use a suitable approximation" means you must choose and name it. Write "X∼B(150,0.02)X \sim B(150, 0.02), approximated by Po(3)\text{Po}(3)".
  • "Justify" or "explain why your approximation is valid": quote the numbers, "n=150>50n = 150 > 50 and np=3<5np = 3 < 5". Words alone ("nn is large and pp is small") are often accepted, but numbers are safer.
  • If you calculate the exact binomial probability when an approximation was asked for, you may lose the method marks even though the answer is close. Follow the instruction.
  • If a question does not ask for an approximation and nn is small enough to compute exactly (for example n=12n = 12), use the exact binomial.
  • Watch for the switch back: once you have found a probability with the Poisson, a follow-up "how many of these boxes" question is an exact binomial with small nn.
Summary
  • B(n,p)≈Po(np)B(n, p) \approx \text{Po}(np) when nn is large and pp is small: n>50n > 50 and np<5np < 5.
  • Means agree exactly; variances npqnpq and npnp are close because q≈1q \approx 1.
  • No continuity correction: both distributions are discrete.
  • If pp is close to 11, count failures: n−X∼B(n,q)≈Po(nq)n - X \sim B(n, q) \approx \text{Po}(nq).
  • Use the normal approximation instead when np>5np > 5 and nq>5nq > 5.
  • Justify with the numbers; follow "use a suitable approximation" literally.

Practice questions

Question
  1. State, with a reason, which approximation (if any) is suitable for each distribution: (a) B(60,0.03)B(60, 0.03), (b) B(40,0.05)B(40, 0.05), (c) B(200,0.45)B(200, 0.45), (d) B(500,0.996)B(500, 0.996).
  2. X∼B(80,0.025)X \sim B(80, 0.025). Use a suitable approximation to find (a) P(X=3)P(X = 3), (b) P(X≤1)P(X \le 1).
  3. A rare condition affects 11 in 20002000 people. Use a suitable approximation to find the probability that, in a town of 60006000 people, more than 44 people have the condition.
  4. 99%99\% of components pass a quality test. A batch of 300300 components is tested. Use a suitable approximation to find the probability that at least 297297 pass.
  5. The probability that an item is defective is 0.010.01. Use a suitable approximation to find the least number of items that must be inspected for the probability of finding at least one defective to be at least 0.90.9.
  6. X∼B(60,0.05)X \sim B(60, 0.05). Calculate P(X=2)P(X = 2) exactly and using a Poisson approximation. Find the percentage error in the approximation.
  7. The probability that a sample of 100100 items from a large batch contains no defective items is 0.1350.135. Use a Poisson approximation to estimate the proportion of defective items in the batch, and hence estimate the probability that a sample of 250250 items contains at most 22 defectives.
  8. The probability that a hen's egg has a double yolk is 0.0040.004. Eggs are packed in crates of 360360. (a) Use a suitable approximation to find the probability that a crate contains at least 33 double-yolked eggs, and justify the approximation. (b) A shop sells 1010 crates in a week. Find the probability that more than one of these crates contains at least 33 double-yolked eggs. (c) Using a Poisson approximation, find the largest number of eggs that can be packed in a box if the probability that the box contains no double-yolked egg is to be greater than 0.950.95.
Answers
  1. (a) Poisson, Po(1.8)\text{Po}(1.8): n=60>50n = 60 > 50 and np=1.8<5np = 1.8 < 5. (b) None needed; n=40n = 40 is not large, so calculate exactly (Poisson would be borderline at best). (c) Normal, N(90,49.5)N(90, 49.5): np=90>5np = 90 > 5 and nq=110>5nq = 110 > 5. (d) Count failures: B(500,0.004)B(500, 0.004) for failures, approximated by Po(2)\text{Po}(2) since n>50n > 50 and nq=2<5nq = 2 < 5.

  2. np=2np = 2, so X≈Po(2)X \approx \text{Po}(2). (a) P(X=3)=e−286=0.180P(X = 3) = e^{-2}\dfrac{8}{6} = 0.180. (b) P(X≤1)=e−2(1+2)=3e−2=0.406P(X \le 1) = e^{-2}(1 + 2) = 3e^{-2} = 0.406.

  3. X∼B(6000,0.0005)≈Po(3)X \sim B(6000, 0.0005) \approx \text{Po}(3). P(X>4)=1−e−3(1+3+4.5+4.5+3.375)=1−16.375e−3=1−0.8153=0.185P(X > 4) = 1 - e^{-3}\left(1 + 3 + 4.5 + 4.5 + 3.375\right) = 1 - 16.375e^{-3} = 1 - 0.8153 = 0.185.

  4. Failures Y∼B(300,0.01)≈Po(3)Y \sim B(300, 0.01) \approx \text{Po}(3). At least 297297 pass means at most 33 fail. P(Y≤3)=e−3(1+3+4.5+4.5)=13e−3=0.647P(Y \le 3) = e^{-3}(1 + 3 + 4.5 + 4.5) = 13e^{-3} = 0.647.

  5. X≈Po(0.01n)X \approx \text{Po}(0.01n). 1−e−0.01n≥0.9⇒e−0.01n≤0.1⇒n≥100ln⁡10=230.31 - e^{-0.01n} \ge 0.9 \Rightarrow e^{-0.01n} \le 0.1 \Rightarrow n \ge 100\ln 10 = 230.3. Least n=231n = 231.

  6. Exact: (602)(0.05)2(0.95)58=0.22588\binom{60}{2}(0.05)^2(0.95)^{58} = 0.22588. Approximation Po(3)\text{Po}(3): e−392=0.22404e^{-3}\dfrac{9}{2} = 0.22404. Percentage error =0.22588−0.224040.22588×100=0.81%= \dfrac{0.22588 - 0.22404}{0.22588} \times 100 = 0.81\%.

  7. e−100p=0.135⇒100p=−ln⁡0.135=2.0025⇒p=0.0200e^{-100p} = 0.135 \Rightarrow 100p = -\ln 0.135 = 2.0025 \Rightarrow p = 0.0200 (3 s.f.). For 250250 items, λ=250×0.020025=5.006\lambda = 250 \times 0.020025 = 5.006, so X≈Po(5.006)X \approx \text{Po}(5.006). P(X≤2)=e−5.006(1+5.006+5.00622)=0.124P(X \le 2) = e^{-5.006}\left(1 + 5.006 + \dfrac{5.006^2}{2}\right) = 0.124. (Using λ=5\lambda = 5 gives 0.1250.125, also acceptable.)

  8. (a) X∼B(360,0.004)X \sim B(360, 0.004); n=360>50n = 360 > 50, np=1.44<5np = 1.44 < 5, so X≈Po(1.44)X \approx \text{Po}(1.44). P(X≥3)=1−e−1.44(1+1.44+1.4422)=1−e−1.44(3.4768)=1−0.82375=0.176P(X \ge 3) = 1 - e^{-1.44}\left(1 + 1.44 + \dfrac{1.44^2}{2}\right) = 1 - e^{-1.44}(3.4768) = 1 - 0.82375 = 0.176. (b) N∼B(10,0.17625)N \sim B(10, 0.17625). P(N>1)=1−0.8237510−10(0.17625)(0.82375)9=1−0.14386−0.30782=0.548P(N > 1) = 1 - 0.82375^{10} - 10(0.17625)(0.82375)^9 = 1 - 0.14386 - 0.30782 = 0.548. (c) For nn eggs, P(none)≈e−0.004n>0.95⇒n<−ln⁡0.950.004=12.82P(\text{none}) \approx e^{-0.004n} > 0.95 \Rightarrow n < \dfrac{-\ln 0.95}{0.004} = 12.82. Largest number =12= 12.

How well do you know this?

Builds on

Where this leads

Console

Search notes, courses and tools, or run an action