physics2u
Depth
PU-106 / L10 · lecture

Discrete Distributions I: Bernoulli and Binomial

~6 min read · Tier toggle changes depth (T1 plain · T2 standard · T3 rigour)

BEFORE YOU START

You should be comfortable with basic combinatorics (counting permutations and combinations) and the definition of expectation and variance for a single random variable. Familiarity with the concept of independent events is also assumed.

Why this matters

The binomial distribution is the fundamental model for any experiment consisting of repeated yes/no trials.

In experimental physics, we frequently encounter situations where a detector either registers a hit or it does not, a particle decays or it does not, or a photon is absorbed or it is not. These are binary outcomes. When we repeat such a trial many times, the total number of successes follows a specific statistical pattern known as the binomial distribution. Understanding this distribution allows us to calculate the probability of observing a specific number of events, determine the most likely outcome, and quantify the uncertainty in our counts. It serves as the discrete counterpart to the Gaussian distribution and is essential for interpreting data from particle detectors, nuclear decay experiments, and digital signal processing.

The Bernoulli Trial

At the heart of the binomial distribution lies the Bernoulli trial. A Bernoulli trial is a single experiment with exactly two possible outcomes, conventionally labeled "success" and "failure." The probability of success is denoted by \( p \), and the probability of failure is \( q = 1 - p \). Crucially, the trials must be independent, meaning the outcome of one trial does not influence the outcome of any other. For example, flipping a fair coin is a Bernoulli trial with \( p = 0.5 \). In a physics context, consider a detector with an efficiency of 90%. Each time a particle passes through, the detector either triggers (success) or fails to trigger (failure). If we assume the detector's state does not change between particles, each passage is an independent Bernoulli trial with \( p = 0.9 \).

While a single Bernoulli trial is simple, real experiments involve many such trials. If we perform \( N \) independent Bernoulli trials, we are interested in the total number of successes, let us call this random variable \( S \). The variable \( S \) can take any integer value from 0 to \( N \). The probability distribution of \( S \) is the binomial distribution.

The Binomial Probability Mass Function

To find the probability that exactly \( k \) successes occur in \( N \) trials, we must consider two factors: the number of ways to arrange these successes and the probability of any specific arrangement. First, the number of distinct sequences of \( N \) trials containing exactly \( k \) successes is given by the binomial coefficient \( \binom{N}{k} \), which counts the combinations of \( k \) successes chosen from \( N \) positions. Second, for any specific sequence with \( k \) successes and \( N-k \) failures, the probability is \( p^k q^{N-k} \), since the trials are independent.

Combining these, the probability mass function (PMF) for the binomial distribution is:

\[ P(S = k) = \binom{N}{k} p^k (1-p)^{N-k} \]

where \( k = 0, 1, \dots, N \). This formula encapsulates the entire statistical behavior of the experiment. It is important to note that the sum of probabilities over all possible values of \( k \) must equal 1, which is guaranteed by the binomial theorem: \( \sum_{k=0}^{N} \binom{N}{k} p^k (1-p)^{N-k} = (p + (1-p))^N = 1 \).

Mean and Variance

For any random variable, the mean (expectation) and variance are key descriptors of its distribution. For the binomial distribution, these have simple closed-form expressions. The mean number of successes is:

\[ \mathbb{E}[S] = Np \]

This result is intuitive: if you flip a coin \( N \) times with a 50% chance of heads, you expect \( N/2 \) heads. More generally, the expected value scales linearly with the number of trials.

The variance, which measures the spread of the distribution around the mean, is:

\[ \mathrm{Var}(S) = Np(1-p) \]

The standard deviation is therefore \( \sigma_S = \sqrt{Np(1-p)} \). Notice that the variance is maximized when \( p = 0.5 \), which corresponds to the most "uncertain" trial (like a fair coin flip). If \( p \) is very close to 0 or 1, the variance is small because the outcome is more predictable.

Connection to the Gaussian Distribution

When \( N \) is large and \( p \) is not too close to 0 or 1, the binomial distribution approximates a Gaussian (normal) distribution. This is a consequence of the Central Limit Theorem. The Gaussian approximation is given by:

\[ P(S = k) \approx \frac{1}{\sqrt{2\pi Np(1-p)}} \exp\left( -\frac{(k - Np)^2}{2Np(1-p)} \right) \]

This approximation is useful for calculations when \( N \) is large, as it allows us to use continuous calculus techniques. However, for small \( N \) or extreme values of \( p \), the exact binomial formula should be used. The transition from discrete to continuous distributions is a central theme in probability theory and has profound implications for how we model physical systems.

The derivation

◆ DERIVATION EMBED · D-115 · not re-derived here

The Binomial Distribution from Bernoulli Trials

Let S = X1 + … + XN be the number of successes in N independent, identically distributed Bernoulli trials, each success carrying probability p ∈ [0,1] and failure probability q = 1 − p. Then S is binomially distributed with probability mass function P(S = k) = NCk pk qN−k for k = 0, 1, …, N, with mean Np and variance Np(1 − p).

P(S=k) = NCk pk(1−p)N−k,    E[S] = Np,    Var(S) = Np(1−p)
Open full derivation →
Tier 3 ยท watch out

A common misconception is that the binomial distribution requires the trials to be physically identical. In fact, independence is the key requirement. If the trials are not independent (e.g., if the detector efficiency changes over time due to heating), the binomial model is invalid, and the variance will be larger than \( Np(1-p) \). This phenomenon is known as over-dispersion. Always check for independence before applying the binomial formula. C

Worked example

Problem. A particle detector has an efficiency of \( p = 0.8 \). A beam of \( N = 100 \) particles is sent through the detector. Calculate the probability that exactly 85 particles are detected. Also, find the mean and standard deviation of the number of detected particles.

  1. Identify the parameters: \( N = 100 \), \( p = 0.8 \), \( q = 1 - p = 0.2 \), and \( k = 85 \).
  2. Use the binomial PMF: \( P(S = 85) = \binom{100}{85} (0.8)^{85} (0.2)^{15} \).
  3. Calculate the binomial coefficient: \( \binom{100}{85} = \frac{100!}{85! 15!} \approx 2.53 \times 10^{17} \).
  4. Calculate the probability terms: \( (0.8)^{85} \approx 5.79 \times 10^{-9} \) and \( (0.2)^{15} \approx 3.28 \times 10^{-11} \).
  5. Multiply the terms: \( P(S = 85) \approx 2.53 \times 10^{17} \times 5.79 \times 10^{-9} \times 3.28 \times 10^{-11} \approx 0.048 \).
  6. Calculate the mean: \( \mathbb{E}[S] = Np = 100 \times 0.8 = 80 \).
  7. Calculate the standard deviation: \( \sigma_S = \sqrt{Np(1-p)} = \sqrt{100 \times 0.8 \times 0.2} = \sqrt{16} = 4 \).
\[ P(S=85) \approx 0.048, \quad \mathbb{E}[S] = 80, \quad \sigma_S = 4 \]

Units check. The probability is dimensionless. The mean and standard deviation are counts of particles, which are dimensionless integers. The results are consistent.

CHECK YOURSELF

Three quick tests

  1. What is the probability of getting exactly 3 heads in 5 fair coin flips?
    Answer\( \binom{5}{3} (0.5)^3 (0.5)^2 = 10 \times 0.125 \times 0.25 = 0.3125 \).
  2. If \( p = 0.1 \) and \( N = 1000 \), what is the mean number of successes?
    Answer\( \mathbb{E}[S] = 1000 \times 0.1 = 100 \).
  3. For which value of \( p \) is the variance of the binomial distribution maximized?
    Answer\( p = 0.5 \), where \( \mathrm{Var}(S) = N/4 \).
PROBLEMS

Problem set for PU-106 →

Each derivation in this unit ends with five graded problems with full solutions.

NOTES

Lecture notes drafted with a local model and checked line by line by the editor; every number in the worked example was recomputed.

Derivations cited: D-115 · Threads: chance