Theory of the Jump-Diffusion Processes

Total Page:16

File Type:pdf, Size:1020Kb

Theory of the Jump-Diffusion Processes OUTLINE OF THE WHOLE LECTURES. PART 1 : JUMP-DIFFUSION PROCESSES. I Theory and simulation. I Stochastic calculus with jump-diffusion processes. I Applications in finance. PART 2 : GENERAL LÉVY PROCESSES. I Theory. I Simulation. I Stochastic calculus with semimartingales. I Applications in finance. A. Popier (ENSTA) Jump-diffusion processes. January 2020 1 / 34 INTRODUCTION. CONTINUOUS TIME MODELS WITH CONTINUOUS TRAJECTORIES. Black-Scholes model. Scale invariance of the Brownian motion. Local volatility models. Possible perfect hedging. Stochastic volatility models. Difficulty to obtain heavy tails, no large sudden moves. CONTINUOUS TIME MODELS WITH DISCONTINUOUS TRAJECTORIES. I Market crash. I Credit risk. I High-frequency trading. Aït-Sahalia & Jacod, High-frequency Financial Econometrics. I Insurance and ruin theory. A. Popier (ENSTA) Jump-diffusion processes. January 2020 2 / 34 INTRODUCTION. CONTINUOUS TIME MODELS WITH CONTINUOUS TRAJECTORIES. Black-Scholes model. Scale invariance of the Brownian motion. Local volatility models. Possible perfect hedging. Stochastic volatility models. Difficulty to obtain heavy tails, no large sudden moves. CONTINUOUS TIME MODELS WITH DISCONTINUOUS TRAJECTORIES. I Market crash. I Credit risk. I High-frequency trading. Aït-Sahalia & Jacod, High-frequency Financial Econometrics. I Insurance and ruin theory. A. Popier (ENSTA) Jump-diffusion processes. January 2020 2 / 34 TWO DIFFERENT CLASSES. 1 JUMP-DIFFUSION MODELS. Prices : diffusion process, with jumps at randon times. Jumps : rare events −! cracks or large losses. 1 Advantages: I price structure : easy to understand, to describe and to simulate ; I then efficient Monte Carlo methods to compute path-depend prices. I Very performant to interpolate the implicit volatility smiles. 2 Inconvenients: I unknown closed formula for the densities, I statistical estimation or moments/quantiles computations : difficult to realize. 2 INFINITE ACTIVITY MODELS (general Lévy processes). A. Popier (ENSTA) Jump-diffusion processes. January 2020 3 / 34 TWO DIFFERENT CLASSES. 1 JUMP-DIFFUSION MODELS. 2 INFINITE ACTIVITY MODELS (general Lévy processes). Models with a infinite number of jumps during any time period. Unnecessary Brownian component. 1 Advantages: I give a more realistic description of the prices at different time scales. I often obtained as subordinator of a Brownian motion (time change), I hence closed formulas or more tractable than for the jump-diffusion models. 2 Inconvenients: I often more complicated to simulate. I Price structure less intuitive. A. Popier (ENSTA) Jump-diffusion processes. January 2020 3 / 34 DEFINITIONS. DEFINITION d A stochastic process (Xt )t≥0 defined on (Ω; F; P), with values in R , is a Lévy process if 1 X0 = 0 a.s. 2 its increments are independent : for every increasing sequence t0;:::; tn, the r.v. Xt0 ; Xt1 − Xt0 ;:::; Xtn − Xtn−1 are independent ; 3 its increments are stationnary : the law of Xt+h − Xt does not depend on t ; 4 X satisfies the property called stochastic continuity : for any " > 0, lim P(jXt+h − Xt j ≥ ") = 0 h!0 5 there exists a subset Ω0 s.t. P(Ω0) = 1 and for every ! 2 Ω0, t 7! Xt (!) is RCLL. A. Popier (ENSTA) Jump-diffusion processes. January 2020 4 / 34 DEFINITIONS. If a filtration (Ft )t≥0 is already given on (Ω; F; P) DEFINITION d A stochastic process (Xt )t≥0 defined on (Ω; F; P), with values in R , is a Lévy process if 1 X0 = 0 a.s. 2 its increments are independent : for any s ≤ t, the r.v. Xt − Xs is independent of Fs ; 3 its increments are stationnary; 4 X satisfies the property called stochastic continuity; 5 a.s. t 7! Xt (!) is RCLL. A. Popier (ENSTA) Jump-diffusion processes. January 2020 4 / 34 REMARKS. REMARKS ON THE DEFINITIONS : X If Ft = Ft , the two definitions are equivalent. X If fFt g is a larger filtration than (Ft ⊂ Ft ) and if Xt − Xs is independent of Fs, then fXt ; 0 ≤ t < +1g is a Lévy process under the large filtration. REMARKS ON THE HYPOTHESES : If we remove Assumption 5, we speak about Lévy process in law. If we remove Assumption 3, we obtain an additive process. Dropping Assumptions 3 and 5, we have an additive process in law. THEOREM A Lévy process (or an additive process) in law has a RCLL modification. We can also prove that 2, 3 and 5 imply 4. A. Popier (ENSTA) Jump-diffusion processes. January 2020 5 / 34 REMARKS. REMARKS ON THE DEFINITIONS : X If Ft = Ft , the two definitions are equivalent. X If fFt g is a larger filtration than (Ft ⊂ Ft ) and if Xt − Xs is independent of Fs, then fXt ; 0 ≤ t < +1g is a Lévy process under the large filtration. REMARKS ON THE HYPOTHESES : If we remove Assumption 5, we speak about Lévy process in law. If we remove Assumption 3, we obtain an additive process. Dropping Assumptions 3 and 5, we have an additive process in law. THEOREM A Lévy process (or an additive process) in law has a RCLL modification. We can also prove that 2, 3 and 5 imply 4. A. Popier (ENSTA) Jump-diffusion processes. January 2020 5 / 34 ENLARGED FILTRATION. For a process X = fXt ; t ≥ 0g, we define N1 = N the set of P-negligible events. X For any 0 ≤ t ≤ 1, augmented filtration : Ft = σ(Ft [N ). THEOREM Let X = fXt ; t ≥ 0g be a Lévy process. Then the augmented filtration fFt g is right-continuous. With respect to the enlarged filtration, fXt ; t ≥ 0g is still a Lévy process. A. Popier (ENSTA) Jump-diffusion processes. January 2020 6 / 34 CHARACTERISTIC FUNCTIONS. PROPOSITION d Let (Xt )t≥0 be a Lévy process in R . Then there exists a function : Rd ! R called characteristic exponent of X s.t. : d ihz;X i t (z) 8z 2 R ; E e t = e : A. Popier (ENSTA) Jump-diffusion processes. January 2020 7 / 34 THEORY OF THE JUMP-DIFFUSION PROCESSES. Alexandre Popier ENSTA, Palaiseau January 2020 A. Popier (ENSTA) Jump-diffusion processes. January 2020 8 / 34 OUTLINE 1 JUMP-DIFFUSION PROCESSES Compound Poisson process Jump-diffusion processes 2 PROPERTIES OF JUMP-DIFFUSION PROCESSES Densities of the jump-diffusion processes. Moments Jump-diffusion processes and martingales. Poisson and Lévy measures A. Popier (ENSTA) Jump-diffusion processes. January 2020 9 / 34 BROWNIAN MOTION. A Brownian motion is a Lévy process satisfying 1 for every t > 0, Xt is Gaussian with mean vector zero and covariance matrix t Id ; 2 the process X has continuous sample paths a.s. Characteristic function : ihz;B i 2 E(e t ) = exp(−tjzj =2): A. Popier (ENSTA) Jump-diffusion processes. January 2020 10 / 34 OUTLINE 1 JUMP-DIFFUSION PROCESSES Compound Poisson process Jump-diffusion processes 2 PROPERTIES OF JUMP-DIFFUSION PROCESSES Densities of the jump-diffusion processes. Moments Jump-diffusion processes and martingales. Poisson and Lévy measures A. Popier (ENSTA) Jump-diffusion processes. January 2020 11 / 34 POISSON PROCESS. DEFINITION A stochastic process (Xt )t≥0, with values in R, is a Poisson process with intensity λ > 0 if it is a Lévy process s.t. for every t > 0,Xt has a Poisson law with parameter λt. PROPOSITION (CONSTRUCTION) If (Tn)n2N is a random walk on R s.t. for every n ≥ 1,Tn − Tn−1 is exponentially distributed with parameter λ (with T0 = 0), then the process (Xt )t≥0 defined by Xt = n () Tn ≤ t < Tn+1 is a Poisson process with intensity λ. A. Popier (ENSTA) Jump-diffusion processes. January 2020 12 / 34 POISSON PROCESS. A. Popier (ENSTA) Jump-diffusion processes. January 2020 13 / 34 JUMP DISTRIBUTION. For any interval I, we denote by N(I) the number of jumps of Xt , t 2 I. PROPOSITION For 0 < s < t and n ≥ 1, the conditional law of Xs knowing Xt = n is binomial with parameters n and s=t. For 0 = t0 < t1 < : : : < tk = t and Ij =]tj−1; tj ], the conditional law of (N(I1);:::; N(Ik )) knowing Xt = n is multinomial with parameters n, (t1 − t0)=t;:::; (tk − tk−1)=t. Ti : jump times of the process. PROPOSITION Let n ≥ 1 and t > 0. The conditional law of T1;:::; Tn knowing Xt = n coincides with the law of the order statistics U(1);:::; u(n) of n independent variables, uniformly distributed on [0; t]. A. Popier (ENSTA) Jump-diffusion processes. January 2020 14 / 34 COMPOUND POISSON PROCESSES. We consider a Poisson process (Pt )t≥0 with intensity λ and jump times d Tn, and a sequence (Yn)n2N∗ of R -valued r.v. s.t. 1 Yn are i.i.d. with distribution measure π ; 2 (Pt )t≥0 and (Yn)n2N∗ are independent. Define P +1 Xt X Xt = Yn = Yn1[0;t](Tn): n=1 n=1 DEFINITION The process (Xt )t≥0 is a compound Poisson processes with intensity λ and jump distribution π. A. Popier (ENSTA) Jump-diffusion processes. January 2020 15 / 34 MERTON MODEL. Compound Poisson processes with Gaussian jumps. A. Popier (ENSTA) Jump-diffusion processes. January 2020 16 / 34 COMPOUND POISSON PROCESSES. PROPOSITION The process (Xt )t≥0 is a Lévy process, with piecewise constant trajectories and characteristic function : Z d ihz;X i ihz;xi 8z 2 R ; E e t = exp tλ (e − 1)π(dx) Rd Z = exp t (eihz;xi − 1)ν(dx) ; Rd DEFINITION ν is a finite measure defined on Rd by : ν(A) = λπ(A); A 2 B(Rd ): ν is called the Lévy measure of the compound Poisson process. Moreover ν(A) = E [#ft 2 [0; 1]; ∆Xt 6= 0; ∆Xt 2 Ag] : A. Popier (ENSTA) Jump-diffusion processes. January 2020 17 / 34 COMPOUND POISSON PROCESSES. DEFINITION The law µ of X1 is called compound Poisson distribution and has a Z characteristic function given by : µ^(z) = exp λ (eihz;xi − 1)π(dx) .
Recommended publications
  • Numerical Solution of Jump-Diffusion LIBOR Market Models
    Finance Stochast. 7, 1–27 (2003) c Springer-Verlag 2003 Numerical solution of jump-diffusion LIBOR market models Paul Glasserman1, Nicolas Merener2 1 403 Uris Hall, Graduate School of Business, Columbia University, New York, NY 10027, USA (e-mail: [email protected]) 2 Department of Applied Physics and Applied Mathematics, Columbia University, New York, NY 10027, USA (e-mail: [email protected]) Abstract. This paper develops, analyzes, and tests computational procedures for the numerical solution of LIBOR market models with jumps. We consider, in par- ticular, a class of models in which jumps are driven by marked point processes with intensities that depend on the LIBOR rates themselves. While this formulation of- fers some attractive modeling features, it presents a challenge for computational work. As a first step, we therefore show how to reformulate a term structure model driven by marked point processes with suitably bounded state-dependent intensities into one driven by a Poisson random measure. This facilitates the development of discretization schemes because the Poisson random measure can be simulated with- out discretization error. Jumps in LIBOR rates are then thinned from the Poisson random measure using state-dependent thinning probabilities. Because of discon- tinuities inherent to the thinning process, this procedure falls outside the scope of existing convergence results; we provide some theoretical support for our method through a result establishing first and second order convergence of schemes that accommodates thinning but imposes stronger conditions on other problem data. The bias and computational efficiency of various schemes are compared through numerical experiments. Key words: Interest rate models, Monte Carlo simulation, market models, marked point processes JEL Classification: G13, E43 Mathematics Subject Classification (1991): 60G55, 60J75, 65C05, 90A09 The authors thank Professor Steve Kou for helpful comments and discussions.
    [Show full text]
  • MERTON JUMP-DIFFUSION MODEL VERSUS the BLACK and SCHOLES APPROACH for the LOG-RETURNS and VOLATILITY SMILE FITTING Nicola Gugole
    International Journal of Pure and Applied Mathematics Volume 109 No. 3 2016, 719-736 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu AP doi: 10.12732/ijpam.v109i3.19 ijpam.eu MERTON JUMP-DIFFUSION MODEL VERSUS THE BLACK AND SCHOLES APPROACH FOR THE LOG-RETURNS AND VOLATILITY SMILE FITTING Nicola Gugole Department of Computer Science University of Verona Strada le Grazie, 15-37134, Verona, ITALY Abstract: In the present paper we perform a comparison between the standard Black and Scholes model and the Merton jump-diffusion one, from the point of view of the study of the leptokurtic feature of log-returns and also concerning the volatility smile fitting. Provided results are obtained by calibrating on market data and by mean of numerical simulations which clearly show how the jump-diffusion model outperforms the classical geometric Brownian motion approach. AMS Subject Classification: 60H15, 60H35, 91B60, 91G20, 91G60 Key Words: Black and Scholes model, Merton model, stochastic differential equations, log-returns, volatility smile 1. Introduction In the early 1970’s the world of option pricing experienced a great contribu- tion given by the work of Fischer Black and Myron Scholes. They developed a new mathematical model to treat certain financial quantities publishing re- lated results in the article The Pricing of Options and Corporate Liabilities, see [1]. The latter work became soon a reference point in the financial scenario. Received: August 3, 2016 c 2016 Academic Publications, Ltd. Revised: September 16, 2016 url: www.acadpubl.eu Published: September 30, 2016 720 N. Gugole Nowadays, many traders still use the Black and Scholes (BS) model to price as well as to hedge various types of contingent claims.
    [Show full text]
  • Sensitivity Estimation of Sabr Model Via Derivative of Random Variables
    Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelspach, K. P. White, and M. Fu, eds. SENSITIVITY ESTIMATION OF SABR MODEL VIA DERIVATIVE OF RANDOM VARIABLES NAN CHEN YANCHU LIU The Chinese University of Hong Kong 709A William Mong Engineering Building Shatin, N. T., HONG KONG ABSTRACT We derive Monte Carlo simulation estimators to compute option price sensitivities under the SABR stochastic volatility model. As a companion to the exact simulation method developed by Cai, Chen and Song (2011), this paper uses the sensitivity of “vol of vol” as a showcase to demonstrate how to use the pathwise method to obtain unbiased estimators to the price sensitivities under SABR. By appropriately conditioning on the path generated by the volatility, the evolution of the forward price can be represented as noncentral chi-square random variables with stochastic parameters. Combined with the technique of derivative of random variables, we can obtain fast and accurate unbiased estimators for the sensitivities. 1 INTRODUCTION The ubiquitous existence of the volatility smile and skew poses a great challenge to the practice of risk management in fixed and foreign exchange trading desks. In foreign currency option markets, the implied volatility is often relatively lower for at-the-money options and becomes gradually higher as the strike price moves either into the money or out of the money. Traders refer to this stylized pattern as the volatility smile. In the equity and interest rate markets, a typical aspect of the implied volatility, which is also known as the volatility skew, is that it decreases as the strike price increases.
    [Show full text]
  • Option Pricing with Jump Diffusion Models
    UNIVERSITY OF PIRAEUS DEPARTMENT OF BANKING AND FINANCIAL MANAGEMENT M. Sc in FINANCIAL ANALYSIS FOR EXECUTIVES Option pricing with jump diffusion models MASTER DISSERTATION BY: SIDERI KALLIOPI: MXAN 1134 SUPERVISOR: SKIADOPOULOS GEORGE EXAMINATION COMMITTEE: CHRISTINA CHRISTOU PITTIS NIKITAS PIRAEUS, FEBRUARY 2013 2 TABLE OF CONTENTS ς ABSTRACT .............................................................................................................................. 3 SECTION 1 INTRODUCTION ................................................................................................ 4 ώ SECTION 2 LITERATURE REVIEW: .................................................................................... 9 ι SECTION 3 DATASET .......................................................................................................... 14 SECTION 4 .............................................................................................................................α 16 4.1 PRICING AND HEDGING IN INCOMPLETE MARKETS ....................................... 16 4.2 CHANGE OF MEASURE ............................................................................................ρ 17 SECTION 5 MERTON’S MODEL .........................................................................................ι 21 5.1 THE FORMULA ........................................................................................................... 21 5.2 PDE APPROACH BY HEDGING ................................................................ε ............... 24 5.3 EFFECT
    [Show full text]
  • A University of Sussex Phd Thesis Available Online Via Sussex
    A University of Sussex PhD thesis Available online via Sussex Research Online: http://sro.sussex.ac.uk/ This thesis is protected by copyright which belongs to the author. This thesis cannot be reproduced or quoted extensively from without first obtaining permission in writing from the Author The content must not be changed in any way or sold commercially in any format or medium without the formal permission of the Author When referring to this work, full bibliographic details including the author, title, awarding institution and date of the thesis must be given Please visit Sussex Research Online for more information and further details NON-STATIONARY PROCESSES AND THEIR APPLICATION TO FINANCIAL HIGH-FREQUENCY DATA Mailan Trinh A thesis submitted for the degree of Doctor of Philosophy University of Sussex March 2018 UNIVERSITY OF SUSSEX MAILAN TRINH A THESIS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY NON-STATIONARY PROCESSES AND THEIR APPLICATION TO FINANCIAL HIGH-FREQUENCY DATA SUMMARY The thesis is devoted to non-stationary point process models as generalizations of the standard homogeneous Poisson process. The work can be divided in two parts. In the first part, we introduce a fractional non-homogeneous Poisson process (FNPP) by applying a random time change to the standard Poisson process. We character- ize the FNPP by deriving its non-local governing equation. We further compute moments and covariance of the process and discuss the distribution of the arrival times. Moreover, we give both finite-dimensional and functional limit theorems for the FNPP and the corresponding fractional non-homogeneous compound Poisson process. The limit theorems are derived by using martingale methods, regular vari- ation properties and Anscombe's theorem.
    [Show full text]
  • Lévy Finance *[0.5Cm] Models and Results
    Stochastic Calculus for L´evyProcesses L´evy-Process Driven Financial Market Models Jump-Diffusion Models General L´evyModels European Style Options Stochastic Calculus for L´evy Processes L´evy-Process Driven Financial Market Models Jump-Diffusion Models Merton-Model Kou-Model General L´evy Models Variance-Gamma model CGMY model GH models Variance-mean mixtures European Style Options Equivalent Martingale Measure Jump-Diffusion Models Variance-Gamma Model NIG Model Professor Dr. R¨udigerKiesel L´evyFinance Stochastic Calculus for L´evyProcesses L´evy-Process Driven Financial Market Models Jump-Diffusion Models General L´evyModels European Style Options Stochastic Integral for L´evyProcesses Let (Xt ) be a L´evy process with L´evy-Khintchine triplet (α, σ, ν(dx)). By the L´evy-It´odecomposition we know X = X (1) + X (2) + X (3), where the X (i) are independent L´evyprocesses. X (1) is a Brownian motion with drift, X (2) is a compound Poisson process with jump (3) distributed concentrated on R/(−1, 1) and X is a square-integrable martingale (which can be viewed as a limit of compensated compound Poisson processes with small jumps). We know how to define the stochastic integral with respect to any of these processes! Professor Dr. R¨udigerKiesel L´evyFinance Stochastic Calculus for L´evyProcesses L´evy-Process Driven Financial Market Models Jump-Diffusion Models General L´evyModels European Style Options Canonical Decomposition From the L´evy-It´odecomposition we deduce the canonical decomposition (useful for applying the general semi-martingale theory) Z t Z X (t) = αt + σW (t) + x µX − νX (ds, dx), 0 R where Z t Z X xµX (ds, dx) = ∆X (s) 0 R 0<s≤t and Z t Z Z t Z Z X X E xµ (ds, dx) = xν (ds, dx) = t xν(dx).
    [Show full text]
  • STOCHASTIC COMPARISONS and AGING PROPERTIES of an EXTENDED GAMMA PROCESS Zeina Al Masry, Sophie Mercier, Ghislain Verdier
    STOCHASTIC COMPARISONS AND AGING PROPERTIES OF AN EXTENDED GAMMA PROCESS Zeina Al Masry, Sophie Mercier, Ghislain Verdier To cite this version: Zeina Al Masry, Sophie Mercier, Ghislain Verdier. STOCHASTIC COMPARISONS AND AGING PROPERTIES OF AN EXTENDED GAMMA PROCESS. Journal of Applied Probability, Cambridge University press, 2021, 58 (1), pp.140-163. 10.1017/jpr.2020.74. hal-02894591 HAL Id: hal-02894591 https://hal.archives-ouvertes.fr/hal-02894591 Submitted on 9 Jul 2020 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Applied Probability Trust (4 July 2020) STOCHASTIC COMPARISONS AND AGING PROPERTIES OF AN EXTENDED GAMMA PROCESS ZEINA AL MASRY,∗ FEMTO-ST, Univ. Bourgogne Franche-Comt´e,CNRS, ENSMM SOPHIE MERCIER & GHISLAIN VERDIER,∗∗ Universite de Pau et des Pays de l'Adour, E2S UPPA, CNRS, LMAP, Pau, France Abstract Extended gamma processes have been seen to be a flexible extension of standard gamma processes in the recent reliability literature, for cumulative deterioration modeling purpose. The probabilistic properties of the standard gamma process have been well explored since the 1970's, whereas those of its extension remain largely unexplored. In particular, stochastic comparisons between degradation levels modeled by standard gamma processes and aging properties for the corresponding level-crossing times are nowadays well understood.
    [Show full text]
  • Lecture Notes
    Lecture Notes Lars Peter Hansen October 8, 2007 2 Contents 1 Approximation Results 5 1.1 One way to Build a Stochastic Process . 5 1.2 Stationary Stochastic Process . 6 1.3 Invariant Events and the Law of Large Numbers . 6 1.4 Another Way to Build a Stochastic Process . 8 1.4.1 Stationarity . 8 1.4.2 Limiting behavior . 9 1.4.3 Ergodicity . 10 1.5 Building Nonstationary Processes . 10 1.6 Martingale Approximation . 11 3 4 CONTENTS Chapter 1 Approximation Results 1.1 One way to Build a Stochastic Process • Consider a probability space (Ω, F, P r) where Ω is a set of sample points, F is an event collection (sigma algebra) and P r assigns proba- bilities to events . • Introduce a function S :Ω → Ω such that for any event Λ, S−1(Λ) = {ω ∈ Ω: S(ω) ∈ Λ} is an event. • Introduce a (Borel measurable) measurement function x :Ω → Rn. x is a random vector. • Construct a stochastic process {Xt : t = 1, 2, ...} via the formula: t Xt(ω) = X[S (ω)] or t Xt = X ◦ S . Example 1.1.1. Let Ω be a collection of infinite sequences of real numbers. Specifically, ω = (r0, r1, ...), S(ω) = (r1, r2, ...) and x(ω) = r0. Then Xt(ω) = rt. 5 6 CHAPTER 1. APPROXIMATION RESULTS 1.2 Stationary Stochastic Process Definition 1.2.1. The transformation S is measure-preserving if: P r(Λ) = P r{S−1(Λ)} for all Λ ∈ F. Proposition 1.2.2. When S is measure-preserving, the process {Xt : t = 1, 2, ...} has identical distributions for every t.
    [Show full text]
  • Arbitrage Free Approximations to Candidate Volatility Surface Quotations
    Journal of Risk and Financial Management Article Arbitrage Free Approximations to Candidate Volatility Surface Quotations Dilip B. Madan 1,* and Wim Schoutens 2,* 1 Robert H. Smith School of Business, University of Maryland, College Park, MD 20742, USA 2 Department of Mathematics, KU Leuven, 3000 Leuven, Belgium * Correspondence: [email protected] (D.B.M.); [email protected] (W.S.) Received: 19 February 2019; Accepted: 10 April 2019; Published: 21 April 2019 Abstract: It is argued that the growth in the breadth of option strikes traded after the financial crisis of 2008 poses difficulties for the use of Fourier inversion methodologies in volatility surface calibration. Continuous time Markov chain approximations are proposed as an alternative. They are shown to be adequate, competitive, and stable though slow for the moment. Further research can be devoted to speed enhancements. The Markov chain approximation is general and not constrained to processes with independent increments. Calibrations are illustrated for data on 2695 options across 28 maturities for SPY as at 8 February 2018. Keywords: bilateral gamma; fast Fourier transform; sato process; matrix exponentials JEL Classification: G10; G12; G13 1. Introduction A variety of arbitrage free models for approximating quoted option prices have been available for some time. The quotations are typically in terms of (Black and Scholes 1973; Merton 1973) implied volatilities for the traded strikes and maturities. Consequently, the set of such quotations is now popularly referred to as the volatility surface. Some of the models are nonparametric with the Dupire(1994) local volatility model being a popular example. The local volatility model is specified in terms of a deterministic function s(K, T) that expresses the local volatility for the stock if the stock were to be at level K at time T.
    [Show full text]
  • Portfolio Optimization and Option Pricing Under Defaultable Lévy
    Universit`adegli Studi di Padova Dipartimento di Matematica Scuola di Dottorato in Scienze Matematiche Indirizzo Matematica Computazionale XXVI◦ ciclo Portfolio optimization and option pricing under defaultable L´evy driven models Direttore della scuola: Ch.mo Prof. Paolo Dai Pra Coordinatore dell’indirizzo: Ch.ma Prof.ssa. Michela Redivo Zaglia Supervisore: Ch.mo Prof. Tiziano Vargiolu Universit`adegli Studi di Padova Corelatori: Ch.mo Prof. Agostino Capponi Ch.mo Prof. Andrea Pascucci Johns Hopkins University Universit`adegli Studi di Bologna Dottorando: Stefano Pagliarani To Carolina and to my family “An old man, who had spent his life looking for a winning formula (martingale), spent the last days of his life putting it into practice, and his last pennies to see it fail. The martingale is as elusive as the soul.” Alexandre Dumas Mille et un fantˆomes, 1849 Acknowledgements First and foremost, I would like to thank my supervisor Prof. Tiziano Vargiolu and my co-advisors Prof. Agostino Capponi and Prof. Andrea Pascucci, for coauthoring with me all the material appearing in this thesis. In particular, thank you to Tiziano Vargiolu for supporting me scientifically and finan- cially throughout my research and all the activities related to it. Thank you to Agostino Capponi for inviting me twice to work with him at Purdue University, and for his hospital- ity and kindness during my double stay in West Lafayette. Thank you to Andrea Pascucci for wisely leading me throughout my research during the last four years, and for giving me the chance to collaborate with him side by side, thus sharing with me his knowledge and his experience.
    [Show full text]
  • A Representation Free Quantum Stochastic Calculus
    JOURNAL OF FUNCTIONAL ANALYSIS 104, 149-197 (1992) A Representation Free Quantum Stochastic Calculus L. ACCARDI Centro Matemarico V. Volterra, Diparfimento di Matematica, Universitci di Roma II, Rome, Italy F. FACNOLA Dipartimento di Matematica, Vniversild di Trento, Povo, Italy AND J. QUAEGEBEUR Department of Mathemarics, Universiteit Leuven, Louvain, Belgium Communicated by L. Gross Received January 2; revised March 1, 1991 We develop a representation free stochastic calculus based on three inequalities (semimartingale inequality, scalar forward derivative inequality, scalar conditional variance inequality). We prove that our scheme includes all the previously developed stochastic calculi and some new examples. The abstract theory is applied to prove a Boson Levy martingale representation theorem in bounded form and a general existence, uniqueness, and unitarity theorem for quantum stochastic differential equations. 0 1992 Academic Press. Inc. 0. INTRODUCTION Stochastic calculus is a powerful tool in classical probability theory. Recently various kinds of stochastic calculi have been introduced in quan- tum probability [18, 12, 10,211. The common features of these calculi are: - One starts from a representation of the canonical commutation (or anticommutation) relations (CCR or CAR) over a space of the form L2(R+, dr; X) (where 3? is a given Hilbert space). - One introduces a family of operator valued measures on R,, related to this representation, e.g., expressed in terms of creation or annihilation or number or field operators. 149 0022-1236192 $3.00 Copyright 0 1992 by Academic Press, Inc. All rights of reproduction in any form reserved. 150 QUANTUM STOCHASTIC CALCULUS -- One shows that it is possible to develop a theory of stochastic integration with respect to these operator valued measures sufficiently rich to allow one to solve some nontrivial stochastic differential equations.
    [Show full text]
  • Normalized Random Measures Driven by Increasing Additive Processes
    The Annals of Statistics 2004, Vol. 32, No. 6, 2343–2360 DOI: 10.1214/009053604000000625 c Institute of Mathematical Statistics, 2004 NORMALIZED RANDOM MEASURES DRIVEN BY INCREASING ADDITIVE PROCESSES By Luis E. Nieto-Barajas1, Igor Prunster¨ 2 and Stephen G. Walker3 ITAM-M´exico, Universit`adegli Studi di Pavia and University of Kent This paper introduces and studies a new class of nonparamet- ric prior distributions. Random probability distribution functions are constructed via normalization of random measures driven by increas- ing additive processes. In particular, we present results for the dis- tribution of means under both prior and posterior conditions and, via the use of strategic latent variables, undertake a full Bayesian analysis. Our class of priors includes the well-known and widely used mixture of a Dirichlet process. 1. Introduction. This paper considers the problem of constructing a stochastic process, defined on the real line, which has sample paths be- having almost surely (a.s.) as a probability distribution function (d.f.). The law governing the process acts as a prior in Bayesian nonparametric prob- lems. One popular idea is to take the random probability d.f. as a normal- ized increasing process F (t)= Z(t)/Z¯, where Z¯ = limt→∞ Z(t) < +∞ (a.s.). For example, the Dirichlet process [Ferguson (1973)] arises when Z is a suitably reparameterized gamma process. We consider the case of normal- ized random d.f.s driven by an increasing additive process (IAP) L, that is, Z(t)= k(t,x) dL(x) and provide regularity conditions on k and L to ensure F is a random probability d.f.
    [Show full text]