Radon-Nikodym derivative and Lebesgue decomposition for kernels #
Let α and γ be two measurable space, where either α is countable or γ is
countably generated. Let κ, η : Kernel α γ be finite kernels.
Then there exists a function Kernel.rnDeriv κ η : α → γ → ℝ≥0∞ jointly measurable on α × γ
and a kernel Kernel.singularPart κ η : Kernel α γ such that
κ = Kernel.withDensity η (Kernel.rnDeriv κ η) + Kernel.singularPart κ η,- for all
a : α,Kernel.singularPart κ η a ⟂ₘ η a, - for all
a : α,Kernel.singularPart κ η a = 0 ↔ κ a ≪ η a, - for all
a : α,Kernel.withDensity η (Kernel.rnDeriv κ η) a = 0 ↔ κ a ⟂ₘ η a.
Furthermore, the sets {a | κ a ≪ η a} and {a | κ a ⟂ₘ η a} are measurable.
When γ is countably generated, the construction of the derivative starts from Kernel.density:
for two finite kernels κ' : Kernel α (γ × β) and η' : Kernel α γ with fst κ' ≤ η',
the function density κ' η' : α → γ → Set β → ℝ is jointly measurable in the first two arguments
and satisfies that for all a : α and all measurable sets s : Set β and A : Set γ,
∫ x in A, density κ' η' a x s ∂(η' a) = (κ' a (A ×ˢ s)).toReal.
We use that definition for β = Unit and κ' = map κ (fun a ↦ (a, ())). We can't choose η' = η
in general because we might not have κ ≤ η, but if we could, we would get a measurable function
f with the property κ = withDensity η f, which is the decomposition we want for κ ≤ η.
To circumvent that difficulty, we take η' = κ + η and thus define rnDerivAux κ η.
Finally, rnDeriv κ η a x is given by
ENNReal.ofReal (rnDerivAux κ (κ + η) a x) / ENNReal.ofReal (1 - rnDerivAux κ (κ + η) a x).
Up to some conversions between ℝ and ℝ≥0, the singular part is
withDensity (κ + η) (rnDerivAux κ (κ + η) - (1 - rnDerivAux κ (κ + η)) * rnDeriv κ η).
The countably generated measurable space assumption is not needed to have a decomposition for
measures, but the additional difficulty with kernels is to obtain joint measurability of the
derivative. This is why we can't simply define rnDeriv κ η by a ↦ (κ a).rnDeriv (ν a)
everywhere unless α is countable (although rnDeriv κ η has that value almost everywhere).
See the construction of Kernel.density for details on how the countably generated hypothesis
is used.
Main definitions #
ProbabilityTheory.Kernel.rnDeriv: a functionα → γ → ℝ≥0∞jointly measurable onα × γProbabilityTheory.Kernel.singularPart: aKernel α γ
Main statements #
ProbabilityTheory.Kernel.mutuallySingular_singularPart: for alla : α,Kernel.singularPart κ η a ⟂ₘ η aProbabilityTheory.Kernel.rnDeriv_add_singularPart:Kernel.withDensity η (Kernel.rnDeriv κ η) + Kernel.singularPart κ η = κProbabilityTheory.Kernel.measurableSet_absolutelyContinuous: the set{a | κ a ≪ η a}is MeasurableProbabilityTheory.Kernel.measurableSet_mutuallySingular: the set{a | κ a ⟂ₘ η a}is Measurable
Uniqueness results: if κ = η.withDensity f + ξ for measurable f and ξ is such that
ξ a ⟂ₘ η a for some a : α then
ProbabilityTheory.Kernel.eq_rnDeriv:f a =ᵐ[η a] Kernel.rnDeriv κ η aProbabilityTheory.Kernel.eq_singularPart:ξ a = Kernel.singularPart κ η a
References #
Theorem 1.28 in [O. Kallenberg, Random Measures, Theory and Applications][kallenberg2017].
Auxiliary function used to define ProbabilityTheory.Kernel.rnDeriv and
ProbabilityTheory.Kernel.singularPart.
This has the properties we want for a Radon-Nikodym derivative only if κ ≪ ν. The definition of
rnDeriv κ η will be built from rnDerivAux κ (κ + η).
Equations
Instances For
A set of points in α × γ related to the absolute continuity / mutual singularity of
κ and η.
Equations
- κ.mutuallySingularSet η = {p : α × γ | 1 ≤ κ.rnDerivAux (κ + η) p.1 p.2}
Instances For
A set of points in α × γ related to the absolute continuity / mutual singularity of
κ and η. That is,
withDensity η (rnDeriv κ η) a (mutuallySingularSetSlice κ η a) = 0,singularPart κ η a (mutuallySingularSetSlice κ η a)ᶜ = 0.
Equations
- κ.mutuallySingularSetSlice η a = {x : γ | 1 ≤ κ.rnDerivAux (κ + η) a x}
Instances For
Alias of ProbabilityTheory.Kernel.notMem_mutuallySingularSetSlice.
Radon-Nikodym derivative of the kernel κ with respect to the kernel η.
Equations
- κ.rnDeriv η a x = ENNReal.ofReal (κ.rnDerivAux (κ + η) a x) / ENNReal.ofReal (1 - κ.rnDerivAux (κ + η) a x)
Instances For
Singular part of the kernel κ with respect to the kernel η.
Equations
- κ.singularPart η = (κ + η).withDensity fun (a : α) (x : γ) => ↑(κ.rnDerivAux (κ + η) a x).toNNReal - ↑(1 - κ.rnDerivAux (κ + η) a x).toNNReal * κ.rnDeriv η a x
Instances For
The singular part of κ with respect to η is mutually singular with η.
Lebesgue decomposition of a finite kernel κ with respect to another one η.
κ is the sum of an absolutely continuous part withDensity η (rnDeriv κ η) and a singular part
singularPart κ η.
The set of points a : α such that κ a ≪ η a is measurable.
The set of points a : α such that κ a ⟂ₘ η a is measurable.
For two kernels κ, η, the singular part of κ a with respect to η a is a measurable
function of a.