Rao–Blackwell theorem
TheoremIntegrability of a convex loss of a conditional expectation#✓fully proved
For a convex loss L, L ∘ μ[W | m] is integrable whenever W and L ∘ W are.
theorem integrable_comp_condExp {m m₀ : MeasurableSpace Ω} (hm : m ≤ m₀)
(μ : Measure[m₀] Ω) [IsProbabilityMeasure μ] {L : ℝ → ℝ} (hL : ConvexOn ℝ Set.univ L)
{W : Ω → ℝ} (hW : Integrable W μ) (hLW : Integrable (fun ω => L (W ω)) μ) :
Integrable (fun ω => L (μ[W | m] ω)) μ
Used by 1
TheoremRao–Blackwell inequality for one law#✓fully proved
For a convex loss L and an integrable W with
integrable loss, ∫ L(μ[W | m]) dμ ≤ ∫ L(W) dμ.
theorem integral_comp_condExp_le {m m₀ : MeasurableSpace Ω} (hm : m ≤ m₀)
(μ : Measure[m₀] Ω) [IsProbabilityMeasure μ] {L : ℝ → ℝ} (hL : ConvexOn ℝ Set.univ L)
{W : Ω → ℝ} (hW : Integrable W μ) (hLW : Integrable (fun ω => L (W ω)) μ) :
∫ ω, L (μ[W | m] ω) ∂μ ≤ ∫ ω, L (W ω) ∂μ
Uses 1 · Used by 1
TheoremSquared-error decomposition#✓fully proved
For W with a finite second moment and any target c,
E (W - c)² = E (μ[W | m] - c)² + E (W - μ[W | m])². This is the Pythagorean identity
MeasureTheory.integral_sq_sub_eq_integral_sq_sub_condExp_add with the constant predictor
fun _ => c, with the two summands commuted; it is not a second Pythagorean proof.
theorem integral_sq_sub_eq_add {m m₀ : MeasurableSpace Ω} (hm : m ≤ m₀)
(μ : Measure[m₀] Ω) [IsProbabilityMeasure μ] {W : Ω → ℝ} (hW : MemLp W 2 μ) (c : ℝ) :
∫ ω, (W ω - c) ^ 2 ∂μ =
∫ ω, (μ[W | m] ω - c) ^ 2 ∂μ + ∫ ω, (W ω - μ[W | m] ω) ^ 2 ∂μ
Uses 1 · Used by 2
TheoremEquality case under squared error#✓fully proved
Conditioning leaves the squared-error risk at c
unchanged if and only if W already equals μ[W | m] almost everywhere.
theorem integral_sq_sub_condExp_eq_iff {m m₀ : MeasurableSpace Ω} (hm : m ≤ m₀)
(μ : Measure[m₀] Ω) [IsProbabilityMeasure μ] {W : Ω → ℝ} (hW : MemLp W 2 μ) (c : ℝ) :
∫ ω, (μ[W | m] ω - c) ^ 2 ∂μ = ∫ ω, (W ω - c) ^ 2 ∂μ ↔ W =ᵐ[μ] μ[W | m]
Uses 1
TheoremThe Rao–Blackwell theorem#✓fully proved
Let m be sufficient for the family of probability measures
P, and let W be a measurable estimator that is integrable under every P θ. There is one
parameter-free m-measurable estimator η, a version of P_θ[W | m] for every θ, with
the same mean as W under every θ and no larger risk under every loss ℓ θ that is convex
in the estimate (for which ℓ θ ∘ W is integrable).
theorem IsSufficient.raoBlackwell {m m₀ : MeasurableSpace Ω}
{P : Θ → Measure[m₀] Ω} [∀ θ, IsProbabilityMeasure (P θ)] (hP : IsSufficient m P)
{W : Ω → ℝ} (hWm : Measurable[m₀] W) (hW : ∀ θ, Integrable W (P θ)) :
∃ η : Ω → ℝ, StronglyMeasurable[m] η ∧ (∀ θ, η =ᵐ[P θ] (P θ)[W | m]) ∧
(∀ θ, ∫ ω, η ω ∂P θ = ∫ ω, W ω ∂P θ) ∧
∀ ℓ : Θ → ℝ → ℝ, ∀ θ, ConvexOn ℝ Set.univ (ℓ θ) →
Integrable (fun ω => ℓ θ (W ω)) (P θ) →
∫ ω, ℓ θ (η ω) ∂P θ ≤ ∫ ω, ℓ θ (W ω) ∂P θ
def IsSufficient (m : MeasurableSpace Ω) {m₀ : MeasurableSpace Ω}
(P : Θ → Measure[m₀] Ω) : Prop :=
m ≤ m₀ ∧ (∀ θ, IsFiniteMeasure (P θ)) ∧
∀ s, MeasurableSet[m₀] s → ∃ f : Ω → ℝ, StronglyMeasurable[m] f ∧
∀ θ, f =ᵐ[P θ] (P θ)[s.indicator (fun _ => (1 : ℝ)) | m]
MathlibExt/Probability/Statistics/Sufficiency/Basic.lean:74
Sources
- Rao (1945); Blackwell (1947); Casella–Berger, Section 7.3.
Uses 3
TheoremThe Rao–Blackwell theorem for a sufficient statistic#✓fully proved
If the statistic T is
sufficient for P and W is a measurable estimator with finite second moment under every
P θ, there is a measurable φ such that φ ∘ T is a version of every E_θ[W | T], has the
same mean as W, and for every target τ θ satisfies
E_θ (W - τ θ)² = E_θ (φ(T) - τ θ)² + E_θ (W - φ(T))².
theorem raoBlackwell_comap {τ : Type*} [MeasurableSpace τ] {m₀ : MeasurableSpace Ω}
{T : Ω → τ} (hT : Measurable T) {P : Θ → Measure[m₀] Ω} [∀ θ, IsProbabilityMeasure (P θ)]
(hP : IsSufficient (MeasurableSpace.comap T ‹MeasurableSpace τ›) P) {W : Ω → ℝ}
(hWm : Measurable[m₀] W) (hW : ∀ θ, MemLp W 2 (P θ)) (tgt : Θ → ℝ) :
∃ φ : τ → ℝ, Measurable φ ∧
(∀ θ, (fun ω => φ (T ω)) =ᵐ[P θ] (P θ)[W | MeasurableSpace.comap T ‹_›]) ∧
(∀ θ, ∫ ω, φ (T ω) ∂P θ = ∫ ω, W ω ∂P θ) ∧
∀ θ, ∫ ω, (W ω - tgt θ) ^ 2 ∂P θ =
∫ ω, (φ (T ω) - tgt θ) ^ 2 ∂P θ + ∫ ω, (W ω - φ (T ω)) ^ 2 ∂P θ
def IsSufficient (m : MeasurableSpace Ω) {m₀ : MeasurableSpace Ω}
(P : Θ → Measure[m₀] Ω) : Prop :=
m ≤ m₀ ∧ (∀ θ, IsFiniteMeasure (P θ)) ∧
∀ s, MeasurableSet[m₀] s → ∃ f : Ω → ℝ, StronglyMeasurable[m] f ∧
∀ θ, f =ᵐ[P θ] (P θ)[s.indicator (fun _ => (1 : ℝ)) | m]
MathlibExt/Probability/Statistics/Sufficiency/Basic.lean:74
Uses 3
Loading dependency explorer…
The dependency explorer needs JavaScript and its graph data.