Tutorials › Statistical Learning (Graduate) › Data Processing and Product-Sample Divergence

Minimax Lower-Bound Methods · Tutorial 77 of 989

Data Processing and Product-Sample Divergence

Learn why common data transformations reduce distinguishability and why independent observations make Kullback–Leibler divergence add across samples.

Intermediate 11 min read

What You'll Learn

  • Define Kullback–Leibler divergence for probability measures, including when it is infinite
  • Prove data processing for a common Markov kernel on arbitrary measurable spaces
  • Apply data processing to deterministic statistics and randomized channels
  • Derive exact divergence additivity for independent product distributions
  • Handle infinite component divergences using coordinate projection
  • Check how coarsening and repeated sampling change distinguishability

Why Transformations and Repeated Samples Matter

In statistical learning, distributions describe what an experiment can reveal about an unknown state. A procedure may discard part of the data, summarize it with a statistic, or randomize it before producing an output. Such processing cannot make two distributions more distinguishable when the same procedure is applied to both. In contrast, independent observations accumulate evidence: their Kullback–Leibler divergence adds.

These facts are useful in lower-bound arguments. Data processing lets us compare the information in an observation with the information retained after a transformation, while product additivity quantifies how distinguishability changes with sample size. The results below apply to probability measures on measurable spaces; no densities with respect to Lebesgue measure are required.

Kullback–Leibler Divergence and Markov Kernels

Definition: Let \(P\) and \(Q\) be probability measures on the same measurable space \((\mathcal X,\mathcal A)\). If \(P\) is absolutely continuous with respect to \(Q\), write \(r=dP/dQ\) and define $$ D(P\|Q)=\int r\log r\,dQ=\int \log\!\left(\frac{dP}{dQ}\right)dP, $$ where \(0\log 0=0\). If \(P\) is not absolutely continuous with respect to \(Q\), set \(D(P\|Q)=+\infty\). The divergence takes values in \([0,+\infty]\).

The negative part of \(r\log r\) is integrable because \(t\log t\geq -1/e\) for \(t\geq0\). Thus the displayed integral is well-defined, possibly with value \(+\infty\). Nonnegativity follows from Jensen’s inequality: \(t\mapsto t\log t\) is convex, and \(\int r\,dQ=1\), so \(\int r\log r\,dQ\geq 1\log 1=0\).

Definition: A Markov kernel \(K\) from \((\mathcal X,\mathcal A)\) to \((\mathcal Y,\mathcal B)\) assigns a probability measure \(K(x,\cdot)\) on \(\mathcal Y\) to each \(x\in\mathcal X\), with \(x\mapsto K(x,B)\) measurable for every \(B\in\mathcal B\). It maps a probability measure \(P\) to the output measure \(PK\), defined by $$ (PK)(B)=\int K(x,B)\,P(dx). $$

A measurable function \(T:\mathcal X\to\mathcal Y\) is a special case: use the deterministic kernel \(K(x,B)=\mathbf 1\{T(x)\in B\}\). A randomized procedure is represented by a general Markov kernel. In both cases, the same kernel must be applied to \(P\) and \(Q\).

Data Processing: The Divergence Cannot Increase

Theorem (Data-processing inequality): For probability measures \(P,Q\) on a measurable space and a Markov kernel \(K\) from that space to another measurable space, $$ D(PK\|QK)\leq D(P\|Q). $$

This statement holds on arbitrary measurable spaces. The proof uses a joint probability measure and conditional expectation, which exist without a standard Borel assumption.

Proof. If \(D(P\|Q)=+\infty\), the inequality holds because every divergence is at most \(+\infty\). Otherwise \(P\ll Q\); set \(r=dP/dQ\). Construct the joint law \(R\) of \((X,Y)\) by first drawing \(X\sim Q\), then drawing \(Y\) from \(K(X,\cdot)\). Its \(Y\)-marginal is \(QK\). Under this joint law, \(E_R[r(X)]=1\).

Because \(P\ll Q\), also \(PK\ll QK\): if \((QK)(B)=0\), then \(K(x,B)=0\) for \(Q\)-almost every \(x\), and hence \((PK)(B)=\int K(x,B)P(dx)=0\). Let \(g=d(PK)/d(QK)\). For each measurable \(B\subseteq\mathcal Y\), $$ E_R\!\left[r(X)\mathbf 1\{Y\in B\}\right] =\int K(x,B)P(dx) =(PK)(B) =\int_B g(y)(QK)(dy). $$ Thus \(g(Y)\) is a version of \(E_R[r(X)\mid \sigma(Y)]\): the equality holds on every event \(\{Y\in B\}\), which generates \(\sigma(Y)\). Conditional Jensen’s inequality for the convex function \(\phi(t)=t\log t\) now gives $$ D(PK\|QK) =E_R[\phi(g(Y))] \leq E_R[\phi(r(X))] =D(P\|Q). $$ The expectations are well-defined because \(\phi\) is bounded below; if the right side is finite, conditional Jensen applies directly. This proves the inequality. \(\square\)

Worked Example: Coarsening a Three-Category Observation

Let \(P=(0.5,0.3,0.2)\) and \(Q=(0.2,0.3,0.5)\) on categories \(a,b,c\). A statistic records only whether the outcome is in \(\{a,b\}\) or is \(c\). Calculate the divergence before and after this coarsening.

Before coarsening, the likelihood ratios are \(0.5/0.2=2.5\), \(0.3/0.3=1\), and \(0.2/0.5=0.4\). Therefore $$ D(P\|Q)=0.5\log(2.5)+0.3\log(1)+0.2\log(0.4) \approx 0.45815+0-0.18326=0.27489. $$

After coarsening, the two distributions are \((0.8,0.2)\) and \((0.5,0.5)\), so $$ D(PK\|QK)=0.8\log(1.6)+0.2\log(0.4) \approx 0.37600-0.18326=0.19274. $$ The divergence falls, as predicted: merging categories removes information that could distinguish \(P\) from \(Q\).

Independent Product Samples

For independent observations, the joint law is a product measure. The next result states exactly how its divergence depends on the component distributions.

Theorem (Product additivity of Kullback–Leibler divergence): Let \(P_i,Q_i\) be probability measures on measurable spaces \((\mathcal X_i,\mathcal A_i)\), for \(i=1,\ldots,n\). Then $$ D\!\left(\bigotimes_{i=1}^n P_i\,\middle\|\,\bigotimes_{i=1}^n Q_i\right) =\sum_{i=1}^n D(P_i\|Q_i), $$ where the sum and equality are interpreted in the extended nonnegative reals. In particular, for identically distributed independent observations, \(D(P^{\otimes n}\|Q^{\otimes n})=nD(P\|Q)\).

Proof when every component divergence is finite. Finiteness implies \(P_i\ll Q_i\); write \(r_i=dP_i/dQ_i\). The density of the product measure \(\bigotimes_i P_i\) with respect to \(\bigotimes_i Q_i\) is \(\prod_i r_i(x_i)\). Indeed, integrating this product density over a measurable rectangle gives the product of the component probabilities, which determines the product measure. Therefore the log density is \(\sum_i\log r_i(x_i)\) wherever the product density is positive.

Each \(\log r_i(X_i)\), under \(X_i\sim P_i\), is integrable when \(D(P_i\|Q_i)<\infty\). To check the negative part, write \(r_i=dP_i/dQ_i\): $$ E_{P_i}\!\left[(-\log r_i(X_i))_+\right] =\int_{\{r_i<1\}} -r_i\log r_i\,dQ_i \leq \frac{1}{e}. $$ The positive part is then finite because the divergence is finite. Thus expectation of the finite sum can be taken term by term under the product law, giving $$ D\!\left(\bigotimes_{i=1}^nP_i\,\middle\|\,\bigotimes_{i=1}^nQ_i\right) =\sum_{i=1}^n E_{P_i}[\log r_i(X_i)] =\sum_{i=1}^nD(P_i\|Q_i). $$

Proof when at least one component divergence is infinite. Choose an index \(j\) with \(D(P_j\|Q_j)=+\infty\). Projection of the product observation onto coordinate \(j\) is a deterministic Markov kernel. The data-processing inequality just proved gives $$ D\!\left(\bigotimes_{i=1}^nP_i\,\middle\|\,\bigotimes_{i=1}^nQ_i\right) \geq D(P_j\|Q_j)=+\infty. $$ Hence the product divergence is infinite, matching the extended sum. This argument is valid on arbitrary measurable spaces because the data-processing result above was proved at that level of generality. \(\square\)

Worked Example: Five Independent Bernoulli Observations

Let one observation have a Bernoulli distribution with success probability \(p=0.6\) under \(P\), and probability \(q=0.4\) under \(Q\). The one-observation divergence is $$ D(P\|Q)=0.6\log\!\left(\frac{0.6}{0.4}\right) +0.4\log\!\left(\frac{0.4}{0.6}\right) =0.2\log(1.5)\approx 0.08109. $$

For five independent observations, the product-additivity theorem gives $$ D(P^{\otimes 5}\|Q^{\otimes 5}) =5(0.08109)\approx 0.40547. $$ The multiplier is the number of independent observations, not the number of possible sample outcomes. Each observation contributes the same expected log-likelihood ratio.

Randomized Processing and a Numerical Check

Worked Example: A Noisy Binary Channel

Suppose the input is binary. Under \(P\), its probability of being 1 is \(0.8\); under \(Q\), it is \(0.3\). A channel flips the input with probability \(0.1\), independently of its value. First calculate the input divergence: $$ D(P\|Q)=0.8\log(0.8/0.3)+0.2\log(0.2/0.7) \approx 0.8(0.98083)+0.2(-1.25276)=0.53411. $$

The output is 1 either when the input is 1 and is not flipped, or when the input is 0 and is flipped. Its probability under \(P\) is \(0.8(0.9)+0.2(0.1)=0.74\); under \(Q\) it is \(0.3(0.9)+0.7(0.1)=0.34\). Thus the output divergence is $$ 0.74\log(0.74/0.34)+0.26\log(0.26/0.66) \approx 0.33330. $$ The calculation verifies that this randomized processing reduces divergence from approximately \(0.53411\) to \(0.33330\).

Common Pitfalls and Why the Results Matter

  • Assuming any processing increases divergence. The inequality goes in the opposite direction when the same kernel is applied to both distributions: \(D(PK\|QK)\leq D(P\|Q)\). The output cannot recover distinctions that the channel has erased.
  • Applying different procedures to the two distributions. Data processing compares \(PK\) with \(QK\) for a common kernel \(K\). If the procedures differ, the theorem does not provide this comparison.
  • Adding divergences without independence. Product additivity concerns product distributions, which encode independence. A dependent joint distribution is not in general the product of its marginals, so the theorem cannot be applied by simply summing marginal divergences.
  • Ignoring infinite component divergences. The finite-density calculation alone does not cover that case. Coordinate projection and data processing show that an infinite divergence in any component forces the product divergence to be infinite.

Together, these properties separate two effects. A common transformation can only discard distinguishability, while independent repetition adds divergence exactly. They allow an analyst to compare the information available in a full experiment with that in a reduced observation, and to track how evidence accumulates across independent samples.

Check Your Understanding

Use the definitions and arguments above to answer these questions.

  1. Why does \(P\ll Q\) imply \(PK\ll QK\) for a Markov kernel \(K\)?
  2. A statistic merges several possible outcomes. Which inequality relates the divergence before and after recording the statistic?
  3. If \(D(P\|Q)=0.12\) is finite, what is \(D(P^{\otimes 8}\|Q^{\otimes 8})\)?
  4. Why does an infinite divergence in one coordinate imply an infinite divergence for the whole product experiment?
  5. Which assumption allows component divergences to be added for a sample of repeated observations?
Clinical Biostats