Session 5: Named Models and Extremes

Movement II: One Random Variable

EE 520: Random Processes (Fall 2026). Instructor: John Lipor.

Launch

Session 4 built random variables from scratch, one ad-hoc count at a time. But the same few shapes keep reappearing, and three questions from a single search make the point. How many of 2525 cells produce a detection on one sweep? Each cell succeeds or fails on its own, a fixed number of independent yes/no trials: that count is binomial. How many passes until this cell finally gives up its target? We repeat one trial until the first success: that waiting time is geometric. If two vehicles each search their own cell, how long until the whole job is done? That is the maximum of two waiting times, and here our tools run out: we have no formula for a maximum. Today we catalog the handful of named models that cover almost every sensor and count in the course, and we forge the one genuinely new tool of the session, the CDF method, which answers every question about a maximum or a minimum.

Development

Bernoulli and binomial

The atom of every count is a single yes/no trial.

Definition 1 (Bernoulli). X∼Bernoulli(p)X \sim \mathrm{Bernoulli}(p) takes the value 11 with probability pp and 00 with probability 1−p1 - p. It is the indicator of a success, so 𝔼[X]=p\mathbb{E}[X] = p and Var⁡(X)=p(1−p)\operatorname{Var}(X) = p(1 - p).

Stacking nn independent copies of the same trial and counting the successes gives the workhorse count model.

Definition 2 (Binomial). X∼Binomial(n,p)X \sim \mathrm{Binomial}(n, p) is the number of successes in nn independent trials, each succeeding with probability pp. Its PMF is pX(k)=(nk)pk(1−p)n−k,k=0,1,…,n,\begin{equation} p_X(k) = \binom{n}{k} p^k (1 - p)^{n - k}, \qquad k = 0, 1, \dots, n, \end{equation} the factor (nk)\binom{n}{k} counting which trials succeeded.

We do not need the PMF to find the mean and variance: Session 4’s machinery delivers them for free. Write X=∑i=1nXiX = \sum_{i=1}^n X_i with Xi∼Bernoulli(p)X_i \sim \mathrm{Bernoulli}(p) the indicator of trial ii. By linearity (no independence needed), 𝔼[X]=∑i=1n𝔼[Xi]=np.\begin{equation} \mathbb{E}[X] = \sum_{i=1}^n \mathbb{E}[X_i] = np. \end{equation} Because the trials are independent, the cross terms in the variance of the sum vanish and the variances add: Var⁡(X)=∑i=1nVar⁡(Xi)=np(1−p).\begin{equation} \operatorname{Var}(X) = \sum_{i=1}^n \operatorname{Var}(X_i) = np(1 - p). \end{equation} This is the whole point of the indicator method: the binomial is not a new object, it is nn Bernoullis added up, and its summary statistics follow from linearity and independence exactly as in the last session.

Geometric: waiting for the first success

Instead of fixing the number of trials, fix the goal (the first success) and let the number of trials be random.

Definition 3 (Geometric). X∼Geometric(q)X \sim \mathrm{Geometric}(q) is the number of independent trials up to and including the first success, where each trial succeeds with probability qq. Its support is {1,2,3,…}\left\{1, 2, 3, \dots\right\} and pX(k)=q(1−q)k−1,k=1,2,3,…,\begin{equation} p_X(k) = q(1 - q)^{k - 1}, \qquad k = 1, 2, 3, \dots, \end{equation} since the first k−1k - 1 trials must fail and the kk-th must succeed.

Convention (fix this once and never deviate). We count the trial on which the first success occurs, so the geometric starts at 11. Many textbooks instead count the number of failures before the first success, which starts at 00; the two conventions differ by exactly 11. Ours is the “number of passes until this cell is searched successfully” count, which is what a search log records.

The survival probability follows directly: the event {X>k}\left\{X > k\right\} is exactly “the first kk trials all failed,” so ℙ(X>k)=(1−q)k,k=0,1,2,….\begin{equation} \label{eq:geomtail} \mathbb{P}(X > k) = (1 - q)^k, \qquad k = 0, 1, 2, \dots. \end{equation} For the mean we introduce a general tool that will pay off repeatedly.

Theorem 1 (Tail-sum formula). For a nonnegative integer-valued random variable XX, 𝔼[X]=∑k≥1ℙ(X≥k).\begin{equation} \mathbb{E}[X] = \sum_{k \ge 1} \mathbb{P}(X \ge k). \end{equation}

Proof: HW2, Problem 1(b). We use only the statement here. Applied to the geometric, ℙ(X≥k)=ℙ(X>k−1)=(1−q)k−1\mathbb{P}(X \ge k) = \mathbb{P}(X > k - 1) = (1 - q)^{k - 1} by [eq:geomtail], so the mean is a geometric series: 𝔼[X]=∑k≥1(1−q)k−1=11−(1−q)=1q.\begin{equation} \mathbb{E}[X] = \sum_{k \ge 1} (1 - q)^{k - 1} = \frac{1}{1 - (1 - q)} = \frac{1}{q}. \end{equation} A cell searched with per-pass success probability qq takes, on average, 1/q1/q passes to give up its target. The rarer the success, the longer the wait, in exact inverse proportion.

Memorylessness

The geometric has a startling property: having waited without success tells you nothing about how much longer you must wait.

Proposition 1 (Memorylessness). If X∼Geometric(q)X \sim \mathrm{Geometric}(q), then for all nonnegative integers m,km, k, ℙ(X>m+k|X>m)=ℙ(X>k).\begin{equation} \mathbb{P}(X > m + k \,|\,X > m) = \mathbb{P}(X > k). \end{equation}

Proof. Using the definition of conditional probability and the survival formula [eq:geomtail], and noting {X>m+k}⊆{X>m}\left\{X > m + k\right\} \subseteq \left\{X > m\right\} so the intersection is {X>m+k}\left\{X > m + k\right\}, ℙ(X>m+k|X>m)=ℙ(X>m+k)ℙ(X>m)=(1−q)m+k(1−q)m=(1−q)k=ℙ(X>k).\begin{equation} \mathbb{P}(X > m + k \,|\,X > m) = \frac{\mathbb{P}(X > m + k)}{\mathbb{P}(X > m)} = \frac{(1 - q)^{m + k}}{(1 - q)^m} = (1 - q)^k = \mathbb{P}(X > k). \end{equation} ◻

In search terms: a cell that has resisted mm passes is, given the target is still there, exactly as hard to clear as it was on the first pass. The trials do not tire and do not learn; each pass faces the same odds.

This is easy to misread, and the misreading collides with everything Movement I taught. Memorylessness is a statement about the conditional law of the waiting time given that the target is present in the cell. It says nothing about whether the target is present. Those are different questions, and Session 2 already answered the second one: each failed pass is evidence, and the failed-search update flows probability out of the searched cell, so unconditionally our belief that the target is even in this cell falls with every empty pass. Both facts hold at once, and they do not conflict. Conditioned on presence, the remaining wait is undiminished; but our confidence in that very conditioning erodes. Persistence that ignores the second fact, and searches forever because “the odds each pass are unchanged,” is exactly the error Session 2 warned against.

Two clocks, not one. Memorylessness governs the wait given presence; Bayes governs belief in presence. A run of failures leaves the first clock untouched and runs the second one down. Never let “the geometric is memoryless” become “the failures told us nothing”: they told us nothing about the wait, and a great deal about the location.

Poisson: the model for rare counts

Definition 4 (Poisson). X∼Poisson(λ)X \sim \mathrm{Poisson}(\lambda) has support {0,1,2,…}\left\{0, 1, 2, \dots\right\} and PMF pX(k)=λke−λk!,k=0,1,2,….\begin{equation} p_X(k) = \frac{\lambda^k e^{-\lambda}}{k!}, \qquad k = 0, 1, 2, \dots. \end{equation} Its mean and variance are both λ\lambda (the mean is a short exercise from the series for eλe^\lambda; the variance is left as one).

The Poisson earns its central place as the limit of a binomial with many trials, each rarely successful. Fix a count kk and let n→∞n \to \infty while holding the mean np=λnp = \lambda fixed, i.e., p=λ/np = \lambda / n: (nk)(λn)k(1−λn)n−k=n(n−1)⋯(n−k+1)nk⏟→1⋅λkk!⋅(1−λn)n⏟→e−λ⋅(1−λn)−k⏟→1→λke−λk!.\begin{equation} \begin{aligned} \binom{n}{k} \left(\frac{\lambda}{n}\right)^k \left(1 - \frac{\lambda}{n}\right)^{n - k} &= \underbrace{\frac{n(n-1)\cdots(n-k+1)}{n^k}}_{\to\, 1} \cdot \frac{\lambda^k}{k!} \cdot \underbrace{\left(1 - \frac{\lambda}{n}\right)^{n}}_{\to\, e^{-\lambda}} \cdot \underbrace{\left(1 - \frac{\lambda}{n}\right)^{-k}}_{\to\, 1} \\ &\longrightarrow \frac{\lambda^k e^{-\lambda}}{k!}. \end{aligned} \end{equation}

The first factor is a product of kk ratios each tending to 11, and the third is the classical limit (1−λ/n)n→e−λ(1 - \lambda/n)^n \to e^{-\lambda}. This is the law of rare events: a great many independent chances, each unlikely, produce a Poisson count.

The limit is not a curiosity. It says the Poisson is what a binomial becomes whenever the batch is large and each success is rare, which is exactly the regime of most sensing counts.

Poisson is the universal rare-count model. Photons striking a detector in a gate, packets arriving at a relay, false alarms across a large sensor array: each is a huge number of independent, individually unlikely opportunities. That is why Poisson is the count model of choice, and why it returns in Session 10 as the noise model for photon-counting detection.

The CDF method for extremes

Now to the tool the Launch promised. When several independent quantities are in play and we ask for the largest or the smallest, neither the PMF nor the PDF is the natural object; the CDF is. The reason is that “max⁡≤x\max \le x” and “min⁡>x\min > x” both factor cleanly across independent variables, whereas “max⁡=x\max = x” does not.

Theorem 2 (Distributions of a maximum and a minimum). Let X1,…,XmX_1, \dots, X_m be independent and identically distributed with common CDF FF. Write M=max⁡iXiM = \max_i X_i and L=min⁡iXiL = \min_i X_i. Then ℙ(M≤x)=F(x)m,ℙ(L>x)=(1−F(x))m.\begin{equation} \mathbb{P}(M \le x) = F(x)^m, \qquad \mathbb{P}(L > x) = \big(1 - F(x)\big)^m. \end{equation}

Proof. The maximum is at most xx exactly when every XiX_i is at most xx: ℙ(M≤x)=ℙ(X1≤x,…,Xm≤x)=∏i=1mℙ(Xi≤x)=F(x)m,\begin{equation} \mathbb{P}(M \le x) = \mathbb{P}(X_1 \le x, \dots, X_m \le x) = \prod_{i=1}^m \mathbb{P}(X_i \le x) = F(x)^m, \end{equation} using independence to factor and identical distribution to collapse each factor to F(x)F(x). Symmetrically, the minimum exceeds xx exactly when every XiX_i exceeds xx: ℙ(L>x)=ℙ(X1>x,…,Xm>x)=∏i=1mℙ(Xi>x)=(1−F(x))m.\begin{equation} \mathbb{P}(L > x) = \mathbb{P}(X_1 > x, \dots, X_m > x) = \prod_{i=1}^m \mathbb{P}(X_i > x) = \big(1 - F(x)\big)^m. \end{equation} In each case independence turned the joint statement about all mm variables into a product, and identical distribution collapsed the factors. ◻

The two formulas share one moral, worth stating on its own before we use them.

When you hear “max” or “min,” reach for the CDF. The event {max≤x}\left\{\max \le x\right\} and the event {min>x}\left\{\min > x\right\} each say something about all the variables at once, and independence turns “all at once” into a product. Never try to attack a maximum or a minimum through its PMF or density directly; write down its CDF (or its survival probability) and differentiate only at the end if you need a density.

Worked examples

(a) When is the job finished? The maximum of two geometrics

The Launch asked when a two-cell job is finished. Each of two cells is known to hold a target, and a dedicated vehicle searches each one, the passes synchronized; on every pass each vehicle independently detects its target with probability q=0.5q = 0.5. Vehicle ii’s time to detection is Xi∼Geometric(0.5)X_i \sim \mathrm{Geometric}(0.5), and the job is done at M=max⁡(X1,X2)M = \max(X_1, X_2).

Distribution. Reach for the CDF. Each XiX_i has F(k)=ℙ(Xi≤k)=1−0.5kF(k) = \mathbb{P}(X_i \le k) = 1 - 0.5^{\,k} by [eq:geomtail], so the maximum half of Theorem 2 gives ℙ(M≤k)=(1−0.5k)2,e.g.ℙ(M≤3)=(78)2=4964≈0.766.\begin{equation} \mathbb{P}(M \le k) = \big(1 - 0.5^{\,k}\big)^2, \qquad\text{e.g.}\qquad \mathbb{P}(M \le 3) = \left(\tfrac 78\right)^2 = \tfrac{49}{64} \approx 0.766. \end{equation} Three passes finish the job only about three times in four, even though each single target is found within three passes with probability 7/87/8.

Mean. The tail-sum formula (Theorem 1) asks for ℙ(M≥k)=ℙ(M>k−1)\mathbb{P}(M \ge k) = \mathbb{P}(M > k - 1), so work with the survival probability ℙ(M>j)=1−(1−0.5j)2=2(0.5)j−(0.25)j,j=0,1,2,….\begin{equation} \mathbb{P}(M > j) = 1 - \big(1 - 0.5^{\,j}\big)^2 = 2\,(0.5)^j - (0.25)^j, \qquad j = 0, 1, 2, \dots. \end{equation} Reindexing with j=k−1j = k - 1 and summing two geometric series, 𝔼[M]=∑j≥0ℙ(M>j)=2∑j≥0(0.5)j−∑j≥0(0.25)j=2⋅2−43=83≈2.67 passes.\begin{equation} \mathbb{E}[M] = \sum_{j \ge 0} \mathbb{P}(M > j) = 2\sum_{j \ge 0} (0.5)^j - \sum_{j \ge 0} (0.25)^j = 2 \cdot 2 - \frac{4}{3} = \frac{8}{3} \approx 2.67 \text{ passes}. \end{equation} A single cell takes 1/q=21/q = 2 passes on average. Two cells searched in parallel take longer, because we must wait for the slower search, but far less than twice as long.

The job has memory. Each search is memoryless, but the job is not. The chance of finishing on the next pass, given that the job is not yet done, is ℙ(M=k+1|M>k)=1−ℙ(M>k+1)/ℙ(M>k)\mathbb{P}(M = k + 1 \,|\,M > k) = 1 - \mathbb{P}(M > k + 1)/\mathbb{P}(M > k), which equals 1/41/4 on the first pass (both targets must be found at once), 5/12≈0.4175/12 \approx 0.417 on the second, and climbs toward 1/21/2 thereafter: the longer the job runs, the more likely one target is already found, leaving a single geometric search. Proposition 1 is a property of one geometric clock, not of a maximum of them.

(b) A binomial that is really a Poisson, and one that is not

A large array has n=1000n = 1000 cells, and each independently raises a false alarm (a detection where nothing is present) with probability 0.0030.003 on a sweep. Let CC be the number of false alarms. Exactly, C∼Binomial(1000,0.003)C \sim \mathrm{Binomial}(1000, 0.003), with mean np=3np = 3. Because nn is large and pp is small with npnp moderate, the law of rare events says CC is very nearly Poisson(3)\mathrm{Poisson}(3).

Compare the two at k=0k = 0, the probability of a clean sweep: ℙ(C=0)exact=(1−0.003)1000=0.9971000≈0.0496,ℙ(C=0)Poisson=e−3≈0.0498.\begin{align} \mathbb{P}(C = 0)_{\text{exact}} &= (1 - 0.003)^{1000} = 0.997^{1000} \approx 0.0496, \\ \mathbb{P}(C = 0)_{\text{Poisson}} &= e^{-3} \approx 0.0498. \end{align} The two agree to three decimal places. Is it the large nn that makes this work? Try a binomial with the same mean but a few likely trials: n=10n = 10 cells with p=0.3p = 0.3, so again np=3np = 3. Now ℙ(C=0)=0.710≈0.0282\mathbb{P}(C = 0) = 0.7^{10} \approx 0.0282, while the Poisson approximation still says e−3≈0.0498e^{-3} \approx 0.0498, almost twice too large. Matching the mean does not make a count Poisson.

The limit computation tells us what does. At k=0k = 0 the exact value is (1−p)n=exp⁡(nln⁡(1−p))(1 - p)^n = \exp\big(n \ln(1 - p)\big), and the series ln⁡(1−p)=−p−p2/2−p3/3−⋯\ln(1 - p) = -p - p^2/2 - p^3/3 - \cdots (valid for 0≤p<10 \le p < 1) gives, with λ=np\lambda = np, (1−p)n=e−λ⋅exp⁡(−λ(p2+p23+⋯)).\begin{equation} (1 - p)^n = e^{-\lambda} \cdot \exp\!\Big(-\lambda\big(\tfrac{p}{2} + \tfrac{p^2}{3} + \cdots\big)\Big). \end{equation} The first factor is the Poisson value; the second is the error, and with λ\lambda fixed it is controlled by pp, not by nn. For p=0.003p = 0.003 it is about e−λp/2=e−0.0045≈0.9955e^{-\lambda p/2} = e^{-0.0045} \approx 0.9955, which is exactly the ratio of the two unrounded values above. For p=0.3p = 0.3 the correction terms are no longer small, and the exact value is only about 0.570.57 of the Poisson one. The law of rare events is a statement about rare: the Poisson model is accurate when each individual chance is unlikely, and then it needs only the single number λ\lambda, no factorials and no thousand-fold product.

Consolidation

Work these in class.

  1. (Mechanical, geometric.) A cell is searched with per-pass success probability q=0.25q = 0.25. What is the expected number of passes to detection, and what is the probability the cell survives more than 44 passes?
    Answer: 𝔼[X]=1/0.25=4\mathbb{E}[X] = 1/0.25 = 4 passes; ℙ(X>4)=(0.75)4=0.3164\mathbb{P}(X > 4) = (0.75)^4 = 0.3164.

  2. (CDF method.) Two independent fair dice are rolled, each with CDF F(k)=k/6F(k) = k/6 for k=1,…,6k = 1, \dots, 6. Using the CDF method, find ℙ(max⁡≤4)\mathbb{P}(\max \le 4) and ℙ(min⁡≥5)\mathbb{P}(\min \ge 5).
    Answer: ℙ(max⁡≤4)=F(4)2=(4/6)2=4/9\mathbb{P}(\max \le 4) = F(4)^2 = (4/6)^2 = 4/9; and ℙ(min⁡≥5)=ℙ(min⁡>4)=(1−F(4))2=(2/6)2=1/9\mathbb{P}(\min \ge 5) = \mathbb{P}(\min > 4) = (1 - F(4))^2 = (2/6)^2 = 1/9.

  3. (Conceptual, memorylessness.) A cell has been searched 1010 times without success. A colleague argues: “the geometric is memoryless, so those failures carry no information; our belief the target is here is unchanged.” Identify the one true clause and the one false clause.
    Answer: True: given the target is present, the remaining wait is unchanged (memorylessness of the conditional law). False: our belief the target is present is unchanged; Session 2’s failed-search update has driven that belief down with every empty pass.

Synthesis

The course thesis asks first, “what is uncertain?” This session’s answer is that most of the time the uncertain count has one of a few recognizable shapes: a fixed batch of independent trials is binomial, a wait for the first success is geometric, and a flood of individually unlikely chances is Poisson. Naming the shape hands you its mean, its spread, and its behavior at a glance. Two ideas reached past the catalog. Memorylessness split one clock into two, the undiminished wait given presence and the eroding belief in presence, and reconciled them with Movement I rather than contradicting it. And the CDF method answered the questions the PMF could not: whenever “max” or “min” is spoken, the CDF factors across independent variables and the PMF does not, so the CDF is the tool. Next session carries the CDF into continuous territory, where variables have no PMF at all but every one of them still has a CDF, and where the geometric’s continuous twin, the memoryless waiting time on a clock rather than a counter, appears as the exponential distribution.

Self-check questions

  1. You need the mean of a Geometric(q)\mathrm{Geometric}(q) waiting time and you know only its survival probability ℙ(X>k)=(1−q)k\mathbb{P}(X > k) = (1-q)^k. Which formula turns the survival probabilities into the mean, and what does it give?
    Answer: the tail-sum formula 𝔼[X]=∑k≥1ℙ(X≥k)=∑k≥1(1−q)k−1=1/q\mathbb{E}[X] = \sum_{k\ge 1}\mathbb{P}(X \ge k) = \sum_{k\ge1}(1-q)^{k-1} = 1/q.

  2. Five independent sensors have iid detection times with common CDF FF. In one line, what is the probability the first detection happens after time xx, and why is this easier than working with the PMF of the minimum?
    Answer: ℙ(min⁡>x)=(1−F(x))5\mathbb{P}(\min > x) = (1 - F(x))^5; the event “all five exceed xx” factors across independent sensors, whereas “the minimum equals xx” does not.

  3. A cell has failed 88 searches, each with per-pass success probability qq. A student says the geometric’s memorylessness means the 88 failures are uninformative. What, precisely, is unchanged, and what has in fact changed?
    Answer: unchanged is the remaining wait given the target is present; changed (fallen) is the belief that the target is present at all, by the Session 2 failed-search update.