Computing via Conditioning — Compute expectations and variances using conditioning - Statistics Graduate Level

SECTION 1

Historical Context & Motivation

The idea of computing expectations by first conditioning on auxiliary information has deep roots in the development of probability theory. Long before the formal measure-theoretic foundations were established, probabilists recognized that breaking a complex calculation into simpler conditional pieces was both natural and powerful. The law of total expectation and the law of total variance emerged as indispensable tools in the probabilist's repertoire, transforming seemingly intractable problems into manageable conditional calculations. Understanding the historical trajectory of these ideas illuminates why conditioning remains one of the most elegant and frequently deployed strategies in modern probability and statistics.

1812 • Laplace's Iterated Expectation

Pierre-Simon Laplace employed iterated expectation arguments in his Théorie analytique des probabilités, decomposing probabilities and expectations over intermediate events to solve problems involving compound random phenomena.

1933 • Kolmogorov's Axiomatization

Andrey Kolmogorov published his Grundbegriffe der Wahrscheinlichkeitsrechnung, providing a rigorous measure-theoretic definition of conditional expectation and establishing the tower property E[E[X|𝒢]|ℋ] = E[X|ℋ] for sub-σ-algebras ℋ ⊆ 𝒢.

1950s • Rao–Blackwell Theorem & Applications

The Rao–Blackwell theorem demonstrated that conditioning an estimator on a sufficient statistic never increases its variance, establishing conditional expectation as a fundamental tool in statistical estimation theory.

1968 • Eve's Law Popularized

The law of total variance, often called the variance decomposition formula or 'Eve's law' (EVVE: E[Var] + Var[E]), became widely taught in graduate curricula as a standard technique for computing variances of compound and mixture distributions.

Modern • Ubiquity in Bayesian & Computational Statistics

Conditioning arguments now underpin Bayesian hierarchical models, Rao–Blackwellized particle filters, Gibbs samplers, and variance reduction techniques in Monte Carlo simulation, making 'computing via conditioning' an everyday workhorse in applied statistics.

The central question that conditioning addresses is straightforward yet profound: given a random variable X whose distribution depends on another random quantity Y, how can we systematically compute E[X] and Var(X) by leveraging what we know about the conditional distribution of X given Y? The two laws we will study provide a complete and elegant answer.

SECTION 2

Core Principles & Definitions

Computing via conditioning rests on two foundational identities—the law of total expectation and the law of total variance—together with the concept of conditional expectation as a random variable. The key conceptual shift is recognizing that E[X|Y] is not a number but a function of the random variable Y, and therefore it is itself a random variable with its own expectation and variance. Mastering this conceptual point is the gateway to deploying the conditioning machinery fluently.

Conditional Expectation E[X|Y]

Given Y = y, the conditional expectation E[X|Y = y] is computed using the conditional distribution of X given Y = y. Viewed as a function of the random Y, E[X|Y] = g(Y) is a random variable.

Law of Total Expectation (Tower Property)

E[X] = E[E[X|Y]]. The outer expectation averages the conditional mean over the distribution of Y. This identity enables two-stage computation: first condition on Y, then average out.

Conditional Variance Var(X|Y)

Var(X|Y) = E[X²|Y] − (E[X|Y])² measures the residual variability in X that remains after observing Y. Like the conditional mean, it is a random variable—a function of Y.

Law of Total Variance (Eve's Law)

Var(X) = E[Var(X|Y)] + Var(E[X|Y]). Total variance decomposes into the mean of within-group variance plus the variance of group means, providing both a computational tool and a conceptual decomposition.

Strategy: Choose Y Wisely

The art of conditioning lies in selecting the right auxiliary variable Y so that E[X|Y] and Var(X|Y) become recognizable distributions. Good choices include mixture parameters, stopping times, or latent variables.

✦

KEY TAKEAWAY

Think of conditioning like a two-phase survey. Suppose you want the average height across an entire university. Rather than measuring everyone at once, you first compute the average height within each department (the conditional mean), then take a weighted average across departments (the outer expectation). The law of total expectation guarantees this two-phase approach yields the correct overall average. The law of total variance goes further: total variability equals the average within-department spread plus the variability of department averages—neatly separating 'within-group' and 'between-group' variation.

SECTION 3

Visual Explanation — The Conditioning Decomposition

The upper half illustrates the law of total expectation: start from the joint distribution, condition on each value of Y to obtain conditional means μₖ, then average over Y to recover E[X]. The lower half shows the total variance decomposition into within-group (purple) and between-group (pink) components.

The diagram above captures the essence of computing via conditioning. In the top half, observe how the joint distribution f(x, y) is partitioned by conditioning on different realized values of Y. Each conditional slice yields its own mean μₖ = E[X|Y = yₖ], and these conditional means are themselves random—they vary as Y varies. The final step applies the ordinary expectation operator to this random function g(Y) = E[X|Y], recovering the unconditional mean. The bottom half extends this logic to variance: the total variance of X is the average within-group variance (how much X spreads around its conditional mean for a fixed Y) plus the between-group variance (how much the conditional means themselves fluctuate as Y changes). This decomposition is both computationally useful and conceptually illuminating.

SECTION 4

Mathematical Framework

Law of Total Expectation — Derivation

We begin with the discrete case for clarity. Let X and Y be jointly distributed discrete random variables. By the definition of marginal expectation and applying the law of total probability, we can decompose E[X] into a double sum that naturally factors via conditional distributions.

LAW OF TOTAL EXPECTATION (DISCRETE)

E[X] = Σ_y E[X | Y = y] · P(Y = y) = E_Y[E[X | Y]]

Proof sketch: E[X] = Σ_x x · P(X = x) = Σ_x x Σ_y P(X = x, Y = y) = Σ_y Σ_x x · P(X = x | Y = y) · P(Y = y) = Σ_y E[X|Y = y] · P(Y = y). The interchange of summation is justified by absolute convergence.

CONTINUOUS ANALOGUE

E[X] = ∫_{−∞}^{∞} E[X | Y = y] · f_Y(y) dy

Here $f_Y$ (y) is the marginal density of Y. The conditional expectation E[X|Y = y] = ∫ x · $f_{X|Y}$ (x|y) dx is evaluated for each fixed y, then integrated against the marginal density of Y.

Law of Total Variance — Derivation

The law of total variance follows from the definition of variance and two applications of the tower property. Write Var(X) = E[X²] − (E[X])². We condition on Y to expand each term and then carefully rearrange.

LAW OF TOTAL VARIANCE (EVE'S LAW)

Var(X) = E[Var(X | Y)] + Var(E[X | Y])

💡

Mnemonic: EVVE

A popular mnemonic for the law of total variance is EVVE: E of the V plus V of the E. Here, the first 'E' and 'V' refer to E[Var(X|Y)] (expected value of the conditional variance), while the second 'V' and 'E' refer to Var(E[X|Y]) (variance of the conditional expectation). The alternating E-V pattern makes it easy to recall.

TOWER PROPERTY (GENERAL FORM)

E[E[X | 𝒢] | ℋ] = E[X | ℋ] whenever ℋ ⊆ 𝒢

In the special case ℋ = {∅, Ω} (the trivial σ-algebra), this reduces to E[E[X|Y]] = E[X]. The tower property is the foundational identity from which both the law of total expectation and the law of total variance are derived.

SECTION 5

Techniques & Common Conditioning Strategies

The power of computing via conditioning lies in the choice of the conditioning variable Y. Different problem structures suggest different natural choices, and recognizing these patterns is a skill honed through practice. Below we catalog the most frequently encountered strategies, along with a visual decision guide.

This decision tree guides the choice of conditioning variable. For mixture/hierarchical problems, condition on the random parameter Θ. For random sums S = X₁ + ⋯ +

X_N

, condition on the random count N. Other patterns include first-step analysis in Markov chains and indicator-based decompositions.

Strategy	Condition on Y =	Typical Problem Class	Example
Mixture parameter	Θ (latent parameter)	Bayesian priors, overdispersed models	X\|Θ ~ Bin(n, Θ), Θ ~ Beta(a, b)
Random count	N (number of terms)	Compound Poisson, random walks	S = X₁ + ⋯ + $X_N$ , N ~ Poisson(λ)
First step	Outcome of first transition	Gambler's ruin, Markov chain hitting times	E[T] = 1 + p·E[T\|up] + q·E[T\|down]
Indicator event	Success/failure of a sub-event	Coupon collector, matching problems	E[X] = E[X\|A]P(A) + E[X\|Aᶜ]P(Aᶜ)

Common conditioning strategies and their typical application domains

SECTION 6

Worked Example — Random Sum with a Poisson Count

An insurance company models the total claims S in a given month as a random sum S = X₁ + X₂ + ⋯ + $X_N$ , where N ~ Poisson(λ = 10) is the number of claims and each claim size Xᵢ is independent and identically distributed with E[Xᵢ] = 500 and Var(Xᵢ) = 40,000. The Xᵢ are also independent of N. Compute E[S] and Var(S).

Computing E[S] and Var(S) via Conditioning on N

1
Step 1 — Identify the Conditioning Variable

The total claims S depends on the random count N. Conditioning on N transforms S into a fixed sum of i.i.d. random variables, which is far easier to analyze. We set Y = N and apply the laws of total expectation and total variance.
2
Step 2 — Compute E[S|N]

Given N = n, the sum S = X₁ + ⋯ + $X_n$ is a sum of n i.i.d. random variables each with mean μ = E[Xᵢ] = 500. Therefore E[S|N = n] = n × 500, or as a random variable, E[S|N] = 500N.

E[S|N] = 500N
3
Step 3 — Apply the Law of Total Expectation

E[S] = E[E[S|N]] = E[500N] = 500 × E[N] = 500 × 10 = 5,000.

E[S] = 5,000
4
Step 4 — Compute Var(S|N)

Given N = n, the sum S has variance n × Var(Xᵢ) = n × 40,000 (by independence of the Xᵢ). As a random variable: Var(S|N) = 40,000 × N.

Var(S|N) = 40,000N
5
Step 5 — Apply the Law of Total Variance (Eve's Law)

Var(S) = E[Var(S|N)] + Var(E[S|N]). The first term: E[40,000N] = 40,000 × E[N] = 40,000 × 10 = 400,000. The second term: Var(500N) = 500² × Var(N) = 250,000 × 10 = 2,500,000. (We used the fact that for N ~ Poisson(λ), Var(N) = λ = 10.) Therefore Var(S) = 400,000 + 2,500,000 = 2,900,000.

Var(S) = 2,900,000
6
Step 6 — Interpretation

The within-group variance E[Var(S|N)] = 400,000 measures the variability due to randomness in individual claim sizes for a fixed number of claims. The between-group variance Var(E[S|N]) = 2,500,000 measures the variability arising from the uncertain number of claims. In this example, the uncertainty in N contributes over 86% of the total variance, highlighting that the random count dominates the overall risk.

SECTION 7

Strengths, Limitations & Comparison with Direct Methods

Computing via conditioning is not the only way to find expectations and variances. Direct integration, moment-generating functions, and characteristic functions offer alternative paths. The table below compares conditioning with these direct approaches so that you can appreciate when conditioning shines and when other methods may be preferable.

Criterion	Conditioning Approach	Direct Integration / MGFs
Ease of setup	Often only requires identifying E[X\|Y] and Var(X\|Y) for a well-known conditional distribution—no explicit joint density needed.	Requires the full joint density or explicit MGF, which may be difficult to derive for compound or hierarchical models.
Interpretability	Naturally decomposes variance into within-group and between-group components, providing deep insight into the sources of variability.	Yields a single number; the decomposition is not automatically provided.
Higher moments	Extending to E[X³] or skewness via conditioning is possible but increasingly cumbersome with each additional moment.	MGFs and characteristic functions handle all moments uniformly through differentiation.
Distributional results	Conditioning yields moments but not the full distribution of X (without further work).	MGFs/CFs can identify the full distribution through uniqueness theorems.
Recursive problems	Excels: first-step analysis, gambler's ruin, and Markov chain expectations reduce to elegant linear systems.	Direct methods struggle with recursive or sequential structure.

Conditioning vs. direct methods for computing moments

🎯

WHEN TO CONDITION

Conditioning is your go-to tool whenever the problem exhibits a natural hierarchical or sequential structure—when X's distribution is determined by some intermediate random quantity Y. If no such structure exists and you have the joint density in closed form, direct integration or MGFs may be more efficient. In practice, many real-world problems (insurance, queuing, Bayesian models) are inherently hierarchical, making conditioning the dominant approach in applied probability.

SECTION 8

Connection to Advanced Theory

The conditioning identities we have explored are special cases of much deeper results in measure-theoretic probability. In the general setting, conditional expectation E[X|𝒢] is defined as a 𝒢-measurable random variable satisfying ∫_A E[X|𝒢] dP = ∫_A X dP for all A ∈ 𝒢. This abstract definition, due to Kolmogorov and refined by Doob, generalizes our formulas from conditioning on a random variable Y to conditioning on an entire σ-algebra 𝒢, enabling applications in filtration-based frameworks such as martingale theory and stochastic calculus.

Concept in This Lesson	Generalization	Application Domain
E[E[X\|Y]] = E[X]	Tower property: E[E[X\|𝒢]\|ℋ] = E[X\|ℋ] for ℋ ⊆ 𝒢	Martingale convergence, Doob decomposition
Var(X) = E[Var(X\|Y)] + Var(E[X\|Y])	ANOVA decomposition into within/between components for multiple nested factors	Hierarchical Bayesian models, mixed-effects models
Conditioning on N in random sums	Wald's identities: E[ $S_N$ ] = E[N]·E[X₁] for stopping times N	Sequential analysis, renewal theory
E[X\|Y] improves estimation of E[X]	Rao–Blackwell theorem: E[T\|S] has ≤ variance than T	Sufficient statistics, UMVUE construction

From conditioning basics to advanced probability and statistics

The Rao–Blackwell theorem deserves special emphasis: it states that if T is an unbiased estimator of θ and S is a sufficient statistic, then T* = E[T|S] is also unbiased and satisfies Var(T*) ≤ Var(T). The proof is a direct application of Eve's law. Since Var(T) = E[Var(T|S)] + Var(E[T|S]) = E[Var(T|S)] + Var(T*), and E[Var(T|S)] ≥ 0, we conclude Var(T*) ≤ Var(T). This beautiful result shows that conditioning on more information never increases mean squared error—a principle that pervades modern statistics.

SECTION 9

Practice Problems

Question 1: CONCEPTUAL

Explain in your own words why E[X|Y] is a random variable rather than a fixed number. What is the source of its randomness? Then state, without proof, how the law of total expectation uses this fact to recover E[X].

Show answer

E[X|Y = y] is a number for each fixed value y, but since Y itself is random, the quantity E[X|Y] = g(Y) inherits randomness from Y. Its distribution is determined by the distribution of Y through the function g. The law of total expectation states that E[X] = E[E[X|Y]] = E[g(Y)], meaning we compute the ordinary expectation of this random variable g(Y) over the distribution of Y. This 'averages out' the conditioning, recovering the unconditional mean.

Question 2: BASIC CALCULATION

Let X|P ~ Binomial(20, P) where P ~ Uniform(0, 1). Use the law of total expectation to find E[X], and then apply Eve's law to find Var(X).

Show answer

First, E[X|P] = 20P, so E[X] = E[20P] = 20 × E[P] = 20 × (1/2) = 10. For the variance, Var(X|P) = 20P(1 − P), so E[Var(X|P)] = 20E[P(1 − P)] = 20(E[P] − E[P²]) = 20(1/2 − 1/3) = 20/6 = 10/3. Also, Var(E[X|P]) = Var(20P) = 400 × Var(P) = 400 × (1/12) = 100/3. By Eve's law, Var(X) = 10/3 + 100/3 = 110/3 ≈ 36.67.

Question 3: INTERMEDIATE

A fair coin is flipped, and the outcome determines a distribution: if Heads, X ~ Exponential(1); if Tails, X ~ Exponential(1/3). Find E[X] and Var(X) using conditioning on the coin outcome Y ∈ {H, T}.

Show answer

Let Y = H or T with P(H) = P(T) = 1/2. Then E[X|Y = H] = 1 and E[X|Y = T] = 3 (means of Exp(1) and Exp(1/3) respectively). So E[X] = (1/2)(1) + (1/2)(3) = 2. For variance: Var(X|Y = H) = 1 and Var(X|Y = T) = 9. E[Var(X|Y)] = (1/2)(1) + (1/2)(9) = 5. Var(E[X|Y]) = E[(E[X|Y])²] − (E[X])² = (1/2)(1) + (1/2)(9) − 4 = 5 − 4 = 1. Thus Var(X) = 5 + 1 = 6. Note that the between-group term (1) and within-group term (5) reveal that most variability comes from the exponential distributions themselves rather than from the coin flip.

Question 4: APPLIED

A call center receives N calls per hour, where N ~ Poisson(λ = 8). Each call lasts a random duration Dᵢ (in minutes), independently drawn from a distribution with E[Dᵢ] = 6 and Var(Dᵢ) = 4. The total agent-time used per hour is T = D₁ + D₂ + ⋯ +

D_N

. Find E[T] and Var(T). If the center has 60 agent-minutes available, find the probability that demand exceeds capacity using a normal approximation.

Show answer

Condition on N. E[T|N] = 6N and Var(T|N) = 4N. By the law of total expectation, E[T] = 6 × E[N] = 6 × 8 = 48 minutes. By Eve's law, Var(T) = E[4N] + Var(6N) = 4 × 8 + 36 × 8 = 32 + 288 = 320. So SD(T) = √320 ≈ 17.89. Using the normal approximation T ≈ N(48, 320), P(T > 60) = P(Z > (60 − 48)/17.89) = P(Z > 0.670) ≈ 1 − 0.749 = 0.251. There is roughly a 25% chance the center exceeds capacity in any given hour—a significant risk that would inform staffing decisions.

Question 5: CRITICAL THINKING

Prove that for any random variable X and any random variable Y (with finite second moments), conditioning on Y can never increase the expected squared error of prediction. Specifically, show that E[(X − E[X])²] ≥ E[(X − E[X|Y])²], and interpret the difference between these two quantities using the law of total variance.

Show answer

We need to show Var(X) ≥ E[(X − E[X|Y])²]. Note that E[(X − E[X|Y])²] = E[E[(X − E[X|Y])² | Y]] = E[Var(X|Y)] by definition of conditional variance. From Eve's law, Var(X) = E[Var(X|Y)] + Var(E[X|Y]). Since Var(E[X|Y]) ≥ 0, we immediately obtain Var(X) ≥ E[Var(X|Y)] = E[(X − E[X|Y])²]. The difference Var(X) − E[(X − E[X|Y])²] = Var(E[X|Y]) is exactly the between-group variance—the amount of variability in X that is 'explained' by Y. This result is the theoretical foundation of R² in regression: R² = Var(E[X|Y]) / Var(X) measures the proportion of variance in X explained by Y. It also underpins the Rao–Blackwell theorem, since replacing a predictor of X by its conditional expectation given more information Y always reduces (or maintains) the expected squared error.

SUMMARY

Summary — Computing via Conditioning

Summary

Computing via conditioning is a cornerstone technique in probability that leverages the law of total expectation E[X] = E[E[X|Y]] and the law of total variance Var(X) = E[Var(X|Y)] + Var(E[X|Y]) to decompose intractable moment calculations into manageable conditional pieces. The central insight is that conditional expectation E[X|Y] is itself a random variable—a function of Y—whose own moments can be computed by averaging over the distribution of Y. The variance decomposition separates total variability into a within-group component E[Var(X|Y)] and a between-group component Var(E[X|Y]), providing both computational efficiency and deep interpretive insight.

The art of the method lies in choosing the right conditioning variable Y so that the conditional distribution of X given Y is a well-known, tractable family. Common strategies include conditioning on mixture parameters in hierarchical models, random counts in compound distributions, and first-step outcomes in sequential and Markov chain problems. These identities generalize to the measure-theoretic tower property and underpin advanced results including the Rao–Blackwell theorem, ANOVA decompositions, and variance reduction in Monte Carlo methods.