We propose a quasi-Monte Carlo (qMC) algorithm to simulate variates from the normal inverse Gaussian (NIG) distribution. The algorithm is based on a Monte Carlo technique found in Rydberg [13], and is based on sampling three independent uniform variables. We apply the algorithm to three problems...

In recent years the Dirichlet process prior has experienced a great success in the context of Bayesian mixture modelling. The idea of overcoming discreteness of its realizations by exploiting it in hierarchical models, combined with the development of suitable sampling techniques, represent one...

It is well known that the distributions of assets returns have heavier tails than the Gaussian's. To capture such a distributional characteristic, the Generalized Hyperbolic(GH) distribution and its subclasses have been applied to assets returns as the distribution with heavier tails. GH...