ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS
|
|
- Dominic Brooks
- 5 years ago
- Views:
Transcription
1 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS MASAHIKO EGAMI AND KAZUTOSHI YAMAZAKI ABSTRACT. We study the scale function of the spectrally negative phase-type Lévy process. Its scale function admits an analytical expression and so do a number of its fluctuation identities. Motivated by the fact that the class of phase-type distributions is dense in the class of all positive-valued distributions, we propose a new approach to approximating the scale function and the associated fluctuation identities for a general spectrally negative Lévy process. Numerical examples are provided to illustrate the effectiveness of the approximation method. The extension to the case with jumps of infinite activity is also discussed. Key words: phase-type models; spectrally negative Lévy processes; scale functions; overshoot and undershoot Mathematics Subject Classification (2) : Primary: 6G51 Secondary: 6J75 1. INTRODUCTION In the last decade, significant progress has been made regarding spectrally negative Lévy processes, and this is mainly due to the scale function. As can be seen in the work of, for example, [7, 25], a number of fluctuation identities concerning spectrally negative Lévy processes can be expressed in terms of scale functions. There are numerous applications in applied probability including optimal stopping, ueuing, branching processes, insurance and credit risk. Despite these advances, a major obstacle still remains in putting these in practice because scale functions are in general known only up to their Laplace transforms, and only a few cases admit explicit expressions. Typically, one needs to rely on numerical Laplace inversion in order to approximate the scale function; see [24, 38]. In this paper, we propose a phase-type (PH)-fitting approach by using the scale function for the class of spectrally negative PH Lévy processes, or Lévy processes with negative PH-distributed jumps. Consider a continuous-time Markov chain with some initial distribution and state space consisting of a single absorbing state and a finite number of transient states. The PH distribution is the distribution of the time to absorption. The class of PH distributions includes, for example, the exponential, hyperexponential, Erlang, hyper-erlang and Coxian distributions; see, e.g. Section 3 of [2]. It is known that the class of PH distributions is dense in the class of all positive-valued distributions, and conseuently the scale function of any spectrally negative Lévy process can be approximated by those of PH Lévy processes. First draft: May 1, 21; this version: July 3, 211. M. Egami is in part supported by Grant-in-Aid for Scientific Research (B) No and (C) No , Japan Society for the Promotion of Science. K. Yamazaki is in part supported by Grant-in-Aid for Young Scientists (B) No , the Ministry of Education, Culture, Sports, Science and Technology, and by Grant-in-Aid for Scientific Research (B) No , Japan Society for the Promotion of Science. 1
2 2 M. EGAMI AND K. YAMAZAKI Our objective is to obtain the scale function of the spectrally negative PH Lévy process, and use it to approximate for a general case. We achieve its closed-form expression by taking advantage of the structure of its Wiener-Hopf factor as obtained in [3]. In order to verify that this is indeed a useful alternative to the Laplace inversion approach, we show convergence results via the continuity theorem and then conduct a series of numerical experiments. One major advantage of this approach is that the approximated scale function is given as a function in a closed form, which enables one to analytically differentiate/integrate to obtain other fluctuation identities explicitly. The derivative of the scale function is used to identify, for example, the fluctuation of its reflected process and is commonly applied in insurance risk literature [5, 27, 3]. The integration of the scale function with respect to the Lévy measure is reuired when overshoots/undershoots are involved. On the other hand, these calculations must be conducted numerically when the scale function is approximated via numerical Laplace inversion; because the scale function increases exponentially and its derivative tends to explode near zero, minimizing numerical errors tends to be a difficult task. Another advantage of the PH-fitting approach is that the Laplace transform of the PH distribution has an explicit expression. The Laplace inversion approach inverts the euality written in terms of the Laplace exponents which admit analytical expressions only in special cases (see (2.1)-(2.2) below). In other words, it contains two types of errors: 1) the approximation error caused while computing the Laplace exponent and 2) the error caused while inverting the Laplace transform. On the other hand, the PH-fitting approach only contains the PH-fitting error thanks to the closed-form Laplace transform of the PH distribution. The PH-fitting approach is particularly efficient when the Lévy density is completely monotone because the fitting can be conducted by hyperexponential distributions (i.e. hyperexponential-fitting) and there exist algorithms that are guaranteed to converge (see, e.g., [15]). The class of Lévy processes with completely monotone Lévy density is rich. It includes compound-poisson processes with long-tailed distributed jumps such as the Pareto, Weibull and Gamma distributions and other processes such as variance gamma [35, 36], CGMY [11], generalized hyperbolic [14] and normal inverse Gaussian [6] processes. The completely monotone assumption of the Lévy density is sometimes reuired, for example, in the optimal dividend problem [33, 34]. In order to give examples of how hyperexponential-fitting can be applied to compute fluctuation identities, we obtain the closed-form expressions of the overshoot/undershoot distributions at the first down-crossing time for the hyperexponential case, and use it to approximate those for the processes with Weibull/Pareto-type jumps. The Lévy measure of the PH Lévy process is only limited to a finite measure. To complement this, we also consider the meromorphic Lévy process [22], whose Lévy measure may be infinite and is given as a countable sum of exponential functions. This generalizes the hyperexponential Lévy process as well as a number of Lévy processes such as Lamperti-stable processes [1, 13], hypergeometric processes [23, 26] and those in the β- and Θ- families [2, 21]. The corresponding scale function can be expressed explicitly as a countable sum of exponential functions, which can be effectively approximated by a finite sum with some analytical error bounds. We show numerically the effectiveness of the approximation procedure using examples with β-processes [2]. The rest of the paper is organized as follows. Section 2 studies the spectrally negative PH Lévy process and its scale function. We show via the continuity theorem that the PH-fitting approach can approximate the scale function for a general spectrally negative Lévy process. We further obtain, for the hyperexponential case, the joint distribution of the overshoot and undershoot at the first down-crossing time. Section 3 obtains the scale function
3 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 3 for the meromorphic Lévy process and its upper and lower bounds. Section 4 verifies the effectiveness of the PH-fitting approach through a series of numerical results. All proofs are given in the appendix. 2. SPECTRALLY NEGATIVE PHASE-TYPE LÉVY PROCESSES 2.1. Scale functions. Let (Ω, F, P) be a probability space hosting a spectrally negative Lévy process X = {X t ; t }, P x be the conditional probability under which X = x, and F := {F t : t } be the filtration generated by X. The process X is uniuely characterized by its Laplace exponent (2.1) ψ(s) := log E [ e sx 1 ] = ˆµs σ2 s 2 + (e sz 1 sz1 {z>1} )Π(dz), s C where Π is a Lévy measure with the support (, ) and satisfies the integrability condition (,) (1 z2 )Π(dz) <. It has paths of bounded variation if and only if σ = and (1 z )Π(dz) < ; (,) see, for example, [25], Lemma In this case, we can rewrite the Laplace exponent (2.1) by with ψ(s) = µs + (e sz 1)Π(dz), s C, µ := ˆµ zπ(dz). 1 Fix and any spectrally negative Lévy process with its Laplace exponent ψ. The scale function W () : R R is a function whose Laplace transform is given by (2.2) where e sx W () (x)dx = 1 ψ(s), s > ζ (2.3) ζ := sup{s : ψ(s) = },. We assume W () (x) = on (, ). Let us define the first down- and up-crossing times, respectively, by (2.4) τ a := inf {t : X t < a} and τ + b := inf {t : X t > b}, a, b R. Then we have for any x < b (2.5) E x [ e τ + b 1{τ + b <τ, τ + b < } ] = W () (x) W () (b) where (2.6) Z () (x) := 1 + and E x [ e τ 1{τ + b >τ, τ < } ] = Z () (x) Z () (b) W () (x) W () (b) x W () (y)dy, x R. Here, we disregard the case when X is a negative subordinator (or decreasing a.s.).
4 4 M. EGAMI AND K. YAMAZAKI We also consider a version of the scale function W ζ : R R that satisfies (2.7) with its Laplace transform W () (x) = e ζx W ζ (x), x R (2.8) e sx W ζ (x)dx = 1 ψ(s + ζ ), s >. Suppose P c, for any given c >, is the probability measure defined by the Esscher transform (2.9) dp c dp = e cxtψ(c)t, t ; Ft see page 78 of [25]. Then W ζ under P ζ is euivalent to W () under P. Furthermore, it is known that W ζ is monotonically increasing and (2.1) W ζ (x) (ψ (ζ )) 1 as x, which also implies that the scale function W () increases exponentially in x; (2.11) W () (x) eζx ψ (ζ ) as x, except for the case = and ψ (ζ ) =. Regarding the smoothness of the scale function, if its jump distribution has no atoms, then W () C 1 (, ); if it has a Gaussian component (σ > ), then W () C 2 (, ); see [12]. In particular, a stronger result holds for the completely monotone jump case. Recall that a density function f is called completely monotone if all the derivatives exist and, for every n 1, where f (n) denotes the n th derivative of f. (1) n f (n) (x), x, Lemma 2.1 (Loeffen [33]). If the Lévy measure has a completely monotone density, W ζ is also completely monotone. Finally, the behavior in the neighborhood of zero is given as follows; see Lemmas of [29]. Lemma 2.2. For every, we have { W (), unbounded variation () = 1 µ, bounded variation } and W () (+) = 2, σ 2 σ >, σ = and Π(, ) = +Π(,), µ 2 compound Poisson.
5 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS Spectrally negative PH Lévy processes. Consider a continuous-time Markov chain Y = {Y t ; t } with finite state space {1,..., m} { } where 1,..., m are transient and is absorbing. Its initial distribution is given by a simplex α = [α 1,..., α m ] such that α i = P {Y = i} for every i = 1,..., m. The intensity matrix Q is partitioned into the m transient states and the absorbing state, and is given by [ ] T t Q :=. Here T is an m m-matrix called the PH-generator, and t = T 1 where 1 = [1,..., 1]. A distribution is called phase-type (PH) with representation (m, α, T ) if it is the distribution of the absorption time to in the Markov chain described above. It is known that T is non-singular and thus invertible; see [1]. Its distribution and density functions are given, respectively, by F (z) = 1 αe T z 1 and f(z) = αe T z t, z. Let X = {X t ; t } be a spectrally negative Lévy process of the form N t (2.12) X t X = µt + σb t Z n, t <, for some µ R and σ. Here B = {B t ; t } is a standard Brownian motion, N = {N t ; t } is a Poisson process with arrival rate λ, and Z = {Z n ; n = 1, 2,...} is an i.i.d. seuence of PH-distributed random variables with representation (m, α, T ). These processes are assumed mutually independent. Its Laplace exponent is then n=1 (2.13) ψ(s) = µs σ2 s 2 + λ ( α(si T ) 1 t 1 ), which is analytic for every s C except for the eigenvalues of T. Disregarding the negative subordinator case, we consider the following two cases: Case 1: when σ > (i.e. X is of unbounded variation), Case 2: when σ = and µ > (i.e. X is a compound Poisson process). Here, in Case 2, down-crossing of a threshold can occur only by jumps; see, for example, Chapter III of [7]. On the other hand, in Case 1, down-crossing can occur also by creeping downward (by the diffusion components). Due to this difference, the form of the scale function differs as we shall see. Fix >. Consider the Cramér-Lundberg euation (2.14) ψ(s) =, and define the set of (the absolute values of) negative roots and the set of poles: I := {i : ψ(ξ i, ) = and R(ξ i, ) > }, { } J := j : ψ(η j ) = and R(η j) >. The elements in I and J may not be distinct, and, in this case, we take each as many times as its multiplicity. By Lemma 1 of [3], we have { J + 1, for Case 1, I = J, for Case 2.
6 6 M. EGAMI AND K. YAMAZAKI In particular, if the representation is minimal (see [3]), we have J = m. Let e be an independent exponential random variable with parameter and denote the running maximum and minimum, respectively, by X t = sup X s and X t = inf X s, t. s t s t The Wiener-Hopf factorization states that /( ψ(s)) = ϕ + (s)ϕ (s) for every s C such that R(s) =, with the Wiener-Hopf factors [ ] (2.15) ϕ (s) := E exp(sx e ) and ϕ + (s) := E [ exp(sx e ) ] that are analytic for s with R(s) > and R(s) <, respectively. By Lemma 1 of [3], we have, for every s such that R(s) >, (2.16) ϕ j J (s) = (s + η j ) i I ξ i, j J η j i I (s + ξ i, ), from which we can obtain the distribution of X e by the Laplace inverse via partial fraction expansion. As in Remark 4 of [3], let n denote the number of different roots in I and m i denote the multiplicity of a root ξ i, for i = 1,..., n. Then we have (2.17) where { } P X e dx = n m i k=1 A (k) i, := 1 (m i k)! A (k) i, ξ (ξ i, x) k1 i, (k 1)! eξ i,x dx, x > m ik s m ik ϕ (s)(s + ξ i, ) m i. s=ξi, Notice that this can be simplified significantly when all the roots in I are distinct Scale functions for spectrally negative PH Lévy processes. We obtain the scale function focusing on the case >. The scale function when = and X drifts to + can be derived by using W () (x) = P x {X } /ψ () and the ruin probability (19) of [3] by taking ; Kyprianou and Palmowski [27] briefly stated the scale function when =, X drifts to + and all the roots in I are distinct. The case = and X drifts to can be obtained indirectly by change of measure (2.9); in this case, W () (x) = e ζ x W ζ (x) with ζ > and X drifts to + under P ζ. Before obtaining the scale function, we shall first represent the positive root ζ (2.3) in terms of the negative roots {ξ i, ; i I }. Let us define (2.18) and by Lemma 2.2 (2.19) ϱ := n θ := ζ W () () + W () (+) = ξ k i, A (1) i, ξ i,, >, { 2, σ 2 for Case 1 ζ µ + +λ, µ 2 for Case 2 }.
7 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 7 Lemma 2.3. For every >, we have ζ = θ ϱ. We now obtain the version of the scale function W ζ ( ). In the lemma below, W ζ () = W () () is either or 1 µ depending on if it is Case 1 or 2; see Lemma 2.2. Lemma 2.4. For every >, we have W ζ (x) W ζ () = ζ n m i k=1 ( ) k A (k) ξi, i, ζ + ξ i, k1 1 e (ζ+ξ ((ζ i,)x + ξ i, )x) j, x. j! Lemmas 2.2 and 2.4 and (2.7) show the following. By Lemma 2.3, these expressions can also be rewritten in a different way. Proposition 2.1. Fix > and x. (1) For Case 1, W () (x) = ζ n m i k=1 ( ) k A (k) ξi, i, ζ + ξ i, j= k1 e ζx e ξ ((ζ i,x + ξ i, )x) j. j! j= (2) For Case 2, W () (x) = ζ n m i k=1 ( ) k A (k) ξi, i, ζ + ξ i, k1 e ζx e ξ ((ζ i,x + ξ i, )x) j + 1 j! µ eζx. j= The scale functions obtained above are infinitely differentiable. In particular, the first derivative becomes W () (x) = ζ n m i ( ) k A (k) ξi, i, ζ e ζx + ξ i, e ξ i,x ((ζ + ξ i, )x) k1 k2 ζ e ξ i,x ζ + ξ i, (k 1)! j! for Case 1 and k=1 j= ((ζ + ξ i, )x) j W () (x) = ζ n m i ( A (k) ξi, i, k=1 ζ + ξ i, ) k ζ e ζx + ξ i, e ξ i,x ((ζ + ξ i, )x) k1 (k 1)! k2 ζ e ξ i,x j= ((ζ + ξ i, )x) j + 1 j! µ ζ e ζx for Case 2. When all the roots in I are distinct, the scale functions above can be simplified and have nice properties as discussed in the following corollary. Corollary 2.1. Suppose all the roots in I are distinct.
8 8 M. EGAMI AND K. YAMAZAKI (1) The scale function can be simplified to W () (x) = i I C i, [e ζx e ξ i,x ], (2.2) W () (x) = i I C i, [e ζx e ξ i,x ] + 1 µ eζx, for Case 1 and 2, respectively where (2) W () is convex. (3) W ζ is completely monotone. C i, := ζ In the lemma above, by (2.1)-(2.11), we must have ξ i, A (1) i, ζ + ξ i,, i I. (2.21) 1 ψ (ζ ) = i I C i, and 1 ψ (ζ ) = i I C i, + 1 µ, respectively for Case 1 and 2. Example 2.1 (Hyperexponential Case). As an important example where all the roots in I are distinct, we consider the case where Z has a hyperexponential distribution with density function f(z) = p j η j e ηjz, z, for some < η 1 < < η m <. Its Laplace exponent (2.1) is then (2.22) ψ(s) = µs σ2 s 2 λ p j s η j + s. Notice in this case that η 1,..., η m are the poles of the Laplace exponent. Furthermore, all the roots in I are distinct and satisfy the following interlacing condition for every > : (1) when σ >, there are m + 1 roots ξ 1,,..., ξ m+1, such that < ξ 1, < η 1 < ξ 2, < < η m < ξ m+1, < ; (2) when σ = and µ >, there are m roots ξ 1,,..., ξ m, such that < ξ 1, < η 1 < ξ 2, < < ξ m, < η m <. The class of hyperexponential distributions is important as it is dense in the class of all positive-valued distributions with completely monotone densities.
9 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS Approximation of the scale function of a general spectrally negative Lévy process. The scale function obtained in Proposition 2.1 can be used to approximate the scale function of a general spectrally negative Lévy process. By Proposition 1 of [3], there exists, for any spectrally negative Lévy process X, a seuence of spectrally negative PH Lévy processes X (n) converging to X in D[, ). This is euivalent to saying that X (n) 1 X 1 in distribution by [18], Corollary VII 3.6; see also [37]. Suppose ψ n (resp. ψ), ζ,n (ζ ) and W n () /W ζ,n (W () /W ζ ) are the Laplace exponent, the positive root (2.3) and the scale function of X (n) (X), respectively. Because W () is an increasing function, the measure W () (dx) associated with the distribution of W () (, x] is well-defined and we obtain as in page 218 of [25], (2.23) e βx W () β (dx) = ψ(β). [, ) Because these processes are spectrally negative and ψ is continuous on [, ), we have, by the continuity theorem, ψ n (β) ψ(β) for every β >. Now in view of (2.23), the convergence of the scale function holds by the continuity of the scale function and the continuity theorem; see [16], Theorem 2a, XIII.1. In other words, W () n (x) W () (x) and W ζ,n(x) W ζ (x) as n for every x. The smoothness and monotonicity properties of the scale function can be additionally used to obtain stronger results. The scale functions in Proposition 2.1 are in C (, ). In addition, when all the roots of I are different, its first derivative W ζ is completely monotone as discussed in Corollary 2.1. Suppose that it is in C 2 (, ) (which holds, for example, when σ > by [12]), W ζ (+) < (i.e., σ > or Π(, ) < ) and W ζ (x) for every x. Because W ζ (x) x, we see that F (x) := (W ζ (+)) 1 x W ζ (y) dy is a probability distribution and e βx F (dx) = (W ζ (+)) [ 1 β 2 ] ψ(β + ζ ) + βw ζ () + W ζ (+), β >. Therefore, noting that W () (x) = ζ W () (x) + e ζx W ζ (x) and assuming that the convergent seuence W ζ,n(x) has the same property, we can obtain by the continuity theorem W ζ n,n (x) W ζ (x) and W () n (x) n W () (x), x. The key assumption of the negativity of W ζ holds, for example, for the completely monotone jump case because W ζ is completely monotone by [33]. We can also choose the seuence W ζ,n completely monotone in view of Corollary 2.1 because approximation can be done via hyperexponential distributions. In fact, it also means that W ζ is C (, ) and the convergence of higher derivatives can be pursued. Even for a general jump distribution, the negativity of W ζ is a reasonable assumption in view of the numerical plots given by [38] Obtaining overshoot and undershoot distributions. We conclude this section by giving an example how we can apply the explicit expression of the scale function obtained above. We consider the joint distribution of
10 1 M. EGAMI AND K. YAMAZAKI overshoot and undershoot (with discounting): (2.24) [ ] h (x; A, B) := E x e τ 1{Xτ B, X τ A, τ < }, B B(, ), A B(, ), focusing on the hyperexponential case addressed in Example 2.1. Our objective is to show an example where the PH (hyperexponential)-fitting can be applied; we focus on this rather simple example because we later evaluate it as an approximation tool by comparing the result with the simulated expectation. The results obtained here can be easily extended to more general fluctuation identities; see e.g. [24]. We assume that X is a hyperexponential Lévy process and >. The joint distribution (2.24) can be easily obtained via the resolvent measure and the compensation formula. Lemma 2.5. For all B B(, ) and A B(, ), h (x; A, B) = { Π(du) W () (x) B (A+u) where Π(du) = λ m p jη j e η ju du for all u (, ). } e ζy dy dyw () (x y), B (A+u) After integrating with respect to Π and assume A and B are open intervals, we obtain the following. Proposition 2.2. Suppose B = (b, b) and A = (a, a) for some a a and b b. Then where, for each 1 j m, h (x; A, B) = λ p j (e ηja e ηja )κ j, (x; B) e ζx κ j, (x; B) := (e ) (η j+ζ )(b x) ψ e (η j+ζ )(b x) (ζ )(η j + ζ ) + [ C i, e ξ 1 ( ) i,x e (η jξ i, )(b x) e (η jξ i, )(b x) 1 ( )] e (η j+ζ )b e (η j+ζ )b. (η j ξ i, ) (η j + ζ ) i I By differentiating the above, we obtain the density: E x [ e τ 1{Xτ da, X τ B, τ < } ] = λ E x [ e τ 1{Xτ db, X τ A, τ < } ] = λ p j η j e ηja κ j, (x; B) p j (e ηja e ηja ) { i I C i, e ξ i,x ( e (η jξ i, )b e (η j+ζ ) )b, 1 ψ (ζ ) eζx(η j+ζ )b i I C i, e (ξ i,x+(η j +ζ )b), b < x b x }.
11 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 11 As a special case, the overshoot (resp. undershoot) density becomes, by setting B = (, ) (A = (, )), (2.25) E x [ e τ 1{Xτ da, τ < } ] = λ p j η j e ηja κ j, (x; (, )), [ ] { λ m E x e τ 1{Xτ db, τ < } = p j i I C i, e ξ i,x ( e (η jξ i, )b e (η j+ζ ) )b, b < x, λ [ m p 1 j ψ (ζ ) eζx(η j+ζ )b ] i I C i, e (ξ i,x+(η j +ζ )b), b x, where κ j, (x; (, )) = 1 1 ψ (ζ ) (η j + ζ ) eη jx + [ C i, i I 1 (η j ξ i, ) Remark 2.1. Based on our numerical experiments, it can be conjectured that 1 1 ψ (ζ ) (η j + ζ ) = 1 C i,, j = 1,..., m, (η j ξ i, ) i I ( ) ] e ξi,x e η jx e ξ 1 i,x. (η j + ζ ) and conseuently κ j, (x; (, )) = ( ) 1 C i, (η j ξ i, ) 1 e ξi,x. (η j + ζ ) i I 3. SCALE FUNCTIONS FOR MEROMORPHIC LÉVY PROCESSES In this section, we consider another class of spectrally negative Lévy processes called meromorphic Lévy processes. We obtain their scale functions and use these as approximation tools for a general spectrally negative Lévy process with a completely monotone Lévy measure. Similarly to the approach applied in the last section, we obtain the scale function using its Wiener-Hopf factorization. It has a form expressed as an countable sum of exponential functions which can be bounded efficiently by finite sums Meromorphic Lévy processes. The following is due to [22], Definition 1. Definition 3.1 (spectrally negative meromorphic Lévy process). A spectrally negative Lévy process X is said to be meromorphic if the following conditions hold. (1) The Laplace exponent ψ(s) (2.1) has a countable set of real negative poles. (2) For every, the Cramér-Lundberg euation (2.14) has a countable set of real negative roots. (3) Let {η k ; k 1} and {ξ k, ; k 1}, respectively, be the sets of the absolute values of the poles and the negative roots of (2.14) for fixed. Then it satisfies the following interlacing conditions: < η k < ξ k, < < η 2 < ξ 2, < η 1 < ξ 1, <. (4) There exists α > 1 2 such that η k ck α as k. (5) The Wiener-Hopf factor (2.15) is expressed as convergent infinite products (3.1) ϕ (s + η k ) ξ k, (s) = η k (s + ξ k, ). k=1
12 12 M. EGAMI AND K. YAMAZAKI The class of Meromorphic Lévy processes complements the class of PH Lévy processes described in the previous section because it also contains those with jumps of infinite activity. As noted by Corollary 3 of [22], the property (3) in Definition 4.1 is euivalent to the condition that the Lévy density has the form π(z) = α j η j e ηj z 1 {z<}, z R. This can be seen as an extension to the hyperexponential case as described in Example 2.1. For more details, see [22]. The Wiener-Hopf factor (3.1) for this process is again a rational function as in (2.16) for the PH case. Therefore, this can be inverted again by partial fraction decomposition, and we have as in Corollary 1 of [22] (3.2) where A i, := s + ξ i, ξ i, { } P X e dx = A i, ξ i, e ξi,x dx, x > ϕ (s) = s=ξi, ( 1 ξ ) i, η i j i 1 ξ i, η j 1 ξ i, ξ j,, i 1. Notice by the interlacing condition that A i, > for every i Scale functions for meromorphic Lévy processes. We now obtain the scale function for the meromorphic Lévy process. We omit the proof because it is similar to the PH case; see Appendix A.1. Lemma 3.1. For every >, we have where W ζ (x) W ζ () = C i, [1 ] e (ζ+ξ i,)x, x (3.3) C i, := ζ ξ i, A i, ζ + ξ i,, i 1. By (2.1) and Lemma 3.1, we have, by taking the limit, (3.4) γ := C i, = (ψ (ζ )) 1 W ζ () <. The scale function can be therefore obtained by Lemma 3.1 and (3.4). Proposition 3.1. For every >, we have (3.5) W () (x) = C i, [e ] ζx e ξ i,x + W ζ ()e ζx = (ψ (ζ )) 1 e ζx C i, e ξi,x, x.
13 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 13 By the dominated convergence theorem thanks to the non-negativity of C i, and convexity of the exponential function, we have, for every > and x >, W () (x) = (ψ (ζ )) 1 ζ e ζx + C i, ξ i, e ξi,x. By adjusting slightly the proof of Lemma 2.3 (given in Appendix A.1), we have the following. Lemma 3.2. (1) The following two statements are euivalent: (a) σ = and Π(, ) =, (b) A i,ξ i, =. (2) Suppose σ > or Π(, ) <. Then, for every >, we have ( ) 1 ζ = θ (3.6) A i, ξ i, (3.7) where θ := ζ W () () + W () (+) = { 2, σ 2 when σ > ζ µ + +Π(,), µ 2 when σ = 3.3. Approximation of the scale functions via finite sum. The scale function obtained in Proposition 3.1 is an countable sum of exponential functions and in reality its exact value cannot be computed. Here, we obtain bounds for W () ( ), W () ( ) and Z () ( ) in terms of finite sums. For every m 1, let ( A (m) i, := 1 {i m} 1 ξ ) i, η i By the interlacing condition, A (m) i, 1 j m,j i 1 ξ i, η j 1 ξ and C (m) i, i, ξ j, := ζ and C (m) i, are all positive and, for every i 1, A (m) i, A i, and C (m) i, C i, as m. Now we define candidates for the upper and lower bounds of W ζ respectively by W (m) ζ (x) := (ψ (ζ )) 1 C (m) i, e(ζ+ξ i,)x, for every m 1 and x, where W (m) ζ (x) := W (m) ζ (x) δ m [e ] ζx + e (ζ+ξ m+1,)x, δ m := γ C (m) i, >, } ξ i, A (m) i, ζ + ξ i,, i 1. which vanishes in the limit as m by (3.4). As candidates for the upper and lower bounds of W (), we also define W (,m) (x) := e ζx W (m) ζ (x) and W (,m) (x) := e ζx W (m) ζ (x), x. The following proposition shows that the scale functions are bounded and approximated by these functions..
14 14 M. EGAMI AND K. YAMAZAKI Proposition 3.2. For every m 1 and x, we have (3.8) Furthermore, we have uniformly on x [, ). W (m) ζ (x) W ζ (x) W (m) ζ (x) and W (,m) (x) W () (x) W (,m) (x). W (m) ζ (x) m W ζ (x) and W (,m) (x) m W () (x), W (m) ζ (x) m W ζ (x) and W (,m) (x) m W () (x), By straightforward calculation, we can bound Z () in (2.6). Let, for every m 1, Z (,m) (x) := 1 + x W (m,) (y)dy and Z (,m) (x) := 1 + Then by Proposition 3.2, we have Z (,m) (x) Z (,m) (x) Z (,m) (x) and Z (,m) (x) Z (,m) (x) = We therefore have the following. x x W (m,) (y)dy, x. ( ) W (,m) (y) W (,m) (y) dy x = δ m [1 ] [ + e ξ m+1,y dy = δ m x + 1 ] (1 e ξm+1,x ). ξ m+1, Corollary 3.1 (Bounds on Z () ). We have Z (,m) (x) Z () (x) and Z (,m) (x) Z () (x) as m pointwise for every x. We now obtain bounds for the derivative. Define, for every x >, w (m) (x) := (ψ (ζ )) 1 ζ e ζx + C (m) i, ξ i,e ξi,x, Here notice that [ ] w (m) (x) := w (m) (x) + max (ξ k,e ξk,x ) + max (ξ k,e ξk,x ) δ m. 1 k m k m+1 max (ξ k,e ξk,x ) = k m+1 Proposition 3.3. For every m 1, we have { 1 ex, ξ m+1, 1 x, ξ m+1, e ξm+1,x, ξ m+1, > 1 x. w (m) (x) W () (x) w (m) (x), x >. Furthermore, we have w (m) (x) W () (x) and w (m) (x) W () (x) uniformly on x x for any x >. A stronger result holds when σ > or Π(, ) <. Recall in this case that θ < by Lemma 3.2-(2) and hence we can define ɛ m := θ ζ ξ i, A (m) i, >, m 1, which vanishes in the limit as m by Lemma 3.2-(2).
15 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 15 Corollary 3.2. When σ > or Π(, ) <, we have where w (m) (x) W () (x) w (m) (x) w (m) (x), x > w (m) (x) := w (m) (x) + max 1 k m (ξ k,e ξ k,x )δ m + e ξ m+1,x ɛ m. j p j η j j p j η j j p j η j E E E-1 6.8E E-9 3.9E E-8 2.2E E-7 1.2E E-6 6.5E E-5 3.5E (i) Weibull(.6,.665) (ii) Pareto(1.2,5) TABLE 1. Parameters of the hyperexponential distributions fitted to (i) Weibull(.6,.665) and (ii) Pareto(1.2, 5) (taken from Tables 3 and 9, respectively, of [15]). 4. NUMERICAL EXAMPLES In this section, we illustrate numerically the effectiveness of the PH-fitting approximation for a general spectrally negative Lévy process. First, we use the classical hyperexponential-fitting algorithm for a completely monotone density function by [15] to approximate the scale function for the case with a (Brownian motion plus) compound Poisson process with Weibull/Pareto-distributed jumps. We approximate (2.5) and (2.24) and compare them with simulated results. Second, we consider, as an example of the meromorphic Lévy process, the β-family introduced by [2] and extend the results to the spectrally negative version of the CGMY process Numerical results on hyperexponential-fitting. As noted earlier, any spectrally negative Lévy process with a completely monotone Lévy density can be approximated arbitrarily closely by hyperexponential-fitting. Here, we use the fitted data computed by [15] to approximate the scale function for a compound Poisson process with i.i.d. Weibull/Pareto-distributed jumps (with or without a Brownian motion component). Recall that the Weibull distribution with parameters c and a (denoted Weibull(c, a)) is give by F (t) = 1 e (t/a)c, t and the Pareto distribution with positive parameters a and b (denoted Pareto(a, b)) is given by F (t) = 1 (1 + bt) a, t. These have long-tails, namely, e δt (1 F (t)) as t for any δ >. See [19] for more details about these distributions. Feldman and Whitt [15] constructed a recursive algorithm to approximate completely monotone
16 16 M. EGAMI AND K. YAMAZAKI σ = 1 σ = x scale function simulation scale function simulation (.313,.3498) (.3797,.3139) (.46517,.4768) (.46792,.47435) (.6337,.63914) (.63372,.63952) (.8123,.81439) (.8131,.81782) Exp(1) σ = 1 σ = x scale function simulation scale function simulation (.4116,.4724) (.41225,.4189) (.55176,.55769) (.568,.5677) (.69767,.7259) (.7229,.7753) (.8474,.85111) (.8517,.85462) Weibull(,6,.665) σ = 1 σ = x scale function simulation scale function simulation (.742,.754) (.7241,.72533) (.8569,.8994) (.81513,.81945) (.883,.88348) (.88491,.88835) (.94248,.9448) (.94635,.94859) Pareto(1.2,5) TABLE 2. Computation of E x [e τ + b 1{τ >τ + b, τ + b < }] via scale function and simulation. The simulation column shows the mean and the 95% confidence interval. densities in terms of hyperexponential densities. We use their results and approximate the scale functions of spectrally negative Lévy processes with Weibull/Pareto-distributed jumps. We consider the Lévy processes X (weibull) and X (pareto) in the form (2.12) where Z is (i) Weibull(.6,.665) and (ii) Pareto(1.2,5), respectively. Table 1 shows the parameters of the hyperexponential distributions obtained by [15] fitted to (i) with m = 6 and to (ii) with m = 14. We use these parameters to construct hyperexponential Lévy processes X (weibull) and X (pareto) (see Example 2.1) that will be used to approximate X (weibull) and X (pareto), respectively. First, we study how accurately the scale function of X (weibull) (resp. X (pareto) ) can be approximated by those of X (weibull) ( X (pareto) ). Toward this end, we use the identity as in (2.5) [ ] (4.1) E x e τ + b 1{τ >τ + b, τ + b < } = W () (x)/w () (b). We compute the right-hand side for X (weibull) ( X (pareto) ) explicitly via (2.2) and approximate the left-hand side for X (weibull) (X (pareto) ) via Monte Carlo simulation based on 1, sample paths. For the simulated results, Brownian motions are approximated by random walks with time step t = T/1 for each interarrival time T
17 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 17 between jumps. In order to confirm the accuracy of the simulated results, we also calculate both sides for the process X (exp) with i.i.d. exponential jumps with parameter 1. For the scale functions of X (exp), X(weibull) and X (pareto), the roots ξ, s and ζ are obtained via the bisection method with error bound 1.E 1. We consider the case σ = and σ = 1 and various values of starting point x with common parameters =.5, λ = 5, b = 5 and µ = 5. Table 2 gives the results. From these, we see that (4.1) for X (weibull) and X (pareto) are approximated very precisely. Based on the accuracy for each x = 1,..., 4, we can also infer that their scale functions are approximated efficiently. Their scale functions and first derivatives for X (exp), X(weibull) and X (pareto) are plotted in Figure 1. It can be easily confirmed that the results are consistent with Lemma 2.2; in particular, W () () = 1/µ =.2 and W () (+) = ( + λ)/(µ 2 ) =.22 for σ = whereas W () () = and W () (+) = 2/(σ 2 ) = 2 for σ = 1. Second, we evaluate the approximation of (2.24) for X (weibull) (X (pareto) ) using those for X (weibull) ( X (pareto) ). Here we compute the overshoot/undershoot density E x [ e τ 1{Xτ da, τ < } ] and E x [ e τ 1{Xτ db, τ < } ]. These values are computed explicitly by (2.25) for X (weibull) ( X (pareto) ). For X (weibull) (X (pareto) ), we simulate [ ] [ ] E x e τ 1{Xτ (a a/2,a+ a/2), τ < } / a and E x e τ 1{Xτ (b b/2,b+ b/2), τ < } / b with a = b =.1 by Monte Carlo simulation with 5, samples. These values are also computed for X (exp). Figures 2-3 show the results for the cases σ = 1 and σ = with common parameters x = 5, µ = 1, λ = 1 and =.5. In Figure 3, the density has a jump at the initial position x = 5 for the case σ = (while it is continuous for the case σ = 1) due to the fact that x = 5 is irregular for (, 5). As can be seen from these figures, the approximation accurately captures the overshoot/undershoot density for X (weibull) and X (pareto). The spike at the initial position in Figure 3 is precisely realized thanks to the closed-form expression (2.25); this would be difficult to realize if the scale function is approximated via numerical Laplace inversion Numerical results on the β-class and CGMY process. We now consider, as an example of meromorphic Lévy processes, the β-class introduced by [2]. The following definition is due to [2], Definition 4. Definition 4.1. A spectrally negative Lévy process is said to be in the β-class if its Lévy density is in the form (4.2) e αβx π(x) = c (1 e βx ) λ 1 {x<}, x R, for some α >, β >, c and λ (, 3). It is euivalent to saying that its Laplace exponent is ψ(z) = ˆµz σ2 z 2 + c {B(α + zβ } β, 1 λ) B(α, 1 λ) where B is the beta function, i.e., B(x, y) := Γ(x)Γ(y)/Γ(x + y). The special case σ = and β = 1 reduces to the class of Lamperti-stable processes, which are obtained by the Lamperti transformation [31] from the stable processes conditioned to stay positive; see [8, 9] and references therein. For the scale function of a related process, see [28].
18 18 M. EGAMI AND K. YAMAZAKI σ σ σ σ Scale function for the Exp(1) case Derivative for the Exp(1) case σ σ σ σ Scale function for the Weibull(.6,.665) case Derivative for the Weibull(.6,.665) case σ σ σ σ Scale function for the Pareto(1.2,5) case Derivative for the Pareto(1.2,5) case FIGURE 1. Scale functions and their derivatives.
19 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 19 Exp(1) with σ = 1 Exp(1) with σ = Weibull(.6,.665) with σ = 1 Weibull(.6,.665) with σ = Pareto(1.2,5) with σ = 1 Pareto(1.2,5) with σ = FIGURE 2. Computation of the overshoot density E x [e τ 1{Xτ da,τ < }].
20 2 M. EGAMI AND K. YAMAZAKI Exp(1) with σ = 1 Exp(1) with σ = Weibull(.6,.665) with σ = 1 Weibull(.6,.665) with σ = Pareto(1.2,5) with σ = 1 Pareto(1.2,5) with σ = FIGURE 3. Computation of the undershoot density E x [e τ 1{Xτ db,τ < }].
21 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 21 It can be also seen that this is a discrete-version of the (spectrally negative) CGMY process, whose Lévy density is given by (4.3) π(x) = c eαx x λ 1 {x<}, x R. Indeed, if we set c = cβ λ and α = αβ 1 in (4.2), we have (4.4) e αβx c (1 e βx ) λ 1 β e αx {x<} c x λ 1 {x<}, x R. See [4] for approximation of (double-sided) CGMY processes using hyperexponential distributions. We use the results in Section 3 to obtain the bounds on the scale functions. We consider the case =.3, σ =.2, ˆµ =.1, λ = 1.5, α = 3, β = 1 and c =.1 in (4.2). It is easy to verify that these have jumps of infinite activity (and bounded variation). First, we plot, in Figure 4, the upper and lower bounds on the scale function and its derivative as obtain in Proposition 3.2 and Corollary 3.2 for m = 1 and m = 1. The difference between the upper and lower bounds indeed converges to zero, although the convergence of the derivative in the neighborhood of zero is relatively slow. Second, motivated by (4.4), we take β to zero and see how the approximation for the CGMY process works. Here we set α = 3, c =.1 and m = 1 and use the same values as the above for the other parameters. Figure 5 shows the upper and lower bounds of scale function and its derivative for various values of β; we can indeed observe the convergence as β. 5. CONCLUDING REMARKS We have studied the scale function for the spectrally negative PH Lévy process and the PH-fitting approach for the approximation of the scale function for a general spectrally negative Lévy process. Because the approximated scale function is given as a function in a closed form, one can analytically differentiate/integrate to obtain other fluctuation identities explicitly. Our numerical results suggest that the PH-fitting is a powerful alternative to the numerical Laplace inversion approach, particularly when the Lévy density is completely monotone. One major challenge for the PH-fitting approach is that, without the completely monotone assumption, there does not exist a fitting algorithm that is guaranteed to work for arbitrary measure. The performance is certainly dependent on the shape of the Lévy measure and it may not be a suitable approach for certain cases. However, there exist a variety of fitting algorithms typically developed in ueueing analysis. Well-known examples are the moment-matching approach (e.g. MEFIT and MEDA) and the maximum-likelihood approach (e.g. MLAPH and EMPHT), and a thorough study of pros and cons of each fitting techniues has been conducted in, for example, [17, 32]. Our next step is, therefore, to apply these existing algorithms for the approximation of the scale function, and analyze its performance for a variety of Lévy measures. Another direction for future research is the PH-fitting construction of scale functions from empirical data as in, for example, [4]. The closed-form expression of the approximated scale function can be used flexibly to identify the fluctuation of the process implied by the empirical data. APPENDIX A. PROOFS
22 22 M. EGAMI AND K. YAMAZAKI Bounds on W () (x) when m = 1 Bounds on W () (x) when m = 1 Bounds on W () (x) when m = 1 Bounds on W () (x) when m = 1 FIGURE 4. Approximation of the scale function and its derivative for the β-class. A.1. Proof of Lemmas By (2.17), it is easy to verify that [ ] E x e τ a 1 {τ a < } = n m i k=1 A (k) i, ξ (ξ i, y) k1 i, xa (k 1)! eξ i,y dy, a < x,
23 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 23 upper bounds for the scale functions upper bounds for the derivatives lower bounds for the scale functions lower bounds for the derivatives FIGURE 5. Convergence to the CGMY model. and hence, because W () C 1 (, ), (A.1) (A.2) [ ] a Ex e τ a 1 {τ a < } = x Ex [ e τ 1{τ < } ] x=+ = n m i A (k) i, ξ (ξ i, (x a)) k1 i, k=1 n A (1) i, ξ i,. (k 1)! e ξ i,(xa), a < x, On the other hand, different representations of (A.1)-(A.2) can be pursued. By Theorem 8.1 of [25] and (2.7), [ ] E x e τ a 1 {τ a < } = Z () (x a) xa W () (x a) = 1 + W () (y)dy e ζ(xa) W ζ (x a) ζ ζ
24 24 M. EGAMI AND K. YAMAZAKI for every a < x. Its derivative with respect to a becomes [ ] a Ex e τ a 1 {τ a < } = W () (x a) + e ζ(xa) W ζ (x a) + e ζ(xa) W ζ ζ (x a) (A.3) = e ζ(xa) W ζ ζ (x a). In particular, when a =, the derivative with respect to x and its limit as x are [ ] x Ex e τ 1{τ < } = ] [ζ W () x + (A.4) (x) + W () (x) = θ. ζ ζ By matching (A.2) and (A.4), Lemma 2.3 is immediate. For the proof of Lemma 2.4, by matching (A.1) and (A.3), we have W ζ (y) = ζ = ζ n m i k=1 m i n k=1 A (k) i, ξ (ξ i, y) k1 i, (k 1)! e(ζ+ξ i,)y Integrating the above and changing variables, we have W ζ (x) W ζ () = ζ ( ) k A (k) i, (ζ ξi, ((ζ + ξ i, )y) k1 + ξ i, ) e (ζ+ξi,)y, y. ζ + ξ i, (k 1)! n m i k=1 ( ) k A (k) ξi, 1 i, ζ + ξ i, (k 1)! (ζ+ξ i, )x z k1 e z dz, x. Lemma 2.4 is now immediate because the integral part is a lower incomplete gamma function. A.2. Proof of Lemma 2.5. Let N(, ) be the Poisson random measure for X and X t := min s t X s for all t. By the compensation formula, we have [ ] [ ] E x e τ 1{Xτ B, X τ A} = E x N(dt, du)e t 1 {Xt u A, X t B, X t >} [ ] = E x e t dt Π(du)1 {Xt u A, X t B, X t >} = = Π(du) Π(du) By using the -resolvent kernel (see, e.g., [25]) dt [ e t P x {X t B (A + u), τ t}] = Substituting this, we obtain the result. B (A+u) = W () (x) dt [ e t P x {X t B (A + u), X t > } ] dt [ e t P x {X t B (A + u), τ t}]. [ ] dy e ζy W () (x) W () (x y) B (A+u) e ζy dy dyw () (x y). B (A+u)
25 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 25 A.3. Proof of Proposition 2.2. Define for all B B(, ), A B(, ) and K R ρ(k; A, B) := Π(du) e Ky dy. Because by Corollary 2.1 and (2.21) B (A+u) Π(du)W () (x) e ζy dy = W () (x)ρ(ζ ; A, B) = B (A+u) 1 ψ (ζ ) eζx C i, e ξ i,x ρ(ζ ; A, B), i I and = Π(du) W () (x y)dy B (,x) (A+u) Π(du) B (,x) (A+u) 1 ψ (ζ ) eζ(xy) C i, e ξ i,(xy) dy i I we obtain via Lemma 2.5 h (x; A, B) = = eζx ψ (ζ ) ρ(ζ ; A, B (, x)) i I C i, e ξ i,x ρ(ξ i, ; A, B (, x)), eζx ψ (ζ ) C i, e ξ i,x ρ(ζ ; A, B) i I eζx ψ (ζ ) ρ(ζ ; A, B (, x)) + i I C i, e ξ i,x ρ(ξ i, ; A, B (, x)) = eζx ψ (ζ ) ρ(ζ ; A, B (x, )) + i I C i, e ξ i,x {ρ(ξ i, ; A, B (, x)) ρ(ζ ; A, B)}. Lemma A.1. Suppose B = (b, b) and A = (a, a). For any K R, we have λp ( ) j ρ(k; A, B) = η j K (eη ja e ηja ) e (ηjk)b e (η jk)b. Proof. For any u > e Ky dy = B (A+u) (b,b) (a+u,a+u) e Ky dy = 1 {a+b<u<a+b} e Ky dy = 1 ((a+u) b,(a+u) b) K For its integral with respect to Π, a+b a+b a+b a+b Π(du)e K((a+u) b) = Π(du)e K((a+u) b) = a+b a+b a+b a+b Π(du)e K(a+u) + Π(du)e Kb + a+b a+b a+b a+b [ e K((a+u) b) e K((a+u) b)] 1 {a+b<u<a+b}. a Π(du)e Kb = e Kb Π(du + b) + a Π(du)e K(a+u) = e Kb Π(du + b) + a a b b b b Π(du + a)e Ku, Π(du + a)e Ku.
26 26 M. EGAMI AND K. YAMAZAKI Therefore ρ(k; A, B) = 1 K ( a e Kb Π(du + b) + a b Using Π(du) = λ m p jη j e ηju du, for any β, β and α, Thus β β Π(du + α)e Ku = λ p j η j e η jα ρ(k; A, B) = λ K λ K = λ b β β a Π(du + a)e Ku e Kb Π(du + b) e (η jk)u = λ a b b Π(du + a)e Ku ) p j η j η j K eη jα ( e (η jk)β e (η jk)β ). p j η j η j K (eη ja e ηja ) (e ) (ηjk)b e (η jk)b p j (e (ηjk)b e (ηjk)b ) ( e ηja e η ja ) We are now ready to prove Proposition 2.2. By Lemma A.1, e ζx ρ(ζ ; A, B (x, )) = e ζx λ p ( ) j η j K (eη ja e ηja ) e (ηjk)b e (η jk)b. p ( ) j (e ηja e ηja ) e (η j+ζ )(b x) e (η j+ζ )(b x), η j + ζ C i, e ξi,x ρ(ζ ; A, B) = λ C i, e ξi,x p ( ) j (e ηja e ηja ) e (η j+ζ )b e (η j+ζ )b, η i I i I j + ζ C i, e ξi,x ρ(ξ i, ; A, B (, x)) i I Hence h (x; A, B) = = λ i I C i, e ξ i,x λ ψ (ζ ) eζx + λ i I C i, e ξ i,x λ i I C i, e ξ i,x = λ p ( ) j (e ηja e ηja ) e (η jξ i, )(b x) e (η jξ i, )(b x). η j ξ i, p j (e ηja e ηja ) (e ) (η j+ζ )(b x) e (η j+ζ )(b x) η j + ζ p j (e ηja e ηja )κ j, (x; B). p j (e ηja e ηja ) (e ) (η jξ i, )(b x) e (η jξ i, )(b x) η j ξ i, p j (e ηja e ηja ) (e ) (η j+ζ )b e (η j+ζ )b η j + ζ.
27 ON SCALE FUNCTIONS OF SPECTRALLY NEGATIVE LÉVY PROCESSES WITH PHASE-TYPE JUMPS 27 A.4. Proof of Proposition 3.2. Notice, for every m 1, that ( ) C i, C (m) i, γ C (m) i, = δ m, (A.5) and hence by (A.5) W (m) ζ (x) W ζ (x) = e ζx i=m+1 ( m ( C i, = γ C i, C (m) i, C i, C (m) i, C i, γ ) e (ζ+ξ i,)x + ) + e (ζ+ξ m+1,)x i=m+1 i=m+1 C (m) i, = δ m, C i, e (ζ+ξ i,)x C i, δ m [e ζx + e (ζ+ξ m+1,)x ]. Therefore we have the bounds for W ζ in (3.8). The bounds for W () are immediate by multiplying e ζx. Finally, the convergence results hold because δ m and e ζx +e (ζ+ξ m+1,)x and 1+e ξ m+1,x are bounded uniformly in x. A.5. Proof of Proposition 3.3. The lower bound is immediate by the fact that C (m) i, C i,. For every x >, we have by (A.5) W () (x) = (ψ (ζ )) 1 ζ e ζx + C i, ξ i, e ξ i,x = w (m) (x) + (C i, C (m) i, )ξ i,e ξ i,x + w (m) (x) + max 1 k m (ξ k,e ξ k,x ) [ w (m) (x) + i=m+1 (C i, C (m) i, max (ξ k,e ξk,x ) + max (ξ k,e ξk,x ) 1 k m k m+1 C i, ξ i, e ξ i,x ) + max (ξ k,e ξk,x ) k m+1 ] δ m = w (m) (x), which shows the first claim. For the second claim, notice that for the given x >, [ ] w (m) (x) w (m) (x) max (ξ k,e ξ k,x ) + max (ξ k,e ξ k,x ) 1 k m k m+1 i=m+1 uniformly on [x, ). Because δ m vanishes as m and sup λ λe λx <, the convergence is immediate. A.6. Proof of Corollary 3.2. For every x >, W () (x) = w (m) (x) + (C i, C (m) i, )ξ i,e ξi,x + w (m) (x) + max 1 k m (ξ k,e ξ k,x ) i=m+1 (C i, C (m) i, w (m) (x) + max 1 k m (ξ k,e ξ k,x )δ m + ɛ m e ξ m+1,x C i, ξ i, e ξ i,x ) + ζ eξ m+1,x δ m i=m+1 C i, A i, ξ i,
28 28 M. EGAMI AND K. YAMAZAKI where the last ineuality holds because ζ k=m+1 ξ i, A i, = θ ζ k=1 This together with Proposition 3.3 shows the claim. ξ i, A i, θ ζ k=1 ξ i, A (m) i, = ɛ m, m 1. REFERENCES [1] S. Asmussen. Fitting phase-type distributions via the EM algorithm. Scand. J. Statist, 23: , [2] S. Asmussen. Applied probability and ueues, volume 51 of Applications of Mathematics (New York). Springer-Verlag, New York, second edition, 23. Stochastic Modelling and Applied Probability. [3] S. Asmussen, F. Avram, and M. R. Pistorius. Russian and American put options under exponential phase-type Lévy models. Stochastic Process. Appl., 19(1):79 111, 24. [4] S. Asmussen, D. B. Madan, and M. R. Pistorius. Pricing euity default swaps under an approximation to the CGMY Lévy model. Journal of Computational Finance, 11(2):79 93, 27. [5] F. Avram, Z. Palmowski, and M. R. Pistorius. On the optimal dividend problem for a spectrally negative Lévy process. Ann. Appl. Probab., 17(1):156 18, 27. [6] O. E. Barndorff-Nielsen. Processes of normal inverse Gaussian type. Finance Stoch., 2(1):41 68, [7] J. Bertoin. Lévy processes, volume 121 of Cambridge Tracts in Mathematics. Cambridge University Press, Cambridge, [8] J. Bertoin and M. Yor. On subordinators, self-similar Markov processes and some factorizations of the exponential variable. Electron. Comm. Probab., 6:95 16 (electronic), 21. [9] M. Caballero, J. Pardo, and J. Pérez. On Lamperti stable processes. Preprint, 28. [1] M. E. Caballero and L. Chaumont. Conditioned stable Lévy processes and the Lamperti representation. J. Appl. Probab., 43(4): , 26. [11] P. Carr, G. Hélyette, D. B. Madan, and M. Yor. The structure of asset returns: an emperical investigation. Journal of Business, 75:35 332, 22. [12] T. Chan, A. E. Kyprianou, and M. Savov. Smoothness of scale functions for spectrally negative Lévy processes. Probability Theory and Related Fields, 29. [13] L. Chaumont, A. E. Kyprianou, and J. C. Pardo. Some explicit identities associated with positive self-similar Markov processes. Stochastic Process. Appl., 119(3):98 1, 29. [14] E. Eberlein, U. Keller, and K. Prause. New insights into smile, mispricing and value at risk: the hyperbolic model. Journal of Business, 71:371 45, [15] A. Feldmann and W. Whitt. Fitting mixtures of exponentials to long-tail distributions to analyze network performance models. Performance Evaluation, (31): , [16] W. Feller. An introduction to probability theory and its applications. Vol. II. Second edition. John Wiley & Sons Inc., New York, [17] A. Horváth and M. Telek. Approximating heavy tailed behaviour with phase type distributions. Advances in Algorithmic Methods for Stochastic Models, 2. [18] J. Jacod and A. N. Shiryaev. Limit theorems for stochastic processes, volume 288 of Grundlehren der Mathematischen Wissenschaften [Fundamental Principles of Mathematical Sciences]. Springer-Verlag, Berlin, second edition, 23. [19] N. L. Johnson and S. Kotz. Distributions in statistics: continuous multivariate distributions. John Wiley & Sons Inc., New York, Wiley Series in Probability and Mathematical Statistics.
Scale functions for spectrally negative Lévy processes and their appearance in economic models
Scale functions for spectrally negative Lévy processes and their appearance in economic models Andreas E. Kyprianou 1 Department of Mathematical Sciences, University of Bath 1 This is a review talk and
More informationAnalytical properties of scale functions for spectrally negative Lévy processes
Analytical properties of scale functions for spectrally negative Lévy processes Andreas E. Kyprianou 1 Department of Mathematical Sciences, University of Bath 1 Joint work with Terence Chan and Mladen
More informationErnesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010
Optimal stopping for Hunt and Lévy processes Ernesto Mordecki 1 Lecture III. PASI - Guanajuato - June 2010 1Joint work with Paavo Salminen (Åbo, Finland) 1 Plan of the talk 1. Motivation: from Finance
More informationBernstein-gamma functions and exponential functionals of Lévy processes
Bernstein-gamma functions and exponential functionals of Lévy processes M. Savov 1 joint work with P. Patie 2 FCPNLO 216, Bilbao November 216 1 Marie Sklodowska Curie Individual Fellowship at IMI, BAS,
More informationPoisson random measure: motivation
: motivation The Lévy measure provides the expected number of jumps by time unit, i.e. in a time interval of the form: [t, t + 1], and of a certain size Example: ν([1, )) is the expected number of jumps
More informationOn Optimal Stopping Problems with Power Function of Lévy Processes
On Optimal Stopping Problems with Power Function of Lévy Processes Budhi Arta Surya Department of Mathematics University of Utrecht 31 August 2006 This talk is based on the joint paper with A.E. Kyprianou:
More informationErik J. Baurdoux Some excursion calculations for reflected Lévy processes
Erik J. Baurdoux Some excursion calculations for reflected Lévy processes Article (Accepted version) (Refereed) Original citation: Baurdoux, Erik J. (29) Some excursion calculations for reflected Lévy
More informationPotentials of stable processes
Potentials of stable processes A. E. Kyprianou A. R. Watson 5th December 23 arxiv:32.222v [math.pr] 4 Dec 23 Abstract. For a stable process, we give an explicit formula for the potential measure of the
More informationA LOGARITHMIC EFFICIENT ESTIMATOR OF THE PROBABILITY OF RUIN WITH RECUPERATION FOR SPECTRALLY NEGATIVE LEVY RISK PROCESSES.
A LOGARITHMIC EFFICIENT ESTIMATOR OF THE PROBABILITY OF RUIN WITH RECUPERATION FOR SPECTRALLY NEGATIVE LEVY RISK PROCESSES Riccardo Gatto Submitted: April 204 Revised: July 204 Abstract This article provides
More informationLeast Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises
Least Squares Estimators for Stochastic Differential Equations Driven by Small Lévy Noises Hongwei Long* Department of Mathematical Sciences, Florida Atlantic University, Boca Raton Florida 33431-991,
More informationarxiv: v2 [math.oc] 5 Feb 2015
OPTIMALITY OF DOUBLY REFLECTED LÉVY PROCESSES IN SINGULAR CONTROL ERIK J. BAURDOUX AND KAZUTOSHI YAMAZAKI arxiv:148.847v2 [math.oc 5 Feb 215 ABSTRACT. We consider a class of two-sided singular control
More informationJump-type Levy Processes
Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic
More informationPrecautionary Measures for Credit Risk Management in. Jump Models
Kyoto University, Graduate School of Economics Research Project Center Discussion Paper Series Precautionary Measures for Credit Risk Management in Jump Models Masahiko Egami and Kazutoshi Yamazaki Discussion
More informationConvergence of price and sensitivities in Carr s randomization approximation globally and near barrier
Convergence of price and sensitivities in Carr s randomization approximation globally and near barrier Sergei Levendorskĭi University of Leicester Toronto, June 23, 2010 Levendorskĭi () Convergence of
More informationStable and extended hypergeometric Lévy processes
Stable and extended hypergeometric Lévy processes Andreas Kyprianou 1 Juan-Carlos Pardo 2 Alex Watson 2 1 University of Bath 2 CIMAT as given at CIMAT, 23 October 2013 A process, and a problem Let X be
More informationABC methods for phase-type distributions with applications in insurance risk problems
ABC methods for phase-type with applications problems Concepcion Ausin, Department of Statistics, Universidad Carlos III de Madrid Joint work with: Pedro Galeano, Universidad Carlos III de Madrid Simon
More informationExponential functionals of Lévy processes
Exponential functionals of Lévy processes Víctor Rivero Centro de Investigación en Matemáticas, México. 1/ 28 Outline of the talk Introduction Exponential functionals of spectrally positive Lévy processes
More informationA NEW PROOF OF THE WIENER HOPF FACTORIZATION VIA BASU S THEOREM
J. Appl. Prob. 49, 876 882 (2012 Printed in England Applied Probability Trust 2012 A NEW PROOF OF THE WIENER HOPF FACTORIZATION VIA BASU S THEOREM BRIAN FRALIX and COLIN GALLAGHER, Clemson University Abstract
More informationA REFINED FACTORIZATION OF THE EXPONENTIAL LAW P. PATIE
A REFINED FACTORIZATION OF THE EXPONENTIAL LAW P. PATIE Abstract. Let ξ be a (possibly killed) subordinator with Laplace exponent φ and denote by I φ e ξs ds, the so-called exponential functional. Consider
More informationThe optimality of a dividend barrier strategy for Lévy insurance risk processes, with a focus on the univariate Erlang mixture
The optimality of a dividend barrier strategy for Lévy insurance risk processes, with a focus on the univariate Erlang mixture by Javid Ali A thesis presented to the University of Waterloo in fulfilment
More informationA Barrier Version of the Russian Option
A Barrier Version of the Russian Option L. A. Shepp, A. N. Shiryaev, A. Sulem Rutgers University; shepp@stat.rutgers.edu Steklov Mathematical Institute; shiryaev@mi.ras.ru INRIA- Rocquencourt; agnes.sulem@inria.fr
More informationThe exact asymptotics for hitting probability of a remote orthant by a multivariate Lévy process: the Cramér case
The exact asymptotics for hitting probability of a remote orthant by a multivariate Lévy process: the Cramér case Konstantin Borovkov and Zbigniew Palmowski Abstract For a multivariate Lévy process satisfying
More informationObstacle problems for nonlocal operators
Obstacle problems for nonlocal operators Camelia Pop School of Mathematics, University of Minnesota Fractional PDEs: Theory, Algorithms and Applications ICERM June 19, 2018 Outline Motivation Optimal regularity
More informationInformation and Credit Risk
Information and Credit Risk M. L. Bedini Université de Bretagne Occidentale, Brest - Friedrich Schiller Universität, Jena Jena, March 2011 M. L. Bedini (Université de Bretagne Occidentale, Brest Information
More informationStochastic Areas and Applications in Risk Theory
Stochastic Areas and Applications in Risk Theory July 16th, 214 Zhenyu Cui Department of Mathematics Brooklyn College, City University of New York Zhenyu Cui 49th Actuarial Research Conference 1 Outline
More informationSolving the Poisson Disorder Problem
Advances in Finance and Stochastics: Essays in Honour of Dieter Sondermann, Springer-Verlag, 22, (295-32) Research Report No. 49, 2, Dept. Theoret. Statist. Aarhus Solving the Poisson Disorder Problem
More informationON ADDITIVE TIME-CHANGES OF FELLER PROCESSES. 1. Introduction
ON ADDITIVE TIME-CHANGES OF FELLER PROCESSES ALEKSANDAR MIJATOVIĆ AND MARTIJN PISTORIUS Abstract. In this note we generalise the Phillips theorem [1] on the subordination of Feller processes by Lévy subordinators
More informationPractical approaches to the estimation of the ruin probability in a risk model with additional funds
Modern Stochastics: Theory and Applications (204) 67 80 DOI: 05559/5-VMSTA8 Practical approaches to the estimation of the ruin probability in a risk model with additional funds Yuliya Mishura a Olena Ragulina
More informationOn Reflecting Brownian Motion with Drift
Proc. Symp. Stoch. Syst. Osaka, 25), ISCIE Kyoto, 26, 1-5) On Reflecting Brownian Motion with Drift Goran Peskir This version: 12 June 26 First version: 1 September 25 Research Report No. 3, 25, Probability
More informationarxiv: v2 [math.pr] 15 Jul 2015
STRIKINGLY SIMPLE IDENTITIES RELATING EXIT PROBLEMS FOR LÉVY PROCESSES UNDER CONTINUOUS AND POISSON OBSERVATIONS arxiv:157.3848v2 [math.pr] 15 Jul 215 HANSJÖRG ALBRECHER AND JEVGENIJS IVANOVS Abstract.
More information1 Basic concepts from probability theory
Basic concepts from probability theory This chapter is devoted to some basic concepts from probability theory.. Random variable Random variables are denoted by capitals, X, Y, etc. The expected value or
More informationNumerical Methods with Lévy Processes
Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:
More informationThe Azéma-Yor Embedding in Non-Singular Diffusions
Stochastic Process. Appl. Vol. 96, No. 2, 2001, 305-312 Research Report No. 406, 1999, Dept. Theoret. Statist. Aarhus The Azéma-Yor Embedding in Non-Singular Diffusions J. L. Pedersen and G. Peskir Let
More informationUniversity Of Calgary Department of Mathematics and Statistics
University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National
More informationAn Introduction to Probability Theory and Its Applications
An Introduction to Probability Theory and Its Applications WILLIAM FELLER (1906-1970) Eugene Higgins Professor of Mathematics Princeton University VOLUME II SECOND EDITION JOHN WILEY & SONS Contents I
More information1 Infinitely Divisible Random Variables
ENSAE, 2004 1 2 1 Infinitely Divisible Random Variables 1.1 Definition A random variable X taking values in IR d is infinitely divisible if its characteristic function ˆµ(u) =E(e i(u X) )=(ˆµ n ) n where
More informationLévy-VaR and the Protection of Banks and Insurance Companies. A joint work by. Olivier Le Courtois, Professor of Finance and Insurance.
Lévy-VaR and the Protection of Banks and Insurance Companies A joint work by Olivier Le Courtois, Professor of Finance and Insurance and Christian Walter, Actuary and Consultant Institutions : EM Lyon
More informationA Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1
Chapter 3 A Change of Variable Formula with Local Time-Space for Bounded Variation Lévy Processes with Application to Solving the American Put Option Problem 1 Abstract We establish a change of variable
More informationarxiv: v1 [math.pr] 30 Mar 2014
Binomial discrete time ruin probability with Parisian delay Irmina Czarna Zbigniew Palmowski Przemys law Świ atek October 8, 2018 arxiv:1403.7761v1 [math.pr] 30 Mar 2014 Abstract. In this paper we analyze
More informationSupermodular ordering of Poisson arrays
Supermodular ordering of Poisson arrays Bünyamin Kızıldemir Nicolas Privault Division of Mathematical Sciences School of Physical and Mathematical Sciences Nanyang Technological University 637371 Singapore
More informationDefinition: Lévy Process. Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes. Theorem
Definition: Lévy Process Lectures on Lévy Processes and Stochastic Calculus, Braunschweig, Lecture 2: Lévy Processes David Applebaum Probability and Statistics Department, University of Sheffield, UK July
More informationOn the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem
On the martingales obtained by an extension due to Saisho, Tanemura and Yor of Pitman s theorem Koichiro TAKAOKA Dept of Applied Physics, Tokyo Institute of Technology Abstract M Yor constructed a family
More informationOptimal stopping, Appell polynomials and Wiener-Hopf factorization
arxiv:1002.3746v1 [math.pr] 19 Feb 2010 Optimal stopping, Appell polynomials and Wiener-Hopf factorization Paavo Salminen Åbo Aademi, Mathematical Department, Fänrisgatan 3 B, FIN-20500 Åbo, Finland, email:
More informationFinite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims
Proceedings of the World Congress on Engineering 29 Vol II WCE 29, July 1-3, 29, London, U.K. Finite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims Tao Jiang Abstract
More informationFunctional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals
Functional Limit theorems for the quadratic variation of a continuous time random walk and for certain stochastic integrals Noèlia Viles Cuadros BCAM- Basque Center of Applied Mathematics with Prof. Enrico
More informationOn the quantiles of the Brownian motion and their hitting times.
On the quantiles of the Brownian motion and their hitting times. Angelos Dassios London School of Economics May 23 Abstract The distribution of the α-quantile of a Brownian motion on an interval [, t]
More informationThreshold dividend strategies for a Markov-additive risk model
European Actuarial Journal manuscript No. will be inserted by the editor Threshold dividend strategies for a Markov-additive risk model Lothar Breuer Received: date / Accepted: date Abstract We consider
More informationIntroduction to self-similar growth-fragmentations
Introduction to self-similar growth-fragmentations Quan Shi CIMAT, 11-15 December, 2017 Quan Shi Growth-Fragmentations CIMAT, 11-15 December, 2017 1 / 34 Literature Jean Bertoin, Compensated fragmentation
More informationThe equivalence of two tax processes
The equivalence of two ta processes Dalal Al Ghanim Ronnie Loeffen Ale Watson 6th November 218 arxiv:1811.1664v1 [math.pr] 5 Nov 218 We introduce two models of taation, the latent and natural ta processes,
More informationOn an Effective Solution of the Optimal Stopping Problem for Random Walks
QUANTITATIVE FINANCE RESEARCH CENTRE QUANTITATIVE FINANCE RESEARCH CENTRE Research Paper 131 September 2004 On an Effective Solution of the Optimal Stopping Problem for Random Walks Alexander Novikov and
More informationAsymptotics for posterior hazards
Asymptotics for posterior hazards Pierpaolo De Blasi University of Turin 10th August 2007, BNR Workshop, Isaac Newton Intitute, Cambridge, UK Joint work with Giovanni Peccati (Université Paris VI) and
More informationTHE VARIANCE CONSTANT FOR THE ACTUAL WAITING TIME OF THE PH/PH/1 QUEUE. By Mogens Bladt National University of Mexico
The Annals of Applied Probability 1996, Vol. 6, No. 3, 766 777 THE VARIANCE CONSTANT FOR THE ACTUAL WAITING TIME OF THE PH/PH/1 QUEUE By Mogens Bladt National University of Mexico In this paper we consider
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationLecture 17 Brownian motion as a Markov process
Lecture 17: Brownian motion as a Markov process 1 of 14 Course: Theory of Probability II Term: Spring 2015 Instructor: Gordan Zitkovic Lecture 17 Brownian motion as a Markov process Brownian motion is
More informationErnesto Mordecki. Talk presented at the. Finnish Mathematical Society
EXACT RUIN PROBABILITIES FOR A CLASS Of LÉVY PROCESSES Ernesto Mordecki http://www.cmat.edu.uy/ mordecki Montevideo, Uruguay visiting Åbo Akademi, Turku Talk presented at the Finnish Mathematical Society
More informationLecture 21 Representations of Martingales
Lecture 21: Representations of Martingales 1 of 11 Course: Theory of Probability II Term: Spring 215 Instructor: Gordan Zitkovic Lecture 21 Representations of Martingales Right-continuous inverses Let
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationStochastic Processes
Stochastic Processes A very simple introduction Péter Medvegyev 2009, January Medvegyev (CEU) Stochastic Processes 2009, January 1 / 54 Summary from measure theory De nition (X, A) is a measurable space
More informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
More informationIntroduction to Rare Event Simulation
Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31
More informationThe Subexponential Product Convolution of Two Weibull-type Distributions
The Subexponential Product Convolution of Two Weibull-type Distributions Yan Liu School of Mathematics and Statistics Wuhan University Wuhan, Hubei 4372, P.R. China E-mail: yanliu@whu.edu.cn Qihe Tang
More information1. Stochastic Processes and filtrations
1. Stochastic Processes and 1. Stoch. pr., A stochastic process (X t ) t T is a collection of random variables on (Ω, F) with values in a measurable space (S, S), i.e., for all t, In our case X t : Ω S
More informationSOME REMARKS ON FIRST PASSAGE OF LÉVY PROCESSES, THE AMERICAN PUT AND PASTING PRINCIPLES
The Annals of Applied Probability 25, Vol. 15, No. 3, 262 28 DOI 1.1214/155165377 Institute of Mathematical Statistics, 25 SOME REMARKS ON FIRST PASSAGE OF LÉVY PROCESSES, THE AMERICAN PUT AND PASTING
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationWiener Measure and Brownian Motion
Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u
More informationPavel V. Gapeev, Neofytos Rodosthenous Perpetual American options in diffusion-type models with running maxima and drawdowns
Pavel V. Gapeev, Neofytos Rodosthenous Perpetual American options in diffusion-type models with running maxima and drawdowns Article (Accepted version) (Refereed) Original citation: Gapeev, Pavel V. and
More informationOn Minimal Entropy Martingale Measures
On Minimal Entropy Martingale Measures Andrii Andrusiv Friedrich Schiller University of Jena, Germany Marie Curie ITN Outline 1. One Period Model Definitions and Notations Minimal Entropy Martingale Measure
More informationAuthor(s) Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 41(1)
Title On the stability of contact Navier-Stokes equations with discont free b Authors Huang, Feimin; Matsumura, Akitaka; Citation Osaka Journal of Mathematics. 4 Issue 4-3 Date Text Version publisher URL
More informationLecture 22 Girsanov s Theorem
Lecture 22: Girsanov s Theorem of 8 Course: Theory of Probability II Term: Spring 25 Instructor: Gordan Zitkovic Lecture 22 Girsanov s Theorem An example Consider a finite Gaussian random walk X n = n
More informationRussian and American put options under exponential phase type Lévy models
Russian and American put options under exponential phase type Lévy models Søren Asmussen a, Florin Avram b, Martijn R. Pistorius c, a Department of Mathematical Sciences, Aarhus University, Ny Munkegade,
More informationBrownian motion. Samy Tindel. Purdue University. Probability Theory 2 - MA 539
Brownian motion Samy Tindel Purdue University Probability Theory 2 - MA 539 Mostly taken from Brownian Motion and Stochastic Calculus by I. Karatzas and S. Shreve Samy T. Brownian motion Probability Theory
More informationStatistical inference on Lévy processes
Alberto Coca Cabrero University of Cambridge - CCA Supervisors: Dr. Richard Nickl and Professor L.C.G.Rogers Funded by Fundación Mutua Madrileña and EPSRC MASDOC/CCA student workshop 2013 26th March Outline
More informationWeb Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.
Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we
More informationReinsurance and ruin problem: asymptotics in the case of heavy-tailed claims
Reinsurance and ruin problem: asymptotics in the case of heavy-tailed claims Serguei Foss Heriot-Watt University, Edinburgh Karlovasi, Samos 3 June, 2010 (Joint work with Tomasz Rolski and Stan Zachary
More informationStrong log-concavity is preserved by convolution
Strong log-concavity is preserved by convolution Jon A. Wellner Abstract. We review and formulate results concerning strong-log-concavity in both discrete and continuous settings. Although four different
More informationAsymptotic Distributions of the Overshoot and Undershoots for the Lévy Insurance Risk Process in the Cramér and Convolution Equivalent Cases
To appear in Insurance Math. Econom. Asymptotic Distributions of the Overshoot and Undershoots for the Lévy Insurance Risk Process in the Cramér and Convolution Euivalent Cases Philip S. Griffin, Ross
More informationUniversal examples. Chapter The Bernoulli process
Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationConvergence of Feller Processes
Chapter 15 Convergence of Feller Processes This chapter looks at the convergence of sequences of Feller processes to a iting process. Section 15.1 lays some ground work concerning weak convergence of processes
More information1 Brownian Local Time
1 Brownian Local Time We first begin by defining the space and variables for Brownian local time. Let W t be a standard 1-D Wiener process. We know that for the set, {t : W t = } P (µ{t : W t = } = ) =
More informationProbability and Measure
Chapter 4 Probability and Measure 4.1 Introduction In this chapter we will examine probability theory from the measure theoretic perspective. The realisation that measure theory is the foundation of probability
More informationRegular Variation and Extreme Events for Stochastic Processes
1 Regular Variation and Extreme Events for Stochastic Processes FILIP LINDSKOG Royal Institute of Technology, Stockholm 2005 based on joint work with Henrik Hult www.math.kth.se/ lindskog 2 Extremes for
More informationRandomization in the First Hitting Time Problem
Randomization in the First Hitting Time Problem Ken Jacson Alexander Kreinin Wanhe Zhang February 1, 29 Abstract In this paper we consider the following inverse problem for the first hitting time distribution:
More informationImportance Sampling Methodology for Multidimensional Heavy-tailed Random Walks
Importance Sampling Methodology for Multidimensional Heavy-tailed Random Walks Jose Blanchet (joint work with Jingchen Liu) Columbia IEOR Department RESIM Blanchet (Columbia) IS for Heavy-tailed Walks
More informationOctavio Arizmendi, Ole E. Barndorff-Nielsen and Víctor Pérez-Abreu
ON FREE AND CLASSICAL TYPE G DISTRIBUTIONS Octavio Arizmendi, Ole E. Barndorff-Nielsen and Víctor Pérez-Abreu Comunicación del CIMAT No I-9-4/3-4-29 (PE /CIMAT) On Free and Classical Type G Distributions
More informationLimit theorems for multipower variation in the presence of jumps
Limit theorems for multipower variation in the presence of jumps Ole E. Barndorff-Nielsen Department of Mathematical Sciences, University of Aarhus, Ny Munkegade, DK-8 Aarhus C, Denmark oebn@imf.au.dk
More informationMonte Carlo Methods. Handbook of. University ofqueensland. Thomas Taimre. Zdravko I. Botev. Dirk P. Kroese. Universite de Montreal
Handbook of Monte Carlo Methods Dirk P. Kroese University ofqueensland Thomas Taimre University ofqueensland Zdravko I. Botev Universite de Montreal A JOHN WILEY & SONS, INC., PUBLICATION Preface Acknowledgments
More informationI forgot to mention last time: in the Ito formula for two standard processes, putting
I forgot to mention last time: in the Ito formula for two standard processes, putting dx t = a t dt + b t db t dy t = α t dt + β t db t, and taking f(x, y = xy, one has f x = y, f y = x, and f xx = f yy
More informationBisection Ideas in End-Point Conditioned Markov Process Simulation
Bisection Ideas in End-Point Conditioned Markov Process Simulation Søren Asmussen and Asger Hobolth Department of Mathematical Sciences, Aarhus University Ny Munkegade, 8000 Aarhus C, Denmark {asmus,asger}@imf.au.dk
More informationFeller Processes and Semigroups
Stat25B: Probability Theory (Spring 23) Lecture: 27 Feller Processes and Semigroups Lecturer: Rui Dong Scribe: Rui Dong ruidong@stat.berkeley.edu For convenience, we can have a look at the list of materials
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationStationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos
N Cufaro Petroni: CYCLOTRONS 2007 Giardini Naxos, 1 5 October, 2007 1 Stationary distributions of non Gaussian Ornstein Uhlenbeck processes for beam halos CYCLOTRONS 2007 Giardini Naxos, 1 5 October Nicola
More informationX n D X lim n F n (x) = F (x) for all x C F. lim n F n(u) = F (u) for all u C F. (2)
14:17 11/16/2 TOPIC. Convergence in distribution and related notions. This section studies the notion of the so-called convergence in distribution of real random variables. This is the kind of convergence
More informationSome SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen
Title Author(s) Some SDEs with distributional drift Part I : General calculus Flandoli, Franco; Russo, Francesco; Wolf, Jochen Citation Osaka Journal of Mathematics. 4() P.493-P.54 Issue Date 3-6 Text
More informationOn optimal stopping of autoregressive sequences
On optimal stopping of autoregressive sequences Sören Christensen (joint work with A. Irle, A. Novikov) Mathematisches Seminar, CAU Kiel September 24, 2010 Sören Christensen (CAU Kiel) Optimal stopping
More informationOther properties of M M 1
Other properties of M M 1 Přemysl Bejda premyslbejda@gmail.com 2012 Contents 1 Reflected Lévy Process 2 Time dependent properties of M M 1 3 Waiting times and queue disciplines in M M 1 Contents 1 Reflected
More informationStability of the Defect Renewal Volterra Integral Equations
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 211 http://mssanz.org.au/modsim211 Stability of the Defect Renewal Volterra Integral Equations R. S. Anderssen,
More informationTesting a Normal Covariance Matrix for Small Samples with Monotone Missing Data
Applied Mathematical Sciences, Vol 3, 009, no 54, 695-70 Testing a Normal Covariance Matrix for Small Samples with Monotone Missing Data Evelina Veleva Rousse University A Kanchev Department of Numerical
More informationAdditive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535
Additive functionals of infinite-variance moving averages Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Departments of Statistics The University of Chicago Chicago, Illinois 60637 June
More informationarxiv: v2 [math.pr] 4 Sep 2017
arxiv:1708.08576v2 [math.pr] 4 Sep 2017 On the Speed of an Excited Asymmetric Random Walk Mike Cinkoske, Joe Jackson, Claire Plunkett September 5, 2017 Abstract An excited random walk is a non-markovian
More information