2025-11-21T10:28:14.581273

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

Gloter, Yoshida

We consider a process X^$Îµ$ solution of a stochastic Volterra equation with an unknown parameter $Î¸$ in the drift function. The Volterra kernel is singular and given by K(u) = cu $Î±$-1 __u>0 with $Î±$ $\in$ (1/2, 1) and it is assumed that the diffusion coefficient is proportional to $Îµ$ $\rightarrow$ 0 Based on the observation of a discrete sampling with mesh h $\rightarrow$ 0 of the Volterra process, we build a Quasi Maximum Likelihood Estimator. The main step is to assess the error arising in the reconstruction of the path of a semi-martingale from the inversion of the Volterra kernel. We show that this error decreases as h^{1/2} whatever is the value of $Î±$. Then, we can introduce an explicit contrast function, which yields an efficient estimator when $Îµ$ $\rightarrow$ 0.

academic

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

基本信息

论文ID: 2510.09028
标题: Drift estimation for rough processes under small noise asymptotic : QMLE approach
作者: Arnaud Gloter (Université d'Evry), Nakahiro Yoshida (University of Tokyo)
分类: math.ST (Statistics Theory)
发表时间: October 13, 2025
论文链接: https://arxiv.org/abs/2510.09028

摘要

本文研究带有未知参数θ⋆的随机Volterra方程解Xε的漂移估计问题。Volterra核是奇异的，形式为K(u) = cuα-1𝟙u>0，其中α ∈ (1/2, 1)，扩散系数与ε → 0成比例。基于Volterra过程的离散采样观测（网格h → 0），构建了拟最大似然估计器(QMLE)。主要步骤是评估从Volterra核反演重构半鞅路径时产生的误差，证明该误差以h1/2的速率递减，与α值无关。进而引入显式对比函数，在ε → 0时得到有效估计器。

研究背景与动机

研究问题: 本文解决的核心问题是在小噪声渐近框架下，基于离散观测数据对粗糙随机过程的漂移参数进行统计推断。
问题重要性:
- 分数阶微分方程在物理学、数学金融、生命科学建模等领域广泛应用，具有记忆性质
- 粗糙过程(α-1/2 < 1/2)在金融资产波动率建模中被广泛研究
- 实际应用中只能获得离散观测数据，连续观测在实践中不可行
现有方法局限性:
- 连续观测下的MLE方法在实际中无法实现
- 现有轨迹拟合估计器(TFE)不是最优的
- 粗糙过程的Hölder正则性(α-1/2 < 1/2)使得直接离散化误差较大
研究动机: 开发基于离散观测的有效估计方法，克服粗糙性带来的技术困难，在小噪声渐近下实现参数的有效估计。

核心贡献

提出了新的QMLE方法: 基于Volterra核反演的拟最大似然估计器，适用于离散观测的粗糙过程
建立了关键的逼近误差界: 证明了重构误差以h1/2速率递减，这一结果对所有α ∈ (1/2, 1)都成立
提供了完整的渐近理论: 包括估计器的一致性和渐近正态性
给出了最优收敛速率: 估计器以ε⁻¹速率收敛，渐近方差达到Fisher信息的逆
进行了数值验证: 通过Monte Carlo实验验证了理论结果的有效性

方法详解

任务定义

考虑d维随机Volterra方程：

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

其中：

K(u) = u^(α-1)/Γ(α) 𝟙_{u>0}，α ∈ (1/2, 1)为奇异核
θ⋆ ∈ Θ ⊂ ℝ^(d_Θ)为未知漂移参数
ε → 0为小噪声参数
观测数据为{X^ε_}_，h → 0

目标：基于离散观测估计参数θ⋆

模型架构

1. 核心思想：半鞅重构

通过Volterra核的反演重构半鞅过程：

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

使用resolvent核L(u) = u^(-α)/Γ(1-α)𝟙_{u>0}，满足L⋆K = 1。

2. 离散重构过程

定义离散观测过程X^(h)t = X{φ_h(t)}，其中φ_h(t) = h⌊t/h⌋，构造：

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

3. 估计器构造

引入子采样方案：设Δ = kh，k → ∞，定义增量：

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

对比函数：

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

估计器：θ̂^ε = argmin_{θ∈Θ} C^(ε,h)(θ)

技术创新点

核反演的误差控制: 通过双指标核的卷积表示，精确控制了g^(h)(t,u) - 1的误差，证明了关键的Lemma 3.5
平滑效应的利用: 虽然X^ε的Hölder正则性为α-1/2 < 1/2，但通过核反演得到的Z^ε重构误差仍能达到h^(1/2)
参数选择条件: 建立了h = o(ε^(1/α²))的条件，平衡了重构误差和统计误差

实验设置

数学模型

考虑线性模型：

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

真实参数：(θ₀*, θ₁*) = (-1, 1)

实验参数

α值：0.6, 0.8, 0.95
时间区间：T = 1, 10, 50
噪声水平：ε = 1/10, 1/20, 1/100
采样步长：h = 10⁻², 2×10⁻²
调节参数：k = 1, 2, 5, 10, 20
Monte Carlo重复次数：1000次

评价指标

经验均值
经验标准差(按ε⁻¹重标准化)

实验结果

主要结果

定理3.1 (重构误差界)：

‖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}

定理4.1 (一致性)：在假设A1-A4下，θ̂^ε →^P θ⋆当ε → 0。

定理4.2 (渐近正态性)：在假设A1-A6下，

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

其中I是Fisher信息矩阵，J是随机项。

数值实验发现

参数k的选择至关重要：
- α = 0.8时，最优k随ε减小而减小
- α = 0.95时，k的选择相对稳定
- α = 0.6时，估计器对k的选择更敏感
时间区间T的影响：
- T越大，估计精度越高
- 长期观测有助于提高估计效率
粗糙度α的影响：
- α较小(更粗糙)时，估计器稳定性下降
- α接近1时，方法表现更稳定