Documentation

TauCeti.Probability.Exchangeability.FiniteDeFinetti

Quantitative finite de Finetti approximation #

Let x : κ → α be a nonempty finite population. Its empirical distribution is the pushforward of the uniform law on κ by x. Sampling ι entries from that distribution independently is therefore the same as choosing a uniform map ι → κ and reading the selected entries of x.

The general sampling laws, their empirical-product characterization, and their collision bounds live in TauCeti.Probability.PopulationSampling.Finite.

For a random population with law ρ, sampleWithReplacement ρ draws the population and, independently, a uniform map ι → κ, and reads the selected entries. It is the same samplePopulation construction as sampleWithoutReplacement ρ, which instead draws a uniform injective map; when κ is nonempty, by the previous paragraph it is the mixture over ρ of the finite product laws of the empirical distributions. The law sampleWithoutReplacement ρ already represents every shorter marginal of a finite exchangeable process. The collision coupling between uniform maps and uniform injective maps consequently gives, for every measurable event A,

prefixLaw μ X m A ≤ sampleWithReplacement (prefixLaw μ X n) A + choose(m, 2) / n
sampleWithReplacement (prefixLaw μ X n) A ≤ prefixLaw μ X m A + choose(m, 2) / n.

Thus, for n > 0, every m-coordinate marginal of an n-exchangeable process is quantitatively approximated by a mixture of m-fold product measures of empirical distributions; for n = m = 0 both sides are the shared sampling construction on the empty index type. The bound is stated eventwise, which fixes its normalization; in terms of Mathlib's total variation distance etvDist, the signed-measure norm ‖P - Q‖ (twice sup_A |P A - Q A|), it gives etvDist P Q ≤ 2 * choose(m, 2) / n.

Counting injective selections instead of using the union bound gives the sharper eventwise error 1 - n.descFactorial m / n^m, and twice that error for the signed-measure norm. This includes m = n = 0, whose error is zero. For a nonempty population the error is the exact index-collision probability; it need not equal the discrepancy of the sampled values, since distinct indices can hold the same value.

Main declarations #

References #

No material is adapted from cameronfreer/exchangeability; that development concerns infinite exchangeable sequences rather than quantitative finite approximation.

theorem TauCeti.Probability.ExchangeableAt.finiteDeFinetti {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) {A : Set (Fin m → α)} (hA : MeasurableSet A) :
(prefixLaw μ X m) A ≤ (sampleWithReplacement (prefixLaw μ X n)) A + ↑(m.choose 2) / ↑n ∧ (sampleWithReplacement (prefixLaw μ X n)) A ≤ (prefixLaw μ X m) A + ↑(m.choose 2) / ↑n

Finite de Finetti theorem. If the first n coordinates of a process are exchangeable and m ≤ n, then its m-prefix law and the with-replacement sample of its n-prefix law differ by at most choose m 2 / n on every measurable event, in both directions. For n > 0 that sample is the empirical-product mixture (sampleWithReplacement_eq_bind_pi_empiricalMeasureOfFintype).

theorem TauCeti.Probability.ExchangeableAt.finiteDeFinetti_exact_collision {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) {A : Set (Fin m → α)} (hA : MeasurableSet A) :
(prefixLaw μ X m) A ≤ (sampleWithReplacement (prefixLaw μ X n)) A + (1 - ↑(n.descFactorial m) / ↑n ^ m) ∧ (sampleWithReplacement (prefixLaw μ X n)) A ≤ (prefixLaw μ X m) A + (1 - ↑(n.descFactorial m) / ↑n ^ m)

Finite de Finetti with the exact collision error. The discrepancy on any measurable event between an exchangeable prefix and its with-replacement sample (the empirical-product mixture when n > 0) is bounded by 1 - (n)_m / n^m, where (n)_m counts ordered injections. This retains the overlap between pairwise collisions lost in the union bound. The empty-sample case m = n = 0 has error zero.

theorem TauCeti.Probability.ExchangeableAt.prefixLaw_le_sampleWithReplacement_add_exact_collision {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) {A : Set (Fin m → α)} (hA : MeasurableSet A) :
(prefixLaw μ X m) A ≤ (sampleWithReplacement (prefixLaw μ X n)) A + (1 - ↑(n.descFactorial m) / ↑n ^ m)

Exact finite de Finetti bound, exchangeable law to empirical mixture. Under finite exchangeability and m ≤ n, the prefix law exceeds its with-replacement sample on a measurable event by at most the exact index-collision error.

theorem TauCeti.Probability.ExchangeableAt.sampleWithReplacement_le_prefixLaw_add_exact_collision {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) {A : Set (Fin m → α)} (hA : MeasurableSet A) :
(sampleWithReplacement (prefixLaw μ X n)) A ≤ (prefixLaw μ X m) A + (1 - ↑(n.descFactorial m) / ↑n ^ m)

Exact finite de Finetti bound, empirical mixture to exchangeable law. Under finite exchangeability and m ≤ n, the with-replacement sample exceeds the prefix law on a measurable event by at most the exact index-collision error.

theorem TauCeti.Probability.ExchangeableAt.prefixLaw_le_sampleWithReplacement_add {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) {A : Set (Fin m → α)} (hA : MeasurableSet A) :
(prefixLaw μ X m) A ≤ (sampleWithReplacement (prefixLaw μ X n)) A + ↑(m.choose 2) / ↑n

Quantitative finite de Finetti bound, exchangeable law to empirical mixture. If the first n coordinates of a process are exchangeable and m ≤ n, then every measurable event under the m-prefix law has mass at most its mass under sampling m entries with replacement from the first n coordinates, plus choose m 2 / n. For n > 0 that sample is the mixture of m-fold products of the empirical distribution of the first n coordinates.

theorem TauCeti.Probability.ExchangeableAt.sampleWithReplacement_le_prefixLaw_add {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) {A : Set (Fin m → α)} (hA : MeasurableSet A) :
(sampleWithReplacement (prefixLaw μ X n)) A ≤ (prefixLaw μ X m) A + ↑(m.choose 2) / ↑n

Quantitative finite de Finetti bound, empirical mixture to exchangeable law. Under the same hypotheses, every measurable event under the with-replacement sample (for n > 0, the empirical-product mixture) has mass at most its mass under the m-prefix law plus choose m 2 / n.

theorem TauCeti.Probability.ExchangeableAt.etvDist_prefixLaw_sampleWithReplacement_le {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) :

Finite de Finetti theorem, total-variation form. If the first n coordinates of a process are exchangeable and m ≤ n, then the total variation distance between the m-prefix law and the with-replacement sample of the n-prefix law (for n > 0, the empirical-product mixture) is at most 2 * choose m 2 / n.

The normalization is Mathlib's: etvDist P Q is the signed-measure norm ‖P - Q‖, twice the distance sup_A |P A - Q A|, so in terms of that distance the bound is choose m 2 / n (etvDist_le_two_mul_iff).

theorem TauCeti.Probability.ExchangeableAt.etvDist_prefixLaw_sampleWithReplacement_le_exact_collision {α : Type u_1} [MeasurableSpace α] {Ω : Type u_2} [MeasurableSpace Ω] {μ : MeasureTheory.Measure Ω} [MeasureTheory.IsProbabilityMeasure μ] {X : ℕ → Ω → α} {m n : ℕ} (h : ExchangeableAt μ X n) (hmn : m ≤ n) (hX : ∀ (i : Fin n), AEMeasurable (X ↑i) μ) :

The finite de Finetti approximation with exact collision error, for Mathlib's total variation distance etvDist. Its signed-measure normalization is twice the probability-distance convention sup_A |P(A) - Q(A)|; no optimality for each individual population law is asserted.