Documentation

Linglib.Core.Probability.Decision.Blackwell

Blackwell comparison of experiments #

A statistical experiment is a Markov kernel P : Kernel Θ 𝓧 generating data in 𝓧 from a parameter in Θ. Experiment P is at least as informative as P' : Kernel Θ 𝓧' when P' can be recovered from P by Markov post-processing ("garbling"): P' = Ξ· βˆ˜β‚– P for some Markov kernel Ξ·. This file develops that order and its characterization through Bayes risk.

[Bla53]'s comparison theorem states that P is at least as informative as P' if and only if, for every decision problem, the Bayes risk under P is no greater than under P'. We state and prove this equivalence over ProbabilityTheory.bayesRisk for finite spaces: the forward direction is the data-processing inequality, the converse the Blackwell–Sherman–Stein separation argument (see the implementation notes).

Main definitions #

Main statements #

Implementation notes #

The development is stated entirely over Mathlib's Kernel and bayesRisk with no further dependencies, so it can serve as a Mathlib.Probability.Decision.Blackwell candidate. On the utility scale, Core.Probability.Decision.ValueOfInformation reads the data-processing direction as the statement that garbling an experiment never raises its value of information, and reads the converse for deterministic experiments as a comparison of the value of observing classifiers.

Kernel.BlackwellDominates quantifies over all decision problems (every measurable action space 𝓨 and loss β„“ : Θ β†’ 𝓨 β†’ ℝβ‰₯0∞) and priors: dominance for a single one does not force garbling. The action-space universe is pinned to u (the universe of 𝓧') because the converse proof instantiates the dominance hypothesis at the action space 𝓨 := 𝓧' (the identity estimator).

The converse proof encodes finite kernels as real vectors of singleton masses (encode) and realizes the garblings of P as a compact convex polytope garblingSet P (the linear image of the product of standard simplices). If encode P' lies outside the polytope, the geometric Hahn–Banach theorem (geometric_hahn_banach_point_closed) yields a separating functional f; its coordinate matrix a ΞΈ x' = f (Pi.single ΞΈ (Pi.single x' 1)), shifted by a constant C to be nonnegative, defines a loss β„“ ΞΈ x' = ENNReal.ofReal (a ΞΈ x' + C) on actions 𝓧'. Under the uniform prior, the Bayes risk of any experiment Q evaluated at the identity estimator equals ENNReal.ofReal (|Θ|⁻¹ Β· f (encode Q) + C): the identity estimator pins P' to f (encode P'), while every estimator for P produces a garbling, whose f-value exceeds the separation level. This realizes bayesRisk β„“ P' Ο€ < bayesRisk β„“ P Ο€, contradicting the hypothesis β€” no minimax theorem is needed, the infimum is bounded below directly. All analytic inputs come from Mathlib (StdSimplex.compactSpace, the geometric_hahn_banach_* lemmas, bayesRisk_fintype). The kernel-to-stochastic-matrix bridge (encode, garblingMap, buildKernel) is currently private proof scaffolding; it is a self-contained finite-kernel ↔ row-stochastic-matrix correspondence that would naturally graduate to its own public file when upstreamed.

References #

def ProbabilityTheory.Kernel.IsGarblingOf {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] (P' : Kernel Θ 𝓧') (P : Kernel Θ 𝓧) :

P' is a garbling of P when a Markov kernel Ξ· post-processes P into P', so that P' = Ξ· βˆ˜β‚– P and P is at least as informative as P'.

Equations
  • P'.IsGarblingOf P = βˆƒ (Ξ· : ProbabilityTheory.Kernel 𝓧 𝓧'), ProbabilityTheory.IsMarkovKernel Ξ· ∧ P' = Ξ·.comp P
Instances For
    theorem ProbabilityTheory.Kernel.IsGarblingOf.refl {Θ : Type u_1} {𝓧 : Type u_2} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] (P : Kernel Θ 𝓧) [IsMarkovKernel P] :
    theorem ProbabilityTheory.Kernel.IsGarblingOf.trans {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] {𝓧'' : Type u_3} [MeasurableSpace 𝓧''] {P'' : Kernel Θ 𝓧''} {P' : Kernel Θ 𝓧'} {P : Kernel Θ 𝓧} (hβ‚‚ : P''.IsGarblingOf P') (h₁ : P'.IsGarblingOf P) :
    def ProbabilityTheory.Kernel.BlackwellDominates {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] (P : Kernel Θ 𝓧) (P' : Kernel Θ 𝓧') :

    P Blackwell-dominates P' when, for every decision problem (an action space 𝓨 with a loss β„“) and every prior Ο€, the Bayes risk under P is at most that under P'.

    Equations
    • One or more equations did not get rendered due to their size.
    Instances For
      theorem ProbabilityTheory.bayesRisk_le_of_isGarblingOf {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] {𝓨 : Type u} [MeasurableSpace 𝓨] (β„“ : Θ β†’ 𝓨 β†’ ENNReal) {P : Kernel Θ 𝓧} {P' : Kernel Θ 𝓧'} (h : P'.IsGarblingOf P) (Ο€ : MeasureTheory.Measure Θ) :
      bayesRisk β„“ P Ο€ ≀ bayesRisk β„“ P' Ο€

      If P' is a garbling of P, the Bayes risk under P is at most that under P' in every decision problem.

      theorem ProbabilityTheory.blackwellDominates_of_isGarblingOf {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] {P : Kernel Θ 𝓧} {P' : Kernel Θ 𝓧'} (h : P'.IsGarblingOf P) :

      A garbling of P is Blackwell-dominated by P.

      The garbling polytope (finite case) #

      Over finite spaces, the Markov garblings {Ξ· βˆ˜β‚– P | Ξ· Markov} of P, encoded by their singleton masses as vectors in Θ β†’ 𝓧' β†’ ℝ, form a compact convex polytope garblingSet P. It is the linear image of the product of standard simplices β€” the stochastic matrices Ξ· β€” under garblingMap P. If encode P' lies outside the polytope, a separating functional gives a decision problem on which P' is strictly worse than P, which proves the converse.

      theorem ProbabilityTheory.isGarblingOf_of_bayesRisk_uniform_le {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] [Fintype Θ] [Fintype 𝓧] [Fintype 𝓧'] [Nonempty Θ] [MeasurableSingletonClass Θ] [MeasurableSingletonClass 𝓧] [MeasurableSingletonClass 𝓧'] {P : Kernel Θ 𝓧} {P' : Kernel Θ 𝓧'} [IsMarkovKernel P] [IsMarkovKernel P'] (h : βˆ€ (β„“ : Θ β†’ 𝓧' β†’ ENNReal), (βˆ€ (ΞΈ : Θ) (x' : 𝓧'), β„“ ΞΈ x' β‰  ⊀) β†’ bayesRisk β„“ P ((↑(Fintype.card Θ))⁻¹ β€’ MeasureTheory.Measure.count) ≀ bayesRisk β„“ P' ((↑(Fintype.card Θ))⁻¹ β€’ MeasureTheory.Measure.count)) :

      If P has Bayes risk at most that of P' under the uniform prior for every everywhere-finite loss, then P' is a garbling of P.

      theorem ProbabilityTheory.isGarblingOf_of_blackwellDominates {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] [Fintype Θ] [Fintype 𝓧] [Fintype 𝓧'] [Nonempty Θ] [MeasurableSingletonClass Θ] [MeasurableSingletonClass 𝓧] [MeasurableSingletonClass 𝓧'] {P : Kernel Θ 𝓧} {P' : Kernel Θ 𝓧'} [IsMarkovKernel P] [IsMarkovKernel P'] (h : P.BlackwellDominates P') :

      On finite spaces with Θ nonempty, if P Blackwell-dominates P' and both are Markov kernels, then P' is a garbling of P.

      Each hypothesis is needed. With Θ empty every Bayes risk is 0, yet no Markov kernel maps a nonempty 𝓧 to an empty 𝓧'. The zero kernel has Bayes risk 0 for every loss and so dominates every P', yet Ξ· βˆ˜β‚– 0 = 0. Over a one-point sample space P' = 2 β€’ P has twice the Bayes risk of P for every loss, yet is no Markov garbling of P. Dominance in a single decision problem does not force garbling either.

      theorem ProbabilityTheory.isGarblingOf_iff_blackwellDominates {Θ : Type u_1} {𝓧 : Type u_2} {𝓧' : Type u} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] [m𝓧' : MeasurableSpace 𝓧'] [Fintype Θ] [Fintype 𝓧] [Fintype 𝓧'] [Nonempty Θ] [MeasurableSingletonClass Θ] [MeasurableSingletonClass 𝓧] [MeasurableSingletonClass 𝓧'] {P : Kernel Θ 𝓧} {P' : Kernel Θ 𝓧'} [IsMarkovKernel P] [IsMarkovKernel P'] :

      On finite spaces with Θ nonempty and Markov P and P', P' is a garbling of P exactly when P Blackwell-dominates P'.

      Deterministic experiments: partitions as kernels #

      A deterministic classifier f : Θ β†’ 𝓧 is the experiment Kernel.deterministic f hf, an error-free observation of the cell of ΞΈ in the partition of Θ into the fibers of f. Between deterministic experiments the garbling order is functional factoring.

      theorem ProbabilityTheory.Kernel.deterministic_isGarblingOf_deterministic_iff {Θ : Type u_1} {𝓧 : Type u_2} [mΘ : MeasurableSpace Θ] [m𝓧 : MeasurableSpace 𝓧] {𝓨 : Type u_3} [MeasurableSpace 𝓨] [Countable 𝓧] [MeasurableSingletonClass 𝓧] [MeasurableSingletonClass 𝓨] [Nonempty 𝓨] {f : Θ β†’ 𝓧} {g : Θ β†’ 𝓨} (hf : Measurable f) (hg : Measurable g) :
      (deterministic g hg).IsGarblingOf (deterministic f hf) ↔ βˆƒ (ψ : 𝓧 β†’ 𝓨), g = ψ ∘ f

      deterministic g is a garbling of deterministic f exactly when g factors through f; randomized post-processing gains nothing between deterministic experiments.