Documentation

Papers.Rockel2026ExactBlest.ExactBlestPaperParams

← Mathematical handbook

The manuscript's parametrization of the extremal families #

In exact-blest-regions.tex every family parameter increases with dependence:

The other modules work with the reflected parameters 1-c, 1-w, and a = 1-b. This module defines the families exactly as printed (paperD, paperA, paperB) and restates every family-dependent claim of the manuscript in that notation: the coefficient formulas, parameter monotonicity and derivatives, graph laws and rank maps, the potentials of Propositions 4.3 and 4.4 with their contact sets, the extremizers in Theorems 1.1 and 1.2 and Corollaries 3.2, 4.5, and 4.6, Remark 4.7, and all values printed above the panels of Figures 2 and 3.

Reflection of the parameter #

The family A_w #

The manuscript's A_w: under its reflected-coordinate coupling the leading fraction w of the first variable is comonotone with the trailing fraction of the second, and the remaining bottom ranks are countermonotone with the top ranks of the second.

Equations
Instances For

    Lemma 4.1: ν(A_w) - η(A_w) = 2(1-w)w³.

    Lemma 4.1: w ↦ η(A_w) is strictly increasing.

    theorem Papers.Rockel2026ExactBlest.paperA_eta_range (w : ↑unitInterval) (hw : 1 / 2 ≤ ↑w) :
    eta (paperA w) ∈ Set.Icc (-3 / 4) 1

    Lemma 4.1: for w ∈ [1/2,1], η(A_w) ranges over [-3/4,1].

    theorem Papers.Rockel2026ExactBlest.paperT_formula (w u : ↑unitInterval) :
    ↑(graphRank (unitInterval.symm w) u) = if ↑u ≤ 1 - ↑w then 1 - ↑u else ↑u + ↑w - 1

    The map T_w of Section 4.1.

    The family B_b #

    noncomputable def Papers.Rockel2026ExactBlest.paperB (b : ↑unitInterval) (hb : ↑b ≤ 1 / 2) :

    The manuscript's B_b for b ∈ [0,1/2]: the leading fraction b of the first variable is split between the two branches of R_b, the rest is countermonotone. The value b = 0 gives W, the continuous endpoint of the family.

    Equations
    Instances For
      theorem Papers.Rockel2026ExactBlest.paperB_symm (a : ↑unitInterval) (ha : 1 / 2 ≤ ↑a) (hb : ↑(unitInterval.symm a) ≤ 1 / 2) :

      e_b in (6) and (17).

      Equations
      Instances For

        ν(B_b) in (17), equivalently Λ(n_b) in (22).

        Equations
        Instances For

          n_b = ν(B_bᵀ) in (22).

          Equations
          Instances For

            Υ(e_b) in (6).

            Equations
            Instances For
              theorem Papers.Rockel2026ExactBlest.eta_paperB (b : ↑unitInterval) (hb : ↑b ≤ 1 / 2) :
              eta (paperB b hb) = paperEtaB ↑b

              Lemma 4.1: ν(B_b) - η(B_b) = Υ(e_b) with the displayed formula (6).

              For b = 1/2, the second branch carries no mass and B_{1/2} = A_{1/2}.

              Lemma 4.1: b ↦ e_b is strictly increasing on [0,1/2].

              theorem Papers.Rockel2026ExactBlest.hasDerivAt_paperEtaB (b : ℝ) (hb : b ≠ 1) :
              HasDerivAt paperEtaB (b ^ 2 * (2 - b) * (3 * b ^ 2 - 8 * b + 6) / (4 * (1 - b) ^ 3)) b

              Lemma 4.1: the displayed derivative ∂_b η(B_b).

              theorem Papers.Rockel2026ExactBlest.deriv_paperEtaB_pos (b : ℝ) (hb0 : 0 < b) (hb : b < 1) :
              0 < b ^ 2 * (2 - b) * (3 * b ^ 2 - 8 * b + 6) / (4 * (1 - b) ^ 3)

              Lemma 4.1: e_{1/2} = -3/4 and e_b → -1 as b → 0.

              theorem Papers.Rockel2026ExactBlest.paperEtaB_range (b : ↑unitInterval) (hb : ↑b ≤ 1 / 2) :
              paperEtaB ↑b ∈ Set.Icc (-1) (-3 / 4)
              noncomputable def Papers.Rockel2026ExactBlest.paperR (b z : ℝ) :

              The map R_b of Section 4.1.

              Equations
              Instances For
                theorem Papers.Rockel2026ExactBlest.cutB_paper (b : ℝ) :
                cutB (1 - b) = b / (2 * (1 - b))
                theorem Papers.Rockel2026ExactBlest.kA_paper (w : ℝ) :
                kA (1 - w) = (3 - 2 * w) / (2 * w)
                theorem Papers.Rockel2026ExactBlest.kB_paper (b : ℝ) :
                kB (1 - b) = 2 * (1 - b) / b

                B_b is the copula associated with the law of (R_b(Z), Z).

                theorem Papers.Rockel2026ExactBlest.paperR_formula (b : ↑unitInterval) (hb0 : 0 < ↑b) (hb : ↑b < 1 / 2) (z : ↑unitInterval) :
                ↑(randomRank ↑(unitInterval.symm b) ⋯ ⋯ z) = paperR ↑b ↑z
                theorem Papers.Rockel2026ExactBlest.paperB_conditional_formula (b : ↑unitInterval) (hb0 : 0 < ↑b) (hb : ↑b ≤ 1 / 2) (x : ↑unitInterval) :
                MeasureTheory.Measure.map (fun (z : ↑unitInterval) => ↑z) ((blockKernel (unitInterval.symm b) (randomSplit (unitInterval.symm b) ⋯)) x) = if ↑x ≤ 1 - ↑b then MeasureTheory.Measure.dirac (1 - ↑x) else ENNReal.ofReal (1 / (2 * (1 - ↑b))) • MeasureTheory.Measure.dirac ((↑x + ↑b - 1) / (2 * (1 - ↑b))) + ENNReal.ofReal ((1 - 2 * ↑b) / (2 * (1 - ↑b))) • MeasureTheory.Measure.dirac ((1 - ↑b - (1 - 2 * ↑b) * ↑x) / (2 * (1 - ↑b)))

                The two conditional atoms (15) and their probabilities for B_b.

                The parameters n_b and Λ #

                theorem Papers.Rockel2026ExactBlest.one_sub_bijOn :
                Set.BijOn (fun (b : ℝ) => 1 - b) (Set.Ioo 0 (1 / 2)) (Set.Ioo (1 / 2) 1)

                Section 4.3: b ↦ n_b increases bijectively from (0,1/2) onto (-1,-7/8).

                theorem Papers.Rockel2026ExactBlest.Lambda_paperNB (b : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) :

                (22): Λ(n_b) = ν(B_b) on the parametric branch.

                Derivatives of Υ along the parametric branch #

                theorem Papers.Rockel2026ExactBlest.hasDerivAt_paperUps (b : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) :
                HasDerivAt (fun (t : ℝ) => etaGap (paperEtaB t)) (b ^ 2 * (2 - 3 * b) * (3 * b ^ 2 - 8 * b + 6) / (4 * (1 - b) ^ 3)) b

                Proof of Corollary 4.5: ∂_b Υ(e_b) = b²(2-3b)(3b²-8b+6)/(4(1-b)³) > 0.

                theorem Papers.Rockel2026ExactBlest.paperUps_deriv_pos (b : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) :
                0 < b ^ 2 * (2 - 3 * b) * (3 * b ^ 2 - 8 * b + 6) / (4 * (1 - b) ^ 3)
                theorem Papers.Rockel2026ExactBlest.hasDerivAt_etaGap_paper (b : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) :
                HasDerivAt etaGap ((2 - 3 * b) / (2 - b)) (paperEtaB b) ∧ (2 - 3 * b) / (2 - b) ∈ Set.Ioo (1 / 3) 1

                Proof of Corollary 4.6: Υ'(e_b) = (2-3b)/(2-b) ∈ (1/3,1).

                theorem Papers.Rockel2026ExactBlest.etaGap_paperEtaB (b : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) :

                (6): the parametric formula for Υ(e_b).

                The family D_c #

                The manuscript's D_c: the copula of (X, ζ_c(X)) for the rank map ζ_c centred at 1-c, so that D_0 = W, D_1 = M, and the support has its kink at u = c.

                Equations
                Instances For
                  theorem Papers.Rockel2026ExactBlest.paperD_values (c : ↑unitInterval) :
                  ((paperD c).spearmanRho = if ↑c ≤ 1 / 2 then 8 * ↑c ^ 3 - 1 else 1 - 8 * (1 - ↑c) ^ 3) ∧ Rockel2026XiBlest.blestNu (paperD c) = if ↑c ≤ 1 / 2 then 16 * ↑c ^ 3 - 12 * ↑c ^ 4 - 1 else 1 - 12 * (1 - ↑c) ^ 4

                  Lemma 3.1, displayed formulas (12).

                  Lemma 3.1: c ↦ ρ(D_c) is strictly increasing.

                  Lemma 3.1: c ↦ ρ(D_c) is a bijection of [0,1] onto [-1,1].

                  theorem Papers.Rockel2026ExactBlest.paperZeta_formula (c x : ↑unitInterval) :
                  ↑(quadraticRank (↑(unitInterval.symm c)) x) = if |↑x + ↑c - 1| ≤ min (↑c) (1 - ↑c) then 2 * |↑x + ↑c - 1| else if 2 - 2 * ↑c ≤ ↑x then ↑x else 1 - ↑x

                  The rank map ζ_c of (11), centred at 1-c.

                  Theorem 1.1 and its proof: D_c is the only copula with ρ = ρ(D_c) on the upper boundary, and its survival copula is the only one on the lower boundary.

                  Corollary 3.2: the maximum ν - ρ = 1/4 is attained only by D_{1/2}.

                  Theorem 1.2 in the manuscript's parametrization #

                  theorem Papers.Rockel2026ExactBlest.cube_root_spec (e : ℝ) (he : e ∈ Set.Icc (-3 / 4) 1) :
                  1 / 2 ≤ ((1 + e) / 2) ^ (1 / 3) ∧ ((1 + e) / 2) ^ (1 / 3) ≤ 1 ∧ (((1 + e) / 2) ^ (1 / 3)) ^ 3 = (1 + e) / 2
                  theorem Papers.Rockel2026ExactBlest.paper_eta_graph_extremizers (e : ℝ) (he : e ∈ Set.Icc (-3 / 4) 1) :
                  ∃ (w : ↑unitInterval), ↑w = ((1 + e) / 2) ^ (1 / 3) ∧ 1 / 2 ≤ ↑w ∧ eta (paperA w) = e ∧ Rockel2026XiBlest.blestNu (paperA w) = e + etaGap e ∧ (∀ (C : ProbabilityTheory.Copula 2), eta C = e → Rockel2026XiBlest.blestNu C = e + etaGap e → C = paperA w) ∧ ∀ (C : ProbabilityTheory.Copula 2), eta C = e → Rockel2026XiBlest.blestNu C = e - etaGap e → C = (paperA w).transpose

                  Theorem 1.2 for η ∈ [-3/4,1]: the unique upper extremizer is A_w with w = ((1+η)/2)^{1/3} ∈ [1/2,1], and the unique lower one is its transpose.

                  theorem Papers.Rockel2026ExactBlest.paper_eta_random_extremizers (e : ℝ) (he : e ∈ Set.Ioo (-1) (-3 / 4)) :
                  ∃ (b : ↑unitInterval) (hb : ↑b < 1 / 2), 0 < ↑b ∧ paperEtaB ↑b = e ∧ eta (paperB b ⋯) = e ∧ Rockel2026XiBlest.blestNu (paperB b ⋯) = e + etaGap e ∧ (∀ (C : ProbabilityTheory.Copula 2), eta C = e → Rockel2026XiBlest.blestNu C = e + etaGap e → C = paperB b ⋯) ∧ ∀ (C : ProbabilityTheory.Copula 2), eta C = e → Rockel2026XiBlest.blestNu C = e - etaGap e → C = (paperB b ⋯).transpose

                  Theorem 1.2 for η ∈ (-1,-3/4): the unique upper extremizer is B_b with e_b = η and b ∈ (0,1/2), and the unique lower one is its transpose.

                  Potentials of Propositions 4.3 and 4.4 #

                  noncomputable def Papers.Rockel2026ExactBlest.paperPhiA (w x : ℝ) :

                  φ_w and ψ_w of (18).

                  Equations
                  Instances For
                    theorem Papers.Rockel2026ExactBlest.hasDerivAt_paperPhiA (w x : ℝ) (hx : x ≠ 1 - w) :
                    HasDerivAt (paperPhiA w) (if x ≤ 1 - w then (1 - x) * (2 * x - (3 - 2 * w) / (2 * w) * (1 - x)) else (x + w - 1) * (2 * x - (3 - 2 * w) / (2 * w) * (x + w - 1))) x
                    theorem Papers.Rockel2026ExactBlest.hasDerivAt_paperPsiA (w z : ℝ) (hz : z ≠ w) :
                    HasDerivAt (paperPsiA w) (if z ≤ w then (z + 1 - w) * (z + 1 - w - 2 * ((3 - 2 * w) / (2 * w)) * z) else (1 - z) * (1 - z - 2 * ((3 - 2 * w) / (2 * w)) * z)) z
                    theorem Papers.Rockel2026ExactBlest.paper_dualA (w x z : ℝ) (hw : 1 / 2 ≤ w) (hw1 : w ≤ 1) (hx : x ∈ Set.Icc 0 1) (hz : z ∈ Set.Icc 0 1) :
                    cost ((3 - 2 * w) / (2 * w)) x z ≤ paperPhiA w x + paperPsiA w z

                    Proposition 4.3 (inequality) with κ_w = (3-2w)/(2w), w ∈ [1/2,1].

                    theorem Papers.Rockel2026ExactBlest.paper_contactA (w x z : ℝ) (hw : 1 / 2 ≤ w) (hw1 : w < 1) (hx : x ∈ Set.Icc 0 1) (hz : z ∈ Set.Icc 0 1) :
                    paperPhiA w x + paperPsiA w z - cost ((3 - 2 * w) / (2 * w)) x z = 0 ↔ x ≤ 1 - w ∧ z = 1 - x ∨ 1 - w ≤ x ∧ z = x + w - 1 ∨ w = 1 / 2 ∧ 1 / 2 ≤ x ∧ z = 1 / 2

                    Proposition 4.3 (equality set): S_w, plus [1/2,1] × {1/2} when w = 1/2.

                    noncomputable def Papers.Rockel2026ExactBlest.paperPhiB (b x : ℝ) :

                    φ_b and ψ_b of (20).

                    Equations
                    Instances For
                      theorem Papers.Rockel2026ExactBlest.paperPotentialsB_normalized (b : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) :
                      paperPhiB b (1 - b) = 0 ∧ paperPsiB b 0 = 0
                      theorem Papers.Rockel2026ExactBlest.hasDerivAt_paperPhiB (b x : ℝ) (hb : b ≠ 1) (hx : x ≠ 1 - b) :
                      HasDerivAt (paperPhiB b) (if x ≤ 1 - b then (1 - x) * (2 * x - 2 * (1 - b) / b * (1 - x)) else (x + b - 1) / (2 * (1 - b)) * (2 * x - 2 * (1 - b) / b * ((x + b - 1) / (2 * (1 - b))))) x
                      theorem Papers.Rockel2026ExactBlest.hasDerivAt_paperPsiB (b z : ℝ) (hzc : z ≠ b / (2 * (1 - b))) (hzb : z ≠ b) :
                      HasDerivAt (paperPsiB b) (paperR b z * (paperR b z - 2 * (2 * (1 - b) / b) * z)) z
                      theorem Papers.Rockel2026ExactBlest.paper_dualB (b x z : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) (hx : x ∈ Set.Icc 0 1) (hz : z ∈ Set.Icc 0 1) :
                      cost (2 * (1 - b) / b) x z ≤ paperPhiB b x + paperPsiB b z

                      Proposition 4.4 (inequality) with κ_b = 2(1-b)/b, b ∈ (0,1/2).

                      theorem Papers.Rockel2026ExactBlest.paper_contactB (b x z : ℝ) (hb0 : 0 < b) (hb : b < 1 / 2) (hx : x ∈ Set.Icc 0 1) (hz : z ∈ Set.Icc 0 1) :
                      paperPhiB b x + paperPsiB b z - cost (2 * (1 - b) / b) x z = 0 ↔ x = paperR b z

                      Proposition 4.4 (equality set): exactly the graph x = R_b(z).

                      Corollaries 4.5 and 4.6 and Remark 4.7 #

                      Corollary 4.5: |ν - η| = 27/128 exactly for A_{3/4} and A_{3/4}ᵀ.

                      Corollary 4.6 for n ≥ -7/8: the maximizer is A_wᵀ with w = ((1+n)/2)^{1/4}.

                      Corollary 4.6 for n < -7/8: the maximizer is B_bᵀ with n_b = n.

                      theorem Papers.Rockel2026ExactBlest.paper_nu_transpose_extremizer_cases (C : ProbabilityTheory.Copula 2) (hn : -1 < Rockel2026XiBlest.blestNu C) :
                      (Rockel2026XiBlest.blestNu C.transpose = Lambda (Rockel2026XiBlest.blestNu C) → (∃ (w : ↑unitInterval), 1 / 2 ≤ ↑w ∧ C = (paperA w).transpose) ∨ ∃ (b : ↑unitInterval) (hb : ↑b < 1 / 2), 0 < ↑b ∧ C = (paperB b ⋯).transpose) ∧ (Rockel2026XiBlest.blestNu C.transpose = LambdaInv (Rockel2026XiBlest.blestNu C) → (∃ (w : ↑unitInterval), 1 / 2 ≤ ↑w ∧ C = paperA w) ∨ ∃ (b : ↑unitInterval) (hb : ↑b < 1 / 2), 0 < ↑b ∧ C = paperB b ⋯)

                      Corollary 4.6: the maximizers of νᵀ at given ν > -1 are the transposes A_wᵀ (w ∈ [1/2,1]) or B_bᵀ (b ∈ (0,1/2)), the minimizers the families themselves.

                      theorem Papers.Rockel2026ExactBlest.paperA_not_extremal (w : ↑unitInterval) (hw0 : 0 < ↑w) (hw : ↑w < 1 / 2) :

                      Remark 4.7(a): for w ∈ (0,1/2), A_w lies strictly below the upper boundary.

                      Remark 4.7(b): the unique maximizer of ρ - η is A_{3/4}^⊥, the shuffle of M that is countermonotone on [0,3/4] and comonotone on [3/4,1]; the minimizer is (A_{3/4}ᵀ)^⊥, its survival copula.

                      Values printed above the panels of Figures 2 and 3 #

                      theorem Papers.Rockel2026ExactBlest.paper_figure_panel_values :
                      ((paperD ⟨1 / 4, ⋯⟩).spearmanRho = -7 / 8 ∧ Rockel2026XiBlest.blestNu (paperD ⟨1 / 4, ⋯⟩) = -51 / 64 ∧ Rockel2026XiBlest.blestNu (paperD ⟨1 / 4, ⋯⟩).survivalCopula = -61 / 64) ∧ ((paperD ⟨1 / 2, ⋯⟩).spearmanRho = 0 ∧ Rockel2026XiBlest.blestNu (paperD ⟨1 / 2, ⋯⟩) = 1 / 4 ∧ Rockel2026XiBlest.blestNu (paperD ⟨1 / 2, ⋯⟩).survivalCopula = -1 / 4) ∧ ((paperD ⟨3 / 4, ⋯⟩).spearmanRho = 7 / 8 ∧ Rockel2026XiBlest.blestNu (paperD ⟨3 / 4, ⋯⟩) = 61 / 64 ∧ Rockel2026XiBlest.blestNu (paperD ⟨3 / 4, ⋯⟩).survivalCopula = 51 / 64) ∧ (eta (paperB ⟨1 / 4, ⋯⟩ ⋯) = -2257 / 2304 ∧ Rockel2026XiBlest.blestNu (paperB ⟨1 / 4, ⋯⟩ ⋯) = -185 / 192 ∧ Rockel2026XiBlest.blestNu (paperB ⟨1 / 4, ⋯⟩ ⋯).transpose = -1147 / 1152) ∧ (eta (paperA ⟨1 / 2, ⋯⟩) = -3 / 4 ∧ Rockel2026XiBlest.blestNu (paperA ⟨1 / 2, ⋯⟩) = -5 / 8 ∧ Rockel2026XiBlest.blestNu (paperA ⟨1 / 2, ⋯⟩).transpose = -7 / 8) ∧ (eta (paperA ⟨7 / 8, ⋯⟩) = 87 / 256 ∧ Rockel2026XiBlest.blestNu (paperA ⟨7 / 8, ⋯⟩) = 1039 / 2048 ∧ Rockel2026XiBlest.blestNu (paperA ⟨7 / 8, ⋯⟩).transpose = 353 / 2048) ∧ 1 / (2 * (1 - 1 / 4)) = 2 / 3 ∧ (1 - 2 * (1 / 4)) / (2 * (1 - 1 / 4)) = 1 / 3

                      Figure 2 (columns D_{1/4}, D_{1/2}, D_{3/4}) and Figure 3 (columns B_{1/4}, A_{1/2} = B_{1/2}, A_{7/8}), each with the value in the bottom row, and the branch probabilities 2/3 and 1/3 of B_{1/4}. The decimals printed for B_{1/4} and A_{7/8} are roundings of these fractions.