English

Chi-square simulation of the CIR process and the Heston model

Computational Finance 2012-07-03 v5 Numerical Analysis Probability Statistics Theory Pricing of Securities Statistics Theory

Abstract

The transition probability of a Cox-Ingersoll-Ross process can be represented by a non-central chi-square density. First we prove a new representation for the central chi-square density based on sums of powers of generalized Gaussian random variables. Second we prove Marsaglia's polar method extends to this distribution, providing a simple, exact, robust and efficient acceptance-rejection method for generalized Gaussian sampling and thus central chi-square sampling. Third we derive a simple, high-accuracy, robust and efficient direct inversion method for generalized Gaussian sampling based on the Beasley-Springer-Moro method. Indeed the accuracy of the approximation to the inverse cumulative distribution function is to the tenth decimal place. We then apply our methods to non-central chi-square variance sampling in the Heston model. We focus on the case when the number of degrees of freedom is small and the zero boundary is attracting and attainable, typical in foreign exchange markets. Using the additivity property of the chi-square distribution, our methods apply in all parameter regimes.

Cite

@article{arxiv.0802.4411,
  title  = {Chi-square simulation of the CIR process and the Heston model},
  author = {Simon J. A. Malham and Anke Wiese},
  journal= {arXiv preprint arXiv:0802.4411},
  year   = {2012}
}

Comments

32 pages, 6 figures, 8 tables, update

R2 v1 2026-06-21T10:17:12.273Z