arxiv: v2 [math.pr] 15 Jul 2015

Transcription

1 STRIKINGLY SIMPLE IDENTITIES RELATING EXIT PROBLEMS FOR LÉVY PROCESSES UNDER CONTINUOUS AND POISSON OBSERVATIONS arxiv: v2 [math.pr] 15 Jul 215 HANSJÖRG ALBRECHER AND JEVGENIJS IVANOVS Abstract. We consider exit problems for general Lévy processes, where the first passage over a threshold is detected either immediately or at an epoch of an independent homogeneous Poisson process. It is shown that the two corresponding one-sided problems are related through a surprisingly simple identity. Moreover, we identify a simple link between two-sided exit problems with one continuous and one Poisson exit. Finally, Poisson exit of a reflected process is connected to the continuous exit of a process reflected at Poisson epochs, and a link between some Parisian type exit problems is established. With the appropriate perspective, the proofs of all these relations turn out to be quite elementary. For spectrally one-sided Lévy processes this approach enables alternative proofs for a number of previously established identities, providing additional insight. 1. Introduction Let X = X t, t be a real-valued Lévy process, and let T i, i 1 be the epochs of an independent Poisson process with intensity λ > ; add T =. The probability law corresponding to X started at u will be denoted by P u with denoting the expectation. When u is not mentioned explicitly we assume that u = and write simply P and E. Define τ = inf{t : X t < }, τ + a = inf{t : X t > a}, τ = min{t i, i N : X Ti < }, τ + a = min{t i, i N : X Ti > a}, which we interpret as the first passage times under continuous and Poisson observations, respectively. Observe that τ < τ and, moreover, τ converges in probability to τ as λ the same is true for τ a + and τ a +. Thus exit theory under Poisson observation can be regarded as a generalization of the classical exit theory. Throughout this paper, however, we keep λ > fixed. Observation at Poisson epochs is both of theoretical and practical interest. Firstly, some exit problems with Poisson observation yield transforms of certain occupation times, e.g. P u τ < τ + a u [exp λ τ 1 {Xt>a}dt ; τ < ]. 21 Mathematics Subject Classification. Primary 6G51; Secondary 91B3. Key words and phrases. Lévy processes, exit problems, Poisson observation, occupation times, Parisian ruin. Financial support by the Swiss National Science Foundation Project is gratefully acknowledged. 1

2 2 H. ALBRECHER AND J. IVANOVS Secondly, Poisson observation is relevant in various applications such as queueing see e.g. [5], reliability and insurance risk theory see e.g. [1, 2]. In particular, in many applications discrete-time observation of stochastic processes would often be considered more natural, but for equidistant discrete time epochs the explicit and tractable analytical structure of continuous-time processes is typically destroyed, so that one is forced towards numerical techniques for the determination of exit probabilities and related quantities. The Poisson observation structure is a bridge between continuous-time and discrete-time observation that still leads to rather explicit, and as will be shown below, also somewhat elegant modifications of the continuous-time formulas Overview and organization. In order to stress the intuition behind the derivation of the identities, we will start with a simple case and gradually generalize the setup. Most of the results are stated in terms of relations between transforms, but can also be understood as relations between the corresponding laws in an obvious way. Some of the wording throughout the manuscript will be in terms of the insurance application, where X is the surplus process of a portfolio of insurance contracts, τ is the time of ruin of the portfolio, {τ = } is the event of infinite-time survival, and τ is the time of observed ruin under Poisson observation of the surplus process in the application the Poisson epochs can for instance be interpreted as the observation times of the regulatory authority. In Section 2 we discuss survival probabilities corresponding to the two observation types, and then proceed to the general one-sided exit problems including the time of exit and the overshoot. In Section 3 we consider more complex problems. Firstly, the two-sided exit problem with one continously observed and one Poisson-discretely observed boundary is related to the one where the observation types at the boundaries are interchanged. Secondly, we provide a link between Poisson exit of a reflected process and continuous exit of the process reflected at Poisson epochs. We also show that a two-sided problem with Poisson exit at both boundaries yields an identity as well, but with a non-standard first passage time. The latter quantity is then linked to a Parisian ruin problem with Erlang-distributed implementation delay. Finally, we establish a link between Parisian ruin problems with continuous and Poisson observations. We conclude with Section 4, where we specialize to the case of spectrally-one sided processes and demonstrate the use of our simple identities, providing simpler proofs and additional insight to some identities established in earlier literature Preliminaries. The Wiener-Hopf factorization plays a crucial role in the derivations below. Define X t = inf{x s, s [, t]}, G t = inf{s [, t] : X s X s = X t }, the infimum and its first time of occurrence up to horizon t. Similarly, the supremum and its last time of occurrence are defined by X t = sup{x s, s [, t]}, G t = sup{s [, t] : X s X s = X t }. Finally, let the pairs D, T D and U, T U be distributed as X T1, G T1 and X T1, G T1 respectively under P, and sampled independently of each other and of everything

3 EXIT PROBLEMS UNDER CONTINUOUS AND POISSON OBSERVATIONS 3 else D and U stand for down and up. Recall that according to the Wiener-Hopf factorization we have X T1, G T1, X T1 X T1, T 1 G T1 d = D, T D, U, T U d = X T1 X T1, T 1 G T1, X T1, G T1, see, e.g., [6, Thm. VI.5], and [7] for applications of factorization embeddings. 2. One-sided exit 2.1. Survival probability. Let us first consider 1 φu := P u τ =, φu := Pu τ =, which in the insurance application are the probabilities of survival with initial capital u under continuous and Poisson observation, respectively. In fact, the two quantities are connected by two very simple relations: Proposition 1. For u it holds that 2 3 φu φu + U, φu φu + D. u D 1 U 1 T 1 T 2 Figure 1. Schematic sample path and embeddings Proof. Survival under Poisson observation is determined by the sequence u + X Ti, whereas survival under continuous observation is determined by the sequence of infima in between the observation epochs black and grey dots in Figure 1, respectively. Let D 1 = X T1, U 1 = X T1 X T1 and define D i+1, U i+1 i 1 in the same way but for the shifted process X Ti+t X Ti and exponential time T i+1 T i. Let Ŝi, i N and S i, i N be the partial sum processes corresponding to, D 1 + U 1, D 2 + U 2, D 3 + U 3,... and D 1, U 1 + D 2, U 2 + D 3,..., respectively; u+ŝi and u+s i are the heights of the black and grey dots in Figure 1. Observe that all U i and D i are independent, because of independence of increments and the Wiener-Hopf factorization. Since the D i s have the law of D we obtain 4 Ŝi + D, i d = S i, i. Similarly, 5 Ŝi, i 1 d = S i + U, i.

4 4 H. ALBRECHER AND J. IVANOVS Hence φu = Pu + min S i = Pu + D + min Ŝ i φu + D, i i and, since u, φu = Pu + min i Ŝ i = Pu + min Ŝ i i 1 = Pu + U + min i S i φu + U. Remark 1. Relation 2 allows to interpret the transition from continuous to discrete Poisson observation simply as a random increase of the starting value initial capital u by U, as far as the survival probability is concerned; that is the structure of φ as a function of u is otherwise completely preserved. Likewise, Relation 3 shows that moving from discrete Poisson to continuous observation preserves the structure, reducing the initial capital by D which has all its probability mass on the negative half-line. Remark 2. Suppose we modify the Poisson observation model, so that there is no observation at time. Then 2 is still valid even for negative u then, whereas 3 does not hold any more. Remark 3. By the same token one can connect the finite-time survival probabilities φu, T i := P u τ > T i = Pu + min j i 1 S j and φu, T i := P u τ > T i = Pu + min j i Ŝ j for i N: φu, T i φu + U, T i, φu, T i φu + D, T i 1. That is, survival under continuous observation up to an independent Erlang distributed time horizon is intimately related to survival under Poisson observation up to a certain arrival epoch The general identities. Proposition 2. For u and α, β it holds that e α τ 6 +βx τ ; τ < ] [e αt U +U e ατ +βx τ ; τ < Ee αt D +βd, e ατ +βx τ ; τ < Ee αt D +βd 7 ] [e αt D +D e α τ +βx τ ; τ <. Observe that the left-hand side of 7 gives the transform of the undershoot of the first grey point below in Figure 1, according to the strong Markov property applied at τ. Now one can establish the relation between black and grey points as in the proof of Proposition 1, additionally taking time into account. Essentially, we just shift the picture so that we start at the first grey point.

5 EXIT PROBLEMS UNDER CONTINUOUS AND POISSON OBSERVATIONS 5 Proof of Proposition 2. Let T1 D = G T1, T1 U = T 1 G T1 and define Ti+1 D, T i+1 U in the same way but for the shifted process X Ti+t X Ti and exponential time T i+1 T i. As in the proof of Proposition 1 we consider the sequences of black and grey dots in Figure 1, but now we also add the time component: Ŝi, T i, i N and S i, G i, i N which are the partial sum processes corresponding to,, D 1 + U 1, T D 1 + T U 1, D 2 + U 2, T D 2 + T U 2,... and D 1, T D 1, U 1 + D 2, T U 1 + T D 2,..., respectively. Similarly to 4 and 5 we observe that Ŝi + D, T i + T D, i d = S i, G i, i Ŝi, T i, i 1 d = S i + U, G i + T U, i. Letting N u = min{i : u + Ŝi < } = min{i 1 : u + Ŝi < } and N u = min{i : u + S i < } be the first passage epochs we can write e α τ +βx τ ; τ < e αt +βu+ŝ Nu Nu ; N u < e αt U +G Nu+U +βu+u+s Nu+U ; N u+u < +U e αt U +τ +βx τ ; τ < Ee αt D +βd, where in the last line we applied the strong Markov property at τ. Identity 7 can be derived analogously. 3. Further exit problems 3.1. Two-sided exit with different observation types. In this section we consider two-sided exit problems with one continuous and one Poisson exit at the boundaries. It turns out that there is a simple relation between the problems when the roles of the continuous and the Poisson exit are interchanged, i.e. problems corresponding to {τ < τ + a } and { τ < τ + a }. Here we extend the ideas of Section 2 to their full potential. Proposition 3. For a u and α, β it holds that e α τ 8 +βx τ ; τ < τ a + ] [e αt U +U e ατ +βx τ ; τ < τ a + Ee αt D +βd, e ατ +βx τ ; τ < τ a + Ee αt D +βd 9 ] [e αt D +D e α τ +βx τ ; τ < τ a +. Proof. The proof is by inspection: For 8, consider the embeddings illustrated in Figure 2. In the left picture the grey dots correspond to the observations and the black to the suprema in between two observations. In the right picture the black dots correspond to observations and the grey to the infima in between observations. Note that the position of a black point with respect to the previous grey point

6 6 H. ALBRECHER AND J. IVANOVS has the same distribution in both cases, namely U, T U. The same is true for the position of the grey points with respect to their previous black points with D, T D. So the patterns of points in each case have the same law up to a certain shifting; we illustrate this by using the same patterns of points in both pictures in Figure 2 and by drawing different sample paths. Now it follows that P u τ < τ a +, see the left picture, must coincide with EP u+u τ < τ a +, see the right picture, because the interpretation of points was reversed. Finally, we include the value of X at first passage and its time using the strong Markov property at τ as in the proof of Proposition 2. The same type of reasoning yields 9. Figure 2. Poisson observations are in grey left picture and in black right picture By considering the negative of X we immediately obtain the following result from Proposition 3. Corollary 1. For a u and α, β it holds that e α τ + a 1 βx τ a + a ; τ a + < τ [e αt D +D e ατ + a βx τ a + a ; τ a + < τ e ατ + a βx τ a + a ; τ a + < τ Ee αt U βu 11 [e αt U +U e α τ + a βx τ a + a ; τ a + < τ ] Ee αt U βu, 3.2. Reflected processes. In this section we consider the process X reflected at a barrier a > in a continuous and Poisson manner, and study its first passage below in Poisson and continuous manner respectively with always opposite manners. Again, these two problems are closely related. Note that in an insurance context reflection at a results when paying out dividends according to a barrier strategy, either continuously or at Poisson epochs see e.g. [1]. Let P a u be the law of X started in u and continuously reflected at a and let R be the corresponding regulator, i.e. X t, R t, t under P a u is ]. X t X t a +, X t a +, t under P u. Similarly, let P a u be the law of X started at u and reflected in Poisson manner at a, i.e. X t, R t, t under P a u is X t max{x Ti : T i t} a +, max{x Ti : T i t} a +, t under P u.

7 EXIT PROBLEMS UNDER CONTINUOUS AND POISSON OBSERVATIONS 7 Proposition 4. For a u and α, β, γ it holds that E a u e α τ +βx τ γr τ ; τ < [ ] e αt U Ê a u+u e ατ +βx τ γr τ ; τ < Ee αt D +βd, Ê a u e ατ +βx τ γr τ ; τ < Ee αt D +βd [ ] e αt D E a u+d e α τ +βx τ γr τ ; τ <. Proof. Again, the proof follows merely by inspection in a similar way as for the previous results. The first relation can be seen from Figure 3, where the left picture depicts continuous reflection at a and Poisson observation at, and the right picture depicts the corresponding shifted Poisson reflection at a and continuous observation at. In the left picture Poisson observations yield the sequence: u, u + U 1 a + D 1, u + U 1 a + D 1 + U 2 a + D 2,.... In the right picture the infima in between Poisson reflection epochs are given by û a + D 1, û a + D 1 + U 1 a + D 2,..., where we choose û to be distributed as u + U. These sequences can be complemented with the respective times as in the proof of Proposition 2. Finally, it is easy to see that R enters the transforms without requiring any changes. The second relation follows accordingly. U 1 D 1 U 2 U 1 D 2 D 1 D 2 Figure 3. Continuous reflection and Poisson exit left, and Poisson reflection and continuous exit right Remark 4. It is easy to see that Proposition 4 can be generalized from reflection to so-called refraction. Concretely, consider the processes X t δx t a + and X t δmax{x Ti : T i t} a + for δ [, 1]. In the insurance context such a refraction has the interpretation of taxation according to a loss-carry-forward scheme and tax rate δ, see e.g. [4] Two-sided Poisson exit. The two-sided exit with Poisson observation at both barriers can be related to a model with another type of exit time. Define the random time τ a + of the first observation such that the process has stayed above a during the entire preceding inter-observation period, i.e. τ a + = TÑ + a with Ñ a + := min{i N : inf{x t, t [T i 1, T i ]} > a}. Similarly, define τ = T Ñ with Ñ := min{i N : sup{x t, t [T i 1, T i ]} < } as the first observation time such that the

8 8 H. ALBRECHER AND J. IVANOVS process has stayed below during the entire preceding inter-observation period. Then for u [, a] it holds that 12 P u τ < τ + a P u+u τ < τ + a P u+d τ < τ + a. To see this, one follows the same ideas as above: for the first equality consider infima in between two observations, see Figure 2, and for the second equality consider suprema in between two observations. Similarly, we also have the reverse identities: 13 P u τ < τ + a P u+d τ < τ + a, P u τ < τ + a P u+u τ < τ + a Parisian ruin. Parisian ruin is defined as the first time when an excursion of X below is longer than some time V sometimes referred to as implementation delay. Whereas the classical definition is in terms of a deterministic V, for analytic tractability it is often assumed that V is a random variable, and that an independent copy of V is assigned to each excursion, see e.g. [11] and [1]. Firstly, from the memoryless property it follows that the time τ of Poisson ruin is also the time of Parisian ruin in the case where V is an exponential random variable with rate λ. Secondly, τ as defined in Section 3.3 is the time of Parisian ruin in the case where V is Erlang2, λ distributed since the latter is the sum of two independent exponential variables. Similarly to 6, Equation 13 can easily be extended to 14 e α τ +βx τ ; τ < [e αt U +U e α τ +βx τ ; τ < ] Ee αt D +βd, which under the present interpretation relates Parisian ruin quantities with exponential and Erlang2-distributed implementation delay here we took a = for simplicity. More generally, consider Parisian ruin with Erlangk, λ implementation delay and let τ k denote the corresponding ruin time. So in particular τ = τ, τ 1 = τ 2 and τ = τ k. On the other hand, define τ as the first epoch T i, i k such that XT i k,..., XT i <, i.e. the process has been observed negative at the last k + 1 Poisson epochs. Then, along the same line of arguments, we can extend 14 cf. Figure 4: τ 3 τ 1 Figure 4. Poisson observations are in grey left picture and in black right picture

9 EXIT PROBLEMS UNDER CONTINUOUS AND POISSON OBSERVATIONS 9 Proposition 5. For k 2 and u we have k ατ +βx k 15 e τ ; τ k < ] [e αt U k 2 α τ k 2 +βx τ +U e ; τ k 2 < Ee αt D +βd. and k 2 α τ k 2 +βx τ 16 e ; τ k 2 < Ee αt D +βd ] [e αt D k ατ +βx k +D e τ ; τ k <. 4. The case of spectrally one-sided Lévy processes If X is a one-sided Lévy process, some of the identities lead to more explicit forms, and this will allow to retrieve a number of results previously obtained in the literature, now with alternative proofs, revealing some more structure of the formulas. Without loss of generality assume that X is a spectrally-negative Lévy process, i.e. it may only have negative jumps and it is not a non-increasing process. Consider its Laplace exponent and put ψ α θ = ψθ α for α. ψθ = log Ee θx1, θ, 4.1. Preliminaries. Let us first recall some basic functions which play a fundamental role in exit theory, see e.g. [8, Ch. 8]. Let Φ α be the largest non-negative zero of ψ α, and let W α u, u be the so-called scale function: a continuous non-negative function determined by its Laplace transform e θu W α udu = 1/ψ α θ, θ > Φ α. In addition, we need a second scale function u Z α u, θ = e 1 θu ψ α θ e θy W α ydy, u, which can be rewritten as 17 Z α u, θ = ψ α θ e θy W α u + ydy for θ > Φ α, see also [3]. The two basic one-sided exit identities under continuous observation are e ατ + a = e Φ αa u, u a e ατ +βxτ, τ < = Z α u, β W α u ψ αβ β Φ α, u, and the Wiener-Hopf factors are given by Ee αt U βu = Φ λ Φ λ+α + β, D Ee αt +βd = λ Φ λ+α β Φ λ λ ψ α β, see e.g. [8, Ch. 8]. In order to apply Formula 6 of Proposition 2 to 19, we first need the following identities:

10 1 H. ALBRECHER AND J. IVANOVS Lemma 1. For u it holds that E e αt U W α u + U = Φ λ λ Z αu, Φ λ+α, E e αt U Z α u + U, β = Proof. Firstly, Φ λ Φ λ+α β 2 Ee αt U ; U dy = Φ λ e Φ λ+αy dy, Z α u, β ψ αβ λ Z αu, Φ λ+α which can be checked by taking transforms and comparing to the Wiener-Hopf factor. Hence E e αt U W α u + U = W α u + yφ λ e Φλ+αy Φ λ dy = ψ α Φ λ+α Z αu, Φ λ+α according to 17. For the first identity it is left to note that ψ α Φ λ+α = λ. Using 17 several times we can write for µ > β > Φ α : = = Z α u + x, βe µx dx = ψ α β. e µ βx e βy+x W α u + y + xdydx e µ βx ψ α β e βy W α u + ydydx x x e Z µ βx α u, β ψ α β e βy W α u + ydy dx 1 = Z α u, β µ β ψ 1 αβ µ β = 1 Z α u, β ψ αβ µ β ψ α µ Z αu, µ e µy W α u + ydy. Plugging in µ = Φ λ+α and multiplying by Φ λ we obtain the second identity. By analytic extension β can be chosen arbitrarily One-sided exit. Assume that EX 1 = ψ > and consider the survival probabilities 1. It is well known that φu = ψ W u, see also 19. According to Proposition 1 and Lemma 1 we have φu φu + U = ψ EW u + U = ψ Φ λ λ Z u, Φ λ, which is Corollary 1 of [9]. Remark 5. Note that due to 2, for the spectrally negative Lévy process the identity 2 simplifies to the pleasant form φu φu + e Φλ, where e Φλ is an exponential random variable with parameter Φ λ. This for instance immediately explains why for a compound Poisson process X with exponential jump sizes the discrete Poisson observation changes the classical ruin probability formula just by a multiplicative factor Φ λ /Φ λ + R, where R is the Lundberg adjustment coefficient cf. [2, Eq.2.18].

11 EXIT PROBLEMS UNDER CONTINUOUS AND POISSON OBSERVATIONS Using the standard identity 19, Proposition 2 and Lemma 1 we obtain e α τ +βx τ, τ < e αt U Z α u + U, β e αt U W α u + U ψ αβ β Φ α = Φ λ Φ λ+α β λ Φ λ Φ λ+α β λ ψ α β = Z α u, β ψ αβ λ λ ψ α β Ee αt D +βd Φ λ λ Z αu, Φ λ+α ψ αβ β Φ α λ Z αu, Φ λ+α Z α u, β Z α u, Φ λ+α ψ αβ λ Φ λ+α Φ α β Φ α Also, by considering Proposition 2 for the negative of X, see also Corollary 1, we arrive at e α τ + a βx τ a + a ; τ a + < ] [e αt D +D e ατ + a βx τ a + a ; τ a + < Ee αt U βu e αt D Φ αa u D Ee αt U βu = e Φαa u λ Φ λ+α Φ α Φ λ λ ψ α Φ α Φ λ+α + β = e Φαa u Φ λ+α Φ α Φ λ+α + β, because of 18 and the fact that X τ + a = a on τ a +. These two identities were obtained in [3, Thm. 3.1] by virtue of a rather technical argument using the expression for the potential density of X Parisian ruin. Finally, we relate our results to previous literature on Parisian ruin. Firstly, taking u = and β = in 21 use Z α, β = 1 one retrieves Corollary 3.2 of [1], which is based on exponential implementation delay. Furthermore, from the form of [1, Eq.49] one can, after some lengthy calculations, obtain the following expression for Erlang2, λ implementation delay: 2 ατ 22 E e ; τ 2 < = λ λ + α α λ + α 2 Φ λ+α Φ α Φ α Φ λ Φ λ+α Φ, λ+α where the derivative is with respect to the subindex. We can alternatively obtain 22 directly using the results of this paper: 15 and 21 imply ] 2 ατ E e ; τ 2 < [e αt U E U e α τ ; τ < Ee αt D = λ [e λ + α E αt U Z α U, Z α U, Φ λ+α α ] Φ λ+α Φ α Ee αt D. λ Φ α Using Lemma 1 and noting again that Z α, β = 1 we get E e αt U Z α U, = Φ λ α + λ Φ λ+α λ, E e αt U Z α U, Φ λ+α = lim h Φ λ Φ λ+α Φ λ+α h h λ = Φ λ λφ λ+α,. which together with the expression for the Wiener-Hopf factor Ee αt D yields 22. readily

12 12 H. ALBRECHER AND J. IVANOVS Acknowledgements The authors would like to thank Ton Dieker for stimulating discussions on the topic. References [1] H. Albrecher, E. C. K Cheung, and S. Thonhauser. Randomized observation periods for the compound poisson risk model: Dividends. Astin Bull., 412: , 211. [2] H. Albrecher, E. C. K. Cheung, and S. Thonhauser. Randomized observation periods for the compound Poisson risk model: the discounted penalty function. Scand. Actuar. J., 6: , 213. [3] H. Albrecher, J. Ivanovs, and X. Zhou. Exit identities for lévy processes observed at poisson arrival times. Bernoulli, 215. in press arxiv: [4] H. Albrecher, J.-F. Renaud, and X. Zhou. A Lévy insurance risk process with tax. J. Appl. Probab., 452: , 28. [5] R. Bekker, O. J. Boxma, and J. A. C. Resing. Lévy processes with adaptable exponent. Adv. in Appl. Probab., 411:177 25, 29. [6] J. Bertoin. Lévy processes, volume 121. Cambridge University Press, [7] A. B. Dieker. Applications of factorization embeddings for Lévy processes. Adv. in Appl. Probab., 383: , 26. [8] A. Kyprianou. Introductory lectures on fluctuations of Lévy processes with applications. Universitext. Springer-Verlag, Berlin, 26. [9] D. Landriault, J.-F. Renaud, and X. Zhou. Occupation times of spectrally negative Lévy processes with applications. Stochastic Process. Appl., 12111: , 211. [1] D. Landriault, J.-F. Renaud, and X. Zhou. An insurance risk model with Parisian implementation delays. Methodol. Comput. Appl. Probab., 163:583 67, 214. [11] R. Loeffen, I. Czarna, and Z. Palmowski. Parisian ruin probability for spectrally negative lévy processes. Bernoulli, 192:599 69, address: hansjoerg.albrecher@unil.ch University of Lausanne and Swiss Finance Institute address: jevgenijs.ivanovs@unil.ch University of Lausanne