Convergence of a splitting method for a general interest rate model

We prove mean-square convergence of a novel numerical method, the tamed-splitting method, for a generalized Ait-Sahalia interest rate model. The method is based on a Lamperti transform, splitting and applying a tamed numerical method for the nonlinearity. The main difficulty in the analysis is caused by the non-globally Lipschitz drift coefficients of the model. We examine the existence, uniqueness of the solution and boundedness of moments for the transformed SDE. We then prove bounded moments and inverses moments for the numerical approximation. The tamed-splitting method is a hybrid method in the sense that a backstop method is invoked to prevent solutions from overshooting zero and becoming negative. We successfully recover the mean-square convergence rate of order one for the tamed-splitting method. In addition we prove that the probability of ever needing the backstop method to prevent a negative value can be made arbitrarily small. In our numerical experiments we compare to other numerical methods in the literature for realistic parameter values.

[1]  Siqing Gan,et al.  First order strong approximation of Ait–Sahalia-type interest rate model with Poisson jumps , 2021, Numerical Algorithms.

[2]  Xuerong Mao,et al.  Truncated EM numerical method for generalised Ait-Sahalia-type interest rate model with delay , 2021, J. Comput. Appl. Math..

[3]  G. Lord,et al.  The role of adaptivity in a numerical method for the Cox-Ingersoll-Ross model , 2020, J. Comput. Appl. Math..

[4]  Yuying Zhao,et al.  On the backward Euler method for a generalized Ait-Sahalia-type rate model with Poisson jumps , 2020, Numer. Algorithms.

[5]  X. Mao,et al.  Generalized Ait-Sahalia-type interest rate model with Poisson jumps and convergence of the numerical approximation , 2018, Physica A: Statistical Mechanics and its Applications.

[6]  G. Karniadakis,et al.  Numerical Methods for Stochastic Partial Differential Equations with White Noise , 2018 .

[7]  Chaman Kumar,et al.  On Milstein approximations with varying coefficients: the case of super-linear diffusion coefficients , 2016, BIT Numerical Mathematics.

[8]  Catherine E. Powell,et al.  An Introduction to Computational Stochastic PDEs , 2014 .

[9]  L. Szpruch,et al.  Numerical simulation of a strongly nonlinear Ait-Sahalia-type interest rate model , 2011 .

[10]  P. Kloeden,et al.  Strong convergence of an explicit numerical method for SDEs with nonglobally Lipschitz continuous coefficients , 2010, 1010.3756.

[11]  Michael B. Giles,et al.  Multilevel Monte Carlo Path Simulation , 2008, Oper. Res..

[12]  Aurélien Alfonsi,et al.  On the discretization schemes for the CIR (and Bessel squared) processes , 2005, Monte Carlo Methods Appl..

[13]  D. Higham,et al.  Convergence of Monte Carlo Simulations involving the Mean-Reverting Square Root Process , 2005 .

[14]  Hengyan Li,et al.  Nonparametric Specification Testing for Continuous-Time Models with Applications to Term Structure of Interest Rates , 2005 .

[15]  A. Gardon.  The Order of Approximations for Solutions of Itô-Type Stochastic Differential Equations with Jumps , 2004 .

[16]  Timothy G. Conley,et al.  Short-term interest rates as subordinated diffusions , 1997 .

[17]  A. Gallant,et al.  ESTIMATION OF CONTINUOUS-TIME MODELS FOR STOCK RETURNS AND INTEREST RATES , 1997, Macroeconomic Dynamics.

[18]  Yacine Ait-Sahalia.  Testing Continuous-Time Models of the Spot Interest Rate , 1995 .

[19]  G. Choe.  Numerical Solution of Stochastic Differential Equations , 2016 .

[20]  Jessica Fuerst,et al.  Stochastic Differential Equations And Applications , 2016 .

[21]  Mireille Bossy,et al.  Euler scheme for SDEs with non-Lipschitz diffusion coefficient: strong convergence , 2006 .

[22]  Mireille Bossy,et al.  An efficient discretisation scheme for one dimensional SDEs with a diffusion coefficient function of the form |x|^a, a in [1/2,1) , 2006 .

[23]  Hongwei Chen,et al.  Tight Bounds for the Normal Distribution: 10611 , 1999 .