σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
σμimplied vol · KKmax(S−K,0)udβρᵢⱼfat tailsPnL · drawdownyield · TdSₜ = μSₜ dt + σSₜ dWₜC = S·N(d₁) − K·e⁻ʳᵀ·N(d₂)σₚ² = wᵀΣwd₁ = [ln(S/K) + (r + σ²⁄2)T] ⁄ σ√T∂V/∂t + ½σ²S²·∂²V/∂S² + rS·∂V/∂S − rV = 0Sharpe = (Rₚ − R_f) ⁄ σₚVaR₉₉ = μ − 2.33σρ = Cov(X,Y) ⁄ σₓσᵧE[Rₚ] = Σ wᵢ·E[Rᵢ]φ(x) = e^(−x²/2) ⁄ √2πE[Rᵢ] = R_f + βᵢ(E[Rₘ] − R_f)dXₜ = θ(μ − Xₜ)dt + σ dWₜσ²ₜ = ω + α·ε²ₜ₋₁ + β·σ²ₜ₋₁C − P = S − K·e⁻ʳᵀΔ = ∂V/∂S Γ = ∂²V/∂S²p = (e^(rΔt) − d) ⁄ (u − d)rₜ = ln(Pₜ ⁄ Pₜ₋₁)E[Xₜ₊₁ | ℱₜ] = Xₜf* = (bp − q) ⁄ bz = (x − μ) ⁄ σP(N=k) = λᵏe^(−λ) ⁄ k!Θ = −∂V/∂Tν = S·φ(d₁)·√Tmin wᵀΣw s.t. wᵀμ = μ*∫Σ∂∞λ
uantCraft AcademyHojas de rutaRutas
Ctrl K
AO
+1,2 k estudiantes
Iniciar sesiónUnirse ↗
Hojas de rutaRutas
AO
+1,2 k estudiantes
Iniciar sesión
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
Ruta de Particle filtering and non-stationary time series›Módulo 2 · Filtros de partículas›Tema 2.4·~1 h 10 minAvanzado

Filtrado de la volatilidad estocástica

Posición en la rutaMódulo 2 de 3 · Tema 4 de 5
2.4 Filtrado de la volatilidad estocástica
—Puntuación del lab— / 100XP total—Racha—

La volatilidad que no se ve

Los rendimientos vienen en períodos tranquilos y en tormentas, así que la dispersión de un rendimiento cambia de un día a otro y nadie la publica. El modelo de volatilidad estocástica trata el logaritmo de la varianza, hth_tht​, como un estado oculto. Se mueve como una autorregresión de primer orden en torno a un nivel de largo plazo μ\muμ: una fracción φ\varphiφ de la brecha de ayer permanece y llega un choque nuevo de tamaño σ\sigmaσ. El rendimiento es una campana con varianza ehte^{h_t}eht​. El logaritmo mantiene la varianza positiva y hace que los choques la multipliquen. Tres números describen el modelo: el nivel μ\muμ, la persistencia φ\varphiφ (cerca de 111, los choques de volatilidad duran) y la volatilidad de la volatilidad σ\sigmaσ.

ht=μ+φ (ht−1−μ)+σ ηt,rt=eht/2 εt,sd(h)=σ1−φ2h_t=\mu+\varphi\,(h_{t-1}-\mu)+\sigma\,\eta_t,\qquad r_t=e^{h_t/2}\,\varepsilon_t,\qquad \mathrm{sd}(h)=\frac{\sigma}{\sqrt{1-\varphi^{2}}}ht​=μ+φ(ht−1​−μ)+σηt​,rt​=eht​/2εt​,sd(h)=1−φ2​σ​(1)

El estado oculto está dentro de una exponencial, así que la observación es una función no lineal de él y el filtro de Kalman no tiene respuesta exacta. Un filtro de partículas bootstrap sí. Se parte de partículas de la distribución de largo plazo de hhh. Cada día se mueve cada partícula con la autorregresión, se pondera con la densidad normal del rendimiento con varianza ehe^{h}eh y se remuestrea cuando el tamaño muestral efectivo es bajo. La volatilidad filtrada es la raíz del promedio ponderado de ehe^{h}eh, y multiplicar por 252\sqrt{252}252​ pasa de diario a anual. El promedio de las verosimilitudes de las partículas da además la log-verosimilitud de los parámetros, que usa el siguiente tema.

Las cifras de un díaTomemos μ=−8,85\mu=-8{,}85μ=−8,85 (volatilidad anual del 19%19\%19%), φ=0,98\varphi=0{,}98φ=0,98 y σ=0,15\sigma=0{,}15σ=0,15. La dispersión de largo plazo de hhh es 0,15/1−0,982=0,7540{,}15/\sqrt{1-0{,}98^2}=0{,}7540,15/1−0,982​=0,754. Una partícula en h=−8h=-8h=−8 significa una varianza diaria e−8e^{-8}e−8 y una volatilidad anual 252 e−8=29,1%\sqrt{252\,e^{-8}}=29{,}1\%252e−8​=29,1%. Su predicción para mañana tiene media −8,85+0,98×0,85=−8,017-8{,}85+0{,}98\times0{,}85=-8{,}017−8,85+0,98×0,85=−8,017. Si el rendimiento es −2%-2\%−2%, su densidad normal vale 12,012{,}012,0. Un rendimiento de −3%-3\%−3% es 505050 veces más probable con h=−8h=-8h=−8 que con una partícula tranquila en h=−9,5h=-9{,}5h=−9,5, así que la partícula tormentosa recibe 505050 veces el peso.

Con los 5 4045\,4045404 rendimientos de SPY desde 2005, con μ=ln⁡(r2‾)=−8,85\mu=\ln(\overline{r^{2}})=-8{,}85μ=ln(r2)=−8,85, φ=0,98\varphi=0{,}98φ=0,98, σ=0,15\sigma=0{,}15σ=0,15 y 100010001000 partículas, la volatilidad anual filtrada promedia 8,7%8{,}7\%8,7% en 2017 y 30,8%30{,}8\%30,8% en 2008, con un máximo de 79,5%79{,}5\%79,5% el 15 de octubre de 2008. La log-verosimilitud es 17 75217\,75217752, frente a 16 25216\,25216252 con volatilidad constante, una ganancia de unos 1 5001\,5001500. El tamaño muestral efectivo tiene mediana 966966966 de 100010001000, pero cae a 3,33{,}33,3 el 10 de octubre de 2018 (−3,2%-3{,}2\%−3,2% tras un período tranquilo, con la volatilidad filtrada saltando de 7,9%7{,}9\%7,9% a 18,9%18{,}9\%18,9%). Tras la calma, casi ninguna partícula tiene una volatilidad lo bastante alta para explicar un choque.

  • El estado oculto es el log-varianza: mantiene la varianza positiva y hace que los choques la multipliquen.
  • Los pesos favorecen la volatilidad que explica el rendimiento: un rendimiento grande eleva las partículas tormentosas y uno pequeño las tranquilas.
  • Vigile el ESS tras la calma: una gran sorpresa puede dejar muy pocas partículas útiles.

Lab: qué se califica

  • La dispersión de largo plazo15 pts
  • Volatilidad anual de una partícula10 pts
  • La media prevista15 pts
  • Una densidad normal10 pts
  • Un cociente de pesos15 pts
  • Por qué el log-varianza10 pts
  • Emparejar el vocabulario10 pts
  • Por qué no usar solo el rendimiento15 pts
0/ 8 respondidos
Ejercicios
1Respuesta numérica15 ptsSin responder

Con φ=0,98\varphi=0{,}98φ=0,98 y σ=0,15\sigma=0{,}15σ=0,15, la desviación típica del log-varianza a largo plazo es σ/1−φ2\sigma/\sqrt{1-\varphi^{2}}σ/1−φ2​. ¿Cuánto vale? Dé tres decimales.

s=s =s=aceptado: ± 0.003
2Respuesta numérica10 ptsSin responder

Una partícula tiene h=−8h=-8h=−8, así que su varianza diaria es e−8=0,000335e^{-8}=0{,}000335e−8=0,000335. ¿Cuál es su volatilidad anual 252×0,000335\sqrt{252\times0{,}000335}252×0,000335​ en porcentaje? Dé un decimal.

v=v =v=aceptado: ± 0.3
3Respuesta numérica15 ptsSin responder

Una partícula tiene h=−8h=-8h=−8 hoy, con μ=−8,85\mu=-8{,}85μ=−8,85 y φ=0,98\varphi=0{,}98φ=0,98. ¿Cuál es la media del hhh de mañana, μ+φ(h−μ)\mu+\varphi(h-\mu)μ+φ(h−μ)? Dé tres decimales.

h=h =h=aceptado: ± 0.003
4Respuesta numérica10 ptsSin responder

Con h=−8h=-8h=−8, el rendimiento tiene desviación típica 0,01830{,}01830,0183. La densidad normal de un rendimiento de −0,02-0{,}02−0,02 es e−z2/20,01832π\dfrac{e^{-z^{2}/2}}{0{,}0183\sqrt{2\pi}}0,01832π​e−z2/2​ con z=−0,02/0,0183z=-0{,}02/0{,}0183z=−0,02/0,0183. Dé su valor con un decimal.

f=f =f=aceptado: ± 0.2
5Respuesta numérica15 ptsSin responder

Un rendimiento de −3%-3\%−3% tiene densidad 5,705{,}705,70 con la partícula tormentosa (h=−8h=-8h=−8) y 0,1130{,}1130,113 con la tranquila (h=−9,5h=-9{,}5h=−9,5). ¿Cuántas veces mayor es el peso de la partícula tormentosa? Dé el cociente al entero más cercano.

n=n =n=aceptado: ± 1.5
6Elección10 ptsSin responder

¿Por qué el modelo pone un proceso tipo paseo aleatorio sobre el logaritmo de la varianza y no sobre la varianza misma?

7Emparejar10 ptsSin responder

Empareje cada parámetro o cantidad con su significado.

Elige la pareja de cada elemento de la izquierda.

Nivel mu
Persistencia phi
Volatilidad de la volatilidad sigma
Log-varianza oculto h
8Explicar con palabras15 ptsSin responder

El 16 de marzo de 2020 SPY cayó un 11,6%11{,}6\%11,6%, pero la volatilidad anual filtrada fue del 69%69\%69% y no del 184%184\%184%. En una o dos frases, explique por qué el filtro no equipara la volatilidad con el tamaño del rendimiento de un solo día.

Modo matemáticasEntrada por teclado o LaTeXEjercicio 1 de 8 · Ctrl ⏎ para comprobar

Inicia sesión para ejecutar tu código.

Todas las clases se pueden leer gratis. Crea una cuenta o inicia sesión para ejecutar código, lanzar las pruebas y enviar este lab.

  • IDE en el navegador con terminal bash — sin instalar nada
  • Corrección instantánea, XP y tu lugar en la ruta
Crear una cuentaIniciar sesión
O
Continuar con GitHub