Exact Simulation of a Truncated Lévy Subordinator

A truncated Lévy subordinator is a Lévy subordinator in R+ with Lévy measure restricted from above by a certain level b. In this article, we study the path and distribution properties of this type of process in detail and set up an exact simulation framework based on a marked renewal process. In particular, we focus on a typical specification of truncated Lévy subordinator, namely the truncated stable process. We establish an exact simulation algorithm for the truncated stable process, which is very accurate and efficient. Compared to the existing algorithm suggested in Chi, our algorithm outperforms over all parameter settings. Using the distributional decomposition technique, we also develop an exact simulation algorithm for the truncated tempered stable process and other related processes. We illustrate an application of our algorithm as a valuation tool for stochastic hyperbolic discounting, and numerical analysis is provided to demonstrate the accuracy and effectiveness of our methods. We also show that variations of the result can also be used to sample two-sided truncated Lévy processes, two-sided Lévy processes via subordinating Brownian motions, and truncated Lévy-driven Ornstein-Uhlenbeck processes.

[1]  StehfestHarald.  Remark on algorithm 368: Numerical inversion of Laplace transforms , 1970 .

[2]  Shibin Zhang,et al.  Exact simulation of tempered stable Ornstein–Uhlenbeck processes , 2011 .

[3]  J. Blanchet,et al.  On exact sampling of stochastic perpetuities , 2011, Journal of Applied Probability.

[4]  N. Shephard,et al.  Non‐Gaussian Ornstein–Uhlenbeck‐based models and some of their uses in financial economics , 2001 .

[5]  Hanspeter Schmidli.  Distribution of the first ladder height of a stationary risk process perturbed by α-stable Lévy motion , 2001 .

[6]  Jean Bertoin,et al.  Cramér's estimate for Lévy processes , 1994 .

[7]  Ward Whitt,et al.  Numerical Inversion of Laplace Transforms of Probability Distributions , 1995, INFORMS J. Comput..

[8]  M. Kanter,et al.  STABLE DENSITIES UNDER CHANGE OF SCALE AND TOTAL VARIATION INEQUALITIES , 1975 .

[9]  V. Zolotarev.  One-dimensional stable distributions , 1986 .

[10]  Xinsheng Zhang,et al.  Exact Simulation of IG-OU Processes , 2008 .

[11]  A. Dassios,et al.  Exact simulation of gamma-driven Ornstein–Uhlenbeck processes with finite and infinite activity jumps , 2019, J. Oper. Res. Soc..

[12]  A. Papachristodoulou,et al.  Behavioural Economics, Hyperbolic Discounting and Environmental Policy , 2010 .

[13]  Neil Shephard,et al.  Integrated OU Processes and Non‐Gaussian OU‐based Stochastic Volatility Models , 2003 .

[14]  Svetlozar T. Rachev,et al.  Tempered stable and tempered infinitely divisible GARCH models , 2010 .

[15]  W. Vervaat.  On a stochastic difference equation and a representation of non–negative infinitely divisible random variables , 1979, Advances in Applied Probability.

[16]  Yan Qu.  Simulations on Lévy subordinators and Lévy driven contagion models , 2019 .

[17]  Claudia Kluppelberg,et al.  Ruin probabilities and overshoots for general Lévy insurance risk processes , 2004 .

[18]  W. Schoutens,et al.  A risk model driven by Lévy processes , 2003 .

[19]  Ken-iti Sato,et al.  Operator-selfdecomposable distributions as limit distributions of processes of Ornstein-Uhlenbeck type , 1984 .

[20]  O. Barndorff-Nielsen,et al.  Some stationary processes in discrete and continuous time , 1998, Advances in Applied Probability.

[21]  E. Nicolato,et al.  Option Pricing in Stochastic Volatility Models of the Ornstein‐Uhlenbeck type , 2003 .

[22]  S. Asmussen,et al.  Approximations of small jumps of Lévy processes with a view towards simulation , 2001, Journal of Applied Probability.

[23]  Peter D. Sozou,et al.  On hyperbolic discounting and uncertain hazard rates , 1998, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[24]  G. Uhlenbeck,et al.  On the Theory of the Brownian Motion , 1930 .

[25]  S. Wolfe.  On a continuous analogue of the stochastic difference equation Xn = ρ X n–1 + Bn , 1980, Advances in Applied Probability.

[26]  D. P. Gaver,et al.  Observing Stochastic Processes, and Approximate Transform Inversion , 1966, Oper. Res..

[27]  N. Shephard,et al.  Realised power variation and stochastic volatility models , 2003 .

[28]  L. Devroye.  Non-Uniform Random Variate Generation , 1986 .

[29]  George Tauchen,et al.  Simulation Methods for Lévy-Driven Continuous-Time Autoregressive Moving Average (CARMA) Stochastic Volatility Models , 2006 .

[30]  Vadim Linetsky,et al.  TIME‐CHANGED ORNSTEIN–UHLENBECK PROCESSES AND THEIR APPLICATIONS IN COMMODITY DERIVATIVE MODELS , 2012, 1204.3679.

[31]  Jan Kallsen,et al.  PRICING OPTIONS ON VARIANCE IN AFFINE STOCHASTIC VOLATILITY MODELS , 2010 .

[32]  Aleksander Weron,et al.  Stable Lévy motion approximation in collective risk theory , 1997 .

[33]  N. Shephard.  Realized power variation and stochastic volatility models , 2003 .

[34]  Ward Whitt,et al.  A Unified Framework for Numerically Inverting Laplace Transforms , 2006, INFORMS J. Comput..

[35]  Hansjörg Furrer,et al.  Risk processes perturbed by α-stable Lévy motion , 1998 .

[36]  Zhiyi Chi.  On Exact Sampling of Nonnegative Infinitely Divisible Random Variables , 2012, Advances in Applied Probability.

[37]  Luc Devroye,et al.  Random variate generation for exponentially and polynomially tilted stable distributions , 2009, TOMC.

[38]  Haijun Wang.  Robust asset pricing with stochastic hyperbolic discounting , 2017 .

[39]  N. Shephard,et al.  Econometric analysis of realized volatility and its use in estimating stochastic volatility models , 2002 .

[40]  S. Wolfe.  On a continuous analogue of the stochastic difference equation Xn=ρXn-1+Bn , 1982 .

[41]  Ward Whitt,et al.  The Fourier-series method for inverting transforms of probability distributions , 1992, Queueing Syst. Theory Appl..

[42]  Jia Wei Lim,et al.  Exact simulation of generalised Vervaat perpetuities , 2019, Journal of Applied Probability.

[43]  Ole E. Barndorff-Nielsen,et al.  Processes of normal inverse Gaussian type , 1997, Finance Stochastics.

[44]  Charles M. Goldie,et al.  Perpetuities and Random Equations , 1994 .