Exact simulation of IG-OU processes

IG-OU processes are a subclass of the non-Gaussian processes of Ornstein-Uhlenbeck type, which are important models appearing in financial mathematics and elsewhere. The simulation of these processes is of interest for its applications in statistical inference. In this paper, a stochastic integral of Ornstein-Uhlenbeck type is represented to be the sum of two independent random variables-one has an inverse Gaussian distribution and the other has a compound Poisson distribution. And in distribution, the compound Poisson random variable is equal to a sum of Poisson-distributed number positive random variables, which are independent identically distributed and have a common specified density function. The exact simulation of the IG-OU processes, proceeding from time 0 and going in steps of time interval Delta, is achieved via the representation of the stochastic integral. Comparing to the approximate method, which is based on Rosinski's infinite series representation of the same stochastic integral, by the quantile-quantile plots, the advantage of the exact simulation method is obvious. In addition, as an application, we provide an estimator of the intensity parameter of the IG-OU processes and validate its superiority to another estimator by our exact simulation method.