2025-11-23T09:58:19.484771

Existence, uniqueness and positivity of solutions to the Guyon-Lekeufack path-dependent volatility model with general kernels

Andrès, Jourdain
We show the existence and uniqueness of a continuous solution to a path-dependent volatility model introduced by Guyon and Lekeufack (2023) to model the price of an equity index and its spot volatility. The considered model for the trend and activity features can be written as a Stochastic Volterra Equation (SVE) with non-convolutional and non-bounded kernels as well as non-Lipschitz coefficients. We first prove the existence and uniqueness of a solution to the SVE under integrability and regularity assumptions on the two kernels and under a condition on the second kernel weighting the past squared returns which ensures that the activity feature is bounded from below by a positive constant. Then, assuming in addition that the kernel weighting the past returns is of exponential type and that an inequality relating the logarithmic derivatives of the two kernels with respect to their second variables is satisfied, we show the positivity of the volatility process which is obtained as a non-linear function of the SVE's solution. We show numerically that the choice of an exponential kernel for the kernel weighting the past returns has little impact on the quality of model calibration compared to other choices and the inequality involving the logarithmic derivatives is satisfied by the calibrated kernels. These results extend those of Nutz and Valdevenito (2023).
academic

Existencia, unicidad y positividad de soluciones al modelo de volatilidad dependiente de la trayectoria de Guyon-Lekeufack con núcleos generales

Información Básica

  • ID del Artículo: 2408.02477
  • Título: Existence, uniqueness and positivity of solutions to the Guyon-Lekeufack path-dependent volatility model with general kernels
  • Autores: Hervé Andrès (Milliman R&D, París & CERMICS, École des Ponts), Benjamin Jourdain (CERMICS, École des Ponts, INRIA)
  • Clasificación: q-fin.CP (Finanzas Computacionales)
  • Fecha de Publicación: 15 de octubre de 2025 (preimpresión arXiv)
  • Enlace del Artículo: https://arxiv.org/abs/2408.02477

Resumen

Este artículo investiga los fundamentos teóricos matemáticos del modelo de volatilidad dependiente de la trayectoria propuesto por Guyon y Lekeufack (2023). El modelo se utiliza para describir los precios de índices bursátiles y su volatilidad al contado, cuyas características de tendencia y actividad pueden expresarse como ecuaciones estocásticas de Volterra (SVE) con núcleos no convolutivos, no acotados y coeficientes no-Lipschitz. Los autores primero demuestran la existencia y unicidad de soluciones de la SVE bajo hipótesis de integrabilidad y regularidad de las funciones núcleo, así como una condición sobre el segundo núcleo que asegura un límite inferior positivo para la característica de actividad. Posteriormente, bajo la hipótesis adicional de que el núcleo que pondera los rendimientos pasados es de tipo exponencial y satisface una condición que involucra una desigualdad de derivadas logarítmicas de ambos núcleos respecto a la segunda variable, demuestran la positividad del proceso de volatilidad como función no lineal de la solución de la SVE. Los experimentos numéricos muestran que la elección de núcleos exponenciales para ponderar los rendimientos pasados tiene un impacto mínimo en la calidad de calibración del modelo, y los núcleos calibrados satisfacen la condición de desigualdad de derivadas logarítmicas.

Antecedentes y Motivación de la Investigación

Contexto del Problema

  1. Desafíos en la Modelización de Volatilidad: Los modelos tradicionales de volatilidad a menudo no capturan adecuadamente las características dependientes de la trayectoria de la volatilidad, en particular el efecto apalancamiento (aumento de volatilidad cuando los precios caen) y el fenómeno de agrupamiento de volatilidad (períodos de alta volatilidad seguidos típicamente por períodos de alta volatilidad).
  2. Éxito Empírico del Modelo Guyon-Lekeufack: El modelo de volatilidad dependiente de la trayectoria propuesto por Guyon y Lekeufack (2023) ha demostrado un desempeño excelente en investigaciones empíricas, siendo capaz de explicar más del 80% de los cambios de volatilidad en varios índices bursátiles principales, demostrando que la volatilidad es principalmente dependiente de la trayectoria.
  3. Ausencia de Teoría Matemática: A pesar del buen desempeño empírico del modelo, sus fundamentos teóricos matemáticos, en particular la garantía de existencia, unicidad y positividad de la volatilidad, aún no han sido establecidos.

Motivación de la Investigación

  1. Perfeccionamiento Teórico: Proporcionar fundamentos matemáticos rigurosos para el modelo Guyon-Lekeufack, asegurando la validez matemática del modelo.
  2. Extensión Generalizada: Generalizar los resultados existentes para núcleos específicos (como núcleos exponenciales) a categorías más generales de funciones núcleo.
  3. Consideraciones Prácticas: Asegurar que los resultados teóricos tengan significado orientador para aplicaciones prácticas, manteniendo simultáneamente el rigor matemático.

Contribuciones Principales

  1. Teorema de Existencia y Unicidad: Se demuestra la existencia y unicidad de soluciones de la ecuación estocástica de Volterra correspondiente al modelo de volatilidad dependiente de la trayectoria bajo condiciones generales de núcleos, generalizando significativamente los resultados existentes.
  2. Garantía de Positividad: Se establecen condiciones suficientes para la positividad del proceso de volatilidad, proporcionando criterios prácticos de verificación de positividad.
  3. Continuidad Hölder: Se demuestra la continuidad local Hölder de las soluciones, proporcionando apoyo teórico para métodos numéricos.
  4. Verificación Empírica: A través de experimentos de calibración con datos de múltiples índices bursátiles, se verifica que las condiciones teóricas son satisfacibles en la práctica.
  5. Contribución Metodológica: Se desarrollan técnicas de localización para tratar ecuaciones estocásticas de Volterra con coeficientes no-Lipschitz.

Explicación Detallada de Métodos

Definición de la Tarea

Considérese el siguiente modelo de volatilidad dependiente de la trayectoria: dStSt=σtdWt\frac{dS_t}{S_t} = \sigma_t dW_t

donde la volatilidad se define como: σt=β0+β1R1,t+β2R2,t\sigma_t = \beta_0 + \beta_1 R_{1,t} + \beta_2 \sqrt{R_{2,t}}

Las características de tendencia y actividad se definen respectivamente como: R1,t=Δ0K1(s,t)σsdBs+0tK1(s,t)σsdWsR_{1,t} = \int_{-\Delta}^0 K_1(s,t)\sigma_s dB_{-s} + \int_0^t K_1(s,t)\sigma_s dW_sR2,t=ΔtK2(s,t)σs2dsR_{2,t} = \int_{-\Delta}^t K_2(s,t)\sigma_s^2 ds

donde β0,β20\beta_0, \beta_2 \geq 0, β10\beta_1 \leq 0, y K1,K2K_1, K_2 son funciones núcleo de ponderación.

Arquitectura del Modelo

El modelo puede reescribirse como la siguiente ecuación estocástica de Volterra: Rt=g(t)+0tK1(s,t)γ(Rs)dWs+0tK2(s,t)b(Rs)dsR_t = g(t) + \int_0^t K_1(s,t)\gamma(R_s)dW_s + \int_0^t K_2(s,t)b(R_s)ds

donde:

  • Rt=(R1,t,R2,t)TR_t = (R_{1,t}, R_{2,t})^T
  • γ(r1,r2)=(β0+β1r1+β2r2,0)T\gamma(r_1, r_2) = (\beta_0 + \beta_1 r_1 + \beta_2\sqrt{r_2}, 0)^T
  • b(r1,r2)=(0,(β0+β1r1+β2r2)2)Tb(r_1, r_2) = (0, (\beta_0 + \beta_1 r_1 + \beta_2\sqrt{r_2})^2)^T
  • g(t)g(t) contiene términos de condiciones iniciales

Puntos de Innovación Técnica

  1. Técnica de Localización: Debido a la naturaleza no-Lipschitz de los coeficientes, los autores construyen soluciones locales utilizando secuencias de funciones de truncamiento {bn,γn}\{b_n, \gamma_n\}, demostrando posteriormente que las soluciones locales son en realidad globales.
  2. Extensión del Principio de Comparación: Se generaliza el teorema de comparación de ecuaciones diferenciales estocásticas a procesos de Itô, utilizándose para demostrar el límite inferior positivo de la volatilidad.
  3. Tratamiento de Núcleos Generales: Mediante condiciones apropiadas de integrabilidad y regularidad, se tratan núcleos no convolutivos que pueden ser no acotados.

Configuración Experimental

Conjunto de Datos

  • Datos de Índices Bursátiles: S&P 500, Stoxx Europe 50, índice FTSE 100
  • Datos de Volatilidad: Índices de volatilidad implícita VIX, VSTOXX, IVI; volatilidad realizada calculada a partir de rendimientos de 5 minutos
  • Período Temporal: Conjunto de entrenamiento del 1 de enero de 2000 al 31 de diciembre de 2018; conjunto de prueba del 1 de enero de 2019 al 15 de mayo de 2022

Selección de Núcleos

Se compararon cuatro combinaciones de funciones núcleo:

  1. Dos núcleos de ley de potencia con cambio temporal (TSPL)
  2. Dos núcleos exponenciales
  3. Combinación convexa de dos núcleos exponenciales
  4. K1K_1 núcleo exponencial, K2K_2 núcleo TSPL

Métricas de Evaluación

Se utiliza el coeficiente de determinación R2R^2 para evaluar la calidad del ajuste del modelo.

Resultados Experimentales

Resultados Principales

Tabla 1: Puntuaciones R2R^2 para diferentes selecciones de núcleos

Selección de NúcleosVIX Entreno/PruebaVSTOXX Entreno/PruebaIVI Entreno/PruebaSPX Entreno/PruebaSTOXX Entreno/PruebaFTSE Entreno/Prueba
Dos núcleos TSPL89.34%/83.43%91.36%/90.91%92.13%/87.97%67.03%/68.84%58.86%/62.48%61.74%/62.72%
Dos núcleos exponenciales89.10%/74.26%90.28%/86.43%90.50%/82.37%66.99%/63.41%58.53%/61.68%61.54%/62.48%
Combinación convexa exponencial90.01%/83.01%91.99%/91.61%92.37%/88.74%67.55%/68.86%59.10%/62.68%62.14%/62.80%
Exponencial + TSPL89.86%/81.33%91.76%/91.03%92.29%/87.67%67.12%/64.61%58.85%/62.78%61.74%/62.88%

Hallazgos Clave

  1. Desempeño Similar: El desempeño de las cuatro selecciones de funciones núcleo es muy similar, indicando que la elección de núcleos exponenciales tiene un impacto limitado en la calidad del ajuste.
  2. Satisfacción de Condiciones de Positividad: Para la combinación exponencial + TSPL, todos los resultados calibrados satisfacen la condición de positividad 2λδζ2\lambda\delta \geq \zeta.
  3. Integración de Teoría y Práctica: La combinación exponencial + TSPL logra un efecto de ajuste comparable al de otras selecciones mientras garantiza teóricamente la positividad.

Trabajo Relacionado

Teoría de Ecuaciones Estocásticas de Volterra

  • Resultados Clásicos: Trabajo pionero de Berger y Mizel (1980), estableciendo existencia y unicidad para SVE con coeficientes Lipschitz
  • Núcleos Singulares: Cochran et al. (1995) y Coutin-Decreusefond (2001) tratan casos con núcleos no acotados
  • Coeficientes No-Lipschitz: Generalización de Wang (2008) bajo condiciones de Yamada-Watanabe

Modelos de Volatilidad Dependiente de la Trayectoria

  • Modelo Guyon-Lekeufack: Propuesta inicial y verificación de la importancia de la dependencia de trayectoria
  • Aproximación Markoviana: Resultados de existencia y positividad de Nutz y Riveros Valdevenito (2024) para el caso de núcleos exponenciales
  • Volatilidad Rugosa: Desarrollo teórico relacionado de SVE, como el trabajo de Abi Jaber et al.

Resultados Teóricos

Teorema 1 (Existencia y Unicidad)

Bajo las hipótesis (I.1)-(I.6), la ecuación (6) posee una única solución continua, y la solución posee continuidad local γ\gamma^*-Hölder, donde γ(0,min(γ,12α1,1α2))\gamma^* \in (0, \min(\gamma, \frac{1}{2\alpha_1^*}, \frac{1}{\alpha_2^*})).

Hipótesis Clave:

  • (I.1)-(I.2): Condiciones de integrabilidad de funciones núcleo
  • (I.3): Acotación de condiciones iniciales
  • (I.4): Integrabilidad de orden superior de funciones núcleo
  • (I.5): Continuidad Hölder de funciones núcleo
  • (I.6): Condición de límite inferior positivo para la característica de actividad

Teorema 2 (Positividad)

Bajo las hipótesis (I.1)-(I.5) e (II.1)-(II.4), el proceso de volatilidad σt\sigma_t posee un límite inferior positivo: σtσ0exp(β10tK1(s,s)dWs+h(t)h(0)12β120tK1(s,s)2ds)\sigma_t \geq \sigma_0 \exp\left(\beta_1\int_0^t K_1(s,s)dW_s + h(t) - h(0) - \frac{1}{2}\beta_1^2\int_0^t K_1(s,s)^2 ds\right)

Hipótesis Adicionales:

  • (II.2): Estructura de tipo exponencial K1(s,t)=f(s)eh(t)K_1(s,t) = f(s)e^{h(t)}
  • (II.3): Relación de derivadas logarítmicas entre núcleos: tK2(s,t)2h(t)K2(s,t)0\partial_t K_2(s,t) - 2h'(t)K_2(s,t) \geq 0

Conclusiones y Discusión

Conclusiones Principales

  1. Perfeccionamiento Teórico: Se establece por primera vez un marco teórico matemático completo para el modelo Guyon-Lekeufack con funciones núcleo generales.
  2. Orientación Práctica: Se proporcionan criterios teóricos para la selección de funciones núcleo, en particular condiciones de garantía de positividad de volatilidad.
  3. Verificación Numérica: Las condiciones teóricas son satisfacibles en calibraciones prácticas, demostrando el valor práctico de la teoría.

Limitaciones

  1. Restricciones de Núcleos: Los resultados de positividad requieren que K1K_1 sea de tipo exponencial, limitando la flexibilidad en la selección de funciones núcleo.
  2. Problema de Truncamiento: La teoría no es aplicable a funciones núcleo con rezagos finitos truncados, aunque estos se utilizan frecuentemente en la práctica.
  3. Complejidad Computacional: La implementación numérica de núcleos generales es más compleja que en el caso Markoviano.

Direcciones Futuras

  1. Relajación de Condiciones de Positividad: Búsqueda de condiciones más generales para garantizar la positividad de volatilidad.
  2. Métodos Numéricos: Desarrollo de algoritmos de solución numérica eficientes.
  3. Extensión Multifactorial: Generalización a modelos de factores de mayor dimensión.

Evaluación Profunda

Fortalezas

  1. Rigor Teórico: El razonamiento matemático es riguroso, las condiciones de hipótesis son claras y razonables, y las técnicas de demostración son avanzadas.
  2. Valor Práctico: Los resultados teóricos guían directamente la aplicación práctica, proporcionando criterios operables para la selección de funciones núcleo.
  3. Innovación Metodológica: Las técnicas de localización y la extensión del principio de comparación proporcionan nuevas herramientas para abordar problemas similares.
  4. Apoyo Empírico: Experimentos numéricos suficientes verifican la efectividad y practicidad de la teoría.

Deficiencias

  1. Limitaciones de Hipótesis: Algunas hipótesis (como funciones núcleo de tipo exponencial) pueden ser demasiado estrictas en la práctica.
  2. Complejidad Técnica: Las técnicas de demostración son relativamente complejas, lo que puede limitar la aplicación generalizada de los resultados.
  3. Tratamiento de Truncamiento: No se aborda adecuadamente el problema de truncamiento de funciones núcleo que es común en la práctica.

Impacto

  1. Contribución Académica: Se realizan contribuciones importantes a la teoría de ecuaciones estocásticas de Volterra y al campo de las finanzas matemáticas.
  2. Orientación Práctica: Se proporciona base teórica para la aplicación práctica de modelos de volatilidad dependiente de trayectoria.
  3. Valor Metodológico: Las técnicas desarrolladas tienen valor de referencia para problemas relacionados.

Escenarios de Aplicabilidad

Esta investigación es particularmente aplicable a:

  • Modelización de volatilidad de acciones e índices
  • Gestión del riesgo de volatilidad en fijación de precios de derivados
  • Modelos de riesgo financiero que requieren considerar dependencia de trayectoria
  • Investigación teórica de ecuaciones estocásticas de Volterra

Referencias Bibliográficas

Las referencias bibliográficas principales incluyen:

  • Guyon, J. and Lekeufack, J. (2023). Volatility is (mostly) path-dependent. Quantitative Finance.
  • Nutz, M. and Riveros Valdevenito, A. (2024). On the Guyon-Lekeufack volatility model. Finance and Stochastics.
  • Zhang, X. (2010). Stochastic Volterra equations in Banach spaces. Journal of Functional Analysis.