Information Theory Inequalities

Jensen's Inequality, Log Sum Inequality, Data Processing Inequality and Sufficient Statistics, and Fano's Inequality

· 32 min read · 6,397 words

Jensen’s inequality

Motivation: Whenever we have a random variable XX, we often need to look at a function of that variable, f(X)f(X). If the function is perfectly linear (like f(x)=2xf(x) = 2x), the average of the function’s outputs exactly equals the function evaluated at the average input: E[f(X)]=f(E[X])\mathbb{E}[f(X)] = f(\mathbb{E}[X]). However, in nearly all branches of engineering and mathematics, systems are non-linear. In information theory, our entire framework is built on the highly non-linear logarithm function (f(x)=logxf(x) = \log x). Because a curve bends, the “average of the curve at points” is not equal to the “curve of the average of those points.” We need a law that tells us exactly how a function’s geometric curvature biases its expected value. This is Jensen’s Inequality. By bounding how expected values behave inside non-linear functions, Jensen’s Inequality allows us to prove some limits of our communication systems.

Convex and concave functions

To be able to reason about functions on random variables in non-linear systems, we need to define the curvature of the functions more rigorously so we can analyze it.

DefinitionConvexity and concavity (geometric)

A function f(x)f(x) is convex over an interval (a,b)(a, b) if you pick any two points on its curve and draw a straight line (a chord) between them, and that chord always lies strictly above or on the curve. A function is concave if the chord always lies strictly below or on the curve. (Note: f(x)f(x) is concave if f(x)-f(x) is convex).

DefinitionConvexity (algebraic)

Algebraically, drawing a chord means taking a weighted average of two inputs x1x_1 and x2x_2, using a weight λ\lambda where 0λ10 \le \lambda \le 1.

A function ff is convex if the function evaluated at the average of the inputs is less than or equal to the average of the function’s outputs:

f(λx1+(1λ)x2)λf(x1)+(1λ)f(x2)f(\lambda x_1 + (1 - \lambda)x_2) \le \lambda f(x_1) + (1 - \lambda)f(x_2)

DefinitionStrict convexity

A function is strictly convex if equality holds only when λ=0\lambda = 0, λ=1\lambda = 1, or x1=x2x_1 = x_2. (The chord never touches the curve except at the endpoints).

TheoremThe second derivative test

If a function ff has a second derivative that is non-negative (f0f'' \ge 0) over an interval, the function is convex. If it is strictly positive (f>0f'' > 0), it is strictly convex.

ProofUsing Taylor series

Expand f(x)f(x) around a point x0x_0:

f(x)=f(x0)+f(x0)(xx0)+f(x)2(xx0)2f(x) = f(x_0) + f'(x_0)(x - x_0) + \frac{f''(x^*)}{2}(x - x_0)^2

Since f(x)0f''(x^*) \ge 0, the last term is always non-negative. Therefore, the curve always sits above its tangent line: f(x)f(x0)+f(x0)(xx0)f(x) \ge f(x_0) + f'(x_0)(x - x_0). If we let x0=λx1+(1λ)x2x_0 = \lambda x_1 + (1-\lambda)x_2 and plug in x=x1x = x_1 and x=x2x = x_2 respectively, multiplying the resulting inequalities by λ\lambda and (1λ)(1-\lambda) and adding them together algebraically yields the exact definition of convexity.

Applications to information theory: f(x)=x2,exf(x) = x^2, e^x, and xlogxx \log x (for x0x \ge 0) have positive second derivatives, making them strictly convex. f(x)=logxf(x) = \log x has a second derivative of 1/x2-1/x^2. Since probability x>0x > 0, this is strictly negative. Therefore, logx\log x is strictly concave and logx-\log x is strictly convex. The use case is usually to move both the summation and the probability inside the log in all the information theoretic terms.

We know a chord between two points on a convex curve sits above the curve. Jensen’s Inequality generalizes this algebraic definition from two points to an entire random variable.

TheoremJensen’s inequality

If ff is a convex function and XX is a random variable, then:

E[f(X)]f(E[X])\mathbb{E}[f(X)] \ge f(\mathbb{E}[X])

(If ff is strictly convex, equality implies X=E[X]X = \mathbb{E}[X] with probability 1, meaning XX is a deterministic constant with no randomness).

ProofBy induction on discrete distributions

Base case (2 mass points): Let XX take values x1,x2x_1, x_2 with probabilities p1,p2p_1, p_2. By the literal definition of convexity (replacing λ\lambda with p1p_1):

p1f(x1)+p2f(x2)f(p1x1+p2x2)p_1 f(x_1) + p_2 f(x_2) \ge f(p_1 x_1 + p_2 x_2)

This is exactly E[f(X)]f(E[X])\mathbb{E}[f(X)] \ge f(\mathbb{E}[X]).

Inductive step: Assume the theorem holds for k1k-1 mass points. We prove it for kk points. Isolate the kk-th term, and re-weight the remaining k1k-1 probabilities by defining pi=pi/(1pk)p_i' = p_i / (1 - p_k):

i=1kpif(xi)=pkf(xk)+(1pk)i=1k1pif(xi)\sum_{i=1}^k p_i f(x_i) = p_k f(x_k) + (1 - p_k) \sum_{i=1}^{k-1} p_i' f(x_i)

Apply the inductive hypothesis to the inner summation (since it has k1k-1 terms), pulling ff outside:

pkf(xk)+(1pk)f(i=1k1pixi)\ge p_k f(x_k) + (1 - p_k) f \left( \sum_{i=1}^{k-1} p_i' x_i \right)

We now have exactly two weighted terms: pkp_k and (1pk)(1-p_k). Apply the base-case convexity definition one final time to pull ff outside everything:

f(pkxk+(1pk)i=1k1pixi)=f(i=1kpixi)=f(E[X])\ge f \left( p_k x_k + (1 - p_k) \sum_{i=1}^{k-1} p_i' x_i \right) = f \left( \sum_{i=1}^k p_i x_i \right) = f(\mathbb{E}[X])

Consequences of Jensen’s inequality

By passing our information measures through Jensen’s Inequality, we establish the fundamental physical bounds of communication theory.

TheoremInformation inequality (Gibbs’ inequality)

Let p(x)p(x) and q(x)q(x) be two probability mass functions.

D(pq)0D(p\vert{}\vert{}q) \ge 0

(Equality holds if and only if p(x)=q(x)p(x) = q(x) for all xx).

Intuition

You can never build an encoding mechanism for a general language and have it be more efficient than a language-specific encoding mechanism by chance.

ProofLet AA be the support set where p(x)>0p(x) > 0

D(pq)=xAp(x)logp(x)q(x)=xAp(x)logq(x)p(x)-D(p\vert{}\vert{}q) = - \sum_{x \in A} p(x) \log \frac{p(x)}{q(x)} = \sum_{x \in A} p(x) \log \frac{q(x)}{p(x)}

Because log(t)\log(t) is strictly concave, E[log(X)]log(E[X])\mathbb{E}[\log(X)] \le \log(\mathbb{E}[X]). We apply this by bringing the logarithm outside the summation:

log(xAp(x)q(x)p(x))=log(xAq(x))\le \log \left( \sum_{x \in A} p(x) \frac{q(x)}{p(x)} \right) = \log \left( \sum_{x \in A} q(x) \right)

Because probabilities must sum to 1 (or less, if AA is a subset), xAq(x)1\sum_{x \in A} q(x) \le 1.

log(1)=0\le \log(1) = 0

Therefore, D(pq)0-D(p\vert{}\vert{}q) \le 0, meaning D(pq)0D(p\vert{}\vert{}q) \ge 0.

CorollaryNon-negativity of mutual information

I(X;Y)0I(X; Y) \ge 0

Intuition

You can never have the wire noise be something where the receiver has less information about the sender’s message than random noise.

Proof

Since I(X;Y)=D(p(x,y)p(x)p(y))I(X;Y) = D(p(x,y) \vert{}\vert{} p(x)p(y)), and all relative entropies are 0\ge 0, mutual information is 0\ge 0. (Equality holds iff XX and YY are independent).

TheoremMaximum entropy bound

H(X)logXH(X) \le \log \vert{}\mathcal{X}\vert{}

(Entropy is maximized when every symbol is equally likely, i.e., a uniform distribution).

Proof

Let u(x)=1/Xu(x) = 1/\vert{}\mathcal{X}\vert{} be the uniform distribution Calculate D(pu)=p(x)logp(x)1/X=p(x)logp(x)+p(x)logX=H(X)+logXD(p\vert{}\vert{}u) = \sum p(x) \log \frac{p(x)}{1/\vert{}\mathcal{X}\vert{}} = \sum p(x) \log p(x) + \sum p(x) \log \vert{}\mathcal{X}\vert{} = -H(X) + \log \vert{}\mathcal{X}\vert{}.

Since D(pu)0D(p\vert{}\vert{}u) \ge 0, we have logXH(X)0\log \vert{}\mathcal{X}\vert{} - H(X) \ge 0.

TheoremConditioning reduces entropy (information can’t hurt)

H(XY)H(X)H(X\vert{}Y) \le H(X)

Intuition

Knowing a second variable YY can only decrease, never increase, your overall uncertainty about XX.

Proof

We know I(X;Y)=H(X)H(XY)I(X;Y) = H(X) - H(X\vert{}Y). We also know I(X;Y)0I(X;Y) \ge 0. Therefore, H(X)H(XY)0H(X) - H(X\vert{}Y) \ge 0.

WarningImportant caveat

This is true on average. A specific, highly unusual observation (e.g., Y=yY=y) might temporarily increase your uncertainty about XX (H(XY=y)>H(X)H(X\vert{}Y=y) > H(X)). However, when weighted across all possible observations of YY, the expected uncertainty mathematically must drop or remain equal.

Medical diagnosis example:

  • Let X{Healthy,Sick}X \in \{\text{Healthy}, \text{Sick}\} be a patient’s true state. A disease is very rare so prior probability is P(X=Sick)=0.01P(X=\text{Sick}) = 0.01 and P(X=Healthy)=0.99P(X=\text{Healthy}) = 0.99.

  • Because you are almost certain the patient is healthy, baseline uncertainty is tiny. H(X)=(0.01log20.01+0.99log20.99)0.081 bitsH(X) = -(0.01 \log_2 0.01 + 0.99 \log_2 0.99) \approx 0.081\text{ bits}

  • You run a lab test Y{Negative,Positive}Y \in \{\text{Negative}, \text{Positive}\}.

  • The test is imperfect. If it comes back Positive (Y=posY = \text{pos}), the posterior probability updates to an ambiguous coin-flip: P(X=SickY=pos)=0.50P(X=\text{Sick} \mid Y=\text{pos}) = 0.50 and P(X=HealthyY=pos)=0.50P(X=\text{Healthy} \mid Y=\text{pos}) = 0.50.

  • For the specific observation (Y=posY=\text{pos}): H(XY=pos)=(0.5log20.5+0.5log20.5)=1.0 bitH(X \mid Y=\text{pos}) = -(0.5 \log_2 0.5 + 0.5 \log_2 0.5) = 1.0\text{ bit} A positive result massively increases your confusion from 0.0810.081 to 1.0 bit1.0\text{ bit}. A calm default of most people being healthy is now a state of total doubt about people after getting a positive on the test.

  • On average across all test outcomes (H(XY)H(X \mid Y)): Most people test negative (Y=negY = \text{neg}), which drives P(X=SickY=neg)0P(X=\text{Sick} \mid Y=\text{neg}) \approx 0, collapsing H(XY=neg)0H(X \mid Y=\text{neg}) \approx 0. When weighted together: H(XY)=yP(Y=y)H(XY=y)H(X)H(X \mid Y) = \sum_y P(Y=y) H(X \mid Y=y) \le H(X) The overwhelming certainty gained on negative tests balances out the rare spike in doubt, ensuring that conditioning reduces entropy on average.

TheoremIndependence bound on entropy

H(X1,X2,,Xn)i=1nH(Xi)H(X_1, X_2, \ldots, X_n) \le \sum_{i=1}^n H(X_i)

This means that the total uncertainty of a system is less than or equal to the sum of its individual parts. Equality holds only if all variables are completely independent.

ProofApply the Chain Rule for Entropy

H(X1,,Xn)=i=1nH(XiXi1,,X1)H(X_1, \ldots, X_n) = \sum_{i=1}^n H(X_i \vert{} X_{i-1}, \ldots, X_1). By theorem above, conditioning reduces entropy, so H(XiXi1,,X1)H(Xi)H(X_i \vert{} X_{i-1}, \ldots, X_1) \le H(X_i) for every term. Substituting this bound into the summation proves the theorem.

The log sum inequality

We have used Jensen’s Inequality to prove the limits of a static communication system (e.g., the penalty of a wrong codebook is always 0\ge 0). But real telegraph networks are dynamic.

Suppose your telegraph operator operates smoothly in the morning (with source probability p1p_1) and gets highly erratic at night due to drinking (with source probability p2p_2). If an engineer is lazy and tries to build a single “daily average” model by mathematically blending the morning and night probabilities together, do they gain or lose precision?

To prove how uncertainty and penalties behave when we mix distributions, we need a tool that handles the sum of blended ratios. This is the Log Sum Inequality. It is the formal proof that blending distinct signals always destroys information.

TheoremThe log sum inequality

For any two sets of non-negative numbers, a1,a2,,ana_1, a_2, \ldots, a_n and b1,b2,,bnb_1, b_2, \ldots, b_n:

i=1nailogaibi(i=1nai)logi=1naii=1nbi\sum_{i=1}^n a_i \log \frac{a_i}{b_i} \ge \left( \sum_{i=1}^n a_i \right) \log \frac{\sum_{i=1}^n a_i}{\sum_{i=1}^n b_i}

Conventionally, 0log0=00 \log 0 = 0, alog(a/0)=a \log (a/0) = \infty, and 0log(0/0)=00 \log (0/0) = 0 based on continuity.

IntuitionWhy this solves the blending problem

The expression alog(a/b)a \log(a / b) measures the mismatch cost between reality (aa) and an assumption (bb). When reality spikes higher than the assumed capacity, the ratio a/b>1a/b > 1 makes the log term positive, increasing the wasted transmission overhead.

Suppose traffic surges irregularly across different telegraph routes or shifts:

  • Shift 1: Actual traffic volume is a1a_1, but your planned capacity was b1b_1.
  • Shift 2: Actual traffic volume is a2a_2, but your planned capacity was b2b_2.

If you evaluate each shift individually and add up the mismatch costs, you get the left side of the inequality:

Isolated Overhead=a1loga1b1+a2loga2b2\text{Isolated Overhead} = a_1 \log \frac{a_1}{b_1} + a_2 \log \frac{a_2}{b_2}

But suppose an engineer gets lazy and looks only at the daily totals: total real traffic A=a1+a2A = a_1 + a_2 versus total planned capacity B=b1+b2B = b_1 + b_2. The mismatch cost calculated on the bulk totals is the right side:

Blended Overhead=AlogAB=(a1+a2)loga1+a2b1+b2\text{Blended Overhead} = A \log \frac{A}{B} = (a_1 + a_2) \log \frac{a_1 + a_2}{b_1 + b_2}

The Log Sum Inequality guarantees:

ailogaibiAlogAB\sum a_i \log \frac{a_i}{b_i} \ge A \log \frac{A}{B}

Why the blended formula hides the problem: The function tlogtt \log t is convex, it curves upward sharply. High ratios blow up the penalty much faster than low ratios reduce it.

When you evaluate shifts in isolation, any extreme bottleneck (a1b1a_1 \gg b_1) triggers a massive penalty that dominates the sum. But when you sum the numbers into bulk totals (AA and BB) before taking the ratio, an overloaded morning (a1/b1=4a_1 / b_1 = 4) algebraically cancels out an empty evening (a2/b2=0.25a_2 / b_2 = 0.25).

Looking at blended totals creates an illusion of balance (A/B1A / B \approx 1), mathematically masking localized congestion. The Log Sum Inequality proves that aggregating data across time or lines strictly hides mismatch peaks, meaning an engineer examining blended averages will always underestimate the actual physical penalties paid during peak imbalances.

ProofVia Jensen’s inequality

To prove this, we first need to build a strictly convex function.

Let f(t)=tlogtf(t) = t \log t.

First derivative: f(t)=logt+1f'(t) = \log t + 1 (if base ee; or 1ln2log2t+1ln2\frac{1}{\ln 2} \log_2 t + \frac{1}{\ln 2}).

Second derivative: f(t)=1/tf''(t) = 1/t (times a positive constant depending on the log base).

Since our inputs tt will be strictly positive, 1/t>01/t > 0. Because the second derivative is strictly positive, the function f(t)=tlogtf(t) = t \log t is strictly convex.

Now, we set up Jensen’s Inequality: λif(ti)f(λiti)\sum \lambda_i f(t_i) \ge f(\sum \lambda_i t_i).

To use Jensen’s, we need our weights λi\lambda_i to act like a valid probability distribution (they must sum to 1).

Let A=aiA = \sum a_i and B=biB = \sum b_i.

We construct our valid probability weights using the bb terms: λi=biB\lambda_i = \frac{b_i}{B}.

We will feed the function the ratio of our terms: ti=aibit_i = \frac{a_i}{b_i}.

Evaluate the left side of Jensen’s (the expected value of the function).

iλif(ti)=ibiB(aibilogaibi)\sum_i \lambda_i f(t_i) = \sum_i \frac{b_i}{B} \left( \frac{a_i}{b_i} \log \frac{a_i}{b_i} \right)

The bib_i in the numerator and denominator cancel out, leaving:

=1Biailogaibi= \frac{1}{B} \sum_i a_i \log \frac{a_i}{b_i}

Evaluate the right side of Jensen’s (the function of the expected value). First, find the expected value inside the function:

iλiti=ibiBaibi=1Biai=AB\sum_i \lambda_i t_i = \sum_i \frac{b_i}{B} \frac{a_i}{b_i} = \frac{1}{B} \sum_i a_i = \frac{A}{B}

Now apply the function f(t)=tlogtf(t) = t \log t to that result:

f(AB)=ABlogABf\left(\frac{A}{B}\right) = \frac{A}{B} \log \frac{A}{B}

Bring them together. Because f(t)f(t) is convex, Left Side \ge Right Side.

1BiailogaibiABlogAB\frac{1}{B} \sum_i a_i \log \frac{a_i}{b_i} \ge \frac{A}{B} \log \frac{A}{B}

Multiply both sides by BB, and you get the Log Sum Inequality:

i=1nailogaibiAlogAB\sum_{i=1}^n a_i \log \frac{a_i}{b_i} \ge A \log \frac{A}{B}

Consequences: applying it to telecom

The Log Sum Inequality looks like abstract algebra, but substituting p(x)p(x) and q(x)q(x) for aa and bb instantly proves two physical properties of communication networks.

Convexity of relative entropy

CorollaryConvexity of relative entropy

If you mix two different communication environments together, the penalty of the blended system is less than or equal to the average penalty of the individual systems.

D(λp1+(1λ)p2λq1+(1λ)q2)λD(p1q1)+(1λ)D(p2q2)D(\lambda p_1 + (1-\lambda)p_2 \vert{}\vert{} \lambda q_1 + (1-\lambda)q_2) \le \lambda D(p_1\vert{}\vert{}q_1) + (1-\lambda) D(p_2\vert{}\vert{}q_2)

Proof

We want to evaluate the blended Relative Entropy for a specific symbol xx. We map our variables to the Log Sum Inequality:

Let a1=λp1(x)a_1 = \lambda p_1(x) and a2=(1λ)p2(x)a_2 = (1-\lambda)p_2(x). (These sum to the blended reality).

Let b1=λq1(x)b_1 = \lambda q_1(x) and b2=(1λ)q2(x)b_2 = (1-\lambda)q_2(x). (These sum to the blended assumption).

Applying the Log Sum Inequality to these two terms:

a1loga1b1+a2loga2b2(a1+a2)loga1+a2b1+b2a_1 \log \frac{a_1}{b_1} + a_2 \log \frac{a_2}{b_2} \ge (a_1 + a_2) \log \frac{a_1 + a_2}{b_1 + b_2}

Substituting our probabilities back in:

λp1(x)logλp1(x)λq1(x)+(1λ)p2(x)log(1λ)p2(x)(1λ)q2(x)(λp1(x)+(1λ)p2(x))logλp1(x)+(1λ)p2(x)λq1(x)+(1λ)q2(x)\lambda p_1(x) \log \frac{\lambda p_1(x)}{\lambda q_1(x)} + (1-\lambda) p_2(x) \log \frac{(1-\lambda) p_2(x)}{(1-\lambda) q_2(x)} \ge \big( \lambda p_1(x) + (1-\lambda)p_2(x) \big) \log \frac{\lambda p_1(x) + (1-\lambda)p_2(x)}{\lambda q_1(x) + (1-\lambda)q_2(x)}

Notice the weights λ\lambda and (1λ)(1-\lambda) cancel out inside the isolated logarithms on the left.

Summing this inequality over all symbols xXx \in \mathcal{X} yields:

λD(p1q1)+(1λ)D(p2q2)D(λp1+(1λ)p2λq1+(1λ)q2)\lambda D(p_1\vert{}\vert{}q_1) + (1-\lambda) D(p_2\vert{}\vert{}q_2) \ge D(\lambda p_1 + (1-\lambda)p_2 \vert{}\vert{} \lambda q_1 + (1-\lambda)q_2)

Intuition

Blending distributions washes out the specific statistical “fingerprints” that set them apart.

If you test two separate wire conditions against their respective models, each setup retains its distinct, localized anomalies, producing a high combined mismatch score. But if you throw all the traffic into a single bin (λp1+(1λ)p2\lambda p_1 + (1-\lambda) p_2) and compare it against a single combined baseline (λq1+(1λ)q2\lambda q_1 + (1-\lambda) q_2), the contrasting spikes average out across the aggregate pool.

Blending never sharpens distinctions, it blurs them. Pooling data across different environments strictly reduces the measurable divergence between reality and your baseline.

Concavity of entropy

Claim: Blending predictable systems creates an unpredictable system. “Averaging destroys certainty.”

CorollaryConcavity of entropy

H(λp1+(1λ)p2)λH(p1)+(1λ)H(p2)H(\lambda p_1 + (1-\lambda)p_2) \ge \lambda H(p_1) + (1-\lambda) H(p_2)

Proof

We can prove this instantly using the convexity of Relative Entropy we just derived. Entropy is simply the maximum possible entropy (logX\log \vert{}\mathcal{X}\vert{}) minus the relative entropy distance to the uniform distribution u(x)=1/Xu(x) = 1/\vert{}\mathcal{X}\vert{}:

H(p)=logXD(pu)H(p) = \log \vert{}\mathcal{X}\vert{} - D(p\vert{}\vert{}u)

We substitute this definition into the convexity inequality, setting both q1q_1 and q2q_2 to the uniform distribution uu:

D(λp1+(1λ)p2u)λD(p1u)+(1λ)D(p2u)D(\lambda p_1 + (1-\lambda)p_2 \vert{}\vert{} u) \le \lambda D(p_1\vert{}\vert{}u) + (1-\lambda) D(p_2\vert{}\vert{}u)

Replace the DD terms with their entropy equivalents:

logXH(λp1+(1λ)p2)λ(logXH(p1))+(1λ)(logXH(p2))\log \vert{}\mathcal{X}\vert{} - H(\lambda p_1 + (1-\lambda)p_2) \le \lambda \big(\log \vert{}\mathcal{X}\vert{} - H(p_1)\big) + (1-\lambda) \big(\log \vert{}\mathcal{X}\vert{} - H(p_2)\big)

The logX\log \vert{}\mathcal{X}\vert{} terms cleanly cancel out from both sides:

H(λp1+(1λ)p2)λH(p1)(1λ)H(p2)-H(\lambda p_1 + (1-\lambda)p_2) \le - \lambda H(p_1) - (1-\lambda) H(p_2)

Multiply by 1-1 and flip the inequality:

H(λp1+(1λ)p2)λH(p1)+(1λ)H(p2)H(\lambda p_1 + (1-\lambda)p_2) \ge \lambda H(p_1) + (1-\lambda) H(p_2)

Intuition

Imagine Operator 1 has a broken telegraph that only types ‘A’ (H(p1)=0H(p_1) = 0 bits). Operator 2 has a broken telegraph that only types ‘B’ (H(p2)=0H(p_2) = 0 bits). The average uncertainty of these two isolated systems is exactly 0. But suppose you walk into the telegraph office and you don’t know which operator is sitting at the desk today (a 50/50 blend, λ=0.5\lambda=0.5). The blended probability of getting an ‘A’ is now 50%, and ‘B’ is 50%. The entropy of this blended system instantly jumps to 1 bit. The math formally proves that when you blur distinct statistical signals together, you permanently destroy structure and inject uncertainty into the universe.

Data processing inequality

Suppose our sender XX transmits a message across a noisy copper wire, and the receiver gets a corrupted signal YY. The receiver decides to plug their machine into a secondary computer to “clean up” the signal, running it through a complex algorithm or a statistical filter to produce a final output ZZ.

We want to know: can this downstream processing step magically recover the original message? Can the final output ZZ have more mutual information with the sender XX than the raw, corrupted signal YY originally did?

To prove this, we must first mathematically define what it means to process a signal in a sequence. We do this using a Markov chain.

Markov chains (XYZX \rightarrow Y \rightarrow Z)

DefinitionMarkov chain

A sequence of three random variables forms a Markov chain (denoted XYZX \rightarrow Y \rightarrow Z) if the future depends on the past only through the present.

IntuitionPhysical intuition

The downstream computer ZZ is only plugged into the receiver YY. It has no secret, secondary wire connected directly to the sender XX. Therefore, everything ZZ knows about XX must have physically passed through YY first.

If Z depends on X only through Y, then if we completely freeze (condition on) Y, X and Z become mathematically independent. The joint probability distribution of a Markov chain factors cleanly into a chain of local conditional probabilities: this means the conditional probability of ZZ given the entire history is just the probability of ZZ given the immediate previous step:

FormulaThe Markov property

p(zx,y)=p(zy)p(z \vert{} x, y) = p(z \vert{} y)

Because knowing XX provides absolutely zero extra predictive power about ZZ once you already know YY, their conditional mutual information is zero: I(X;ZY)=0I(X ; Z \vert{} Y) = 0.

TheoremThe data processing inequality

If X→Y→Z forms a Markov chain, then the information that Z shares with X can never exceed the information that Y shares with X.

I(X;Y)I(X;Z)I(X;Y) \ge I(X;Z)

IntuitionPhysical meaning

No algorithm, function, or secondary physical process applied to a received signal can ever increase the amount of information it contains about the original source. Processing can only preserve information or destroy it.

ProofThrough the chain rule for mutual information

We want to calculate the total mutual information between the original sender XX and the entire receiver system (Y,Z)(Y, Z). We can expand this joint mutual information I(X;Y,Z)I(X; Y, Z) in two different directions.

Expansion 1 (start with YY). Using the chain rule, the information XX shares with the whole system is the information it shares with YY, plus the new information it shares with ZZ (given that we already know YY):

I(X;Y,Z)=I(X;Y)+I(X;ZY)I(X; Y, Z) = I(X; Y) + I(X; Z \vert{} Y)

Expansion 2 (start with ZZ). Alternatively, we can expand it backwards. The information XX shares with the whole system is the information it shares with ZZ, plus the new information it shares with YY (given that we already know ZZ):

I(X;Y,Z)=I(X;Z)+I(X;YZ)I(X; Y, Z) = I(X; Z) + I(X; Y \vert{} Z)

Equating and simplifying. Because both expansions equal the exact same total I(X;Y,Z)I(X;Y,Z), we set them equal to each other:

I(X;Y)+I(X;ZY)=I(X;Z)+I(X;YZ)I(X; Y) + I(X; Z \vert{} Y) = I(X; Z) + I(X; Y \vert{} Z)

Now, we apply our physical constraint. Because X → Y → Z is a Markov chain, we established that I(X;ZY)=0I(X; Z \vert{} Y) = 0. (Once you have the raw signal Y, the processed signal Z tells you absolutely nothing new about X). Substitute 0 into the left side of the equation:

I(X;Y)=I(X;Z)+I(X;YZ)I(X; Y) = I(X; Z) + I(X; Y \vert{} Z)

Finally, we proved via Jensen’s Inequality that all mutual information quantities are non-negative. Therefore, I(X;YZ)0I(X; Y \vert{} Z) \ge 0. If I(X;Y)I(X;Y) is equal to I(X;Z)I(X;Z) plus some positive number, then I(X;Y)I(X;Y) must be greater than or equal to I(X;Z)I(X;Z).

I(X;Y)I(X;Z)I(X; Y) \ge I(X; Z)

Consequences and corollaries

CorollaryDeterministic processing

What if the downstream processing isn’t a random variable, but a strict deterministic function (like a software algorithm)? Let Z=g(Y)Z = g(Y). Because ZZ is strictly calculated from YY, it perfectly forms the Markov chain XYg(Y)X \rightarrow Y \rightarrow g(Y). Therefore:

I(X;Y)I(X;g(Y))I(X; Y) \ge I(X; g(Y))

This proves that you cannot gain information by applying a function to your data. If you take a noisy audio file and run it through a complex equalizer function, the resulting file might sound better to a human, but mathematically, it contains less (or equal) information about the original studio recording than the raw file did.

When does equality hold?

The only way for I(X;Y)=I(X;Z)I(X; Y) = I(X; Z) to be true is if the leftover term I(X;YZ)I(X; Y \vert{} Z) is exactly 00. This means XZYX \rightarrow Z \rightarrow Y must also form a Markov chain. In engineering terms, this means the processing step ZZ perfectly preserved every single detail of YY that was relevant to XX. (In statistics, a function that achieves this equality is called a Sufficient Statistic).

Sufficient statistics

In our communication systems, we don’t always know the exact physical reality of the wire. Often, there is a hidden parameter θ\theta governing the system (e.g., the exact probability θ\theta that a static burst flips a bit). We cannot observe θ\theta directly; we can only observe a raw sequence of received data XX.

To save bandwidth and storage, an engineer might want to compress that massive raw dataset XX into a small, calculated summary value, T(X)T(X) — such as the sample mean or the total number of errors.

Because the parameter generates the data, and the data generates the summary, this forms a strict Markov chain:

θXT(X)\theta \rightarrow X \rightarrow T(X)

By the Data Processing Inequality, we know that processing the raw data into a summary must destroy (or at best, preserve) information about the hidden parameter:

I(θ;T(X))I(θ;X)I(\theta; T(X)) \le I(\theta; X)

We need to know: is it mathematically possible to compress the raw data XX into a smaller footprint T(X)T(X) such that we hit exact equality, losing absolutely zero information about the hidden parameter θ\theta?

DefinitionSufficient statistic

A statistic T(X)T(X) is sufficient for θ\theta if it contains all the information in XX about θ\theta. Mathematically, this means it achieves perfect equality in the Data Processing Inequality:

I(θ;T(X))=I(θ;X)I(\theta; T(X)) = I(\theta; X)

By the equality condition of the Data Processing Inequality, for the information to remain perfectly intact, the Markov chain must be reversible. The sequence θT(X)X\theta \rightarrow T(X) \rightarrow X must also form a valid Markov chain. This gives us the formal definition: T(X)T(X) is a sufficient statistic if the raw data XX is conditionally independent of the parameter θ\theta, given T(X)T(X).

If I give you the summary T(X)T(X), you have captured the “essence” of the system so perfectly that the actual parameter θ\theta provides no additional help in guessing which specific raw sequence XX actually occurred.

Examples of sufficient statistics

  1. The coin toss (discrete binary channel). Suppose a noisy wire flips bits with an unknown probability θ\theta. We send a sequence of nn zeroes, and receive a raw sequence of bits X1,X2,,XnX_1, X_2, \ldots, X_n. Instead of storing the exact sequence of bits (e.g., 0, 1, 1, 0, 1), the engineer simply stores the total number of flipped bits: T(X)=Xi=3T(X) = \sum X_i = 3. Is TT a sufficient statistic? Yes. If you know exactly 3 bits flipped out of 5, every single permutation that contains three 1s (like 1,1,1,0,0 or 0,1,0,1,1) is equally likely. The specific arrangement is purely random combinatorial noise. Knowing the actual physical parameter θ\theta does not give you any extra predictive power to guess which specific arrangement happened. The sum TT captured 100% of the useful information about θ\theta.
  2. The normal distribution (continuous). If XX is a sequence of normally distributed voltage readings with an unknown mean θ\theta and a known variance of 1, the sample mean Xˉ=1nXi\bar{X} = \frac{1}{n}\sum X_i is a sufficient statistic. Once you know the sample mean, the exact distribution of the individual voltage spikes around that mean is independent of θ\theta.
  3. The uniform distribution (continuous). If noise is distributed uniformly between θ\theta and θ+1\theta + 1, the sample mean is actually not sufficient. The boundaries dictate the reality of the parameter. A sufficient statistic here is the pair of the absolute highest and lowest spikes: T(X)=(max{Xi},min{Xi})T(X) = (\max\{X_i\}, \min\{X_i\}).

Minimal sufficient statistics

While a sufficient statistic guarantees no information is lost, it does not guarantee maximum compression. For example, if we are tracking the normal distribution of voltages, storing the pair of values (mean of even-indexed readings, mean of odd-indexed readings) is technically a sufficient statistic, it preserves all information about θ\theta. But it is redundant. It is storing two numbers when one would do the job.

DefinitionMinimal sufficient statistic

A statistic T(X)T(X) is a minimal sufficient statistic if it maximally compresses the information about θ\theta. Mathematically, T(X)T(X) is minimal if it is a direct function of every other possible sufficient statistic U(X)U(X). In terms of the Data Processing Inequality, this forces the minimal statistic to sit at the absolute end of the Markov chain bottleneck, stripping out every last bit of irrelevant combinatorial noise while preserving θ\theta:

θT(X)U(X)X\theta \rightarrow T(X) \rightarrow U(X) \rightarrow X

In the previous examples (the sum of coin tosses, the sample mean, and the min/max pair), those sufficient statistics are all minimal.

Fano’s inequality: the physical limits of guessing

The problem: can a smart algorithm guess its way out of noise? Until now, we have measured the noise in the wire as a continuous amount of leftover uncertainty, H(XY)H(X|Y). But in the real world, the receiver’s telegraph machine doesn’t just output a “cloud of uncertainty.” It must eventually print a specific, physical letter. It runs the noisy signal YY through a decoding algorithm (a deterministic function) to produce a final, best guess: X^=g(Y)\hat{X} = g(Y).

Because the guess X^\hat{X} is calculated directly from the received signal YY, the system forms a Markov chain:

XYX^X \rightarrow Y \rightarrow \hat{X}

We define a communication error as the event where the machine’s final guess does not match the sender’s original keystroke. The probability of error is Pe=Pr(X^X)P_e = \Pr(\hat{X} \neq X).

We know intuitively that if the wire’s noise H(XY)H(X|Y) is exactly zero, we can estimate XX from YY with a zero probability of error. Extending this, we expect to be able to achieve a low probability of error only if the conditional entropy H(XY)H(X|Y) is small.

Fano’s Inequality quantifies this exact relationship. It bridges the gap between information theory (entropy) and statistics (probability of error), establishing an absolute, unbreakable mathematical floor on how low the error rate can go based on the physical noise in the wire.

The theorem

TheoremFano’s inequality

For any estimator X^\hat{X} such that XYX^X \rightarrow Y \rightarrow \hat{X} forms a Markov chain, with error probability Pe=Pr(X^X)P_e = \Pr(\hat{X} \neq X), we have:

H(Pe)+PelogXH(XX^)H(XY)H(P_e) + P_e \log |\mathcal{X}| \ge H(X|\hat{X}) \ge H(X|Y)

(Here, H(Pe)H(P_e) is the binary entropy of the error event itself: PelogPe(1Pe)log(1Pe)-P_e \log P_e - (1-P_e) \log (1-P_e), and X|\mathcal{X}| is the size of the sender’s alphabet).

IntuitionPhysical intuition

Read the inequality from right to left. The right side, H(XY)H(X|Y), is the sheer physical noise introduced by the wire. The inequality states that this physical noise must be “absorbed” by the error rate on the left side. If Pe=0P_e = 0, the left side collapses to 0, which correctly implies that H(XY)H(X|Y) must also be 0.

The proof

ProofVia an error indicator variable

We prove this by tracking the machine’s success using an artificial binary random variable, EE. Let E=1E = 1 if the machine failed (X^X\hat{X} \neq X), and E=0E = 0 if the machine succeeded (X^=X\hat{X} = X).

We want to calculate the joint conditional entropy H(E,XX^)H(E, X | \hat{X}). We expand this using the Chain Rule for Entropy in two different directions.

Expansion 1 (condition on X first). H(E,XX^)=H(XX^)+H(EX,X^)H(E, X | \hat{X}) = H(X | \hat{X}) + H(E | X, \hat{X}) If we know both the sender’s true symbol XX and the machine’s guess X^\hat{X}, there is absolutely zero uncertainty about whether an error occurred. Therefore, H(EX,X^)=0H(E | X, \hat{X}) = 0. H(E,XX^)=H(XX^)H(E, X | \hat{X}) = H(X | \hat{X})

Expansion 2 (condition on E first). H(E,XX^)=H(EX^)+H(XE,X^)H(E, X | \hat{X}) = H(E | \hat{X}) + H(X | E, \hat{X}) By Theorem 2.6.5 (Conditioning Reduces Entropy), we can drop X^\hat{X} from the first term to upper-bound it: H(EX^)H(E)H(E | \hat{X}) \le H(E). Since EE is a simple binary coin flip (error or no error), its entropy is just H(Pe)H(P_e). Now we bound the second term, H(XE,X^)H(X | E, \hat{X}), by averaging over the two possible states of the error indicator EE: H(XE,X^)=Pr(E=0)H(XX^,E=0)+Pr(E=1)H(XX^,E=1)H(X | E, \hat{X}) = \Pr(E = 0)H(X | \hat{X}, E=0) + \Pr(E = 1)H(X | \hat{X}, E=1)

  • Case E=0 (no error). If there is no error, XX is exactly X^\hat{X}. Since we are conditioning on X^\hat{X}, there is 0 uncertainty about XX.
  • Case E=1 (error). If an error occurred, we know XX must be one of the remaining symbols in the alphabet. The absolute maximum uncertainty of any variable is the log of its total possible outcomes. We can safely upper bound this by logX\log |\mathcal{X}|.

Substituting these bounds back into the expansion gives: H(XE,X^)(1Pe)(0)+PelogXH(X | E, \hat{X}) \le (1 - P_e)(0) + P_e \log |\mathcal{X}| Therefore, Expansion 2 is bounded by: H(E,XX^)H(Pe)+PelogXH(E, X | \hat{X}) \le H(P_e) + P_e \log |\mathcal{X}|

Equating and applying data processing. Equating Expansion 1 and our bound for Expansion 2, we get: H(Pe)+PelogXH(XX^)H(P_e) + P_e \log |\mathcal{X}| \ge H(X|\hat{X}) Finally, because XYX^X \rightarrow Y \rightarrow \hat{X} is a Markov chain, the Data Processing Inequality guarantees that I(X;X^)I(X;Y)I(X; \hat{X}) \le I(X; Y). Expanding this into entropies (H(X)H(XX^)H(X)H(XY)H(X) - H(X|\hat{X}) \le H(X) - H(X|Y)) proves that H(XX^)H(XY)H(X|\hat{X}) \ge H(X|Y). We chain this inequality onto the end of our proof to finish it: H(Pe)+PelogXH(XX^)H(XY)H(P_e) + P_e \log |\mathcal{X}| \ge H(X|\hat{X}) \ge H(X|Y)


Useful corollaries and bounds

1. Weakening the bound for simplicity. Because the binary entropy of the error event, H(Pe)H(P_e), can never exceed 1 bit (which occurs when the machine is right 50% of the time), we can replace H(Pe)H(P_e) with 1.

1+PelogXH(XY)1 + P_e \log |\mathcal{X}| \ge H(X|Y)

Algebraically rearranging this gives us a direct formula for the absolute lowest possible error rate:

PeH(XY)1logXP_e \ge \frac{H(X|Y) - 1}{\log |\mathcal{X}|}

2. A tighter bound (if the guess must be in the alphabet). If the decoding algorithm is restricted to only guessing symbols that actually exist in the sender’s alphabet (X^X\hat{X} \in \mathcal{X}), then in the event of an error (E=1E=1), we know the true symbol XX must be one of the other (X1)(|\mathcal{X}| - 1) symbols. We can replace logX\log |\mathcal{X}| with log(X1)\log(|\mathcal{X}| - 1):

H(Pe)+Pelog(X1)H(XY)H(P_e) + P_e \log(|\mathcal{X}| - 1) \ge H(X|Y)

3. Direct observation bound. If we don’t use a decoding algorithm at all, and just let our “guess” be the raw received signal (X^=Y\hat{X} = Y), then PeP_e is simply the probability that the physical transmission was flipped: p=Pr(XY)p = \Pr(X \neq Y).

H(p)+plogXH(XY)H(p) + p \log |\mathcal{X}| \ge H(X|Y)

RemarkBlind guessing (the sharpness of Fano)

What if the wire is completely cut, and there is no YY at all? You just have to guess XX blindly. Your best mathematical guess is just picking the single most common symbol in the language (probability p1p_1). Your probability of error is Pe=1p1P_e = 1 - p_1. Fano’s inequality becomes: H(Pe)+Pelog(m1)H(X)H(P_e) + P_e \log(m - 1) \ge H(X). If the rest of the symbols p2pmp_2 \dots p_m are all perfectly uniform, this inequality evaluates to exactly equal. This proves that Fano’s Inequality is “sharp” (the bound cannot be mathematically tightened any further, because a physical scenario exists where the bound is perfectly reached).


A sidelight: probability of collision between independent sources

While discussing the probability of guessing and matching strings, we can introduce a related, beautiful inequality regarding two independent sources.

Motivation. Suppose two isolated telegraph operators, XX and XX', are completely independent of each other but both type following the same language statistics (they are i.i.d. with entropy H(X)H(X)). What is the probability that they both happen to press the exact same key at the exact same time?

Pr(X=X)=xPr(X=x)Pr(X=x)=xp2(x)\Pr(X = X') = \sum_x \Pr(X = x)\Pr(X' = x) = \sum_x p^2(x)

LemmaCollision probability

If XX and XX' are i.i.d. with entropy H(X)H(X), the probability they match is bounded by their entropy:

Pr(X=X)2H(X)\Pr(X = X') \ge 2^{-H(X)}

(Equality holds if and only if the operators type completely at random — a uniform distribution).

ProofVia Jensen’s inequality

We want to bound p2(x)\sum p^2(x). Let’s rewrite this cleverly to look like an Expected Value. The expected value of the probability itself is: E[p(X)]=p(x)p(x)=p2(x)\mathbb{E}[p(X)] = \sum p(x) p(x) = \sum p^2(x). Now, notice that p(x)p(x) can be rewritten using logs as 2logp(x)2^{\log p(x)}. Therefore, our collision probability is E[2logp(X)]\mathbb{E}[2^{\log p(X)}].

The function f(y)=2yf(y) = 2^y is strictly convex (its second derivative is positive). Therefore, by Jensen’s Inequality, the Expected Value of the function is \ge the function of the Expected Value:

E[2logp(X)]2E[logp(X)]\mathbb{E}[2^{\log p(X)}] \ge 2^{\mathbb{E}[\log p(X)]}

Look at the exponent on the right side: E[logp(X)]=p(x)logp(x)=H(X)\mathbb{E}[\log p(X)] = \sum p(x) \log p(x) = -H(X). Substituting this back in proves the lemma:

p2(x)2H(X)\sum p^2(x) \ge 2^{-H(X)}

CorollaryCollision of mismatched operators

What if operator XX follows English p(x)p(x), but operator XX' follows French r(x)r(x)? They are independent, but not identically distributed. The probability they collide is p(x)r(x)\sum p(x)r(x). We apply the exact same Jensen’s Inequality proof, but taking the Expected Value over pp:

Ep[2logr(X)]2Ep[logr(X)]\mathbb{E}_p[2^{\log r(X)}] \ge 2^{\mathbb{E}_p[\log r(X)]}

The right side exponent is p(x)logr(x)\sum p(x) \log r(x). We can add and subtract p(x)logp(x)p(x) \log p(x) to force it into recognizable terms:

p(x)logr(x)=p(x)logp(x)+p(x)logr(x)p(x)=H(p)D(pr)\sum p(x) \log r(x) = \sum p(x) \log p(x) + \sum p(x) \log \frac{r(x)}{p(x)} = -H(p) - D(p||r)

Therefore, the collision probability is bounded by both their baseline uncertainty and their Relative Entropy distance:

Pr(X=X)2H(p)D(pr)\Pr(X = X') \ge 2^{-H(p) - D(p||r)}

(By symmetry, you can evaluate the expectation over rr to get the alternate bound 2H(r)D(rp)\ge 2^{-H(r) - D(r||p)}).

← All notes