English

Numerical Computation of First-Passage Times of Increasing Levy Processes

Probability 2009-04-28 v1

Abstract

Let {D(s),s0}\{D(s), s \geq 0\} be a non-decreasing L\'evy process. The first-hitting time process {E(t)t0}\{E(t) t \geq 0\} (which is sometimes referred to as an inverse subordinator) defined by E(t)=inf{s:D(s)>t}E(t) = \inf \{s: D(s) > t \} is a process which has arisen in many applications. Of particular interest is the mean first-hitting time U(t)=EE(t)U(t)=\mathbb{E}E(t). This function characterizes all finite-dimensional distributions of the process EE. The function UU can be calculated by inverting the Laplace transform of the function U~(λ)=(λϕ(λ))1\widetilde{U}(\lambda) = (\lambda \phi(\lambda))^{-1}, where ϕ\phi is the L\'evy exponent of the subordinator DD. In this paper, we give two methods for computing numerically the inverse of this Laplace transform. The first is based on the Bromwich integral and the second is based on the Post-Widder inversion formula. The software written to support this work is available from the authors and we illustrate its use at the end of the paper.

Keywords

Cite

@article{arxiv.0904.4232,
  title  = {Numerical Computation of First-Passage Times of Increasing Levy Processes},
  author = {Mark S. Veillette and Murad S. Taqqu},
  journal= {arXiv preprint arXiv:0904.4232},
  year   = {2009}
}

Comments

31 Pages, 7 sections, 11 figures, 2 tables

R2 v1 2026-06-21T12:55:32.090Z