Point process simulation of Generalised inverse Gaussian processes and estimation of the Jaeger Integral

by   Simon Godsill, et al.

In this paper novel simulation methods are provided for the Generalised inverse Gaussian (GIG) Lévy process. Such processes are intractable for simulation except in certain special edge cases, since the Lévy density associated with the GIG process is expressed as an integral involving certain Bessel Functions, known as the Jaeger Integral in diffusive transport applications. We here show for the first time how to solve the problem indirectly, using generalised shot-noise methods to simulate the underlying point processes and constructing an auxiliary variables approach that avoids any direct calculation of the integrals involved. The augmented bivariate process is still intractable and so we propose a novel thinning method based on upper bounds on the intractable integrand. Moreover our approach leads to lower and upper bounds on the Jaeger integral itself, which may be compared with other approximation methods. We note that the GIG process is the required Brownian motion subordinator for the generalised hyperbolic (GH) Lévy process and so our simulation approach will straightforwardly extend also to the simulation of these intractable proceses. Our new methods will find application in forward simulation of processes of GIG and GH type, in financial and engineering data, for example, as well as inference for states and parameters of stochastic processes driven by GIG and GH Lévy processes.


page 1

page 2

page 3

page 4


The Berry-Esséen Upper Bounds of Vasicek Model Estimators

The Berry-Esséen upper bounds of moment estimators and least squares est...

Determinantal shot noise Cox processes

We present a new class of cluster point process models, which we call de...

Lower and upper bounds of the superposition of renewal processes and extensions

Consider a system consisting of multiple sockets into each of which a co...

Mean-Field Variational Inference for Gradient Matching with Gaussian Processes

Gradient matching with Gaussian processes is a promising tool for learni...

Stochastic Process Bandits: Upper Confidence Bounds Algorithms via Generic Chaining

The paper considers the problem of global optimization in the setup of s...

Composite Inference for Gaussian Processes

Large-scale Gaussian process models are becoming increasingly important ...

Hawkes processes as competing hazards models and a simulation algorithm

In this paper, we establish a competing hazards interpretation of Hawkes...

Code Repositories


Point process simulation of Generalised inverse-Gaussian processes

view repo