The Normal Inverse Gaussian Lévy Process

The demonstration shows a path of the Normal Inverse Gaussian (NIG) Lévy process and the graph of the probability density of the process at various moments in time. The NIG process is a pure jump Lévy process with infinite variation, which has been used successfully in modeling the distribution of stock returns on the German and Danish exchanges. The version of the model represented here is controlled by three parameters, which arise from the realization of the process as a time changed Brownian motion with drift. The parameters are the drift and the volatility of the Brownian process and the variance of the (inverse Gaussian) subordinator (whose expectation is assumed to be 1). In the limiting case when the variance of the subordinator is set to zero the NIG process coincides with Brownian motion and the probability density is normal. For other values of the variance the NIG probability density has non-zero excess kurtosis and skewness, displayed beneath the graph on the RHS.
Note that a part of the path represented on the LHS graph may sometimes disappear from view, in which case one should adjust the range of values parameter.
  • Contributed by: Piotr Teklinski, Piotr Wroblewski, Andrzej Kozlowski

The Normal Inverse Gaussian distribution and associated stochastic processes was introduced by Barndorff-Nielsen in 1) and 2). The name derives from its representation as the distribution of Brownian motion with drift time changed by the inverse Gaussian Lévy process. Since the distribution is infinitely divisible it gives rise to a corresponding Lévy process, which has been shown capable of accurately modelling the returns on a number of assets on German, Danish and US exchanges, see 4).
The Normal Inverse Gaussian Lévy process is in many ways similar to the Variance Gamma process due to Madan and Seneta. Both belong to the family of Lévy processes of the Generalized Hyperbolic type, however they posses unique properties that make them particularly tractable and convenient for option pricing. In particular alone among hyperbolic family they have the property of being closed under convolution (i.e. the sum of independently distributed variables of the given type has the same type). Both processes can be represented as a time change of a Brownian motion with drift by a Lévy process with increasing increments. In the case of the Variance Gamma process the time change process is the Gamma process, in the case of a NIG process it is an Inverse Gamma Process. Both processes are pure jump Lévy processes (they have no continuous Brownian component) but they differ in the nature of jumps: the Variance Gamma process has jumps of finite variation while the variation of the NIG process is infinite.
There are several common parametrization of the NIG process. Most of them use 4 parameters. Here, however, we use a 3-parameter representation given in 4). The parameters are the drift and volatility of the time changed Brownian process and the variance of the time change process, which is assumed to have expectation 1. (This does not loose generality due to scaling properties for the NIG process). This representation brings out clearly the similarity between the NIG process and the Variance Gamma process; in particular when the variance of the time change tends to 0 both processes approximate the Brownian motion process.
1.) O. E. Barndorff-Nielsen, "Normal inverse Gaussian distributions and stochastic
volatility modelling". Scandinavian Journal of Statistics 24,1997 pp. 1–13.
2). O.E. Barndorff-Nielsen. "Processes of normal inverse Gaussian type", Finance &
2, 1998 pp. 41–68.
3). R. Cont, P.Tankov, Financial Modelling with Jump Processes, CRC Press, 2004
4). T. H. Rydberg, : The Normal Inverse Gaussian Lévy Process: Simulation and Approximation, Comm. Stat.: Stoch. Models 13, 1997 pp. 887–910.
Mathematica »
The #1 tool for creating Demonstrations
and anything technical.
Wolfram|Alpha »
Explore anything with the first
computational knowledge engine.
MathWorld »
The web's most extensive
mathematics resource.
Course Assistant Apps »
An app for every course—
right in the palm of your hand.
Wolfram Blog »
Read our views on math,
science, and technology.
Computable Document Format »
The format that makes Demonstrations
(and any information) easy to share and
interact with.
STEM Initiative »
Programs & resources for
educators, schools & students.
Computerbasedmath.org »
Join the initiative for modernizing
math education.
Step-by-step Solutions »
Walk through homework problems one step at a time, with hints to help along the way.
Wolfram Problem Generator »
Unlimited random practice problems and answers with built-in Step-by-step solutions. Practice online or make a printable study sheet.
Wolfram Language »
Knowledge-based programming for everyone.
Powered by Wolfram Mathematica © 2014 Wolfram Demonstrations Project & Contributors  |  Terms of Use  |  Privacy Policy  |  RSS Give us your feedback
Note: To run this Demonstration you need Mathematica 7+ or the free Mathematica Player 7EX
Download or upgrade to Mathematica Player 7EX
I already have Mathematica Player or Mathematica 7+