Estimation de la dérive pour les processus rugueux sous asymptotique de petit bruit : approche QMLE

Informations de base

• ID de l'article: 2510.09028
• Titre: Drift estimation for rough processes under small noise asymptotic : QMLE approach
• Auteurs: Arnaud Gloter (Université d'Evry), Nakahiro Yoshida (University of Tokyo)
• Classification: math.ST (Théorie statistique)
• Date de publication: 13 octobre 2025
• Lien de l'article: https://arxiv.org/abs/2510.09028

Résumé

Cet article étudie le problème d'estimation de la dérive pour la solution X^ε d'équations stochastiques de Volterra avec paramètre inconnu θ⋆. Le noyau de Volterra est singulier, de la forme K(u) = cu^(α-1)𝟙_{u>0}, où α ∈ (1/2, 1), et le coefficient de diffusion est proportionnel à ε → 0. Sur la base d'observations discrètes du processus de Volterra (grille h → 0), un estimateur de quasi-vraisemblance maximale (QMLE) est construit. L'étape principale consiste à évaluer l'erreur générée lors de la reconstruction de la trajectoire semi-martingale par inversion du noyau de Volterra, en prouvant que cette erreur décroît au taux h^(1/2), indépendamment de la valeur de α. Une fonction de contraste explicite est ensuite introduite pour obtenir un estimateur efficace lorsque ε → 0.

Contexte et motivation de la recherche

1. Problème de recherche: L'article résout le problème fondamental de l'inférence statistique des paramètres de dérive pour les processus stochastiques rugueux basée sur des données d'observations discrètes dans le cadre de l'asymptotique de petit bruit.
2. Importance du problème:
   • Les équations différentielles fractionnaires sont largement appliquées en physique, en mathématiques financières, en modélisation des sciences de la vie, etc., avec des propriétés de mémoire
   • Les processus rugueux (α-1/2 < 1/2) sont largement étudiés dans la modélisation de la volatilité des actifs financiers
   • En pratique, seules des données d'observations discrètes sont disponibles; l'observation continue n'est pas réalisable
3. Limitations des méthodes existantes:
   • Les méthodes MLE sous observation continue ne peuvent pas être implémentées en pratique
   • Les estimateurs d'ajustement de trajectoire existants (TFE) ne sont pas optimaux
   • La régularité Höldérienne des processus rugueux (α-1/2 < 1/2) entraîne une erreur de discrétisation directe importante
4. Motivation de la recherche: Développer des méthodes d'estimation efficaces basées sur des observations discrètes, surmonter les difficultés techniques causées par la rugosité, et réaliser une estimation efficace des paramètres sous asymptotique de petit bruit.

Contributions principales

1. Proposition d'une nouvelle méthode QMLE: Estimateur de quasi-vraisemblance maximale basé sur l'inversion du noyau de Volterra, applicable aux processus rugueux avec observations discrètes
2. Établissement de bornes d'erreur d'approximation clés: Preuve que l'erreur de reconstruction décroît au taux h^(1/2), résultat valide pour tout α ∈ (1/2, 1)
3. Fourniture d'une théorie asymptotique complète: Incluant la cohérence et la normalité asymptotique de l'estimateur
4. Obtention du taux de convergence optimal: L'estimateur converge au taux ε⁻¹, avec variance asymptotique égale à l'inverse de l'information de Fisher
5. Vérification numérique: Validation des résultats théoriques par expériences de Monte Carlo

Détails de la méthode

Définition de la tâche

Considérons l'équation stochastique de Volterra d-dimensionnelle:

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

Où:

• K(u) = u^(α-1)/Γ(α) 𝟙_{u>0}, α ∈ (1/2, 1) est le noyau singulier
• θ⋆ ∈ Θ ⊂ ℝ^(d_Θ) est le paramètre de dérive inconnu
• ε → 0 est le paramètre de petit bruit
• Les données observées sont {X^ε_}_, h → 0

Objectif: Estimer le paramètre θ⋆ basé sur les observations discrètes

Architecture du modèle

1. Idée centrale: reconstruction semi-martingale

Reconstruction du processus semi-martingale par inversion du noyau de Volterra:

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

Utilisation du noyau résolvant L(u) = u^(-α)/Γ(1-α)𝟙_{u>0}, satisfaisant L⋆K = 1.

2. Processus de reconstruction discrète

Définition du processus d'observation discrète X^(h)t = X{φ_h(t)}, où φ_h(t) = h⌊t/h⌋, construction:

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

3. Construction de l'estimateur

Introduction d'un schéma de sous-échantillonnage: soit Δ = kh, k → ∞, définition des accroissements:

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

Fonction de contraste:

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

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

Points d'innovation technique

1. Contrôle d'erreur de l'inversion du noyau: Par représentation de convolution du noyau à double indice, contrôle précis de l'erreur g^(h)(t,u) - 1, preuve du Lemme 3.5 clé
2. Utilisation de l'effet de lissage: Bien que la régularité Höldérienne de X^ε soit α-1/2 < 1/2, l'erreur de reconstruction de Z^ε obtenue par inversion du noyau peut atteindre h^(1/2)
3. Conditions de sélection des paramètres: Établissement de la condition h = o(ε^(1/α²)), équilibrant l'erreur de reconstruction et l'erreur statistique

Configuration expérimentale

Modèle mathématique

Considérons le modèle linéaire:

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

Paramètres vrais: (θ₀*, θ₁*) = (-1, 1)

Paramètres expérimentaux

• Valeurs de α: 0.6, 0.8, 0.95
• Intervalle de temps: T = 1, 10, 50
• Niveaux de bruit: ε = 1/10, 1/20, 1/100
• Pas d'échantillonnage: h = 10⁻², 2×10⁻²
• Paramètre d'ajustement: k = 1, 2, 5, 10, 20
• Nombre de répétitions de Monte Carlo: 1000

Indicateurs d'évaluation

• Moyenne empirique
• Écart-type empirique (renormalisé par ε⁻¹)

Résultats expérimentaux

Résultats principaux

Théorème 3.1 (Borne d'erreur de reconstruction):

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

Théorème 4.1 (Cohérence): Sous les hypothèses A1-A4, θ̂^ε →^P θ⋆ lorsque ε → 0.

Théorème 4.2 (Normalité asymptotique): Sous les hypothèses A1-A6,

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

où I est la matrice d'information de Fisher, J est le terme aléatoire.

Découvertes des expériences numériques

1. Le choix du paramètre k est crucial:
   • Pour α = 0.8, le k optimal diminue avec ε
   • Pour α = 0.95, le choix de k est relativement stable
   • Pour α = 0.6, l'estimateur est plus sensible au choix de k
2. Influence de l'intervalle de temps T:
   • Plus T est grand, meilleure est la précision d'estimation
   • L'observation à long terme aide à améliorer l'efficacité d'estimation
3. Influence de la rugosité α:
   • Lorsque α est petit (plus rugueux), la stabilité de l'estimateur diminue
   • Lorsque α est proche de 1, la méthode fonctionne plus stabilement

Comparaison avec la méthode TFE

• La méthode QMLE a une variance asymptotique plus petite, particulièrement dans le cas de grand T
• La méthode TFE ne nécessite pas d'ajustement de paramètre, mais son efficacité est inférieure

Travaux connexes

1. Théorie des équations de Volterra: Travaux fondateurs de Berger & Mizel (1980), extensions ultérieures aux noyaux singuliers, coefficients non-Lipschitz, etc.
2. Statistique des processus rugueux:
   • Modèle de volatilité rugueuse de Gatheral et al. (2018)
   • Théorie de l'inférence statistique de Chong et al. (2024)
3. Asymptotique de petit bruit:
   • Théorie classique de Kutoyants (1994)
   • Méthode d'observation discrète de Sørensen & Uchida (2003)
4. Estimation des processus de Volterra:
   • MLE d'observation continue de Zuniga (2021)
   • Méthode d'ajustement de trajectoire de Gloter & Yoshida (2025)

Conclusions et discussion

Conclusions principales

1. Construction réussie d'un estimateur QMLE basé sur des observations discrètes, réalisant le taux de convergence optimal ε⁻¹
2. L'erreur de reconstruction par inversion du noyau atteint h^(1/2), résultat potentiellement optimal
3. Établissement d'une théorie asymptotique complète, incluant la cohérence et la normalité asymptotique

Limitations

1. Conditions d'échantillonnage strictes: Nécessite h = o(ε^(1/α²)), condition très stricte lorsque α approche 1/2
2. Ajustement des paramètres: Nécessite le choix d'une valeur k appropriée, pouvant nécessiter une validation croisée en pratique
3. Hypothèses du modèle: Exige que les coefficients satisfassent les conditions de Lipschitz et de croissance polynomiale

Directions futures

1. Étude de conditions de coefficient plus générales
2. Développement de méthodes de sélection adaptative du paramètre k
3. Extension aux paramètres multidimensionnels et aux cadres non-paramétriques
4. Étude des propriétés sur échantillon fini

Évaluation approfondie

Avantages

1. Rigueur théorique: Fourniture de preuves mathématiques complètes, traitement technique raffiné
2. Innovation méthodologique: Utilisation ingénieuse de l'effet de lissage de l'inversion du noyau, surmonter les difficultés de rugosité
3. Valeur pratique: Méthode applicable aux données d'observations discrètes réelles
4. Vérification numérique: Expériences de Monte Carlo suffisantes soutenant les résultats théoriques

Insuffisances

1. Limitations des conditions: La condition d'échantillonnage h = o(ε^(1/α²)) peut être trop stricte pour les applications pratiques
2. Sélection des paramètres: Le choix de la valeur k manque de guidance théorique, dépendant de l'ajustement empirique
3. Limitations du modèle: Considère uniquement des formes spécifiques du noyau de Volterra

Impact

1. Contribution théorique: Fourniture d'un nouveau cadre théorique pour l'inférence statistique des processus rugueux
2. Valeur méthodologique: La méthode QMLE peut être généralisée à d'autres problèmes similaires
3. Perspectives d'application: Valeur d'application importante dans la modélisation de la volatilité financière et autres domaines

Scénarios d'application

1. Modélisation de la volatilité rugueuse sur les marchés financiers
2. Estimation des paramètres de systèmes physiques avec effets de mémoire
3. Modélisation de la dynamique fractionnaire en sciences de la vie
4. Applications de calcul scientifique nécessitant une estimation de paramètres haute précision

Références

L'article cite 37 travaux importants, couvrant plusieurs domaines incluant la théorie des équations de Volterra, la statistique des processus rugueux, l'asymptotique de petit bruit, etc., fournissant une base théorique solide pour la recherche.