Strong convergence rates for Cox-Ingersoll-Ross processes - Full parameter range

We study strong (pathwise) approximation of Cox-Ingersoll-Ross processes. We propose a Milstein-type scheme that is suitably truncated close to zero, where the diffusion coefficient fails to be locally Lipschitz continuous. For this scheme we prove positive convergence rates for the full parameter range including the accessible boundary regime. The error criterion is given by the maximal Lp-distance of the solution and its approximation on a compact interval. In the particular case of a squared Bessel process of dimension δ>0 the convergence rate is given by min⁡(1,δ)/(2p).