Drift estimation for rough processes under small noise asymptotic: QMLE approach

Basic Information

• Paper ID: 2510.09028
• Title: Drift estimation for rough processes under small noise asymptotic: QMLE approach
• Authors: Arnaud Gloter (Université d'Evry), Nakahiro Yoshida (University of Tokyo)
• Classification: math.ST (Statistics Theory)
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.09028

Abstract

This paper investigates the drift estimation problem for solutions X^ε of stochastic Volterra equations with unknown parameter θ⋆. The Volterra kernel is singular with the form K(u) = cu^(α-1)𝟙_{u>0}, where α ∈ (1/2, 1), and the diffusion coefficient is proportional to ε → 0. Based on discrete sampling observations of the Volterra process (grid h → 0), a quasi-maximum likelihood estimator (QMLE) is constructed. The main step involves evaluating the error generated when reconstructing a semimartingale path through inversion of the Volterra kernel, proving that this error decreases at rate h^(1/2), independent of the value of α. An explicit contrast function is then introduced to obtain an efficient estimator as ε → 0.

Research Background and Motivation

1. Core Problem: This paper addresses statistical inference for drift parameters of rough stochastic processes based on discrete observation data within the small noise asymptotic framework.
2. Problem Significance:
   • Fractional differential equations are widely applied in physics, mathematical finance, and life sciences modeling with memory properties
   • Rough processes (α - 1/2 < 1/2) are extensively studied in financial asset volatility modeling
   • In practical applications, only discrete observation data are available; continuous observation is infeasible in practice
3. Limitations of Existing Methods:
   • MLE methods under continuous observation cannot be implemented in practice
   • Existing trajectory fitting estimators (TFE) are suboptimal
   • The Hölder regularity of rough processes (α - 1/2 < 1/2) leads to large direct discretization errors
4. Research Motivation: Develop effective estimation methods based on discrete observations, overcome technical difficulties arising from roughness, and achieve efficient parameter estimation under small noise asymptotics.

Core Contributions

1. Proposed a novel QMLE method: A quasi-maximum likelihood estimator based on Volterra kernel inversion, applicable to rough processes with discrete observations
2. Established key approximation error bounds: Proved that reconstruction error decreases at rate h^(1/2), a result valid for all α ∈ (1/2, 1)
3. Provided complete asymptotic theory: Including consistency and asymptotic normality of the estimator
4. Derived optimal convergence rates: The estimator converges at rate ε^(-1), with asymptotic variance achieving the inverse of Fisher information
5. Conducted numerical verification: Monte Carlo experiments validate the theoretical results

Methodology Details

Problem Formulation

Consider a d-dimensional stochastic Volterra equation:

X^ε_t = X^ε_0 + ε∫₀ᵗ K(t-s)a(X^ε_s)dBs + ∫₀ᵗ K(t-s)b(X^ε_s, θ⋆)ds

where:

• K(u) = u^(α-1)/Γ(α) 𝟙_{u>0}, α ∈ (1/2, 1) is a singular kernel
• θ⋆ ∈ Θ ⊂ ℝ^(d_Θ) is the unknown drift parameter
• ε → 0 is the small noise parameter
• Observation data: {X^ε_}_, h → 0

Objective: Estimate parameter θ⋆ based on discrete observations

Model Architecture

1. Core Idea: Semimartingale Reconstruction

Reconstruct the semimartingale process through Volterra kernel inversion:

Z^ε_t := ε∫₀ᵗ a(X^ε_s)dBs + ∫₀ᵗ b(X^ε_s, θ⋆)ds

Using the resolvent kernel L(u) = u^(-α)/Γ(1-α)𝟙_{u>0}, satisfying L⋆K = 1.

2. Discrete Reconstruction Process

Define the discrete observation process X^(h)t = X{φ_h(t)}, where φ_h(t) = h⌊t/h⌋, and construct:

Z^(ε,h)_t = ∫₀ᵗ L(t-s)(X^(h)_s - X^(h)_0)ds

3. Estimator Construction

Introduce a subsampling scheme: Let Δ = kh, k → ∞, and define increments:

Ξ^(ε,h)_j(θ) = Z^(ε,h)_{(j+1)Δ} - Z^(ε,h)_{jΔ} - Δb(X^ε_{jΔ}, θ)

Contrast function:

C^(ε,h)(θ) = ∑_{j=0}^{N-1} Ξ^(ε,h)_j(θ)* H(X^ε_{jΔ}) Ξ^(ε,h)_j(θ)

Estimator: θ̂^ε = argmin_{θ∈Θ} C^(ε,h)(θ)

Technical Innovations

1. Error Control in Kernel Inversion: Through convolution representation of double-indexed kernels, precisely controlled the error of g^(h)(t,u) - 1, proving the key Lemma 3.5
2. Exploitation of Smoothing Effects: Although X^ε has Hölder regularity α - 1/2 < 1/2, the reconstruction error of Z^ε obtained through kernel inversion still achieves h^(1/2)
3. Parameter Selection Conditions: Established the condition h = o(ε^(1/α²)), balancing reconstruction error and statistical error

Experimental Setup

Mathematical Model

Consider a linear model:

X^ε_t = X^ε_0 + ε∫₀ᵗ K(t-s)dBs + ∫₀ᵗ K(t-s)(θ₀X^ε_s + θ₁)ds

True parameters: (θ₀*, θ₁*) = (-1, 1)

Experimental Parameters

• α values: 0.6, 0.8, 0.95
• Time intervals: T = 1, 10, 50
• Noise levels: ε = 1/10, 1/20, 1/100
• Sampling step sizes: h = 10⁻², 2×10⁻²
• Tuning parameters: k = 1, 2, 5, 10, 20
• Monte Carlo repetitions: 1000 times

Evaluation Metrics

• Empirical mean
• Empirical standard deviation (rescaled by ε⁻¹)

Experimental Results

Main Results

Theorem 3.1 (Reconstruction Error Bound):

‖Z^(h)_t - Z_t‖_{L^p} ≤ c(p)h^(1/2) sup_{s∈[0,T]} ‖a(X_s)‖_{L^p} + c(p)h^α sup_{s∈[0,T]} ‖b(X_s)‖_{L^p}

Theorem 4.1 (Consistency): Under assumptions A1-A4, θ̂^ε →^P θ⋆ as ε → 0.

Theorem 4.2 (Asymptotic Normality): Under assumptions A1-A6,

ε⁻¹(θ̂^ε - θ⋆) →^P I⁻¹J

where I is the Fisher information matrix and J is the random term.

Numerical Experiment Findings

1. Parameter k Selection is Critical:
   • When α = 0.8, optimal k decreases as ε decreases
   • When α = 0.95, k selection is relatively stable
   • When α = 0.6, the estimator is more sensitive to k selection
2. Impact of Time Interval T:
   • Larger T yields higher estimation accuracy
   • Longer observation periods improve estimation efficiency
3. Impact of Roughness α:
   • When α is smaller (rougher), estimator stability decreases
   • When α approaches 1, the method performs more stably

Comparison with TFE Method

• QMLE method has smaller asymptotic variance, especially for large T
• TFE method requires no parameter tuning but has lower efficiency

Related Work

1. Volterra Equation Theory: Pioneering work by Berger & Mizel (1980), subsequently extended to singular kernels and non-Lipschitz coefficients
2. Statistics of Rough Processes:
   • Rough volatility model by Gatheral et al. (2018)
   • Statistical inference theory by Chong et al. (2024)
3. Small Noise Asymptotics:
   • Classical theory by Kutoyants (1994)
   • Discrete observation methods by Sørensen & Uchida (2003)
4. Volterra Process Estimation:
   • Continuous observation MLE by Zuniga (2021)
   • Trajectory fitting method by Gloter & Yoshida (2025)

Conclusions and Discussion

Main Conclusions

1. Successfully constructed a QMLE estimator based on discrete observations, achieving optimal convergence rate ε^(-1)
2. Reconstruction error from kernel inversion achieves h^(1/2), which may be optimal
3. Established complete asymptotic theory, including consistency and asymptotic normality

Limitations

1. Strict Sampling Conditions: Requires h = o(ε^(1/α²)), which becomes very restrictive when α approaches 1/2
2. Parameter Tuning: Requires selecting appropriate k value, which may require cross-validation in practice
3. Model Assumptions: Requires coefficients to satisfy Lipschitz conditions and polynomial growth conditions

Future Directions

1. Study more general coefficient conditions
2. Develop adaptive selection methods for parameter k
3. Extend to multidimensional parameters and nonparametric settings
4. Investigate finite sample properties

In-Depth Evaluation

Strengths

1. Rigorous Theory: Provides complete mathematical proofs with refined technical treatment
2. Methodological Innovation: Cleverly exploits smoothing effects of kernel inversion, overcoming roughness difficulties
3. Practical Value: Method is applicable to real discrete observation data
4. Numerical Verification: Sufficient Monte Carlo experiments support theoretical results

Weaknesses

1. Restrictive Conditions: The sampling condition h = o(ε^(1/α²)) may be overly restrictive for practical applications
2. Parameter Selection: k value selection lacks theoretical guidance, relying on empirical tuning
3. Model Limitations: Only considers specific forms of Volterra kernels

Impact

1. Theoretical Contribution: Provides new theoretical framework for statistical inference of rough processes
2. Methodological Value: QMLE method can be generalized to similar problems
3. Application Prospects: Important application value in financial volatility modeling and related fields

Applicable Scenarios

1. Rough volatility modeling in financial markets
2. Parameter estimation for physical systems with memory effects
3. Fractional-order dynamics modeling in life sciences
4. Scientific computing applications requiring high-precision parameter estimation

References

The paper cites 36 important references covering classical and frontier work in Volterra equation theory, rough process statistics, small noise asymptotics, and other related fields, providing a solid theoretical foundation for the research.