Abstract:The inverse Gaussian distribution is a useful distribution with important applications. But there is less discussion in the literature on sampling of this distribution. The method given in [Atkinson, A.C., 1982. The simulation of generalized inverse Gaussian and hyperbolic random variables. SIAM Journal on Scientific and Statistical Computing 3(4), 502-515] is based on rejection method where some (uniform) random numbers from the sample are discarded. This feature makes it difficult to take advantage of the low discrepancy sequences which have important applications. In [Michael, J., Schucany, W., Haas, R., 1976. Generating random variates using transformations with multiple roots. The American Statistician 30(2), 88-90], Michael et al. give a method to generate random variables with inverse Gaussian distribution. In their method, two pseudorandom numbers uniformly distributed on (0, 1) are needed in order to generate one inverse Gaussian random variable. In this short paper, we present a new method, based on direct approximate inversion, to generate the inverse Gaussian random variables. In this method, only one pseudorandom number is needed in generating one inverse Gaussian random variate. This method enables us to make use of the better convergence of low discrepancy sequence than the pseudorandom sequence. Numerical results show the superiority of low discrepancy sequence than the pseudorandom sequence in simulating the mean of the inverse Gaussian distribution by using our sampling method. Further application of this method in exotic option pricing under the normal inverse Gaussian model is under investigation.

Inverse Gaussian Process regression for likelihood-free inference

Likelihood approximations via Gaussian approximate inference

Likelihood-Free Inference with Deep Gaussian Processes

Constrained Gaussian Process for Signal Integrity applications using Variational Inference

Introduction To Gaussian Process Regression In Bayesian Inverse Problems, With New ResultsOn Experimental Design For Weighted Error Measures

Continuous Gaussian mixture solution for linear Bayesian inversion with application to Laplace priors

Improving Hyperparameter Learning under Approximate Inference in Gaussian Process Models

An adaptive Gaussian process method for multi-modal Bayesian inverse problems

Adaptive Gaussian Process Regression for Bayesian inverse problems

Laplace approximation for logistic Gaussian process density estimation and regression

Saddlepoint approximation for the generalized inverse Gaussian Lévy process

Fast, accurate and lightweight sequential simulation-based inference using Gaussian locally linear mappings

Pseudo-marginal Bayesian inference for Gaussian process latent variable models

Composite Inference for Gaussian Processes

Generating Inverse Gaussian Random Variates by Approximation

Likelihood Based Inference for the Shape Parameter of the Inverse Gaussian Distribution

An efficient likelihood-free Bayesian inference method based on sequential neural posterior estimation

An adaptive mixture-population Monte Carlo method for likelihood-free inference

Amortized Variational Inference for Deep Gaussian Processes

Preconditioned Neural Posterior Estimation for Likelihood-free Inference