This paper deals with a Skorokhod's integral based projection type estimator $\widehat b_m$ of the drift function $b_0$ computed from $N\in\mathbb N^*$ independent copies $X^1,\dots,X^N$ of the solution $X$ of $dX_t = b_0(X_t)dt +ÏdB_t$, where $B$ is a fractional Brownian motion of Hurst index $H\in (1/2,1)$. Skorokhod's integral based estimators cannot be calculated directly from $X^1,\dots,X^N$, but in this paper an $\mathbb L^2$-error bound is established on a calculable approximation of $\widehat b_m$.
論文ID : 2307.04949タイトル : On a Calculable Skorokhod's Integral Based Projection Estimator of the Drift Function in Fractional SDE著者 : Nicolas Marie (Laboratoire Modal'X, Université Paris Nanterre)分類 : math.ST（統計理論）、math.PR（確率論）、stat.TH（統計理論）発表日時 : 2024年1月15日（arXiv v2）論文リンク : https://arxiv.org/abs/2307.04949 本論文は、分数ブラウン運動に駆動される確率微分方程式における漂流関数 b 0 b_0 b 0 b ^ m \hat{b}_m b ^ m X X X N N N X 1 , … , X N X^1,\ldots,X^N X 1 , … , X N d X t = b 0 ( X t ) d t + σ d B t dX_t = b_0(X_t)dt + \sigma dB_t d X t = b 0 ( X t ) d t + σ d B t B B B H ∈ ( 1 / 2 , 1 ) H \in (1/2,1) H ∈ ( 1/2 , 1 ) L 2 L^2 L 2
本論文が解決する中核的問題は、分数確率微分方程式における漂流関数の非パラメトリック推定問題である。以下の分数SDE を考察する：
X t = X 0 + ∫ 0 t b 0 ( X s ) d s + σ B t , t ∈ [ 0 , T ] X_t = X_0 + \int_0^t b_0(X_s)ds + \sigma B_t, \quad t \in [0,T] X t = X 0 + ∫ 0 t b 0 ( X s ) d s + σ B t , t ∈ [ 0 , T ]
ここで B B B H ∈ ( 1 / 2 , 1 ) H \in (1/2,1) H ∈ ( 1/2 , 1 ) b 0 ∈ C 1 ( R ) b_0 \in C^1(\mathbb{R}) b 0 ∈ C 1 ( R )
実用的応用の必要性 ：薬物動力学などの分野では、単一標本の長時間挙動に依存するのではなく、複数の独立標本（臨床試験における複数患者など）に基づいてパラメータを推定する必要がある。計算上の課題 ：従来のスコロホッド積分ベースの推定器は理論的には実行可能だが、観測データから直接計算できず、その実用性を著しく制限している。理論的ギャップ ：既存研究は主にパラメトリック推定または長時間挙動ベースの推定器に焦点を当てており、独立複製ベースの非パラメトリック推定に関する研究は相対的に少ない。計算可能な推定器の提案 ：不動点戦略を通じて、計算不可能なスコロホッド積分推定器を計算可能な不動点推定器 b ~ m \tilde{b}_m b ~ m 理論的誤差界の確立 ：補助推定器 b ^ m \hat{b}_m b ^ m b ~ m \tilde{b}_m b ~ m L 2 L^2 L 2 収束率分析 ：定常状況下で、収束率が N − ( 2 H − 1 ) / ( 10 H ) N^{-(2H-1)/(10H)} N − ( 2 H − 1 ) / ( 10 H ) N N N T T T m m m 非パラメトリック枠組みの拡張 ：Marie 12 のパラメトリック推定不動点戦略を非パラメトリック推定枠組みに成功裏に拡張する。分数SDE の N N N X 1 , … , X N X^1,\ldots,X^N X 1 , … , X N [ 0 , T ] [0,T] [ 0 , T ] b 0 b_0 b 0
H = { h : ∥ h ∥ H < ∞ } H = \{h: \|h\|_H < \infty\} H = { h : ∥ h ∥ H < ∞ } B T B_T B T ⟨ h , η ⟩ H = α H ∫ 0 T ∫ 0 T h ( s ) η ( t ) ∣ t − s ∣ 2 H − 2 d s d t \langle h,\eta \rangle_H = \alpha_H \int_0^T \int_0^T h(s)\eta(t)|t-s|^{2H-2}dsdt ⟨ h , η ⟩ H = α H ∫ 0 T ∫ 0 T h ( s ) η ( t ) ∣ t − s ∣ 2 H − 2 d s d t α H = H ( 2 H − 1 ) \alpha_H = H(2H-1) α H = H ( 2 H − 1 )
スコロホッド積分はマリアヴァン導数の随伴作用素として定義される。過程 Y Y Y ∫ 0 t Y s δ X s = ∫ 0 t Y s b 0 ( X s ) d s + σ ∫ 0 t Y s δ B s \int_0^t Y_s \delta X_s = \int_0^t Y_s b_0(X_s)ds + \sigma \int_0^t Y_s \delta B_s ∫ 0 t Y s δ X s = ∫ 0 t Y s b 0 ( X s ) d s + σ ∫ 0 t Y s δ B s
区間 I I I ( ϕ 1 , … , ϕ m ) (\phi_1,\ldots,\phi_m) ( ϕ 1 , … , ϕ m ) b ^ m ( x ) = 1 f ( x ) ∑ j = 1 m ( 1 N T ∑ i = 1 N ∫ 0 T ϕ j ( X s i ) δ X s i ) ϕ j ( x ) \hat{b}_m(x) = \frac{1}{f(x)} \sum_{j=1}^m \left(\frac{1}{NT} \sum_{i=1}^N \int_0^T \phi_j(X^i_s)\delta X^i_s\right) \phi_j(x) b ^ m ( x ) = f ( x ) 1 ∑ j = 1 m ( NT 1 ∑ i = 1 N ∫ 0 T ϕ j ( X s i ) δ X s i ) ϕ j ( x )
ここで f ( x ) = 1 T ∫ 0 T f s ( x ) d s f(x) = \frac{1}{T}\int_0^T f_s(x)ds f ( x ) = T 1 ∫ 0 T f s ( x ) d s
スコロホッド積分は直接計算できないため、Nualart 15 の結果を利用する：
∫ 0 T ϕ j ( X s i ) δ X s i = I i , j − a ∫ 0 T ∫ 0 t ϕ j ′ ( X t i ) exp ( ∫ s t b 0 ′ ( X u i ) d u ) ∣ t − s ∣ 2 H − 2 d s d t \int_0^T \phi_j(X^i_s)\delta X^i_s = I_{i,j} - a\int_0^T \int_0^t \phi'_j(X^i_t)\exp\left(\int_s^t b'_0(X^i_u)du\right)|t-s|^{2H-2}dsdt ∫ 0 T ϕ j ( X s i ) δ X s i = I i , j − a ∫ 0 T ∫ 0 t ϕ j ′ ( X t i ) exp ( ∫ s t b 0 ′ ( X u i ) d u ) ∣ t − s ∣ 2 H − 2 d s d t
ここで I i , j = ∫ 0 T ϕ j ( X s i ) d X s i I_{i,j} = \int_0^T \phi_j(X^i_s)dX^i_s I i , j = ∫ 0 T ϕ j ( X s i ) d X s i a = σ 2 H ( 2 H − 1 ) a = \sigma^2 H(2H-1) a = σ 2 H ( 2 H − 1 )
これは関係式 b ^ m = Φ m ( b 0 ) ≈ Φ m ( b ^ m ) \hat{b}_m = \Phi_m(b_0) \approx \Phi_m(\hat{b}_m) b ^ m = Φ m ( b 0 ) ≈ Φ m ( b ^ m ) Φ m \Phi_m Φ m Φ m ( ϕ ) ( x ) = 1 f ( x ) ∑ j = 1 m [ 1 N T ∑ i = 1 N [ I i , j − a ∫ 0 T ∫ 0 t ϕ j ′ ( X t i ) exp ( ∫ s t ϕ ′ ( X u i ) d u ) ∣ t − s ∣ 2 H − 2 d s d t ] ] ϕ j ( x ) \Phi_m(\phi)(x) = \frac{1}{f(x)} \sum_{j=1}^m \left[\frac{1}{NT}\sum_{i=1}^N \left[I_{i,j} - a\int_0^T \int_0^t \phi'_j(X^i_t)\exp\left(\int_s^t \phi'(X^i_u)du\right)|t-s|^{2H-2}dsdt\right]\right]\phi_j(x) Φ m ( ϕ ) ( x ) = f ( x ) 1 ∑ j = 1 m [ NT 1 ∑ i = 1 N [ I i , j − a ∫ 0 T ∫ 0 t ϕ j ′ ( X t i ) exp ( ∫ s t ϕ ′ ( X u i ) d u ) ∣ t − s ∣ 2 H − 2 d s d t ] ] ϕ j ( x )
不動点推定器 b ~ m \tilde{b}_m b ~ m S m , c = { ϕ ∈ S m : ∥ ϕ ′ ∥ ∞ ≤ c } S_{m,c} = \{\phi \in S_m: \|\phi'\|_\infty \leq c\} S m , c = { ϕ ∈ S m : ∥ ϕ ′ ∥ ∞ ≤ c } Φ m \Phi_m Φ m c > ∥ b 0 ′ ∥ ∞ c > \|b'_0\|_\infty c > ∥ b 0 ′ ∥ ∞
論文は主に理論分析を行い、2つの主要な状況を考察する：
決定論的初期値の場合 ：X 0 ( ⋅ ) = x 0 ∈ R X_0(\cdot) = x_0 \in \mathbb{R} X 0 ( ⋅ ) = x 0 ∈ R 定常状況 ：b 0 b_0 b 0 ∃ m > 0 : ∀ x ∈ R , b 0 ′ ( x ) ≤ − m \exists m > 0: \forall x \in \mathbb{R}, b'_0(x) \leq -m ∃ m > 0 : ∀ x ∈ R , b 0 ′ ( x ) ≤ − m 主に三角基底関数の場合を考察する：コンパクト区間 I = [ ℓ , r ] I = [\ell,r] I = [ ℓ , r ] ϕ 1 ( x ) = 1 r − ℓ 1 I ( x ) \phi_1(x) = \sqrt{\frac{1}{r-\ell}}1_I(x) ϕ 1 ( x ) = r − ℓ 1 1 I ( x ) ϕ 2 j + 1 ( x ) = 2 r − ℓ sin ( 2 π j ( x − ℓ ) r − ℓ ) 1 I ( x ) \phi_{2j+1}(x) = \sqrt{\frac{2}{r-\ell}}\sin\left(\frac{2\pi j(x-\ell)}{r-\ell}\right)1_I(x) ϕ 2 j + 1 ( x ) = r − ℓ 2 sin ( r − ℓ 2 πj ( x − ℓ ) ) 1 I ( x ) ϕ 2 j ( x ) = 2 r − ℓ cos ( 2 π j ( x − ℓ ) r − ℓ ) 1 I ( x ) \phi_{2j}(x) = \sqrt{\frac{2}{r-\ell}}\cos\left(\frac{2\pi j(x-\ell)}{r-\ell}\right)1_I(x) ϕ 2 j ( x ) = r − ℓ 2 cos ( r − ℓ 2 πj ( x − ℓ ) ) 1 I ( x )
仮定3.1 ：X t X_t X t f t f_t f t t ↦ f t ( x ) ∈ L 1 ( [ 0 , T ] ) t \mapsto f_t(x) \in L^1([0,T]) t ↦ f t ( x ) ∈ L 1 ([ 0 , T ]) 仮定4.1 ：m f > 0 m_f > 0 m f > 0 f ( x ) ≥ m f , ∀ x ∈ I f(x) \geq m_f, \forall x \in I f ( x ) ≥ m f , ∀ x ∈ I 仮定4.4 ：定数 c ϕ > 0 c_\phi > 0 c ϕ > 0 I ( m ) ∨ L ( m ) ≤ c ϕ R ( m ) I(m) \vee L(m) \leq c_\phi R(m) I ( m ) ∨ L ( m ) ≤ c ϕ R ( m ) ここで L ( m ) = ∑ j = 1 m ∥ ϕ j ∥ ∞ 2 L(m) = \sum_{j=1}^m \|\phi_j\|_\infty^2 L ( m ) = ∑ j = 1 m ∥ ϕ j ∥ ∞ 2 R ( m ) = ∑ j = 1 m ∥ ϕ j ′ ∥ ∞ 2 R(m) = \sum_{j=1}^m \|\phi'_j\|_\infty^2 R ( m ) = ∑ j = 1 m ∥ ϕ j ′ ∥ ∞ 2
定理3.3（補助推定器誤差界） ：
E ( ∥ b ^ m − b 0 ∥ f 2 2 ) ≤ ∥ b m − b 0 ∥ f 2 2 + 2 N ( ∥ b 0 ∥ f 2 L ( m ) + c 2.5 σ 2 m T T 2 − 2 H ( L ( m ) + R ( m ) ) ) E(\|\hat{b}_m - b_0\|_{f^2}^2) \leq \|b_m - b_0\|_{f^2}^2 + \frac{2}{N}\left(\|b_0\|_f^2 L(m) + c_{2.5}\sigma^2 \frac{m_T}{T^{2-2H}}(L(m) + R(m))\right) E ( ∥ b ^ m − b 0 ∥ f 2 2 ) ≤ ∥ b m − b 0 ∥ f 2 2 + N 2 ( ∥ b 0 ∥ f 2 L ( m ) + c 2.5 σ 2 T 2 − 2 H m T ( L ( m ) + R ( m )) )
定理4.7（不動点推定器誤差界） ：定常状況下で、
E ( ∥ b ~ m c , l − b 0 ∥ f 2 2 ) ≤ 2 ∥ b m − b 0 ∥ f 2 2 + c 4.7 R ( m ) V ( N , T ) E(\|\tilde{b}_m^{c,l} - b_0\|_{f^2}^2) \leq 2\|b_m - b_0\|_{f^2}^2 + c_{4.7}R(m)V(N,T) E ( ∥ b ~ m c , l − b 0 ∥ f 2 2 ) ≤ 2∥ b m − b 0 ∥ f 2 2 + c 4.7 R ( m ) V ( N , T ) V ( N , T ) = N − 1 / 2 T − 1 + T 2 H − 1 V(N,T) = N^{-1/2}T^{-1} + T^{2H-1} V ( N , T ) = N − 1/2 T − 1 + T 2 H − 1
三角基底関数に対して、( b 0 f ) ∣ I ∈ W 2 β ( I ) (b_0f)|_I \in W_2^\beta(I) ( b 0 f ) ∣ I ∈ W 2 β ( I ) β \beta β E ( ∥ b ~ m c , l − b 0 ∥ f 2 2 ) ≤ c 4.7 ( m − 2 β + m 3 N − ( 2 H − 1 ) / ( 4 H ) ) E(\|\tilde{b}_m^{c,l} - b_0\|_{f^2}^2) \leq c_{4.7}(m^{-2\beta} + m^3 N^{-(2H-1)/(4H)}) E ( ∥ b ~ m c , l − b 0 ∥ f 2 2 ) ≤ c 4.7 ( m − 2 β + m 3 N − ( 2 H − 1 ) / ( 4 H ) )
最適選択：
時間窓：T ( N ) = N − 1 / ( 4 H ) T(N) = N^{-1/(4H)} T ( N ) = N − 1/ ( 4 H ) 射影次元：m ≍ N ( 2 H − 1 ) / ( 4 H ( 3 + 2 β ) ) m \asymp N^{(2H-1)/(4H(3+2\beta))} m ≍ N ( 2 H − 1 ) / ( 4 H ( 3 + 2 β )) 収束率：N − ( 2 H − 1 ) / ( 10 H ) N^{-(2H-1)/(10H)} N − ( 2 H − 1 ) / ( 10 H ) 補題4.6 ：適切な条件下で、
P ( Δ m c ) ≤ c 4.6 , 2 R ( m ) ( N − 1 / 2 T − 1 + T 2 H − 1 ) P(\Delta_m^c) \leq c_{4.6,2}R(m)(N^{-1/2}T^{-1} + T^{2H-1}) P ( Δ m c ) ≤ c 4.6 , 2 R ( m ) ( N − 1/2 T − 1 + T 2 H − 1 ) Δ m \Delta_m Δ m Φ m \Phi_m Φ m
長時間挙動推定器 ：Kleptsyna & Le Breton 9 、Tudor & Viens 17 、Hu & Nualart 7 など、単一標本の長時間観測に基づくパラメトリック推定を研究した。非パラメトリック推定 ：Saussereau 16 、Comte & Marie 2 は分数SDE の非パラメトリック推定を研究した。独立複製法 ：Comte & Genon-Catalot 1 は古典的拡散過程における独立複製ベースの推定を初めて研究し、Marie 12 はこれを分数情形のパラメトリック推定に拡張した。本論文は初めて独立複製ベースの方法を分数SDE の非パラメトリック推定に拡張し、スコロホッド積分の計算問題を解決した。
分数SDE 漂流関数の計算可能な非パラメトリック推定器の構成に成功した 存在性、一意性、収束率を含む完全な理論分析枠組みを確立した 標本数、観測時間、射影次元間の複雑なバランス関係を明らかにした 計算複雑性 ：不動点法は反復計算を必要とし、実際の実装はより複雑である密度推定への依存 ：実用的応用では未知の密度関数 f f f 仮定条件 ：より強い正則性仮定と有界導数条件を必要とする密度関数が未知の場合の結合推定問題の研究 より効率的な数値アルゴリズムの開発 多次元情形と一般的な分数過程への拡張 理論的革新性 ：分数SDE におけるスコロホッド積分推定器の計算問題を初めて解決した方法の厳密性 ：完全な数学的証明と誤差分析実用的価値 ：分数過程の統計的推論に新しい計算枠組みを提供する技術的深さ ：マリアヴァン分析と不動点理論を巧みに結合している実際的検証の欠如 ：論文は純理論的であり、数値実験検証を欠く計算効率 ：不動点反復の収束速度と計算複雑度が深く分析されていない仮定の制限 ：いくつかの技術的仮定は実際の応用では検証が困難な可能性がある学術的貢献 ：分数確率過程統計学に新しい方向を開く方法論的価値 ：不動点戦略は他の類似問題に適用可能である応用前景 ：金融、生物学などの分野で潜在的な応用価値を有する金融市場における分数ブラウン運動モデルのパラメータ推定 生物医学における複数個体データに基づくパラメータ推論 物理システムにおける長記憶過程のモデリング 論文は17篇の関連文献を引用しており、主に以下を含む：
マリアヴァン分析基礎理論（Nualart 15 、Decreusefond 4 ） 分数SDE 理論（Hairer & Ohashi 6 ） 統計的推論方法（Comte & Marie 2,3 、Marie 12 ） 近似理論（DeVore & Lorentz 5 ） 総合評価 ：これは分数確率微分方程式の非パラメトリック推定分野における重要な貢献をなす高質量の理論統計学論文である。数値実験を欠いているが、理論分析は深く厳密であり、当該分野のさらなる発展のための堅固な基礎を築いている。
