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

المعلومات الأساسية

  • معرّف الورقة: 2510.09028
  • العنوان: Drift estimation for rough processes under small noise asymptotic : QMLE approach
  • المؤلفون: Arnaud Gloter (جامعة Evry)، Nakahiro Yoshida (جامعة طوكيو)
  • التصنيف: math.ST (نظرية الإحصاء)
  • تاريخ النشر: 13 أكتوبر 2025
  • رابط الورقة: https://arxiv.org/abs/2510.09028

الملخص

تدرس هذه الورقة مشكلة تقدير الانجراف للحل Xε من معادلة Volterra العشوائية مع معامل غير معروف θ⋆. نواة Volterra متفردة بالشكل K(u) = cuα-1𝟙u>0، حيث α ∈ (1/2, 1)، ومعامل الانتشار يتناسب مع ε → 0. بناءً على الملاحظات المأخوذة بعينات منفصلة من عملية Volterra (الشبكة h → 0)، تم بناء مقدّر الاحتمالية العظمى شبه الكلاسيكي (QMLE). الخطوة الرئيسية هي تقييم الخطأ الناشئ من إعادة بناء المسار شبه الرياضي من خلال عكس نواة Volterra، مما يثبت أن هذا الخطأ يتناقص بمعدل h1/2، بغض النظر عن قيمة α. ثم يتم إدخال دالة تباين صريحة للحصول على مقدّر فعال عندما ε → 0.

الخلفية البحثية والدافع

  1. المشكلة البحثية: تحل هذه الورقة المشكلة الأساسية المتمثلة في الاستدلال الإحصائي على معاملات الانجراف للعمليات العشوائية الخشنة بناءً على بيانات الملاحظات المنفصلة في إطار التقارب المقارب للضوضاء الصغيرة.
  2. أهمية المشكلة:
    • تُستخدم المعادلات التفاضلية الكسرية على نطاق واسع في الفيزياء والرياضيات المالية والنمذجة العلمية الحيوية، وتتمتع بخصائص الذاكرة
    • تُدرس العمليات الخشنة (α-1/2 < 1/2) على نطاق واسع في نمذجة تقلبات الأصول المالية
    • في التطبيقات العملية، يمكن الحصول فقط على بيانات الملاحظات المنفصلة، والملاحظة المستمرة غير قابلة للتطبيق عملياً
  3. قيود الطرق الموجودة:
    • طرق MLE تحت الملاحظة المستمرة غير قابلة للتطبيق عملياً
    • مقدّرات ملاءمة المسار الموجودة (TFE) ليست مثالية
    • الانتظام Hölder للعمليات الخشنة (α-1/2 < 1/2) يؤدي إلى خطأ تقدير مباشر أكبر
  4. الدافع البحثي: تطوير طرق تقدير فعالة بناءً على الملاحظات المنفصلة، والتغلب على الصعوبات التقنية الناجمة عن الخشونة، وتحقيق التقدير الفعال للمعاملات تحت التقارب المقارب للضوضاء الصغيرة.

المساهمات الأساسية

  1. اقتراح طريقة QMLE جديدة: مقدّر الاحتمالية العظمى شبه الكلاسيكي بناءً على عكس نواة Volterra، مناسب للعمليات الخشنة مع الملاحظات المنفصلة
  2. إنشاء حدود خطأ التقريب الرئيسية: إثبات أن خطأ إعادة البناء يتناقص بمعدل h1/2، وهذه النتيجة صحيحة لجميع α ∈ (1/2, 1)
  3. توفير نظرية تقارب مقارب كاملة: تشمل الاتساق والحالة الطبيعية المقاربة للمقدّر
  4. إعطاء معدل تقارب أمثل: يتقارب المقدّر بمعدل ε⁻¹، والتباين المقارب يصل إلى معكوس معلومات Fisher
  5. إجراء التحقق العددي: التحقق من صحة النتائج النظرية من خلال تجارب Monte Carlo

شرح الطريقة

تعريف المهمة

ضع في الاعتبار معادلة Volterra العشوائية ذات البعد d:

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)(θ)

نقاط الابتكار التقني

  1. التحكم في خطأ عكس النواة: من خلال تمثيل الالتفاف النواة ثنائي الفهرس، يتم التحكم الدقيق في خطأ g^(h)(t,u) - 1، مما يثبت Lemma 3.5 الرئيسي
  2. استخدام تأثير التمويه: على الرغم من أن الانتظام Hölder لـ X^ε هو α-1/2 < 1/2، فإن خطأ إعادة البناء Z^ε الذي تم الحصول عليه من خلال عكس النواة يمكن أن يصل إلى h^(1/2)
  3. شروط اختيار المعاملات: إنشاء شرط 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 هي الحد العشوائي.

اكتشافات التجارب العددية

  1. اختيار المعامل k حاسم:
    • عندما α = 0.8، يقل k الأمثل مع تناقص ε
    • عندما α = 0.95، يكون اختيار k مستقراً نسبياً
    • عندما α = 0.6، يكون المقدّر أكثر حساسية لاختيار k
  2. تأثير الفترة الزمنية T:
    • كلما زاد T، زادت دقة التقدير
    • تساعد الملاحظات طويلة الأجل على تحسين كفاءة التقدير
  3. تأثير درجة الخشونة α:
    • عندما تكون α أصغر (أكثر خشونة)، تنخفض استقرارية المقدّر
    • عندما تقترب α من 1، يكون أداء الطريقة أكثر استقراراً

المقارنة مع طريقة TFE

  • طريقة QMLE لها تباين مقارب أصغر، خاصة في حالة T الكبيرة
  • طريقة TFE لا تتطلب معاملات تعديل، لكن الكفاءة أقل

الأعمال ذات الصلة

  1. نظرية معادلات Volterra: العمل الرائد لـ Berger & Mizel (1980)، مع امتدادات لاحقة إلى النوى المتفردة والمعاملات غير Lipschitz
  2. إحصائيات العمليات الخشنة:
    • نموذج التقلبات الخشنة لـ Gatheral et al. (2018)
    • نظرية الاستدلال الإحصائي لـ Chong et al. (2024)
  3. التقارب المقارب للضوضاء الصغيرة:
    • النظرية الكلاسيكية لـ Kutoyants (1994)
    • طريقة الملاحظات المنفصلة لـ Sørensen & Uchida (2003)
  4. تقدير عمليات Volterra:
    • MLE للملاحظة المستمرة لـ Zuniga (2021)
    • طريقة ملاءمة المسار لـ Gloter & Yoshida (2025)

الخلاصة والمناقشة

الاستنتاجات الرئيسية

  1. تم بناء مقدّر QMLE بنجاح بناءً على الملاحظات المنفصلة، مما حقق معدل التقارب الأمثل ε⁻¹
  2. خطأ إعادة بناء عكس النواة يصل إلى h^(1/2)، وقد تكون هذه النتيجة مثالية
  3. تم إنشاء نظرية تقارب مقارب كاملة، تشمل الاتساق والحالة الطبيعية المقاربة

القيود

  1. شروط أخذ العينات صارمة: تتطلب h = o(ε^(1/α²))، وعندما تقترب α من 1/2 تكون الشروط صارمة جداً
  2. تعديل المعاملات: يتطلب اختيار قيمة k مناسبة، وقد يتطلب التحقق المتقاطع في الممارسة العملية
  3. افتراضات النموذج: يتطلب أن تحقق المعاملات شروط Lipschitz وشروط النمو متعدد الحدود

الاتجاهات المستقبلية

  1. دراسة شروط المعاملات الأكثر عمومية
  2. تطوير طرق اختيار تكيفية لقيمة المعامل k
  3. التوسع إلى الإعدادات متعددة الأبعاد وغير البارامترية
  4. دراسة الخصائص ذات العينات المحدودة

التقييم المتعمق

المزايا

  1. صرامة نظرية: توفير إثبات رياضي كامل، معالجة تقنية دقيقة
  2. ابتكار الطريقة: استخدام ذكي لتأثير التمويه من عكس النواة، التغلب على صعوبات الخشونة
  3. القيمة العملية: الطريقة مناسبة لبيانات الملاحظات المنفصلة الفعلية
  4. التحقق العددي: تجارب Monte Carlo كافية تدعم النتائج النظرية

أوجه القصور

  1. تقييد الشروط: شرط أخذ العينات h = o(ε^(1/α²)) قد يكون صارماً جداً في التطبيقات العملية
  2. اختيار المعاملات: اختيار قيمة k يفتقر إلى التوجيه النظري، ويعتمد على التعديل التجريبي
  3. قيود النموذج: يأخذ في الاعتبار فقط شكلاً محدداً من نوى Volterra

التأثير

  1. المساهمة النظرية: توفير إطار نظري جديد للاستدلال الإحصائي للعمليات الخشنة
  2. قيمة الطريقة: يمكن تعميم طريقة QMLE على مشاكل مماثلة أخرى
  3. آفاق التطبيق: قيمة تطبيقية مهمة في مجالات مثل نمذجة تقلبات الأسواق المالية

السيناريوهات المعمول بها

  1. نمذجة التقلبات الخشنة في الأسواق المالية
  2. تقدير معاملات الأنظمة الفيزيائية ذات تأثيرات الذاكرة
  3. نمذجة الديناميكيات الكسرية في العلوم الحيوية
  4. تطبيقات الحسابات العلمية التي تتطلب تقديراً دقيقاً للمعاملات

المراجع

تستشهد الورقة بـ 36 مرجعاً مهماً، تغطي نظرية معادلات Volterra وإحصائيات العمليات الخشنة والتقارب المقارب للضوضاء الصغيرة وغيرها من المجالات الكلاسيكية والمتقدمة، مما يوفر أساساً نظرياً متيناً للبحث.