InformationTheory.Shannon.Chernoff.Basic
Chernoff information and the Hoeffding tradeoff exponent #
The Chernoff exponent C(P₁, P₂) and the Hoeffding tradeoff exponent E₂(α)
(Cover–Thomas) for finite-alphabet pmfs P₁, P₂ : α → ℝ, together with the
convexity of λ ↦ log Z(λ) and the achievability side of the Chernoff bound on the Bayes
error of an n-sample binary hypothesis test.
Main definitions #
chernoffZSum P₁ P₂ lam— the Chernoff partition functionZ(λ) := ∑ a, P₁(a)^(1-λ) · P₂(a)^λ.chernoffInfo P₁ P₂— the Chernoff information-sInf ((log ∘ Z) '' Icc 0 1).chernoffMediator P₁ P₂ lam— the tilted mediator pmfP₁(a)^(1-λ) · P₂(a)^λ / Z(λ).hoeffdingE2 P₁ P₂ alpha— the Hoeffding tradeoff exponent at Type I levelalpha.hoeffdingConstraintSet P₁ alpha— the constraint set{Q ∈ stdSimplex | klDivPmf Q P₁ ≤ alpha}.bayesErrorMinPmf P₁ P₂ n— the optimaln-sample Bayes error with equal priors.
Main statements #
chernoffInfo_attained,chernoffInfo_nonneg— the Chernoff information is attained and nonnegative.convexOn_chernoffLogZ—λ ↦ log Z(λ)is convex onIcc 0 1.hoeffdingE2_attained,hoeffdingE2_nonneg,hoeffdingE2_unique— existence, nonnegativity, and uniqueness of the Hoeffding minimizer.bayesErrorMinPmf_le_half_Z_pow— the Chernoff boundbayesErrorMinPmf ≤ (1/2) · Z(λ)^n.chernoff_lemma_achievability— the achievability rateliminf ≥ chernoffInfo.
Implementation notes #
Everything is stated in finite-alphabet pmf form (α → ℝ) so that CsiszarProjection.klDivPmf
can be reused directly. chernoffInfo and hoeffdingE2 are stated as sInf over an image of a
compact set, so attainment follows from IsCompact.exists_sInf_image_eq. The convexity of
log Z is obtained from the multiplicative Hölder bound Z(αλ₁ + βλ₂) ≤ Z(λ₁)^α · Z(λ₂)^β.
References #
- T. M. Cover and J. A. Thomas, Elements of Information Theory (2nd ed.), Wiley, 2006.
Chernoff exponent: definition and basic properties #
InformationTheory.Shannon.Chernoff.chernoffZSum
sourceThe Chernoff partition function Z(λ) := ∑_a P₁(a)^(1-λ) · P₂(a)^λ.
Equations
Instances For
Used by
InformationTheory.Shannon.Chernoff.chernoffInfo
sourceThe Chernoff information C(P₁, P₂) := -min_{λ ∈ [0,1]} log Z(λ), stated as -sInf of the
image; chernoffInfo_attained shows the minimum is realized.
Equations
- InformationTheory.Shannon.Chernoff.chernoffInfo P₁ P₂ = -sInf ((fun (lam : ℝ) => Real.log (InformationTheory.Shannon.Chernoff.chernoffZSum P₁ P₂ lam)) '' Set.Icc 0 1)
Instances For
Used by
Endpoint values Z(0) = Z(1) = 1 #
InformationTheory.Shannon.Chernoff.chernoffZSum_lam_zero
sourceZ(0) = ∑ P₁(a) = 1 for pmf P₁.
Used by
InformationTheory.Shannon.Chernoff.chernoffZSum_lam_one
sourceZ(1) = ∑ P₂(a) = 1 for pmf P₂.
Used by
Strict positivity #
InformationTheory.Shannon.Chernoff.chernoffZSum_term_pos
sourceEach summand (P₁ a)^(1-λ) · (P₂ a)^λ is strictly positive under full support.
Used by
InformationTheory.Shannon.Chernoff.chernoffZSum_pos
sourceZ(λ) > 0 under full support (any λ ∈ ℝ).
Used by
Continuity #
InformationTheory.Shannon.Chernoff.chernoffZSum_continuous
sourceλ ↦ Z(λ) is continuous on ℝ.
Used by
InformationTheory.Shannon.Chernoff.chernoffLogZ_continuous
sourceλ ↦ log Z(λ) is continuous on ℝ (under full-support strict positivity of Z).
Used by
Attainment and nonnegativity of chernoffInfo #
InformationTheory.Shannon.Chernoff.chernoffInfo_attained
sourceChernoff information is attained: ∃ λ* ∈ Icc 0 1, chernoffInfo = -log Z(λ*).
Used by
InformationTheory.Shannon.Chernoff.chernoffInfo_nonneg
sourcechernoffInfo P₁ P₂ ≥ 0.
Used by
Symmetry Z_{P₁,P₂}(λ) = Z_{P₂,P₁}(1 - λ) #
InformationTheory.Shannon.Chernoff.chernoffZSum_swap
sourceSymmetry of chernoffZSum under λ ↔ 1 - λ: Z_{P₁,P₂}(λ) = Z_{P₂,P₁}(1 - λ).
Used by
Hoeffding tradeoff exponent #
InformationTheory.Shannon.Chernoff.klDivPmf_self_eq_zero
sourceklDivPmf P P = 0: the Kullback–Leibler divergence of a positive pmf with itself is zero.
Used by
InformationTheory.Shannon.Chernoff.hoeffdingE2
sourceHoeffding tradeoff exponent at Type I level alpha:
E₂(α) := min { klDivPmf Q P₂ | Q ∈ stdSimplex ∧ klDivPmf Q P₁ ≤ α }.
Equations
- One or more equations did not get rendered due to their size.
Instances For
Used by
InformationTheory.Shannon.Chernoff.hoeffdingConstraintSet
sourceThe Hoeffding constraint set
K := {Q ∈ stdSimplex | klDivPmf Q P₁ ≤ α}.
Equations
- InformationTheory.Shannon.Chernoff.hoeffdingConstraintSet P₁ alpha = {Q : α → ℝ | Q ∈ stdSimplex ℝ α ∧ InformationTheory.Shannon.CsiszarProjection.klDivPmf Q P₁ ≤ alpha}
Instances For
Used by
InformationTheory.Shannon.Chernoff.hoeffdingConstraintSet_nonempty
sourceThe Hoeffding constraint set is non-empty when α ≥ 0 and P₁ is a positive pmf:
P₁ itself satisfies klDivPmf P₁ P₁ = 0 ≤ α.
Used by
InformationTheory.Shannon.Chernoff.hoeffdingConstraintSet_subset_stdSimplex
sourceThe Hoeffding constraint set is a subset of the simplex.
Used by
InformationTheory.Shannon.Chernoff.hoeffdingConstraintSet_isClosed
sourceThe Hoeffding constraint set is closed (intersection of the closed simplex with
the closed sublevel set of the continuous function Q ↦ klDivPmf Q P₁).
Used by
InformationTheory.Shannon.Chernoff.hoeffdingE2_attained
sourceThe Hoeffding infimum is attained: there exists Q* ∈ K with
hoeffdingE2 P₁ P₂ alpha = klDivPmf Q* P₂.
Used by
InformationTheory.Shannon.Chernoff.hoeffdingE2_nonneg
sourcehoeffdingE2 P₁ P₂ alpha ≥ 0.
Used by
Convexity of log Z(λ) via Hölder #
The convexity of λ ↦ log Z(λ) on Icc 0 1 (Cover–Thomas). From Hölder's inequality one
obtains the multiplicative form Z(αλ₁ + βλ₂) ≤ Z(λ₁)^α · Z(λ₂)^β; taking logarithms gives the
convexity. The endpoints α = 0 and α = 1 are handled separately, since they fall outside the
1 < p requirement of Real.HolderConjugate.
InformationTheory.Shannon.Chernoff.chernoffZSum_holder_mul
sourceHölder multiplicative form for the Chernoff partition function:
Z(αλ₁ + βλ₂) ≤ Z(λ₁)^α · Z(λ₂)^β for α, β ∈ (0, 1), α + β = 1.
Stated under full-support P₁, P₂ > 0 to avoid 0^x corner cases.
Used by
InformationTheory.Shannon.Chernoff.convexOn_chernoffLogZ
sourceλ ↦ log Z(λ) is convex on Icc 0 1.
Used by
Chernoff mediator T_λ (pmf form) #
InformationTheory.Shannon.Chernoff.chernoffMediator
sourceThe Chernoff mediator pmf T_λ(a) := P₁(a)^(1-λ) · P₂(a)^λ / Z(λ).
Equations
- InformationTheory.Shannon.Chernoff.chernoffMediator P₁ P₂ lam a = P₁ a ^ (1 - lam) * P₂ a ^ lam / InformationTheory.Shannon.Chernoff.chernoffZSum P₁ P₂ lam
Instances For
Used by
InformationTheory.Shannon.Chernoff.chernoffMediator_pos
sourceT_λ(a) > 0 under full-support P₁, P₂ > 0.
Used by
InformationTheory.Shannon.Chernoff.chernoffMediator_sum_eq_one
sourceT_λ is a pmf: ∑ a, T_λ(a) = 1 (by construction).
Used by
InformationTheory.Shannon.Chernoff.chernoffMediator_lam_zero
sourceMediator pmf bound: at endpoints lam = 0 (resp lam = 1), T_λ = P₁ (resp P₂)
under full support + probability conditions.
Used by
InformationTheory.Shannon.Chernoff.chernoffMediator_lam_one
sourceUsed by
Uniqueness of the Hoeffding minimizer #
InformationTheory.Shannon.Chernoff.hoeffdingE2_unique
sourceThe Hoeffding minimizer is unique: klDivPmf · P₂ is strictly convex on the convex constraint
set K, so the minimizer of hoeffdingE2 is unique.
Used by
Chernoff achievability (Bayes error upper bound) #
The achievability side of Cover–Thomas: the Chernoff bound
bayesErrorMinPmf P₁ P₂ n ≤ (1/2) · Z(λ)^n for every λ ∈ Icc 0 1 yields
liminf_n -(1/n) log bayesErrorMinPmf ≥ chernoffInfo P₁ P₂. The n-fold IID structure is written
directly with a finite sum ∑_{x : Fin n → α} of finite products ∏ i, P (x i).
InformationTheory.Shannon.Chernoff.bayesErrorMinPmf
sourcen-IID Bayes error in pmf form:
bayesErrorMinPmf P₁ P₂ n := (1/2) · ∑_{x : Fin n → α} min(∏ P₁(x_i), ∏ P₂(x_i)).
This is the optimal Bayes error for the 2-class hypothesis test with equal priors
1/2 : 1/2 on n IID samples (Bayes-optimal rule decides i := argmax_i P_i^n(x),
giving error contribution (1/2) · min(P₁^n(x), P₂^n(x)) per x).
Equations
Instances For
Used by
Per-point bound min(a, b) ≤ a^(1-λ) · b^λ #
InformationTheory.Shannon.Chernoff.min_le_rpow_mul_rpow
sourceGeometric mean inequality (degenerate Hölder form):
min(a, b) ≤ a^{1-λ} · b^λ for a, b ≥ 0, λ ∈ [0, 1].
Used by
bayesErrorMinPmf ≤ (1/2) Z(λ)^n #
InformationTheory.Shannon.Chernoff.prod_rpow_mul_rpow
sourceAuxiliary: n-IID per-point factorization:
∏ i, (P₁ (x i)) ^ (1-lam) * (P₂ (x i)) ^ lam = (∏ i, P₁ (x i)) ^ (1-lam) * (∏ i, P₂ (x i)) ^ lam
under P₁, P₂ ≥ 0.
Used by
InformationTheory.Shannon.Chernoff.sum_prod_rpow_eq_Z_pow
sourcen-IID Chernoff partition function via product factorization:
∑_{x : Fin n → α} (∏ i, P₁ (x i))^(1-lam) · (∏ i, P₂ (x i))^lam = Z(λ)^n.
Used by
InformationTheory.Shannon.Chernoff.bayesErrorMinPmf_le_half_Z_pow
sourceThe Chernoff bound bayesErrorMinPmf ≤ (1/2) · Z(λ)^n for each λ ∈ Icc 0 1.
Used by
Positivity of bayesErrorMinPmf #
InformationTheory.Shannon.Chernoff.bayesErrorMinPmf_pos
sourcebayesErrorMinPmf > 0 under full support P₁, P₂ > 0.
Used by
Per-λ rate lower bound #
InformationTheory.Shannon.Chernoff.chernoff_rate_ge_neg_log_Z_per_lam
sourceFor each fixed λ ∈ Icc 0 1,
-(1/n) log bayesErrorMinPmf ≥ -log Z(λ) + (log 2)/n (eventually for n ≥ 1).
The + log 2 / n slack term vanishes as n → ∞, leaving -log Z(λ) (and after min
over λ, chernoffInfo).
Used by
InformationTheory.Shannon.Chernoff.chernoff_rate_ge_chernoffInfo_eventually
sourceFor each fixed λ* ∈ Icc 0 1 attaining chernoffInfo (= -log Z(λ*)),
-(1/n) log bayesErrorMinPmf ≥ chernoffInfo + (log 2)/n (eventually for n ≥ 1).
Used by
InformationTheory.Shannon.Chernoff.chernoff_achievability
sourceChernoff achievability (rate-side lower bound):
liminf_n -(1/n) log bayesErrorMinPmf ≥ chernoffInfo P₁ P₂.
Used by
Achievability main statement #
InformationTheory.Shannon.Chernoff.chernoff_lemma_achievability
sourceChernoff bound (achievability half): the exponential convergence rate of
bayesErrorMinPmf is at least chernoffInfo P₁ P₂.
See also chernoff_achievability.