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

小ノイズ漸近下の粗過程のドリフト推定：QMLE手法

基本情報

論文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 (統計理論)
発表日: 2025年10月13日
論文リンク: 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に近い場合、手法の性能がより安定している