σμ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 statsmodels›Módulo 2 · Modelos de series temporales›Tema 2.5·~55 minAvanzado

Cointegración y pares

Posición en la rutaMódulo 2 de 3 · Tema 5 de 5
2.5 Cointegración y pares
—Puntuación del lab— / 100XP total—Racha—

Dos precios que permanecen juntos

Los precios vagan: se comportan como paseos aleatorios, así que cada uno es no estacionario. Dos series que vagan están cointegradas cuando alguna combinación de ellas no vaga: el diferencial st=yt−βxts_t=y_t-\beta x_tst​=yt​−βxt​ vuelve siempre hacia una media. Esa es la idea del trading de pares. La prueba de Engle–Granger tiene dos pasos. Primero se regresa yyy sobre xxx por mínimos cuadrados para obtener la razón de cobertura β^\hat\betaβ^​ y el diferencial. Luego se hace una prueba de Dickey–Fuller aumentada sobre el diferencial: Δst=γst−1+εt\Delta s_t=\gamma s_{t-1}+\varepsilon_tΔst​=γst−1​+εt​, y se pregunta si γ<0\gamma<0γ<0. El estadístico t debe estar por debajo de valores críticos especiales, cerca de −3,34-3{,}34−3,34 al 5 % para dos series, más estrictos que los habituales porque el diferencial fue estimado. En statsmodels es coint(y, x).

st=yt−β^xt,Δst=γst−1+εt,ϕ=1+γ,t1/2=ln⁡0.5ln⁡ϕs_t=y_t-\hat\beta x_t,\qquad \Delta s_t=\gamma s_{t-1}+\varepsilon_t,\qquad \phi=1+\gamma,\qquad t_{1/2}=\frac{\ln 0.5}{\ln\phi}st​=yt​−β^​xt​,Δst​=γst−1​+εt​,ϕ=1+γ,t1/2​=lnϕln0.5​(1)
Un diferencial de cinco días a manoTomemos x=1,2,3,4,5x=1,2,3,4,5x=1,2,3,4,5 e y=3,5,6,9,10y=3,5,6,9,10y=3,5,6,9,10. La regresión da β^=18/10=1,8\hat\beta=18/10=1{,}8β^​=18/10=1,8, así que el diferencial y−1,8xy-1{,}8xy−1,8x es 1,2;1,4;0,6;1,8;1,01{,}2;1{,}4;0{,}6;1{,}8;1{,}01,2;1,4;0,6;1,8;1,0, con media 1,21{,}21,2 y desviación típica 0,4470{,}4470,447. El diferencial de hoy 1,01{,}01,0 tiene z-score (1,0−1,2)/0,447=−0,45(1{,}0-1{,}2)/0{,}447=-0{,}45(1,0−1,2)/0,447=−0,45: algo por debajo de lo normal. Si el diferencial sigue st=0,9 st−1+εs_t=0{,}9\,s_{t-1}+\varepsilonst​=0,9st−1​+ε, una brecha se cierra a la mitad en ln⁡0,5/ln⁡0,9=6,58\ln0{,}5/\ln0{,}9=6{,}58ln0,5/ln0,9=6,58 días. Una regresión de prueba Δs=γst−1\Delta s=\gamma s_{t-1}Δs=γst−1​ con γ^=−0,1\hat\gamma=-0{,}1γ^​=−0,1 y error estándar 0,040{,}040,04 tiene t=−2,5t=-2{,}5t=−2,5, por encima del crítico −3,34-3{,}34−3,34, así que no se puede llamar cointegrado al par. Con datos diarios reales desde 2000 (6 6966\,6966696 días), regresar un log-precio sobre otro da un R2R^2R2 de 0,850{,}850,85 a 0,930{,}930,93 en pares como KO y PEP, XOM y CVX o HD y LOW, y sin embargo los estadísticos de Engle–Granger son −2,69-2{,}69−2,69, −1,80-1{,}80−1,80 y −3,97-3{,}97−3,97: solo HD y LOW pasa al 5 % (p=0,008p=0{,}008p=0,008). Incluso eso es frágil: cada mitad de la muestra por separado falla (t=−1,58t=-1{,}58t=−1,58 y −2,63-2{,}63−2,63) y la razón de cobertura pasa de 0,130{,}130,13 a 0,930{,}930,93. El diferencial de toda la muestra tiene ϕ=0,998\phi=0{,}998ϕ=0,998, una vida media de unos 332332332 días, demasiado lenta para operar. Como comparación, en 100010001000 pares simulados de paseos aleatorios independientes, el 85%85\%85% muestra una pendiente con ∣t∣>1,96|t|>1{,}96∣t∣>1,96.
Un R cuadrado alto no es evidenciaDos precios que vagan siempre parecen relacionados al regresar uno sobre otro. Solo la estacionariedad del diferencial, probada con los valores críticos de Engle–Granger, dice que el vínculo es real. Compruebe la estabilidad en subperiodos antes de operarlo.
  • Dos pasos: regresar para la razón de cobertura y luego probar una raíz unitaria en el diferencial con los valores de Engle–Granger.
  • La vida media dice si una reversión real a la media es lo bastante rápida para operarse.

Lab: qué se califica

  • La razón de cobertura10 pts
  • El diferencial medio10 pts
  • Un z-score15 pts
  • Una vida media15 pts
  • Un estadístico de Dickey–Fuller15 pts
  • Leer la prueba10 pts
  • Pasar no basta25 pts
0/ 7 respondidos
Ejercicios
1Respuesta numérica10 ptsSin responder

Para x=1,2,3,4,5x=1,2,3,4,5x=1,2,3,4,5 e y=3,5,6,9,10y=3,5,6,9,10y=3,5,6,9,10 las sumas son Sxx=10S_{xx}=10Sxx​=10 y Sxy=18S_{xy}=18Sxy​=18. ¿Cuánto vale la razón de cobertura β^=Sxy/Sxx\hat\beta=S_{xy}/S_{xx}β^​=Sxy​/Sxx​?

b=b =b=aceptado: ± 0.01
2Respuesta numérica10 ptsSin responder

Con β^=1,8\hat\beta=1{,}8β^​=1,8 el diferencial y−1,8xy-1{,}8xy−1,8x es 1,2;1,4;0,6;1,8;1,01{,}2;1{,}4;0{,}6;1{,}8;1{,}01,2;1,4;0,6;1,8;1,0. ¿Cuál es su media?

m=m =m=aceptado: ± 0.01
3Respuesta numérica15 ptsSin responder

El diferencial de hoy es 1,01{,}01,0, la media 1,21{,}21,2 y la desviación típica 0,4470{,}4470,447. ¿Cuánto vale el z-score (1,0−1,2)/0,447(1{,}0-1{,}2)/0{,}447(1,0−1,2)/0,447? Dé dos decimales.

z=z =z=aceptado: ± 0.01
4Respuesta numérica15 ptsSin responder

Un diferencial sigue st=0,9 st−1+εts_t=0{,}9\,s_{t-1}+\varepsilon_tst​=0,9st−1​+εt​. ¿Después de cuántos días se reduce una brecha a la mitad, ln⁡0,5/ln⁡0,9\ln0{,}5/\ln0{,}9ln0,5/ln0,9? Use ln⁡0,5=−0,6931\ln0{,}5=-0{,}6931ln0,5=−0,6931, ln⁡0,9=−0,1054\ln0{,}9=-0{,}1054ln0,9=−0,1054 y dé un decimal.

hl=hl =hl=aceptado: ± 0.1
5Respuesta numérica15 ptsSin responder

En Δst=γst−1+εt\Delta s_t=\gamma s_{t-1}+\varepsilon_tΔst​=γst−1​+εt​ la estimación es γ^=−0,1\hat\gamma=-0{,}1γ^​=−0,1 con error estándar 0,040{,}040,04. ¿Cuánto vale el estadístico t? Dé un decimal.

t=t =t=aceptado: ± 0.05
6Elección10 ptsSin responder

El estadístico es -2,5 y el valor crítico de Engle–Granger al 5 % es -3,34. ¿Qué concluye?

7Explicar con palabras25 ptsSin responder

HD y LOW pasan la prueba de Engle–Granger en toda la muestra, pero cada mitad falla y la razón de cobertura pasa de 0,130{,}130,13 a 0,930{,}930,93. En una o dos frases, explique por qué eso vuelve arriesgado operar el par.

Modo matemáticasEntrada por teclado o LaTeXEjercicio 1 de 7 · 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