Proposition 7 —
OpenModelRiskOT.Duality.proposition_7Let be a Polish space, a Borel probability measure on , a lower semicontinuous cost with if and only if (Assumption (A1)), upper semicontinuous and -integrable (Assumption (A2)), and . Let be a probability measure on such that
- (a) ;
- (b) ;
- (c) for every Borel set .
Let and let be the set of pairs with , universally measurable and for all (29). Then
Note that is not required to satisfy the budget . Proposition 7 carries the compact duality of Proposition 6 to the -compact set .
Formalization Note The space is a Polish space with its Borel σ-algebra; the cost is real-valued and written curried, ; (A1) is the structure AssumptionA1; (A2) is the pair of hypotheses UpperSemicontinuous f and Integrable f μ. Values that can be infinite (, , , , ) live in EReal; an integral of an extended-real function is with lower Lebesgue integrals, and evaluates to , so a coupling with never raises the primal supremum (the paper's footnote 2 reading).
import Mathlib import Definitions.Def_ModelRiskOT_Duality_AssumptionA1 import Definitions.Def_ModelRiskOT_Duality_primalValue import Definitions.Def_ModelRiskOT_Duality_dualValue import Definitions.Def_ModelRiskOT_Duality_admissibleSet open MeasureTheory
namespace ModelRiskOT.Duality
/-- **Proposition 7** (Blanchet & Murthy, arXiv:1604.01446v2, §4.3, p. 20). Under (A1) and (A2),
with `δ > 0`: let `π` be a probability measure on `S × S` with (a) `∫ c dπ < ∞`,
(b) `∫ f(y) dπ(x, y) ∈ (−∞, ∞)` and (c) `π(A × S) = μ(A)` for every Borel `A`. Then
`inf_{(λ, φ) ∈ Λ(S_π × S_π)} J(λ, φ) ≤ I`, where `S_π = Spt(π_X) ∪ Spt(π_Y)` and `Λ(K × K)` is
(29). No budget constraint `∫ c dπ ≤ δ` is assumed. -/
theorem proposition_7 {S : Type*} [TopologicalSpace S] [PolishSpace S] [MeasurableSpace S] [BorelSpace S]
(μ : Measure S) [IsProbabilityMeasure μ] (c : S → S → ℝ) (hc : AssumptionA1 c)
(f : S → ℝ) (hf_usc : UpperSemicontinuous f) (hf_int : Integrable f μ)
(δ : ℝ) (hδ : 0 < δ)
(π : Measure (S × S)) [IsProbabilityMeasure π]
(ha : ∫⁻ p, ENNReal.ofReal (c p.1 p.2) ∂π < ⊤)
(hb : Integrable (fun p => f p.2) π)
(hmarg : ∀ A : Set S, MeasurableSet A → π (A ×ˢ Set.univ) = μ A) :
(⨅ p ∈ dualFeasible c f (supportUnion π), dualObj μ δ p.1 p.2) ≤ primalValue c f μ δ := by sorry
end ModelRiskOT.Duality