The Esscher Transform of the Densities of a Symmetric NIG Lévy Process
Requires a Wolfram Notebook System
Interact on desktop, mobile and cloud with the free Wolfram Player or other Wolfram Language products.
This Demonstration illustrates the effect of applying the Esscher transform to one-dimensional densities when the logarithm of the price process is the so-called Normal Inverse Gaussian (NIG) Lévy process. The distribution of this process is determined by four parameters , , , and : steepness, asymmetry, scale, and location. The location parameter is set to 0 (in other words, the process is centered), which makes the implementation of the Esscher transfer particularly simple (while at the same time being quite acceptable, on empirical grounds, for modeling stock prices). We display on the same graph the density of the logarithm of the price process at a given time , and the density of the process obtained by applying to it the Esscher transform. It is the latter density that is used in computing option prices.[more]
In addition to the three model parameters, you can vary the interest rate, the moment in time for which the density of the log of the price is shown, and the range of values over which the densities are shown. You can identify the graphs by hovering over them.[less]
Contributed by: Andrzej Kozlowski (March 2011)
Open content licensed under CC BY-NC-SA
The fundamental theorem of asset pricing (e.g., ) states (loosely speaking) that the arbitrage-free prices of "contingent claims" (e.g., options) on an asset (e.g., a stock) whose price follows a (semi-martingale) stochastic process can be computed as expected values of the discounted payoff, but not (as was once thought) under the "real world measure", but under a measure equivalent to it and having the property that the discounted price process of the asset is a local martingale with respect to it. In the case of the classical Black–Scholes model, this "equivalent martingale measure" is unique and is obtained by changing the "drift" of the logarithm of the price process (the Girsanov theorem). However, in more realistic models, such as those involving jumps or stochastic volatility, there are usually infinitely many equivalent martingale measures and one has to choose one based on economic considerations. One such measure that is frequently used is obtained by means of the so-called Esscher transform, which had a long history of being used in actuarial mathematics before being introduced to option pricing by Gerber and Shiu . For a Lévy process, the Esscher transform can be defined by transforming the one-dimensional marginal densities in the manner of . The transformed marginal densities depend on a single parameter and are related to the original marginal densities by the relation . By writing down the condition that the transformed process is a martingale, you can obtain an equation for the parameter . In the case of the centered (i.e., when the location parameter ) NIG process, considered in this Demonstration, the equation takes a particularly simple form (see , p. 28) and the marginal densities of the transformed process are also centered NIG densities. We note that, as shown on page 4 of , for certain parameter values an exponential NIG process may not have an Esscher transform. In such a case only the graph of the real-world marginal PDF will be displayed.
 F. Delbaen and W. Schachermayer, The Mathematics of Arbitrage, Berlin-Heidelberg-New York: Springer, 2006.
 H. U. Gerber and E. S. W. Shiu, "Option Pricing by Esscher Transforms," Transactions of the Society of Actuaries, 46, 1994 pp. 99–191.
 K. Prause, "The Generalized Hyperbolic Model: Estimation, Financial Derivatives, and Risk Measures," PhD thesis, University of Freiburg, 1999.
 S. Raible, "Lévy Processes in Finance: Theory, Numerics, and Empirical Facts," PhD thesis, University of Freiburg, 2000.