Documentation

LeanPool.RareSubtreePatterns.Bridge.RootMicroscopic

Inherited microscopicity for declumped marks: Lemma lem:root_declumping(iii) #

Work package WP3.5 of the bridge plan. Source lines 2197–2299 (Lemma 5.1, proof step P16). Parts (i), (ii) and the parent formula are in RootDeclumping/ParentFormula.lean and RootDeclumping/Lemma51.lean; this file proves part (iii).

Notation: Ghat = declumpedMark B G (root degree in B and some marked child), D the root degree, 𝒒_β„“ a forest of β„“ independent BienaymΓ© trees, c_B = cB ΞΌ B.

Main statements #

Departure from the paper (periodic case, correction F3): the bridge-weighted estimate is used in numerator form with the time shift j = 1 of consequenceM_weighted_sup, instead of "replacing n by n - 1" in a ratio with Ξ¨_{n-1}(-1) (which vanishes for d β‰₯ 2 at admissible n); the denominator Ξ¨_n(-1) β‰₯ c / a_n enters only in microscopic_of_large.

The branching decomposition at a marked child #

E[f(𝒯)] = βˆ‘_T w(T) f(T).

theorem PoissonDeclumping.Bridge.tsum_weight_mul_eq_sigma (ΞΌ : PMF β„•) (F : PlaneTree β†’ ENNReal) :
βˆ‘' (T : PlaneTree), PlaneTree.weight ΞΌ T * F T = βˆ‘' (d : β„•) (f : Fin d β†’ PlaneTree), (ΞΌ d * ∏ i : Fin d, PlaneTree.weight ΞΌ (f i)) * F (PlaneTree.node (List.ofFn f))

Decomposition of βˆ‘_T w(T) F(T) according to the root degree and the children.

theorem PoissonDeclumping.Bridge.tsum_pi_succAbove (d : β„•) (i : Fin (d + 1)) (F : PlaneTree β†’ (Fin d β†’ PlaneTree) β†’ ENNReal) :
(βˆ‘' (f : Fin (d + 1) β†’ PlaneTree), F (f i) fun (j : Fin d) => f (i.succAbove j)) = βˆ‘' (t : PlaneTree) (g : Fin d β†’ PlaneTree), F t g

Separating the i-th child: for F depending on the i-th tree and the others, βˆ‘_{f : Fin (d+1) β†’ 𝕋} F(f_i, (f_j)_{j β‰  i}) = βˆ‘_t βˆ‘_g F(t, g).

theorem PoissonDeclumping.Bridge.tsum_child_mark (ΞΌ : PMF β„•) (d : β„•) (i : Fin (d + 1)) (G : PlaneTree β†’ Prop) (Ξ¦ : β„• β†’ ENNReal) :
βˆ‘' (f : Fin (d + 1) β†’ PlaneTree), (∏ j : Fin (d + 1), PlaneTree.weight ΞΌ (f j)) * ({t : PlaneTree | G t}.indicator 1 (f i) * Ξ¦ (1 + βˆ‘ j : Fin (d + 1), (f j).vertexCount)) = ∫⁻ (t : PlaneTree), {t : PlaneTree | G t}.indicator (fun (t : PlaneTree) => βˆ‘' (q : β„•), PlaneTree.forestSizeWeight ΞΌ d q * Ξ¦ (1 + t.vertexCount + q)) t βˆ‚PlaneTree.bienayme ΞΌ

Branching property at a marked child (source lines 2236–2240, 2247–2251): for the i-th child among d + 1, βˆ‘_f ∏_j w(f_j) 1{G(f_i)} Ξ¦(1 + βˆ‘_j |f_j|) = E[1{G(𝒯)} βˆ‘_q P(|𝒒_d| = q) Ξ¦(1 + |𝒯| + q)]: the other d children form an independent forest.

The declumped mark and its marked children (source lines 2236–2241, union bound over the marked child): E[1{Ghat(𝒯)} Ξ¦(|𝒯|)] ≀ βˆ‘_{d ∈ B} d ΞΌ_d E[1{G(𝒯)} βˆ‘_q P(|𝒒_{d-1}| = q) Ξ¦(1 + |𝒯| + q)].

∫⁻ 1_s f dν = ν(s) ∫⁻ f dν(· | s) for a finite measure on plane trees.

theorem PoissonDeclumping.Bridge.lintegral_declumped_le_of_bound (ΞΌ : PMF β„•) (B : Set β„•) (G : PlaneTree β†’ Prop) (Ξ¦ : β„• β†’ ENNReal) (A : PlaneTree β†’ ENNReal) (ΞΊ : β„• β†’ ENNReal) (hb : βˆ€ (d : β„•) (t : PlaneTree), 1 ≀ d β†’ G t β†’ βˆ‘' (q : β„•), PlaneTree.forestSizeWeight ΞΌ (d - 1) q * Ξ¦ (1 + t.vertexCount + q) ≀ A t + ΞΊ d) :

Union bound over the marked child, with a pointwise forest bound (source lines 2236–2241, 2247–2270): if for every d β‰₯ 1 and every marked t, βˆ‘_q P(|𝒒_{d-1}| = q) Ξ¦(1 + |t| + q) ≀ A(t) + ΞΊ(d), then E[1{Ghat(𝒯)} Ξ¦(|𝒯|)] ≀ Ο€ (c_B E[A(𝒯) | G(𝒯) = 1] + βˆ‘_{d ∈ B} d ΞΌ_d ΞΊ(d)).

Forest convolution and the min(1, D / a_{M_n}) domination #

theorem PoissonDeclumping.Bridge.tsum_forestSizeWeight_mul_walkProb (ΞΌ : PMF β„•) (β„“ : β„•) (N : β„€) :
βˆ‘' (q : β„•), PlaneTree.forestSizeWeight ΞΌ β„“ q * ENNReal.ofReal (walkProb ΞΌ (N - ↑q) 0) = ENNReal.ofReal (walkProb ΞΌ N (-↑ℓ))

Forest first-passage decomposition, tsum form (source lines 2250–2255, eq:parent_convolution): for every N : β„€, βˆ‘_q P(|𝒒_β„“| = q) Ξ¨_{N-q}(0) = Ξ¨_N(-β„“) (both sides vanish for N < 0).

theorem PoissonDeclumping.Bridge.tendsto_tsum_min_div {ΞΌ : PMF β„•} {Ξ± : ℝ} {a : β„• β†’ ℝ} (hT : StableToolkit ΞΌ Ξ± a) (B : Set β„•) {l : Filter β„•} {m : β„• β†’ β„•} (hm : Filter.Tendsto m l Filter.atTop) :
Filter.Tendsto (fun (n : β„•) => βˆ‘' (d : β„•), B.indicator (fun (d : β„•) => ↑d * LeafHeight.mass ΞΌ d * min 1 (↑d / a (m n))) d) l (nhds 0)

The dominated-convergence step of source lines 2242–2245: βˆ‘_{d ∈ B} d ΞΌ_d min(1, d / a_{m_n}) β†’ 0 whenever m_n β†’ ∞.

theorem PoissonDeclumping.Bridge.min_le_mul_min_one {u₁ uβ‚‚ K y : ℝ} (h1 : u₁ ≀ K) (h2 : uβ‚‚ ≀ K * y) :
min u₁ uβ‚‚ ≀ K * min 1 y

min u₁ uβ‚‚ ≀ K min(1, y) when u₁ ≀ K and uβ‚‚ ≀ K y.

theorem PoissonDeclumping.Bridge.forest_bound_unweighted {ΞΌ : PMF β„•} {Ξ± : ℝ} {a : β„• β†’ ℝ} (hT : StableToolkit ΞΌ Ξ± a) :
βˆƒ (K : ℝ), 1 ≀ K ∧ βˆ€ (M : β„•), 8 ≀ M β†’ βˆ€ (d : β„•) (t : PlaneTree), 1 ≀ d β†’ βˆ‘' (q : β„•), PlaneTree.forestSizeWeight ΞΌ (d - 1) q * {s : β„• | M < s}.indicator 1 (1 + t.vertexCount + q) ≀ {t : PlaneTree | ↑M / 4 < ↑t.vertexCount}.indicator 1 t + ENNReal.ofReal (K * min 1 (↑d / a (M / 2)))

Unweighted forest bound at a marked child (source lines 2236–2245): for M β‰₯ 8, d β‰₯ 1 and every tree t, P(1 + |t| + |𝒒_{d-1}| > M) ≀ 1{|t| > M/4} + K min(1, d / a_{⌊M/2βŒ‹}): if |t| ≀ M/4 the other forest has more than M/2 vertices.

theorem PoissonDeclumping.Bridge.forest_bound_weighted {ΞΌ : PMF β„•} {Ξ± : ℝ} {a : β„• β†’ ℝ} (hT : StableToolkit ΞΌ Ξ± a) :
βˆƒ (K : ℝ), 0 ≀ K ∧ βˆ€ (M n : β„•), 8 ≀ M β†’ 2 * M ≀ n β†’ βˆ€ (d : β„•) (t : PlaneTree), 1 ≀ d β†’ βˆ‘' (q : β„•), PlaneTree.forestSizeWeight ΞΌ (d - 1) q * ENNReal.ofReal ({s : β„• | M < s}.indicator (fun (s : β„•) => walkProb ΞΌ (↑n - ↑s) 0) (1 + t.vertexCount + q)) ≀ ENNReal.ofReal ({t : PlaneTree | ↑M / 4 < ↑t.vertexCount}.indicator (fun (t : PlaneTree) => supWalk ΞΌ (↑n - 1 - ↑t.vertexCount)) t) + ENNReal.ofReal (K / a n * min 1 (↑d / a (M / 2)))

Bridge-weighted forest bound at a marked child (source lines 2247–2264, eq:parent_convolution, eq:parent_forest, numerator form): for M β‰₯ 8, n β‰₯ 2M, d β‰₯ 1 and every tree t, E[Ξ¨_{n-1-|t|-|𝒒_{d-1}|}(0) 1{1 + |t| + |𝒒_{d-1}| > M}] ≀ 1{|t| > M/4} sup_x Ξ¨_{n-1-|t|}(x) + (K / a_n) min(1, d / a_{⌊M/2βŒ‹}). The first term is the convolution bound, the second combines it with the weighted forest tail (the other forest has more than M/2 vertices).

Lemma lem:root_declumping(iii) #

theorem PoissonDeclumping.Bridge.tsum_indicator_mul_ofReal_min {ΞΌ : PMF β„•} (hc : IsCritical ΞΌ) (B : Set β„•) {K x : ℝ} (hK : 0 ≀ K) (hx : 0 < x) :
βˆ‘' (d : β„•), B.indicator (fun (d : β„•) => ↑d * ΞΌ d * ENNReal.ofReal (K * min 1 (↑d / x))) d = ENNReal.ofReal (K * βˆ‘' (d : β„•), B.indicator (fun (d : β„•) => ↑d * LeafHeight.mass ΞΌ d * min 1 (↑d / x)) d)

βˆ‘_{d ∈ B} d ΞΌ_d K min(1, d / x) = K βˆ‘_{d ∈ B} d ΞΌ_d min(1, d / x) (from ℝβ‰₯0∞ to ℝ).

theorem PoissonDeclumping.Bridge.toReal_le_mul_add {Y P C X : ENNReal} {y : ℝ} (hP : P β‰  ⊀) (hC : C β‰  ⊀) (hX : X β‰  ⊀) (hy : 0 ≀ y) (h : Y ≀ P * (C * X + ENNReal.ofReal y)) :

From Y ≀ P (C X + y) in ℝβ‰₯0∞ to ℝ.

theorem PoissonDeclumping.Bridge.indicator_indicator_vertexCount {Ξ² : Type u_1} [Zero Ξ²] (S : Set PlaneTree) (P : β„• β†’ Prop) (f : β„• β†’ Ξ²) (T : PlaneTree) :
S.indicator (fun (T : PlaneTree) => {s : β„• | P s}.indicator f T.vertexCount) T = {T : PlaneTree | T ∈ S ∧ P T.vertexCount}.indicator (fun (T : PlaneTree) => f T.vertexCount) T
theorem PoissonDeclumping.Bridge.tsum_indicator_min_nonneg {ΞΌ : PMF β„•} (B : Set β„•) (x : ℝ) (hx : 0 ≀ x) :
0 ≀ βˆ‘' (d : β„•), B.indicator (fun (d : β„•) => ↑d * LeafHeight.mass ΞΌ d * min 1 (↑d / x)) d
theorem PoissonDeclumping.Bridge.microscopic_declumped {ΞΌ : PMF β„•} {Ξ± : ℝ} {a : β„• β†’ ℝ} (hT : StableToolkit ΞΌ Ξ± a) {G : β„• β†’ PlaneTree β†’ Prop} {k M : β„• β†’ β„•} {B : Set β„•} (hS : StructuralMicroscopic ΞΌ a G k M) (hM : IsCutoff M) (hB0 : 0 < cB ΞΌ B) (hB1 : cB ΞΌ B < ⊀) :
Microscopic ΞΌ (fun (j : β„•) => PlaneTree.declumpedMark B (G j)) k M

Lemma lem:root_declumping(iii) (source lines 2212–2214, 2231–2298, P16): if the raw marks (G_k) satisfy Assumption (S) along (k_n) with respect to the cutoff (M_n) and c_B ∈ (0, ∞), then the declumped marks Ghat_k = declumpedMark B G_k satisfy Assumption (M) along (k_n) with respect to the same cutoff. Periodic offspring laws are included: the paper's "replace n by n - 1" (line 2286) is replaced by the numerator form of Lemma lem:consequenceM with time shift j = 1 (consequenceM_weighted_sup), so no Ξ¨_{n-1}(-1) appears (correction F3).