σμ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.3·~55 minAvanzado

VAR y causalidad de Granger

Posición en la rutaMódulo 2 de 3 · Tema 3 de 5
2.3 VAR y causalidad de Granger
—Puntuación del lab— / 100XP total—Racha—

Variables que se predicen entre sí

Una autorregresión vectorial, VAR, modela varias series juntas. Cada variable se regresa sobre los valores pasados de todas, así que para dos series xxx e yyy y un rezago, yt=c+A yt−1+uty_t=c+A\,y_{t-1}+u_tyt​=c+Ayt−1​+ut​ con yt=(xt,yt)⊤y_t=(x_t,y_t)^\topyt​=(xt​,yt​)⊤ y una matriz AAA de 2×22\times22×2. Cada ecuación se ajusta por mínimos cuadrados ordinarios, y en statsmodels VAR(data).fit(p) las hace todas a la vez. El modelo es estable cuando todos los valores propios de AAA están dentro del círculo unidad, lo que significa que los choques se apagan en lugar de explotar. Los pronósticos salen por multiplicaciones repetidas. Para elegir el orden ppp se comparan criterios de información como AIC y BIC, que premian el ajuste y penalizan los coeficientes extra.

La causalidad de Granger hace una pregunta estrecha: ¿ayudan los valores pasados de yyy a predecir xxx una vez conocidos los valores pasados del propio xxx? Se ajusta la ecuación de xxx dos veces, con los rezagos de yyy (sin restricción, suma de cuadrados de residuos RSSu\mathrm{RSS}_uRSSu​) y sin ellos (restringida, RSSr\mathrm{RSS}_rRSSr​), y se compara con un estadístico F. Un resultado significativo dice que yyy precede a xxx en el tiempo. No prueba que yyy cause xxx: una tercera variable, o inversores que anticipan, pueden producir el mismo patrón.

yt=c+A1yt−1+⋯+Apyt−p+ut,F=(RSSr−RSSu)/qRSSu/(n−k),∣λi(A)∣<1y_t=c+A_1y_{t-1}+\cdots+A_py_{t-p}+u_t,\qquad F=\frac{(\mathrm{RSS}_r-\mathrm{RSS}_u)/q}{\mathrm{RSS}_u/(n-k)},\qquad |\lambda_i(A)|<1yt​=c+A1​yt−1​+⋯+Ap​yt−p​+ut​,F=RSSu​/(n−k)(RSSr​−RSSu​)/q​,∣λi​(A)∣<1(1)
Un VAR de dos variables a manoSea A=(0,50,20,10,6)A=\begin{pmatrix}0{,}5&0{,}2\\0{,}1&0{,}6\end{pmatrix}A=(0,50,1​0,20,6​) y partamos de (x,y)=(1,2)(x,y)=(1,2)(x,y)=(1,2). Un paso adelante es (0,5⋅1+0,2⋅2;  0,1⋅1+0,6⋅2)=(0,9;1,3)(0{,}5\cdot1+0{,}2\cdot2;\;0{,}1\cdot1+0{,}6\cdot2)=(0{,}9;1{,}3)(0,5⋅1+0,2⋅2;0,1⋅1+0,6⋅2)=(0,9;1,3) y dos pasos (0,5⋅0,9+0,2⋅1,3;  0,1⋅0,9+0,6⋅1,3)=(0,71;0,87)(0{,}5\cdot0{,}9+0{,}2\cdot1{,}3;\;0{,}1\cdot0{,}9+0{,}6\cdot1{,}3)=(0{,}71;0{,}87)(0,5⋅0,9+0,2⋅1,3;0,1⋅0,9+0,6⋅1,3)=(0,71;0,87). La traza es 1,11{,}11,1 y el determinante 0,5⋅0,6−0,2⋅0,1=0,280{,}5\cdot0{,}6-0{,}2\cdot0{,}1=0{,}280,5⋅0,6−0,2⋅0,1=0,28, así que los valores propios resuelven λ2−1,1λ+0,28=0\lambda^2-1{,}1\lambda+0{,}28=0λ2−1,1λ+0,28=0 y valen 0,70{,}70,7 y 0,40{,}40,4: ambos dentro del círculo unidad, el sistema es estable. Para la prueba F, si quitar los rezagos de yyy sube la suma de cuadrados de RSSu=10\mathrm{RSS}_u=10RSSu​=10 a RSSr=12\mathrm{RSS}_r=12RSSr​=12, con q=2q=2q=2 restricciones y n−k=40n-k=40n−k=40, entonces F=(12−10)/210/40=4,0F=\frac{(12-10)/2}{10/40}=4{,}0F=10/40(12−10)/2​=4,0, por encima del valor crítico del 5 % de 3,233{,}233,23.

Con datos mensuales reales, 258258258 meses desde marzo de 2005, tome el rendimiento de SPY en porcentaje y el diferencial del Tesoro a 10 años menos 2 años. El AIC elige 333 rezagos y el BIC 111, y un VAR(2) es estable, con la raíz más cercana en 1,021{,}021,02 porque el diferencial es muy persistente (sus propios rezagos suman 0,980{,}980,98). La prueba de Granger de si el diferencial ayuda a predecir los rendimientos da F=0,22F=0{,}22F=0,22 (p=0,81p=0{,}81p=0,81), y de si los rendimientos ayudan a predecir el diferencial da F=1,88F=1{,}88F=1,88 (p=0,15p=0{,}15p=0,15). Ninguna dirección es significativa: conocer la pendiente de la curva no aporta nada a un pronóstico mensual de SPY.

Lab: qué se califica

  • Un paso adelante10 pts
  • Dos pasos adelante15 pts
  • El mayor valor propio15 pts
  • Un estadístico F de Granger15 pts
  • Contar coeficientes15 pts
  • Leer una prueba de Granger10 pts
  • Predecir no es causar20 pts
0/ 7 respondidos
Ejercicios
1Respuesta numérica10 ptsSin responder

Con A=(0,50,20,10,6)A=\begin{pmatrix}0{,}5&0{,}2\\0{,}1&0{,}6\end{pmatrix}A=(0,50,1​0,20,6​) y valores actuales (x,y)=(1,2)(x,y)=(1,2)(x,y)=(1,2), ¿cuánto vale el pronóstico a un paso de xxx, 0,5x+0,2y0{,}5x+0{,}2y0,5x+0,2y?

h1=h1 =h1=aceptado: ± 0.01
2Respuesta numérica15 ptsSin responder

El pronóstico a un paso es (0,9;1,3)(0{,}9;1{,}3)(0,9;1,3). ¿Cuánto vale el pronóstico a dos pasos de yyy, 0,1⋅0,9+0,6⋅1,30{,}1\cdot0{,}9+0{,}6\cdot1{,}30,1⋅0,9+0,6⋅1,3? Dé dos decimales.

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

Los valores propios resuelven λ2−1,1λ+0,28=0\lambda^2-1{,}1\lambda+0{,}28=0λ2−1,1λ+0,28=0, con traza 1,11{,}11,1 y determinante 0,280{,}280,28. ¿Cuál es la raíz mayor? Use 1,12−4⋅0,28=0,3\sqrt{1{,}1^2-4\cdot0{,}28}=0{,}31,12−4⋅0,28​=0,3 y dé un decimal.

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

Quitar los rezagos de yyy sube la suma de cuadrados de residuos de RSSu=10\mathrm{RSS}_u=10RSSu​=10 a RSSr=12\mathrm{RSS}_r=12RSSr​=12. Con q=2q=2q=2 restricciones y n−k=40n-k=40n−k=40, ¿cuánto vale F=(RSSr−RSSu)/qRSSu/(n−k)F=\frac{(\mathrm{RSS}_r-\mathrm{RSS}_u)/q}{\mathrm{RSS}_u/(n-k)}F=RSSu​/(n−k)(RSSr​−RSSu​)/q​?

F=F =F=aceptado: ± 0.05
5Respuesta numérica15 ptsSin responder

Un VAR(2) tiene 333 variables. Cada una de sus 333 ecuaciones tiene una constante más 222 rezagos de cada una de las 333 variables. ¿Cuántos coeficientes tiene el modelo entero?

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

La prueba de si el diferencial de rendimiento causa en el sentido de Granger los rendimientos mensuales de SPY da F = 0,22 con p = 0,81. ¿Qué dice?

7Explicar con palabras20 ptsSin responder

Suponga que yyy causa xxx en el sentido de Granger con un valor p muy pequeño. En una o dos frases, explique por qué eso sigue sin probar que yyy mueve a xxx.

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