Estimación de deriva para procesos rugosos bajo asintótica de ruido pequeño: enfoque QMLE

Información Básica

• ID del Artículo: 2510.09028
• Título: Drift estimation for rough processes under small noise asymptotic : QMLE approach
• Autores: Arnaud Gloter (Université d'Evry), Nakahiro Yoshida (University of Tokyo)
• Clasificación: math.ST (Teoría de Estadística)
• Fecha de Publicación: 13 de octubre de 2025
• Enlace del Artículo: https://arxiv.org/abs/2510.09028

Resumen

Este artículo estudia el problema de estimación de deriva para la solución X^ε de ecuaciones estocásticas de Volterra con parámetro desconocido θ⋆. El núcleo de Volterra es singular con la forma K(u) = cu^(α-1)𝟙_{u>0}, donde α ∈ (1/2, 1), y el coeficiente de difusión es proporcional a ε → 0. Basándose en observaciones de muestreo discreto de procesos de Volterra (malla h → 0), se construye un estimador de cuasi-máxima verosimilitud (QMLE). El paso principal es evaluar el error generado al reconstruir la trayectoria de semimartingala mediante la inversión del núcleo de Volterra, demostrando que este error disminuye a una tasa de h^(1/2), independientemente del valor de α. Posteriormente, se introduce una función de contraste explícita para obtener un estimador eficiente cuando ε → 0.

Antecedentes y Motivación de la Investigación

1. Problema de Investigación: El problema central que aborda este artículo es la inferencia estadística de parámetros de deriva para procesos estocásticos rugosos basada en datos de observación discreta bajo el marco de asintótica de ruido pequeño.
2. Importancia del Problema:
   • Las ecuaciones diferenciales fraccionarias se aplican ampliamente en física, finanzas matemáticas, modelado de ciencias de la vida, etc., con propiedades de memoria
   • Los procesos rugosos (α-1/2 < 1/2) se estudian ampliamente en el modelado de volatilidad de activos financieros
   • En aplicaciones prácticas solo se pueden obtener datos de observación discreta; la observación continua no es viable en la práctica
3. Limitaciones de Métodos Existentes:
   • Los métodos MLE bajo observación continua no se pueden implementar en la práctica
   • Los estimadores existentes de ajuste de trayectoria (TFE) no son óptimos
   • La regularidad de Hölder de procesos rugosos (α-1/2 < 1/2) produce errores de discretización directa más grandes
4. Motivación de la Investigación: Desarrollar métodos de estimación eficientes basados en observación discreta, superando las dificultades técnicas causadas por la rugosidad, y lograr estimación eficiente de parámetros bajo asintótica de ruido pequeño.

Contribuciones Principales

1. Se propone un nuevo método QMLE: Estimador de cuasi-máxima verosimilitud basado en inversión del núcleo de Volterra, aplicable a procesos rugosos con observación discreta
2. Se establecen cotas de error de aproximación clave: Se demuestra que el error de reconstrucción disminuye a una tasa de h^(1/2), resultado válido para todo α ∈ (1/2, 1)
3. Se proporciona teoría asintótica completa: Incluyendo consistencia y normalidad asintótica del estimador
4. Se proporcionan tasas de convergencia óptimas: El estimador converge a una tasa de ε⁻¹, con varianza asintótica igual a la inversa de la información de Fisher
5. Se realiza verificación numérica: Se validan los resultados teóricos mediante experimentos de Monte Carlo

Explicación Detallada del Método

Definición de la Tarea

Considérese la ecuación estocástica de Volterra d-dimensional:

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

Donde:

• K(u) = u^(α-1)/Γ(α) 𝟙_{u>0}, α ∈ (1/2, 1) es el núcleo singular
• θ⋆ ∈ Θ ⊂ ℝ^(d_Θ) es el parámetro de deriva desconocido
• ε → 0 es el parámetro de ruido pequeño
• Los datos observados son {X^ε_}_, h → 0

Objetivo: Estimar el parámetro θ⋆ basándose en observación discreta

Arquitectura del Modelo

1. Idea Principal: Reconstrucción de Semimartingala

Reconstruir el proceso de semimartingala mediante inversión del núcleo de Volterra:

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

Utilizando el núcleo resolvente L(u) = u^(-α)/Γ(1-α)𝟙_{u>0}, que satisface L⋆K = 1.

2. Proceso de Reconstrucción Discreta

Defínase el proceso de observación discreta X^(h)t = X{φ_h(t)}, donde φ_h(t) = h⌊t/h⌋, y constrúyase:

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

3. Construcción del Estimador

Introdúzcase un esquema de submuestreo: Sea Δ = kh, k → ∞, defínanse los incrementos:

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

Función de contraste:

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

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

Puntos de Innovación Técnica

1. Control de Error de Inversión del Núcleo: Mediante representación de convolución de núcleos de doble índice, se controla precisamente el error de g^(h)(t,u) - 1, demostrando el Lema 3.5 clave
2. Utilización del Efecto de Suavizado: Aunque la regularidad de Hölder de X^ε es α-1/2 < 1/2, el error de reconstrucción de Z^ε obtenido mediante inversión del núcleo aún puede alcanzar h^(1/2)
3. Condiciones de Selección de Parámetros: Se establece la condición h = o(ε^(1/α²)), equilibrando el error de reconstrucción y el error estadístico

Configuración Experimental

Modelo Matemático

Considérese el modelo lineal:

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

Parámetros verdaderos: (θ₀*, θ₁*) = (-1, 1)

Parámetros Experimentales

• Valores de α: 0.6, 0.8, 0.95
• Intervalo de tiempo: T = 1, 10, 50
• Niveles de ruido: ε = 1/10, 1/20, 1/100
• Pasos de muestreo: h = 10⁻², 2×10⁻²
• Parámetro de ajuste: k = 1, 2, 5, 10, 20
• Repeticiones de Monte Carlo: 1000 veces

Indicadores de Evaluación

• Media empírica
• Desviación estándar empírica (reescalada por ε⁻¹)

Resultados Experimentales

Resultados Principales

Teorema 3.1 (Cota de Error de Reconstrucción):

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

Teorema 4.1 (Consistencia): Bajo los supuestos A1-A4, θ̂^ε →^P θ⋆ cuando ε → 0.

Teorema 4.2 (Normalidad Asintótica): Bajo los supuestos A1-A6,

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

donde I es la matriz de información de Fisher, J es el término aleatorio.

Hallazgos de Experimentos Numéricos

1. La Selección del Parámetro k es Crítica:
   • Cuando α = 0.8, el k óptimo disminuye a medida que ε disminuye
   • Cuando α = 0.95, la selección de k es relativamente estable
   • Cuando α = 0.6, el estimador es más sensible a la selección de k
2. Influencia del Intervalo de Tiempo T:
   • Cuanto mayor sea T, mayor será la precisión de estimación
   • La observación a largo plazo ayuda a mejorar la eficiencia de estimación
3. Influencia de la Rugosidad α:
   • Cuando α es más pequeño (más rugoso), la estabilidad del estimador disminuye
   • Cuando α se aproxima a 1, el método funciona más estable

Comparación con el Método TFE

• El método QMLE tiene varianza asintótica más pequeña, especialmente en casos de T grande
• El método TFE no requiere ajuste de parámetros, pero tiene menor eficiencia

Trabajo Relacionado

1. Teoría de Ecuaciones de Volterra: Trabajo pionero de Berger & Mizel (1980), con extensiones posteriores a núcleos singulares, coeficientes no-Lipschitz, etc.
2. Estadística de Procesos Rugosos:
   • Modelo de volatilidad rugosa de Gatheral et al. (2018)
   • Teoría de inferencia estadística de Chong et al. (2024)
3. Asintótica de Ruido Pequeño:
   • Teoría clásica de Kutoyants (1994)
   • Método de observación discreta de Sørensen & Uchida (2003)
4. Estimación de Procesos de Volterra:
   • MLE de observación continua de Zuniga (2021)
   • Método de ajuste de trayectoria de Gloter & Yoshida (2025)

Conclusiones y Discusión

Conclusiones Principales

1. Se construye exitosamente un estimador QMLE basado en observación discreta, logrando una tasa de convergencia óptima de ε⁻¹
2. El error de reconstrucción de inversión del núcleo alcanza h^(1/2), resultado que posiblemente sea óptimo
3. Se establece teoría asintótica completa, incluyendo consistencia y normalidad asintótica

Limitaciones

1. Condiciones de Muestreo Estrictas: Se requiere h = o(ε^(1/α²)), condición muy estricta cuando α se aproxima a 1/2
2. Ajuste de Parámetros: Es necesario seleccionar un valor apropiado de k, que en la práctica puede requerir validación cruzada
3. Supuestos del Modelo: Se requiere que los coeficientes satisfagan condiciones de Lipschitz y crecimiento polinomial

Direcciones Futuras

1. Investigar condiciones de coeficientes más generales
2. Desarrollar métodos de selección adaptativa para el parámetro k
3. Extender a configuraciones de parámetros multidimensionales y no paramétricos
4. Investigar propiedades de muestras finitas

Evaluación Profunda

Ventajas

1. Teoría Rigurosa: Proporciona pruebas matemáticas completas con tratamiento técnico fino
2. Innovación Metodológica: Utiliza ingeniosamente el efecto de suavizado de inversión del núcleo, superando dificultades de rugosidad
3. Valor Práctico: El método es aplicable a datos de observación discreta reales
4. Verificación Numérica: Experimentos de Monte Carlo suficientes apoyan los resultados teóricos

Insuficiencias

1. Restricciones de Condiciones: La condición de muestreo h = o(ε^(1/α²)) puede ser demasiado estricta en aplicaciones prácticas
2. Selección de Parámetros: La selección del valor k carece de orientación teórica, dependiendo del ajuste empírico
3. Limitaciones del Modelo: Solo se considera una forma específica del núcleo de Volterra

Impacto

1. Contribución Teórica: Proporciona un nuevo marco teórico para inferencia estadística de procesos rugosos
2. Valor del Método: El método QMLE puede generalizarse a otros problemas similares
3. Perspectivas de Aplicación: Tiene valor de aplicación importante en modelado de volatilidad financiera y otros campos

Escenarios Aplicables

1. Modelado de volatilidad rugosa en mercados financieros
2. Estimación de parámetros de sistemas físicos con efectos de memoria
3. Modelado de dinámicas de orden fraccionario en ciencias de la vida
4. Aplicaciones de computación científica que requieren estimación de parámetros de alta precisión

Referencias Bibliográficas

El artículo cita 36 referencias importantes, abarcando múltiples campos incluyendo teoría de ecuaciones de Volterra, estadística de procesos rugosos, asintótica de ruido pequeño, etc., proporcionando una base teórica sólida para la investigación.