Direct Inversion for the Squared Bessel Process and Applications

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Malham, Simon J. A., Wiese, Anke, Xu, Yifan
Format: Preprint
Published: 2024
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866916538150289408
author Malham, Simon J. A.
Wiese, Anke
Xu, Yifan
author_facet Malham, Simon J. A.
Wiese, Anke
Xu, Yifan
contents In this paper we derive a new direct inversion method to simulate squared Bessel processes. Since the transition probability of these processes can be represented by a non-central chi-square distribution, we construct an efficient and accurate algorithm to simulate non-central chi-square variables. In this method, the dimension of the squared Bessel process, equivalently the degrees of freedom of the chi-square distribution, is treated as a variable. We therefore use a two-dimensional Chebyshev expansion to approximate the inverse function of the central chi-square distribution with one variable being the degrees of freedom. The method is accurate and efficient for any value of degrees of freedom including the computationally challenging case of small values. One advantage of the method is that noncentral chi-square samples can be generated for a whole range of values of degrees of freedom using the same Chebyshev coefficients. The squared Bessel process is a building block for the well-known Cox-Ingersoll-Ross (CIR) processes, which can be generated from squared Bessel processes through time change and linear transformation. Our direct inversion method thus allows the efficient and accurate simulation of these processes, which are used as models in a wide variety of applications.
format Preprint
id arxiv_https___arxiv_org_abs_2412_16655
institution arXiv
publishDate 2024
record_format arxiv
spellingShingle Direct Inversion for the Squared Bessel Process and Applications
Malham, Simon J. A.
Wiese, Anke
Xu, Yifan
Computation
Probability
Computational Finance
Methodology
In this paper we derive a new direct inversion method to simulate squared Bessel processes. Since the transition probability of these processes can be represented by a non-central chi-square distribution, we construct an efficient and accurate algorithm to simulate non-central chi-square variables. In this method, the dimension of the squared Bessel process, equivalently the degrees of freedom of the chi-square distribution, is treated as a variable. We therefore use a two-dimensional Chebyshev expansion to approximate the inverse function of the central chi-square distribution with one variable being the degrees of freedom. The method is accurate and efficient for any value of degrees of freedom including the computationally challenging case of small values. One advantage of the method is that noncentral chi-square samples can be generated for a whole range of values of degrees of freedom using the same Chebyshev coefficients. The squared Bessel process is a building block for the well-known Cox-Ingersoll-Ross (CIR) processes, which can be generated from squared Bessel processes through time change and linear transformation. Our direct inversion method thus allows the efficient and accurate simulation of these processes, which are used as models in a wide variety of applications.
title Direct Inversion for the Squared Bessel Process and Applications
topic Computation
Probability
Computational Finance
Methodology
url https://arxiv.org/abs/2412.16655