1The Poisson Distribution
▶
The Poisson distribution is one of the most important and widely applied probability distributions in statistics. It models the number of times a particular event occurs within a fixed interval — whether that interval is measured in time, distance, area, volume, or any other continuous dimension. Unlike the binomial distribution, which counts successes across a fixed number of discrete trials, the Poisson distribution applies to situations where events can happen at any point within a continuous interval, and the total number of possible occurrences is not bounded in advance. Classic examples include the number of phone calls arriving at a call centre per hour, the number of typographical errors per page of text, the number of radioactive decay events per second, or the number of potholes per kilometre of road. In each of these cases, we observe a count of events over a fixed interval, and the Poisson distribution gives us a precise mathematical tool for calculating the probability of any particular count.
The Poisson distribution rests on three foundational conditions that must hold for the model to be appropriate. First, the events must occur independently of one another — the occurrence of one event gives no information about whether or when the next event will occur. A phone call arriving at 2:03 p.m. does not make another call at 2:04 p.m. more or less likely. Second, the average rate at which events occur must be constant throughout the interval being observed. If calls arrive at an average rate of 10 per hour, that rate must apply uniformly across the hour, not spiking during some minutes and falling during others. Third, theoretically any number of events can occur within the interval — the count can be 0, 1, 2, 3, and so on without an upper ceiling. Additionally, two events cannot occur at precisely the same instant; the probability of simultaneous events is effectively zero. When all of these conditions are satisfied, the Poisson distribution provides an accurate and elegant model.
The foundation of the Poisson distribution is its rate parameter, denoted by the Greek letter λ (lambda). Lambda represents the average number of events expected to occur within the specified interval. It is always a positive real number — λ must be strictly greater than zero, because a rate of zero would mean events never happen and there would be nothing to model. Lambda can be a whole number (e.g., λ = 5 calls per hour) or a decimal (e.g., λ = 2.3 defects per metre), depending on the context.
A particularly useful property of λ is that it scales proportionally with the size of the interval. If the average number of events in one hour is λ = 12, then the average number of events in half an hour is λ = 6, in two hours it is λ = 24, and in ten minutes (one-sixth of an hour) it is λ = 2. This scaling principle is essential: before applying the Poisson formula, you must always ensure that the value of λ you use corresponds exactly to the interval specified in the problem. For example, suppose a bakery receives an average of 30 customer orders per hour. If you want to find the probability of receiving exactly 5 orders in a 10-minute window, you must first convert: 10 minutes is one-sixth of an hour, so the appropriate λ for that window is 30 × (10/60) = 5. Failing to adjust λ to match the interval in the question is one of the most common errors when applying this distribution.
The Poisson probability formula calculates the probability of observing exactly x events in an interval, given that the average rate is λ:
P(X = x) = (e^(−λ) × λ^x) / x!
Each component of this formula has a specific role:
- e is Euler's number, the base of the natural logarithm, approximately equal to 2.71828. It appears naturally in many growth and decay processes and is central to the mathematics underlying the Poisson distribution.
- e^(−λ) is the normalising factor. It ensures that when you sum P(X = x) over all possible values of x (that is, x = 0, 1, 2, 3, …), the total equals exactly 1. In other words, it guarantees that the formula produces a valid probability distribution.
- λ^x captures how the probability of observing x events changes with the rate — higher rates make higher counts more probable.
- x! (x factorial) is the product of all positive integers up to x. For example, 4! = 4 × 3 × 2 × 1 = 24. Factorials grow extremely rapidly, which is why the probability of very large counts quickly becomes very small even when λ is moderate. By convention, 0! = 1.
- x can take any non-negative integer value: 0, 1, 2, 3, and so on with no upper limit.
To make this concrete, suppose a hospital emergency department receives an average of 4 patients per hour. What is the probability that exactly 6 patients arrive in the next hour? Here λ = 4 and x = 6:
P(X = 6) = (e^(−4) × 4^6) / 6!
= (0.018316 × 4096) / 720
= 75.02 / 720
≈ 0.1042
So there is approximately a 10.4% chance of exactly 6 arrivals. Notice how each component is evaluated step by step: first compute e^(−λ), then λ^x, then x!, and finally combine them. This methodical approach minimises arithmetic errors.
The mean and variance of the Poisson distribution share a remarkable and defining property: they are both equal to λ.
- Expected value (mean): E(X) = λ
- Variance: Var(X) = λ
- Standard deviation: SD(X) = √λ
The fact that the mean equals the variance is unique to the Poisson distribution and is often used as a diagnostic check — if you observe a dataset where the sample mean and sample variance are approximately equal, this is evidence that a Poisson model may be appropriate. The practical implication of this property is also important: as λ increases, both the average count and the spread of the distribution grow together. A process with λ = 1 has a tight distribution concentrated near 0 and 1, while a process with λ = 20 has a much broader distribution centred around 20. The standard deviation in that case is √20 ≈ 4.47, meaning values between roughly 15 and 25 would be quite typical.
The following table illustrates how the Poisson probabilities are distributed for λ = 3, giving a sense of the shape of the distribution:
| x (number of events) | Calculation | P(X = x) |
|---|---|---|
| 0 | (e^(−3) × 3^0) / 0! = 0.049787 × 1 / 1 | 0.0498 |
| 1 | (e^(−3) × 3^1) / 1! = 0.049787 × 3 / 1 | 0.1494 |
| 2 | (e^(−3) × 3^2) / 2! = 0.049787 × 9 / 2 | 0.2240 |
| 3 | (e^(−3) × 3^3) / 3! = 0.049787 × 27 / 6 | 0.2240 |
| 4 | (e^(−3) × 3^4) / 4! = 0.049787 × 81 / 24 | 0.1680 |
| 5 | (e^(−3) × 3^5) / 5! = 0.049787 × 243 / 120 | 0.1008 |
| 6 | (e^(−3) × 3^6) / 6! = 0.049787 × 729 / 720 | 0.0504 |
| 7 | (e^(−3) × 3^7) / 7! = 0.049787 × 2187 / 5040 | 0.0216 |
Notice that the probabilities peak at x = 2 and x = 3 (both equal to 0.2240 here), which is consistent with the mean of λ = 3. The probabilities taper off symmetrically for smaller and larger counts, and the sum of all probabilities across every non-negative integer value of x equals 1.
In practice, many Poisson problems ask for cumulative probabilities — questions phrased as "at most," "fewer than," "at least," or "more than" a certain number of events. For cumulative probabilities, you simply sum the individual Poisson probabilities across the relevant range of x values.
- "At most k" events: P(X ≤ k) = P(X=0) + P(X=1) + … + P(X=k)
- "Fewer than k" events: P(X < k) = P(X ≤ k−1)
- "At least k" events: P(X ≥ k) = 1 − P(X ≤ k−1)
- "More than k" events: P(X > k) = 1 − P(X ≤ k)
For example, returning to the hospital scenario with λ = 4, suppose we want the probability that at most 2 patients arrive in the next hour. We compute:
P(X = 0) = (e^(−4) × 4^0) / 0! = 0.018316 × 1 / 1 = 0.018316
P(X = 1) = (e^(−4) × 4^1) / 1! = 0.018316 × 4 / 1 = 0.073263
P(X = 2) = (e^(−4) × 4^2) / 2! = 0.018316 × 16 / 2 = 0.146525
P(X ≤ 2) = 0.018316 + 0.073263 + 0.146525 = 0.2381
There is approximately a 23.8% chance of 2 or fewer arrivals in the hour. Now suppose instead we need the probability of more than 2 patients: P(X > 2) = 1 − P(X ≤ 2) = 1 − 0.2381 = 0.7619, or about 76.2%.
To summarise the step-by-step approach for solving any Poisson problem:
- Step 1 — Identify λ: Read the problem to find the stated average rate of events. Check whether the interval in the problem matches the interval for the given rate. If they differ, scale λ proportionally to match the interval of interest.
- Step 2 — Identify x: Determine the specific count (or counts) for which you need to calculate a probability. Make sure x is a non-negative integer.
- Step 3 — Apply the formula: Compute P(X = x) = (e^(−λ) × λ^x) / x! carefully, evaluating e^(−λ), λ^x, and x! as separate components before combining them.
- Step 4 — Handle cumulative questions: If the problem asks for a range of values (e.g., "at most 3"), sum the individual probabilities for each relevant value of x, or use the complement rule to simplify the calculation where appropriate.
- Step 5 — Interpret the result: Express your answer as a probability between 0 and 1, or as a percentage, and interpret it clearly in the context of the original problem.
The Poisson distribution is a cornerstone of applied statistics precisely because it captures something very real about the natural world: many processes generate events randomly in time or space at a roughly constant rate, with each event independent of the others. Mastering the conditions, the rate parameter, the formula, and the cumulative probability technique gives you a powerful and versatile tool for analysing and predicting such processes across a vast range of fields, from engineering and medicine to ecology, finance, and telecommunications.