If you see this, something is wrong
We investigate the problem of estimating the drift parameter from \( N\) independent copies of the solution of a stochastic differential equation driven by a multiplicative fractional Brownian noise with Hurst parameter \( H\in (1/3,1)\) . Building on a least-squares-type object involving the Skorokhod integral, a key challenge consists in approximating this unobservable quantity with a computable fixed-point estimator, which requires addressing the correction induced by replacing the Skorokhod integral with its pathwise counterpart. To this end, a crucial technical contribution of this work is the reformulation of the Malliavin derivative of the process in a way that does not depend explicitly on the driving noise, enabling control of the approximation error in the multiplicative setting. For the case \( H\in (1/3,1/2]\) , we further exploit results on two-dimensional Young integrals to manage the more intricate correction term that appears. As a result, we establish the well-posedness of a fixed-point estimator for any \( H\in (1/3,1)\) , together with both an asymptotic confidence interval and a non-asymptotic risk bound. Finally, a numerical study illustrates the good practical performance of the proposed estimator.
Assumption 1
The function \( \sigma\) is bounded and \( \inf_{\mathbb R}|\sigma| > 0\) .
Proposition 1
Under Assumption 1, for every \( s,t\in [0,T]\) ,
Proposition 2
Consider
Assume that \( H\in (1/3,1/2]\) , and let \( x : [0,T]^2\rightarrow\mathbb R\) be a function such that
(4)
and
(5)
where \( \alpha\in (1 - 2H,H)\) and \( \mathfrak c_x\) is a positive constant. Then, the 2D Young integral of \( x\) with respect to the covariance function \( R\) of the fBm is well-defined, and
Proposition 3
Assume that \( H\in (1/3,1/2]\) . Under Assumption 1, if \( b\) is bounded, then
(6)
and
(7)
Remark 1
Let us add a few remarks on Proposition 3.
It is worth noting that, after application of Proposition 2, the relationships between Skorokhod and pathwise integrals in the two cases \( H > 1/2\) and \( H\in (1/3,1/2]\) , summarized in Equations (3) and (7) respectively, appear quite similar. However, there is a crucial difference: in Equation (7), the right-hand side cannot be simplified in such a way that the second integral cancels the \( -1\) appearing in the third term, as happens in Equation (3). This is because the map
does not satisfy the assumptions of Proposition 2.
Proposition 4
Assume that \( b'\) , \( \varphi\) and \( \psi\) are nonpositive. Under Assumption 1, if
(10)
where \( \mathfrak c\) is a deterministic constant arbitrarily chosen in \( (0,1)\) ,
then \( \Theta_N\) is a contraction from \( \mathbb R_+\) into itself. Therefore, \( R_N\) exists and is unique.
Remark 2
Note that, in particular, the conditions of Proposition 4 on \( (\varphi,\psi)\) imply that it is bounded. Indeed,
leading to
Since \( b'\) and \( \sigma\) are bounded, the ratio \( \sigma'b/\sigma\) is also bounded, and therefore so are \( (\varphi,\psi)\) .
Example 1
Let us provide examples of drift and volatility functions satisfying the condition \( {\rm (A)}\) :
Assume that \( b(x) = -x\) , which is quite common. Then,
For instance,
If \( \sigma(x) =\pi +\arctan(x)\) , then \( \sigma'(x) = (1 + x^2)^{-1}\) , leading to
If \( \sigma(x) = 1 + e^{-x^2}\) , then \( \sigma'(x) = -2xe^{-x^2}\) , leading to
Proposition 5
Assume that \( b'\) , \( \varphi\) and \( \psi\) are nonpositive, \( \theta_0 > 0\) and that
(12)
Under Assumption 1, there exists a constant \( \mathfrak c_{5} > 0\) , not depending on \( N\) , such that
Proposition 6
Assume that \( b'\) , \( \varphi\) and \( \psi\) are nonpositive, and that \( \theta_0 > 0\) . Under Assumption 1, if \( T\) satisfies the condition (12) with \( \mathfrak c = 1/2\) , then
for every \( \alpha\in (0,1)\) , where \( u_{\cdot} =\phi^{-1}(\cdot)\) , \( \phi\) is the standard normal distribution function, and
Proposition 7
Assume that \( b'\) , \( \varphi\) and \( \psi\) are nonpositive, and that \( \theta_0 > 0\) . Under Assumption 1, if \( T\) satisfies the condition (12), then there exists a constant \( \mathfrak c_{7} > 0\) , not depending on \( N\) , such that
Proposition 8
Assume that \( b\) is bounded, and that \( b'\) , \( \varphi\) and \( \psi\) are nonpositive. Under Assumption 1, if \( M_N/D_N\) satisfies the condition (10), then \( \widetilde\Theta_N\) is a contraction from \( \mathbb R_+\) into itself. Therefore, \( R_N\) exists and is unique.
Proposition 9
Assume that \( b\) is bounded, \( b'\) , \( \varphi\) and \( \psi\) are nonpositive, and that \( \theta_0\leqslant\theta_{\max}\) with a known \( \theta_{\max} > 0\) . Under Assumption 1, if \( T\) satisfies the condition (12) with \( \mathfrak c = 1/2\) , then
where
Remark 3
Comparing the confidence interval in Proposition 6 with that of Proposition 9, one can notice a difference in the random variables appearing in the widths of the intervals: \( Y_N\) is now replaced by \( \mathfrak Y_N\) . This discrepancy arises because, in the case \( H > 1/2\) , when analyzing the variance of the Skorokhod integral, we could invoke Theorem 3.11.1 in [67], which allows us to propose a Lebesgue integral bounded by \( Y_N\) and converging in probability to the variance of the Skorokhod integral.
For \( H\in (1/3,1/2]\) , this theorem is no longer applicable. Instead, we handle the variance of the Skorokhod integral by directly using the relationship between the Skorokhod and pathwise integrals, as established through our probabilistic tools in Equality (7).
Proposition 10
Assume that \( b\) is bounded, and that \( b'\) , \( \varphi\) and \( \psi\) are nonpositive. Under Assumption 1, if \( T\) satisfies the condition (12), then there exists a constant \( \mathfrak c_{10} > 0\) , not depending on \( N\) , such that
Lemma 1
For every \( s,t\in [0,T]\) , consider
Under the assumptions of Proposition 3,
There exists a deterministic constant \( \mathfrak c_{1} > 0\) such that
and for every \( \alpha\in (0,H)\) ,
where \( \|X\|_{\alpha,T}\) is the \( \alpha\) -H
specialChar{34}older norm of \( X\) over \( [0,T]\) .
[1] Marie, N. (2025). On a Computable Skorokhod's Integral Based Estimator of the Drift Parameter in Fractional SDE. Scandinavian Journal of Statistics 52(1), 1-37.
[2] Biagini, F., Hu, Y., Oksendal, B. & Zhang, T. (2008). Stochastic Calculus for Fractional Brownian Motion and Applications. Springer, London.