AN INTRODUCTION to LÉVY PROCESSES Contents 1. Introduction 2 2. Definition of Lévy Processes 2 3. Toy Example: a Lévy Jump-Di

Total Page:16

File Type:pdf, Size:1020Kb

Load more

AN INTRODUCTION TO LEVY´ PROCESSES ERIK BAURDOUX AND ANTONIS PAPAPANTOLEON Abstract. Lecture notes from courses at TU Berlin, LSE and CIMAT Contents 1. Introduction 2 2. Definition of L´evy processes 2 3. Toy example: a L´evyjump-diffusion 3 4. Infinitely Divisible distributions 6 5. The L´evy{Khintchine representation 11 6. The L´evy{It^odecomposition 19 7. The L´evymeasure and path properties 32 8. Elementary operations 38 9. Moments and Martingales 41 10. Popular examples 46 11. Simulation of L´evyprocesses 52 12. Stochastic integration 53 References 60 Date: February 6, 2015. 1 2 ERIK BAURDOUX AND ANTONIS PAPAPANTOLEON 1. Introduction These lecture notes aim at providing a (mostly) self-contained introduc- tion to L´evyprocesses. We start by defining L´evyprocesses in Section 2 and study the simple but very important example of a L´evyjump-diffusion in Section 3. In Section 4 we discuss infinitely divisible distributions, as the distribution of a L´evyprocess at any fixed time has this property. In the sub- sequent section we prove the L´evy{Khintchine formula, which characterises all infinitely divisible distributions. We then prove the L´evy{It^odecomposi- tion of a L´evyprocess, which is an explicit existence result explaining how to construct a L´evyprocess based on a given infinitely divisible distribution. In Section 7, we make use of the L´evy{It^odecomposition to derive necesarry and sufficient conditions for a L´evyprocess to be of (in)finite variation and a subordinator, respectively. In Section 8 we study elementary operations such as linear transformations, projections and subordination. We then move on to moments and martingales followed by a brief section on simulation fol- lowed by a number of popular models used mostly in mathematical finance. Then there is a section on simulation and we finish with a brief treatment of stochastic integration. 2. Definition of Levy´ processes 2.1. Notation and auxiliary definitions. Let (Ω; F; F; IP) denote a sto- chastic basis, or filtered probability space, i.e. a probability space (Ω; F; IP) endowed with a filtration F = (Ft)t≥0. A filtration is an increasing family of sub-σ-algebras of F, i.e. Fs ⊂ Ft W for s ≤ t. By convention F1 = F and F∞− = s≥0 Fs. A stochastic basis satisfies the usual conditions if it is right-continuous, T i.e. Ft = Ft+, where Ft+ = s>t Fs, and is complete, i.e. the σ-algebra F is IP-complete and every Ft contains all IP-null sets of F. Definition 2.1. A stochastic process X = (Xt)t≥0 has independent incre- ments if, for any n ≥ 1 and 0 ≤ t0 ≤ t1 ≤ · · · ≤ tn, the random variables Xt0 ;Xt1 − Xt0 ;:::;Xtn − Xtn−1 are independent. Alternatively, we say that X has independent increments if, for any 0 ≤ s < t, Xt − Xs is independent of Fs. Definition 2.2. A stochastic process X = (Xt)t≥0 has stationary incre- ments if, for any s; t ≥ 0, the distribution of Xt+s − Xs does not depend on s. Alternatively, we say that X has stationary increments if, for any 0 ≤ s ≤ t, Xt − Xs is equal in distribution to Xt−s. Definition 2.3. A stochastic process X = (Xt)t≥0 is stochastically contin- uous if, for every t ≥ 0 and > 0 lim IP(jXs − Xtj > ) = 0: s!t 2.2. Definition of L´evyprocesses. We will now define L´evyprocesses and then present some well-known examples, like the Brownian motion and the Poisson process. d Definition 2.4 (L´evyprocess). An adapted, R -valued stochastic process X = (Xt)t≥0 with X0 = 0 a.s. is called a L´evyprocess if: AN INTRODUCTION TO LEVY´ PROCESSES 3 (L1) X has independent increments, (L2) X has stationary increments, (L3) X is stochastically continuous. In the sequel, we will always assume that X has c`adl`agpaths. The next two results provide the justification. Lemma 2.5. If X is a L´evyprocess and Y is a modification of X (i.e. IP(Xt 6= Yt) = 0 a.s. for each t ≥ 0), then Y is a L´evyprocess and has the same characteristics as X. Proof. [App09, Lemma 1.4.8]. Theorem 2.6. Every L´evyprocess has a unique c`adl`agmodification that is itself a L´evyprocess. Proof. [App09, Theorem 2.1.8] or [Pro04, Theorem I.30]. 2.3. Examples. The following are some well-known examples of L´evypro- cesses: • The linear drift is the simplest L´evyprocess, a deterministic process; see Figure 5.3 for a sample path. • The Brownian motion is the only non-deterministic L´evyprocess with continuous sample paths; see Figure 5.3 for a sample path. • The Poisson, the compound Poisson and the compensated (com- pound) Poisson processes are also examples of L´evyprocesses; see Figure 5.3 for a sample path of a compound Poisson process. The sum of a linear drift, a Brownian motion and a (compound or com- pensated) Poisson process is again a L´evyprocess. It is often called a \jump- diffusion” process. We shall call it a L´evyjump-diffusion process, since there exist jump-diffusion processes which are not L´evyprocesses. See Figure 5.3 for a sample path of a L´evyjump-diffusion process. 3. Toy example: a Levy´ jump-diffusion Let us study the L´evyjump-diffusion process more closely; it is the sim- plest L´evyprocess we have encountered so far that contains both a diffusive part and a jump part. We will calculate the characteristic function of the L´evyjump-diffusion, since it offers significant insight into the structure of the characteristic function of general L´evyprocesses. Assume that the process X = (Xt)t≥0 is a L´evyjump-diffusion, i.e. a linear deterministic process, plus a Brownian motion, plus a compensated compound Poisson process. The paths of this process are described by N Xt Xt = bt + σWt + Jk − tλβ ; (3.1) k=1 where b 2 R, σ 2 R>0, W = (Wt)t≥0 is a standard Brownian motion, N = (Nt)t≥0 is a Poisson process with intensity λ 2 R>0 (i.e. IE[Nt] = λt), and J = (Jk)k≥1 is an i.i.d. sequence of random variables with probability distribution F and IE[Jk] = β < 1. Here F describes the distribution of the jumps, which arrive according to the Poisson process N. All sources of randomness are assumed mutually independent. 4 ERIK BAURDOUX AND ANTONIS PAPAPANTOLEON 0.2 0.05 0.04 0.1 0.03 0.0 0.02 0.01 −0.1 0.0 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 0.4 0.8 0.3 0.6 0.4 0.2 0.2 0.1 0.0 −0.2 0.0 −0.4 −0.1 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 Figure 2.1. Sample paths of a linear drift processs (top- left), a Brownian motion (top-right), a compound Poisson process (bottom-left) and a L´evyjump-diffusion. The characteristic function of Xt, taking into account that all sources of randomness are independent, is h Nt i iuXt X IE e = IE exp iu bt + σWt + Jk − tλβ k=1 Nt h i h X i = exp iubt IE exp iuσWt IE exp iu Jk − iutλβ ; k=1 recalling that the characteristic functions of the normal and the compound Poisson distributions are 1 2 2 iuσWt − σ u t IE[e ] = e 2 ;Wt ∼ N (0; t) PNt iuJ iu Jk λt(IE[e k −1]) IE[e k=1 ] = e ;Nt ∼ Poi(λt) (cf. Example 4.14 and Exercise 1), we get h 1 i h i = exp iubt exp − u2σ2t exp λtIE[eiuJk − 1] − iuIE[J ] 2 k h 1 i h i = exp iubt exp − u2σ2t exp λtIE[eiuJk − 1 − iuJ ] ; 2 k AN INTRODUCTION TO LEVY´ PROCESSES 5 and since the distribution of Jk is F we have h 1 i h Z i = exp iubt exp − u2σ2t exp λt eiux − 1 − iuxF (dx) : 2 R Finally, since t is a common factor, we can rewrite the above equation as u2σ2 Z IEeiuXt = exp t iub − + (eiux − 1 − iux)λF (dx) : (3.2) 2 R We can make the following observations based on the structure of the char- acteristic function of the random variable Xt from the L´evyjump-diffusion: (O1) time and space factorize; (O2) the drift, the diffusion and the jump parts are separated; (O3) the jump part decomposes to λ × F , where λ is the expected number of jumps and F is the distribution of jump size. One would naturally ask if these observations are true for any L´evyprocess. The answer for (O1) and (O2) is yes, because L´evy processes have stationary and independent increments. The answer for (O3) is no, because there exist L´evyprocesses with infinitely many jumps (on any compact time interval), thus their expected number of jumps is also infinite. Since the characteristic function of a random variable determines its dis- tribution, (3.2) provides a characterization of the distribution of the random variables Xt from the L´evyjump-diffusion X. We will soon see that this dis- tribution belongs to the class of infinitely divisible distributions and that equation (3.2) is a special case of the celebrated L´evy-Khintchineformula. 3.1. The basic connections. The next sections will be devoted to estab- lishing the connection between the following mathematical objects: • L´evyprocesses X = (Xt)t≥0 • infinitely divisible distributions ρ = L(X1) • L´evytriplets (b; c; ν).
Recommended publications
  • Levy Processes

    Levy Processes

    LÉVY PROCESSES, STABLE PROCESSES, AND SUBORDINATORS STEVEN P.LALLEY 1. DEFINITIONSAND EXAMPLES d Definition 1.1. A continuous–time process Xt = X(t ) t 0 with values in R (or, more generally, in an abelian topological groupG ) isf called a Lévyg ≥ process if (1) its sample paths are right-continuous and have left limits at every time point t , and (2) it has stationary, independent increments, that is: (a) For all 0 = t0 < t1 < < tk , the increments X(ti ) X(ti 1) are independent. − (b) For all 0 s t the··· random variables X(t ) X−(s ) and X(t s ) X(0) have the same distribution.≤ ≤ − − − The default initial condition is X0 = 0. A subordinator is a real-valued Lévy process with nondecreasing sample paths. A stable process is a real-valued Lévy process Xt t 0 with ≥ initial value X0 = 0 that satisfies the self-similarity property f g 1/α (1.1) Xt =t =D X1 t > 0. 8 The parameter α is called the exponent of the process. Example 1.1. The most fundamental Lévy processes are the Wiener process and the Poisson process. The Poisson process is a subordinator, but is not stable; the Wiener process is stable, with exponent α = 2. Any linear combination of independent Lévy processes is again a Lévy process, so, for instance, if the Wiener process Wt and the Poisson process Nt are independent then Wt Nt is a Lévy process. More important, linear combinations of independent Poisson− processes are Lévy processes: these are special cases of what are called compound Poisson processes: see sec.
  • Random Measures, Theory and Applications

    Random Measures, Theory and Applications

    Olav Kallenberg Random Measures, Theory and Applications 123 Olav Kallenberg Department of Mathematics and Statistics Auburn University Auburn, Alabama USA ISSN 2199-3130 ISSN 2199-3149 (electronic) Probability Theory and Stochastic Modelling ISBN 978-3-319-41596-3 ISBN 978-3-319-41598-7 (eBook) DOI 10.1007/978-3-319-41598-7 Library of Congress Control Number: 2017936349 Mathematics Subject Classification (2010): 60G55, 60G57 © Springer International Publishing Switzerland 2017 This work is subject to copyright. All rights are reserved by the Publisher, whether the whole or part of the material is concerned, specifically the rights of translation, reprinting, reuse of illustrations, recitation, broadcasting, reproduction on microfilms or in any other physical way, and transmission or information storage and retrieval, electronic adaptation, computer software, or by similar or dissimilar methodology now known or hereafter developed. The use of general descriptive names, registered names, trademarks, service marks, etc. in this publication does not imply, even in the absence of a specific statement, that such names are exempt from the relevant protective laws and regulations and therefore free for general use. The publisher, the authors and the editors are safe to assume that the advice and information in this book are believed to be true and accurate at the date of publication. Neither the publisher nor the authors or the editors give a warranty, express or implied, with respect to the material contained herein or for any errors or omissions
  • Dependent Normalized Random Measures

    Dependent Normalized Random Measures

    Dependent Normalized Random Measures Changyou Chen A thesis submitted for the degree of Doctor of Philosophy The Australian National University February 2014 c Changyou Chen 2014 To all my family, on earth and in heaven, for their selfless care and love. Acknowledgments It is time to say goodbye to my PhD research in the Australian National University. It has been an exciting and memorable experience and I would like to take this opportunity to thank everyone who has helped me during my PhD. My utmost gratitude goes to my supervisor, Professor Wray Buntine, for both his extraordinary supervision in my research and selfless care in my daily life. Wray brought me into the Bayesian nonparametrics world, he has taught me how to per- form quality research from the beginning of my PhD. He taught me how to do critical thinking, how to design professional experiments, how to write quality papers, and how to present my research work. I benefited greatly from his experience which has helped me develop as a machine learning researcher. I would like to express my gratitude and respect to Professor Yee Whye Teh, who hosted me as a visiting student at UCL for a short time in 2012. My special thanks also goes to Vinayak Rao, who explained to me his spatial normalized Gamma pro- cesses, a precursor of our joint work together with Yee Whye. I was impressed by Yee Whye and Vinayak’s knowledge in Bayesian nonparametrics, and their sophisticated ideas for Bayesian posterior inference. I was lucky to work with them and had our joint work published in ICML.
  • Levy Processes, Stable Processes, and Subordinators

    Levy Processes, Stable Processes, and Subordinators

    LÉVY PROCESSES, STABLE PROCESSES, AND SUBORDINATORS STEVEN P.LALLEY 1. DEFINITIONS AND EXAMPLES d Definition 1.1. A continuous–time process X t = X (t ) t 0 with values in R (or, more generally, in an abelian topological group G{) is called} a≥Lévy process if (1) its sample paths are right-continuous and have left limits at every time point t , and (2) it has stationary, independent increments, that is: (a) For all 0 = t0 < t1 < < tk , the increments X (ti ) X (ti 1) are independent. (b) For all 0 s t the··· random variables X (t ) X (−s ) and−X (t s ) X (0) have the same distribution. − − − The default initial condition is X0 = 0. It can be shown that every process with sta- tionary, independent increments has a version with paths that are right-continuous with left limits, but this fact won’t be needed in this course. See BREIMAN, Probability, chs. 9, 14 for the proof. Definition 1.2. A subordinator is a real-valued Lévy process with nondecreasing sam- ple paths. A stable process is a real-valued Lévy process X t t 0 with initial value X0 = 0 that satisfies the self-similarity property { } ≥ 1/↵ (1.1) X t /t =D X1 t > 0. 8 The parameter ↵ is called the exponent of the process. Example 1.1. The fundamental Lévy processes are the Wiener process and the Poisson process. The Poisson process is a subordinator, but is not stable; the Wiener process is stable, with exponent ↵ = 2. Any linear combination of independent Lévy processes is again a Lévy process, so, for instance, if the Wiener process Wt and the Poisson process Nt are independent then Wt Nt is a Lévy process.
  • Poisson Point Processes and Subordinators

    Poisson Point Processes and Subordinators

    CHAPTER 1 Poisson point processes and subordinators. In this chapter, we introduce basic notions on Poisson point processes and subordina- tors. Poisson processes and two remarkable families of related martingales are studied. We also introduce the notion of Poisson random measures in order to define the Poisson point process. The last part of this chapter concerns to subordinators and their connection with the Levy-Kinthine´ formula. 1. Poisson point processes 1.1. Poisson processes. DEFINITION 1. A Poisson process with parameter c > 0 is a renewal process where the time between occurrences is exponentially distributed with parameter c. More precisely take a sequence (τn, n ≥ 1) of independent exponential random variables with parameter c and introduce the partial sums Sn = τ1 + ··· + τn, n ∈ IN. The right-continuous inverse n o Nt = sup n ∈ IN : Sn ≤ t , t ≥ 0, is called a Poisson process of parameter c. Let us explain some details of the above definition. We first recall that Sn has the same law as a Gamma distribution with parameters c and n. Therefore, for any fixed k ∈ IN and t ∈ IR+ P(Nt = k) = E 1I{Sk≤t<Sk+1} = E 1I{Sk≤t}P(τn+1 ≥ t − Sk) 1 Z t cktk = e−c(t−x)cke−cxxk−1dx = e−ct. Γ(k) 0 k! This implies that for any fixed t > 0, Nt is a Poisson r.v. with parameter tc, from where this process takes his name. The lack of memory property of the exponential law implies that for every 0 ≤ s ≤ t, the increment Nt+s − Nt has the Poisson distribution with parameter cs and is independent of the σ-field generated by (Nu, 0 ≤ u ≤ t), i.e.
  • Notes on the Poisson Point Process

    Notes on the Poisson Point Process

    Notes on the Poisson point process Paul Keeler March 20, 2018 This work is licensed under a “CC BY-SA 3.0” license. Abstract The Poisson point process is a type of random object in mathematics known as a point process. For over a century this point process has been the focus of much study and application. This survey aims to give an accessible but detailed account of the Poisson point process by covering its history, mathematical defini- tions in a number of settings, and key properties as well detailing various termi- nology and applications of the point process, with recommendations for further reading. 1 Introduction In probability, statistics and related fields, a Poisson point process or a Pois- son process (also called a Poisson random measure or a Poisson point field) is a type of random object known as a point process that consists of randomly po- sitioned points located on some underlying mathematical space [51, 54]. The Poisson point process has interesting mathematical properties [51][54, pages xv- xvi], which has led to it being frequently defined in Euclidean space and used as a mathematical model for seemingly random objects in numerous disciplines such as astronomy [4], biology [62], ecology [84], geology [17], physics [74], im- age processing [10], and telecommunications [33]. For example, defined on the real line, it can be interpreted as a stochastic process 1 called the Poisson pro- cess [83, pages 7–8][2][70, pages 94–95], forming an example of a Levy´ process and a Markov process [73, page ix] [2]. In this setting, the Poisson process plays 1Point processes and stochastic processes are usually considered as different objects, with point processes being associated or related to certain stochastic processes [20, page 194][18, page 3][83, pages 7–8][25, page 580].
  • An Introduction to L\'{E} Vy Processes with Applications in Finance

    An Introduction to L\'{E} Vy Processes with Applications in Finance

    AN INTRODUCTION TO LEVY´ PROCESSES WITH APPLICATIONS IN FINANCE ANTONIS PAPAPANTOLEON Abstract. These lectures notes aim at introducing L´evyprocesses in an informal and intuitive way, accessible to non-specialists in the field. In the first part, we focus on the theory of L´evyprocesses. We analyze a `toy' example of a L´evyprocess, viz. a L´evyjump-diffusion, which yet offers significant insight into the distributional and path structure of a L´evyprocess. Then, we present several important results about L´evy processes, such as infinite divisibility and the L´evy-Khintchine formula, the L´evy-It^odecomposition, the It^oformula for L´evyprocesses and Gir- sanov's transformation. Some (sketches of) proofs are presented, still the majority of proofs is omitted and the reader is referred to textbooks instead. In the second part, we turn our attention to the applications of L´evyprocesses in financial modeling and option pricing. We discuss how the price process of an asset can be modeled using L´evyprocesses and give a brief account of market incompleteness. Popular models in the literature are presented and revisited from the point of view of L´evy processes, and we also discuss three methods for pricing financial deriva- tives. Finally, some indicative evidence from applications to market data is presented. Contents Part 1. Theory 2 1. Introduction2 2. Definition5 3. `Toy' example: a L´evyjump-diffusion6 4. Infinitely divisible distributions and the L´evy-Khintchine formula8 5. Analysis of jumps and Poisson random measures 11 6. The L´evy-It^odecomposition 12 7. The L´evymeasure, path and moment properties 14 8.
  • Looking for Poisson-Like Random Measures

    Looking for Poisson-Like Random Measures

    Received: 18 March 2019 DOI: 10.1002/mma.6224 RESEARCH ARTICLE Throwing stones and collecting bones: Looking for Poisson-like random measures Caleb Deen Bastian1 Grzegorz A. Rempala2 1Program in Applied and Computational Mathematics, Princeton University, We show that in a broad class of random counting measures, one may identify Princeton, New Jersey only three that are rescaled versions of themselves when restricted to a sub- 2 Division of Biostatistics and Department space. These are Poisson, binomial, and negative binomial random measures. of Mathematics, The Ohio State University, Columbus, Ohio We provide some simple examples of possible applications of such measures. Correspondence KEYWORDS Caleb Deen Bastian, Program in Applied Laplace functional, Poisson-type (PT) distributions, random counting measure, stone throwing and Computational Mathematics, construction, strong invariance, thinning Washington Road, Fine Hall, Princeton University, Princeton, NJ. 08544. Email: [email protected] MSC CLASSIFICATION 60G57; 60G18; 60G55 Communicatedby: S. Wise Funding information National Science Foundation, Grant/Award Number: DMS 1853587 1 INTRODUCTION Random counting measures, also known as point processes, are the central objects of this note. For general introduc- tion, see, for instance, monographs by Olav1,2 or Erhan.3 Random counting measures have numerous uses in statistics and applied probability, including representation and construction of stochastic processes, Monte Carlo schemes, etc. For example, the Poisson random measure is a fundamental random counting measure that is related to the structure of Lévy processes, Markov jump processes, or the excursions of Brownian motion, and is prototypical to the class of completely random (additive) random measures.3 In particular, it is also well known that the Poisson random measure is self-similar in the sense of being invariant under restriction to a subspace (invariant under thinning).
  • Poisson Random Measures and Their Numerical Approximation

    Poisson Random Measures and Their Numerical Approximation

    SPDES driven by Poisson Random Measures and their numerical Approximation Erika Hausenblas joint work with Thomas Dunst and Andreas Prohl University of Salzburg, Austria SPDES driven by Poisson Random Measures and their numerical Approximation – p.1 Outline Lévy processes - Poisson Random Measure Stochastic Integration in Banach spaces SPDEs driven by Poisson Random Measure Their Numerical Approximation SPDES driven by Poisson Random Measures and their numerical Approximation – p.2 A typical Example Let be a bounded domain in Rd with smooth boundary. O The Equation: du(t,ξ) ∂2 = 2 u(t, ξ)+ α u(t, ξ)+ g(u(t, ξ))L˙ (t) dt ∂ξ ∇ + f(u(t, ξ)), ξ , t> 0; ∈ O (⋆) u(0, ξ) = u (ξ) ξ ; 0 ∈ O u(t, ξ) = u(t, ξ) = 0, t 0, ξ ∂ ; ≥ ∈ O where u Lp( ), p 1, g a certain mapping and L(t) is a space time 0 ∈ O ≥ Lévy noise specified later. Problem: To find a function u : [0, ) R ∞ ×O −→ solving Equation (⋆). SPDES driven by Poisson Random Measures and their numerical Approximation – p.3 The Lévy Process L Definition 1 Let E be a Banach space. A stochastic process L = L(t), 0 t< is an E–valued Lévy process over (Ω; ; P) if the { ≤ ∞} F following conditions are satisfied: L(0) = 0; L has independent and stationary increments; L is stochastically continuous, i.e. for φ bounded and measurable, the function t Eφ(L(t)) is continuous on R+; 7→ L has a.s. càdlàg paths; SPDES driven by Poisson Random Measures and their numerical Approximation – p.4 The Lévy Process L E denotes a separable Banach space and E′ the dual on E.