1. Introduction and main results
1.1. Introduction
Consider the following ultraspherical Jacobi process
where $W=\{W_t\colon t\geq0\}$ is a standard Brownian motion with unknown drift parameters $b \in (\!-\!\infty,-1)$ . This process is well defined under the assumption $b<-1$ , i.e. we have $|X_t|<1$ for all $t\geq 0$ . The Jacobi process, also called Wright–Fisher diffusion, was originally used to model gene frequencies [Reference Ethier and Kurtz9, Reference Karlin and Taylor16]. More recently, it has also been applied to describe financial factors. For example, [Reference Delbane and Shirakawa7] models interest rates by the Jacobi process and studies moment-based techniques for pricing bonds. Moreover, this process has also been applied to model stochastic correlation matrices [Reference Ahdida and Alfonsi2] and credit default swap indexes [Reference Bernis and Scotti5]. For the multivariate case, see [Reference Ackerer, Filipović and Pulido1, Reference Gouriéroux and Jasiak13].
For $b \in (\!-\!\infty,-1)$ , the Jacobi process (1.1) has stationary distribution $\textrm{Beta}(\!-\!b,-b)$ , i.e. the Beta distribution with shape parameter $-b$ and scale parameter $-b$ . Here, b also represents the mean-reverting parameter. For practical applications, it is crucial to construct the estimator for b and study the corresponding asymptotic properties. Let $\mathbb{P}_b$ be the probability distribution of the solution of (1.1) on $C(\mathbb{R}_+,\mathbb{R})$ . Under $\mathbb{P}_b$ , the maximum likelihood estimator of b is given by
The Jacobi process was subordinated in [Reference Demni and Zani8] by the method of random time change, and the corresponding semi-group density was obtained. Together with the Girsanov formula technique [Reference Bercu and Richou4, Reference Florens-Landais and Pham12], they gave the large deviations for $\widehat b_T$ . For the Girsanov formula technique, see also [Reference Zang and Zhang20]. Moreover, [Reference Jiang, Gao and Zhao14] studied the moderate deviations for $\widehat b_T$ .
From [Reference Kutoyants17], it follows that
which immediately implies that
as $T\to\infty$ . The Berry–Esseen bound for (1.3) bounds the global convergence speed which is uniform in all $x\in\mathbb{R}$ . In this paper, our motivation is to study the nonuniform convergence rate of the difference
which will depend on x as well as T. This local bound provides more accurate information than the Berry–Esseen bound; see Corollaries 1.1 and 1.2 and Remark 1.1. Our approach is based on the change of measure method [Reference Wang and Jing19] and precise asymptotic analysis techniques [Reference Jiang and Zhou15]. The explicit calculations and estimations for the Laplace functionals of $\int_{0}^{T} {X_t}/({1-X_t^2}) \,\textrm{d}{X_t}$ and $\int_{0}^{T} {X_t^2}/({1-X_t^2}) \,\textrm{d}{t}$ play crucial roles. Moreover, the techniques used here could potentially be used for a wider scope of diffusions that have a structure similar to our model.
1.2. Main results
We now state the main results of this paper, i.e. the exponential nonuniform Berry–Esseen bound and its application for $\widehat b_T$ .
Theorem 1.1. There exists some positive constant $C>0$ depending only on b such that, for T large enough and any $\rho>0$ with $|x|<\rho T^{{1}/{6}}$ ,
where the O(1) term only depends on b.
As applications of the above nonuniform Berry–Esseen bound, we can obtain the following optimal uniform Berry–Esseen bound and optimal Cramér-type moderate deviations for $\widehat b_T$ .
Corollary 1.1. There exists some positive constant $C>0$ depending only on b such that, for T large enough,
Corollary 1.2. Let $\{\lambda_T,T>0\}$ be a family of positive numbers satisfying $\lambda_T/T^{{1}/{6}}\rightarrow 0$ , $T\rightarrow \infty$ . Then, for any $\rho\geq0$ ,
Remark 1.1. Note that, for any $x>0$ ,
Hence, when $x\geq(\log T)^{1/2}$ , $1-\Phi(x)=o\big(T^{-1/2}\big)$ . Moreover, using [Reference Jiang, Gao and Zhao14, Theorem 1.1], for $x\geq(\log T)^{1/2}$ ,
Then,
Therefore, the uniform Berry–Esseen bound cannot characterize the error
for large x depending on T. Our nonuniform Berry–Esseen bound (1.4) obtained in Theorem 1.1) can fill this gap.
Remark 1.2. When T takes values in positive integers, letting
we have
where $\mathcal{F}_{i}=\sigma(W_t, t\leq i)$ . Then $\{(\xi_i,\mathcal{F}_i), i\in\mathbb{N}\}$ is a sequence of martingale differences, and $\{\widehat{b}_T-b, T>0\}$ can be viewed as self-normalized martingales. The exponential nonuniform Berry–Esseen bound in Theorem 1.1 is a parallel result to the self-normalized sum of independent variables [Reference Wang and Jing19].
Moreover, it is worth noting that [Reference Fan and Shao11] obtained an exponential nonuniform Berry–Esseen bound and Cramér-type moderate deviation for self-normalized martingales. However, for the martingale differences $\{(\xi_i,\mathcal{F}_i), i\in\mathbb{N}\}$ it is difficult to verify the Bernstein condition in [Reference Fan and Shao11, (A1)]. So, our results cannot be covered or deduced directly by [Reference Fan and Shao11]. For more details on this topic, see [Reference Fan, Ion, Liu and Shao10] and the references therein.
The rest of this paper is organized as follows. In Section 2, we first give the explicit calculation of Laplace functionals related to a Jacobi process, which plays an important role in our asymptotic analysis. The proofs of the main results and their corollaries then follow in Section 3 and Appendix A.
2. Explicit calculation of Laplace functionals related to a Jacobi process
Recall from (1.2) that we can write
Observe that, for $x>0$ ,
and similarly for $x<0$ :
Set
then,
Write
and define a new probability measure $\mathbb{Q}$ which is absolutely continuous with respect to $\mathbb{P}_b$ , with Radon–Nikodym density
Now, we can obtain, for $x>0$ ,
and, for $x<0$ ,
To prove Theorem 1.1, we need to analyze the Laplace functionals of $G_{T,x}$ under $\mathbb{P}_b$ and $\mathbb{Q}$ explicitly. Recall that if $\{X_t, t\geq 0\}$ is defined by (1.1), then its transition semi-group density is given by [Reference Demni and Zani8, p. 526]
where $*$ is the convolution of two functions and
with
The following lemma plays an important role in our analysis.
Lemma 2.1. For any $u\in\mathbb{R}$ and $h_{T,x}$ as defined in (2.2),
holds for all $T>0$ . Here, $\textrm{i}=\sqrt{-1}$ and
Proof. From the definition of $G_{T,x}$ (2.1), it follows that
where $\lambda_{T,x}(u)$ and $\mu_{T,x}(u)$ are defined by
Let us define the joint log-Laplace transformation
Denote the Brownian filtration by $\mathcal{F}_t=\sigma(W_s,s \leq t)$ ; then, by virtue of the Girsanov theorem, we have, for any $\textrm{Re}(b_0)<-1$ ,
Now, choose $b_0= -1 -\sqrt{(b+1)^2 +2(\lambda_{T,x}(u) +\mu_{T,x}(u))}\;:\!=\;c_{T,x}(u)\;:\!=\;c$ ; then
Applying Itô’s formula,
and, as a consequence,
Using the transition density function of $X_{T}$ [Reference Demni and Zani8], see also (2.6), we obtain
where B stands for the Beta function and
Recalling the definition of $\mathcal{R}_{T,x}(u)$ in the above lemma we can write
where
In order to derive the precise estimation of $\mathcal{R}_{T,x}(u)$ , the following expansion of Gamma functions plays a crucial role.
Lemma 2.2. Let $\Gamma(z)$ be the Gamma function defined on the complex plane $\mathbb{C}$ . Then
where $z_1,z_2\in \mathbb{C}$ , $\textrm{Re}(z_1), \textrm{Re}(z_2)>0$ , and $\gamma=0.5772\ldots$ is the Euler constant. Moreover, when $z_1-z_2=c\in\mathbb{R_+}$ and $\textrm{Re}(z_2)>1$ ,
Proof. By Euler’s formula [Reference Ahlfors3, p. 199] we have, for all $z\in\mathbb{C}$ ,
where $\gamma=0.5772\ldots$ is the Euler constant. As a consequence,
Moreover, when $z_1-z_2=c\in\mathbb{R_+}$ and $\textrm{Re}(z_2)>0$ , we have
Thus, we obtain the desired estimate,
Define the Laplace functional of $G_{T,x}$ under $\mathbb{Q}$ as
In the following propositions, we give explicit asymptotic expansions of $\varphi_{\mathbb{Q},G_{T,x}}({u}/{\sqrt{T}})$ for different ranges of u.
Proposition 2.1. For any constant $\rho>0$ , and $|u|\leq\rho T^{{1}/{6}}$ , $|x|\leq\rho T^{{1}/{6}}$ ,
where
Proof. By using (2.3) and (2.7), we can write
and
On the one hand, by Taylor’s expansion, straightforward but tedious calculations will give
Consequently, we find that
On the other hand, by Proposition A.1,
Therefore, combining (2.10), (2.15), and (A.1), we get (2.8). Moreover, (2.9) can be obtained via (2.11), (2.14), and (A.1).
Proposition 2.2. For any constant $\rho>0$ , and $|u|\leq({|b|}/{32})T^{{1}/{2}}$ , $|x|\leq\rho T^{{1}/{6}}$ ,
where $\gamma$ is the Euler constant.
Proof. Recall the calculation in (2.11); it suffices to analyze the real part of $-(c_{T,x}(0)-c_{T,x}(u))T/2$ and give an upper bound for $|\exp(\mathcal{R}_{T,x}(u))|$ . Using
we have
Now, for T large enough and $|u|\leq({|b|}/{32})T^{{1}/{2}}$ , $|x|\leq\rho T^{{1}/{6}}$ ,
Consequently,
which immediately implies that
Next, we turn to estimating $|\exp(\mathcal{R}_{T,x}(u))|$ . From (2.17),
Together with (2.12) and (A.2)–(A.4), we have
By Lemma 2.2, we have
which, together with (A.8), implies that
Finally, by (2.11), (2.18), and (2.19), we can complete the proof of this proposition.
3. Proofs of the main results
In this section we give the proof of our main result, Theorem 1.1, and its two corollaries. Recall that the definitions of $G_{T,x}$ and the probability measure $\mathbb{Q}$ are given in (2.1) and (2.3) respectively. Consider the normalized version of $G_{T,x}$ ,
Using the simple identity $\int_{0}^{\infty}\exp\!{\{-xu\}}\,\textrm{d}{\Phi(u)} = \textrm{e}^{{x^2}/{2}}(1-\Phi(x))$ and (2.4), we can write, for $x>0$ ,
and for $x<0$ , using (2.5),
where $h_{T,x}=x\sqrt{-{2(b+1)}/{T}}$ and
Next, we need to give the estimation of $\mathcal{J}_{T,x}$ and $\widehat{\mathcal{J}}_{T,x}$ . The following lemma [Reference Petrov18, Theorem 2, p. 109] will play a crucial role.
Lemma 3.1. Assume F and G are probability distribution functions, and that G has bounded derivative. Define $\varphi(u)=\int_{\mathbb{R}}\textrm{e}^{\textrm{i}ux}\,\textrm{d}{F(x)}$ and $\psi(u)=\int_{\mathbb{R}}\textrm{e}^{\textrm{i}ux}\,\textrm{d}{G(x)}$ . Then, for any $M>0$ ,
Proof of Theorem 1.1. Without of loss of generality, we assume that $x>0$ ; the proof for $x<0$ is similar. $\mathcal{J}_{T,x}$ , as defined by (3.2), can be estimated as follows. By using Lemma 3.1, we have
Moreover, we have
From Proposition 2.1, it follows that, for $|u|\leq T^{{1}/{6}}$ ,
Together with Proposition 2.2, we have
Therefore, we obtain
Applying (2.7), (2.8), (3.1), and (3.3),
Therefore, we get
which completes the proof of (1.4).
Note that
Together with (1.4), we can get (1.5), i.e.
Proof of Corollary 1.1. According to Theorem 1.1, for any $\rho>0$ ,
Moreover, using [Reference Jiang, Gao and Zhao14, Theorem 1.1] and (3.4), we have
By (3.5) and (3.6), we can complete the proof of Corollary 1.1.
Appendix A. Estimations of $\mathcal{R}_{\textbf{\textit{T,x}}}(\textbf{u})$ and $\mathcal{R}_{\textbf{\textit{T,x}}}(\textbf{0})$
In this appendix, we give precise estimates of $\mathcal{R}_{T,x}(u)$ and $\mathcal{R}_{T,x}(0)$ , which play crucial roles in the proof of Proposition 2.1.
Proposition A.1. For any constant $\rho>0$ , and $|u|\leq\rho T^{{1}/{6}}$ , $|x|\leq\rho T^{{1}/{6}}$ ,
Proof. From (2.12) and (2.13), we have
By Lemma 2.2, we have the following asymptotic results:
Similarly, we also have
Recall the Jacobi Theta function [Reference Biane, Pitman and Yor6]:
We have
Moreover, from the fact that
it follows that
Recall the fact that $\Theta(x)=(1/\sqrt{x})\Theta (1/x)$ , which yields
By the asymptotic expansion of the Jacobi Theta function,
we obtain
Consequently, we have
Using Lemma 2.2,
so we obtain
and the small $o(\cdot)$ term is hold uniformly in n; this further implies that
Similarly, we have
Applying (A.5) and (A.7), we get
Acknowledgements
We would like to express our great gratitude to the two anonymous referees and the AE for their careful reading and insightful comments, which surely led to an improved presentation of this paper.
Funding information
Hui Jiang is supported by the Natural Science Foundation of Jiangsu Province of China (No. BK20231435) and Fundamental Research Funds for the Central Universities (No. NS2022069). Shaochen Wang is partially supported by Guangdong Basic and Applied Basic Research Foundation (No. 2023A1515012125).
Competing interests
There were no competing interests to declare which arose during the preparation or publication process of this article.