Documentation

BGW.StableTrees.StablePassageSubordinator

Passage subordinator and local time at the minimum #

The strict downward passage times of the spectrally positive stable process form a càdlàg subordinator. Its right-continuous inverse is the negative running infimum, whose Stieltjes measure is carried by the zero set of the reflected path.

Reference: J. Bertoin, Lévy Processes, Cambridge University Press, 1996, Chapter VII, Section 1 and Chapter IV, for the local time at the infimum and its inverse.

The ℝ≥0-valued passage time (junk value 0 on {T_x = ⊤}, a null set).

Equations
Instances For

    First hitting time of (-∞, -x], H_x = inf{t | X_t ≤ -x} (weak hitting, NOT strict).

    Equations
    Instances For

      The law of the passage time T_x under the local stable law, as a measure on ℝ (finite reading passageTime).

      Equations
      Instances For

        Joint measurability of (x, f) ↦ T_x(f) (finite reading).

        passageLaw θ x is a probability measure (x = 0 uses T_0 = 0).

        Pathwise structure: finiteness, monotonicity, strict right-continuity, T_0 = 0 #

        Almost surely, for every x ≥ 0, T_x is finite and equals its finite reading, and T_0 = 0.

        theorem BGW.StableTrees.ae_isCadlag_monotone_passage {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) :
        ∀ᵐ (f : Skorokhod.DLoc) ∂↑(stableLocalLaw θ), (IsCadlag fun (x : NNReal) => ↑(passageTime x f)) ∧ (Monotone fun (x : NNReal) => passageTime x f) ∧ passageTime 0 f = 0

        Monotone, càdlàg, starts at 0 Right-continuity here is the one for the strict passage time T_x = inf{t | X_t < -x}.

        Almost surely T is unbounded: T_x → ∞ as x → ∞.

        Inverse relation with the running infimum (L = -runningInf is the local time) #

        T is the strict right inverse of L = -runningInf: T_x = inf{t | L_t > x}. This is a deterministic identity, valid for every path (bounded-below on compacts), every x ≥ 0, and the inf over the empty set is ⊤.

        Event sandwich for the running infimum, strict side. For every path, level and time: {T_x < t} ⊆ {I_t < -x} ⊆ {T_x ≤ t}.

        Weak side of the sandwich. For every path, H_x ≤ t → I_t ≤ -x; also T_x ≤ t → I_t ≤ -x (deterministic; stated a.s. for uniformity with the other items). The converses are NOT claimed: at a plateau level x of L (an excursion from the minimum straddling t) one has I_t = -x but T_x > t.

        Strict passage versus first hitting (deterministic sandwich). For every path and level x > 0, sup_{y<x} T_y ≤ H_x ≤ T_x. The left limit T_{x-} of the passage process is at most the first hitting time of -x; a plateau of -runningInf at level x is exactly where these are <.

        Local time L = -I: measurability, adaptedness, Stieltjes concentration, straddling #

        noncomputable def BGW.StableTrees.localTime (f : Skorokhod.DLoc) (t : NNReal) :

        The local time at the minimum, L_t = -I_t.

        Equations
        Instances For

          L is jointly measurable in (t, f).

          L is adapted to the right-continuous filtration of DLoc.

          Almost surely L_0 = 0 (from f 0 = 0), L is monotone, continuous, and unbounded.

          The Stieltjes function s ↦ L_{s⁺} on ℝ (s⁺ = max s 0) when L is monotone and continuous; junk value the zero function otherwise (a null set of paths).

          Equations
          • One or more equations did not get rendered due to their size.
          Instances For

            The ℝ≥0-valued left limit T_{x-} = sup_{y<x} T_y of the passage process, with T_{0-} = 0; junk 0 where some T_y = ⊤ (null set).

            Equations
            Instances For

              L_t as a level in ℝ≥0 (L_t ≥ 0 almost surely).

              Equations
              Instances For

                Joint measurability of (x, f) ↦ T_{x-}(f).

                First hitting agrees with strict passage at each FIXED level (a.s., NOT simultaneously for all levels).

                Straddling. Almost surely for all t: T_{L_t -} ≤ t ≤ T_{L_t}.

                theorem BGW.StableTrees.ae_localTime_passageTime {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) :

                Generalized inverse identities. Almost surely for all x ≥ 0: L (T_x) = x and L (T_{x-}) = x.

                Strong Markov at T_x and independent stationary increments #

                Additivity from creeping (deterministic on the creeping event): a.s., for all x y ≥ 0, T_{x+y}(f) = T_x(f) + T_y(shiftAt T_x f).

                Reflected process after T_x: on the creeping event, X_{T_x + t} - I_{T_x + t} = X'_t - I'_t for the centered shift X' = shiftAt T_x.

                theorem BGW.StableTrees.measure_inter_passageIncrement {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) (x y : NNReal) {A : Set Skorokhod.DLoc} (hA : MeasurableSet A) {S : Set ℝ} (hS : MeasurableSet S) :
                ↑(stableLocalLaw θ) (A ∩ {f : Skorokhod.DLoc | ↑(passageTime (x + y) f) - ↑(passageTime x f) ∈ S}) = ↑(stableLocalLaw θ) A * (passageLaw θ y) S

                Fresh increment: independence of F_{T_x} and law. For x, y ≥ 0, any event A of the stopped σ-algebra F_{T_x} and Borel S ⊆ ℝ: P(A ∩ {T_{x+y} - T_x ∈ S}) = P(A) · P(T_y ∈ S).

                theorem BGW.StableTrees.map_passageIncrement {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) (x y : NNReal) :

                Stationary increments. T_{x+y} - T_x has the law of T_y.

                One-dimensional law: Laplace transform, positivity, scaling #

                theorem BGW.StableTrees.integral_exp_neg_passageLaw {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) (x : NNReal) {q : ℝ} (hq : 0 ≤ q) :
                ∫ (t : ℝ), Real.exp (-q * t) ∂passageLaw θ x = Real.exp (-↑x * q ^ (1 / θ))

                Laplace transform: ∫ e^{-q t} d passageLaw θ x = e^{-x q^{1/θ}} for x ≥ 0, q ≥ 0.

                theorem BGW.StableTrees.passageLaw_Iic_zero {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) {x : NNReal} (hx : 0 < x) :
                (passageLaw θ x) (Set.Iic 0) = 0

                For x > 0 the law of T_x lives on (0, ∞).

                theorem BGW.StableTrees.passageLaw_scaling {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) (x : NNReal) {c : NNReal} (hc : 0 < c) :
                passageLaw θ (c * x) = MeasureTheory.Measure.map (fun (t : ℝ) => ↑c ^ θ * t) (passageLaw θ x)

                Scaling. T_{c x} =d c^θ T_x, for c > 0.

                theorem BGW.StableTrees.eq_passageLaw_of_integral_exp_neg {θ : ℝ} (hθ : 1 < θ ∧ θ ≤ 2) (x : NNReal) {ρ : MeasureTheory.Measure ℝ} [MeasureTheory.IsProbabilityMeasure ρ] (h0 : ρ (Set.Iio 0) = 0) (hρ : ∀ (q : ℝ), 0 < q → ∫ (t : ℝ), Real.exp (-q * t) ∂ρ = Real.exp (-↑x * q ^ (1 / θ))) :
                ρ = passageLaw θ x

                Identification of the law by its Laplace transform: any probability law on ℝ with the transform e^{-x q^{1/θ}} is passageLaw θ x (the positive (1/θ)-stable law).