σμ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 AcademyRoadmapsSkills
Ctrl K
AO
+1,2 Tsd. Lernende
AnmeldenBeitreten ↗
RoadmapsSkills
AO
+1,2 Tsd. Lernende
Anmelden
  • 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 AcademyRoadmapsSkills
Ctrl K
AO
+1,2 Tsd. Lernende
AnmeldenBeitreten ↗
RoadmapsSkills
AO
+1,2 Tsd. Lernende
Anmelden
  • English
  • Français
  • Español
  • العربية
  • Português
  • Deutsch
  • 日本語
  • 中文
statsmodels-Skill›Modul 2 · Modelle für Zeitreihen›Thema 2.3·~55 Min.Fortgeschritten

VAR und Granger-Kausalität

Position in der SkillModul 2 von 3 · Thema 3 von 5
2.3 VAR und Granger-Kausalität
—Lab-Punkte— / 100XP gesamt—Serie—

Variablen, die sich gegenseitig vorhersagen

Eine Vektorautoregression, VAR, modelliert mehrere Reihen gemeinsam. Jede Variable wird auf die vergangenen Werte aller Variablen regressiert, für zwei Reihen xxx und yyy und eine Verzögerung also yt=c+A yt−1+uty_t=c+A\,y_{t-1}+u_tyt​=c+Ayt−1​+ut​ mit yt=(xt,yt)⊤y_t=(x_t,y_t)^\topyt​=(xt​,yt​)⊤ und einer 2×22\times22×2-Matrix AAA. Jede Gleichung wird mit der gewöhnlichen Kleinste-Quadrate-Methode angepasst, und in statsmodels erledigt VAR(data).fit(p) alle auf einmal. Das Modell ist stabil, wenn alle Eigenwerte von AAA im Einheitskreis liegen, Schocks also abklingen statt zu explodieren. Prognosen entstehen durch wiederholte Multiplikation. Die Ordnung ppp wählt man durch Vergleich von Informationskriterien wie AIC und BIC, die Anpassung belohnen und zusätzliche Koeffizienten bestrafen.

Granger-Kausalität stellt eine enge Frage: Helfen vergangene Werte von yyy, xxx vorherzusagen, wenn die vergangenen Werte von xxx selbst bekannt sind? Man passt die Gleichung für xxx zweimal an, einmal mit den Verzögerungen von yyy (uneingeschränkt, Residuenquadratsumme RSSu\mathrm{RSS}_uRSSu​) und einmal ohne sie (eingeschränkt, RSSr\mathrm{RSS}_rRSSr​), und vergleicht mit einer F-Statistik. Ein signifikantes Ergebnis besagt, dass yyy zeitlich vor xxx liegt. Es beweist nicht, dass yyy xxx verursacht: Eine dritte Variable oder vorausschauende Anleger können dasselbe Muster erzeugen.

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)
Ein Zwei-Variablen-VAR von HandSei 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​) und der Start (x,y)=(1,2)(x,y)=(1,2)(x,y)=(1,2). Ein Schritt voraus ergibt (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), zwei Schritte (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). Die Spur ist 1,11{,}11,1, die 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, die Eigenwerte lösen also λ2−1,1λ+0,28=0\lambda^2-1{,}1\lambda+0{,}28=0λ2−1,1λ+0,28=0 und sind 0,70{,}70,7 und 0,40{,}40,4: beide im Einheitskreis, das System ist stabil. Für den F-Test: Steigt die Residuenquadratsumme beim Weglassen der Verzögerungen von yyy von RSSu=10\mathrm{RSS}_u=10RSSu​=10 auf RSSr=12\mathrm{RSS}_r=12RSSr​=12, mit q=2q=2q=2 Restriktionen und n−k=40n-k=40n−k=40, so ist F=(12−10)/210/40=4,0F=\frac{(12-10)/2}{10/40}=4{,}0F=10/40(12−10)/2​=4,0, über dem 5-%-Wert 3,233{,}233,23.

Mit echten Monatsdaten, 258258258 Monate seit März 2005, nehmen wir die SPY-Rendite in Prozent und den Spread zehnjährige minus zweijährige Staatsanleihe. Das AIC wählt 333 Verzögerungen, das BIC 111, und ein VAR(2) ist stabil, wobei die nächste Wurzel bei 1,021{,}021,02 liegt, weil der Spread sehr persistent ist (seine eigenen Verzögerungen summieren sich auf 0,980{,}980,98). Der Granger-Test, ob der Spread bei der Vorhersage der Renditen hilft, ergibt F=0,22F=0{,}22F=0,22 (p=0,81p=0{,}81p=0,81), ob Renditen bei der Vorhersage des Spreads helfen, F=1,88F=1{,}88F=1,88 (p=0,15p=0{,}15p=0,15). Keine Richtung ist signifikant: Die Steigung der Zinskurve bringt einer monatlichen SPY-Prognose nichts.

Lab: Was bewertet wird

  • Ein Schritt voraus10 Pkt.
  • Zwei Schritte voraus15 Pkt.
  • Der größte Eigenwert15 Pkt.
  • Eine Granger-F-Statistik15 Pkt.
  • Koeffizienten zählen15 Pkt.
  • Einen Granger-Test lesen10 Pkt.
  • Vorhersage ist keine Verursachung20 Pkt.
0/ 7 beantwortet
Aufgaben
1Numerische Antwort10 Pkt.Nicht beantwortet

Mit 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​) und den aktuellen Werten (x,y)=(1,2)(x,y)=(1,2)(x,y)=(1,2): Wie groß ist die Ein-Schritt-Prognose von xxx, 0,5x+0,2y0{,}5x+0{,}2y0,5x+0,2y?

h1=h1 =h1=akzeptiert: ± 0.01
2Numerische Antwort15 Pkt.Nicht beantwortet

Die Ein-Schritt-Prognose ist (0,9;1,3)(0{,}9;1{,}3)(0,9;1,3). Wie groß ist die Zwei-Schritt-Prognose von yyy, 0,1⋅0,9+0,6⋅1,30{,}1\cdot0{,}9+0{,}6\cdot1{,}30,1⋅0,9+0,6⋅1,3? Zwei Dezimalstellen.

h2=h2 =h2=akzeptiert: ± 0.01
3Numerische Antwort15 Pkt.Nicht beantwortet

Die Eigenwerte lösen λ2−1,1λ+0,28=0\lambda^2-1{,}1\lambda+0{,}28=0λ2−1,1λ+0,28=0 mit Spur 1,11{,}11,1 und Determinante 0,280{,}280,28. Wie groß ist die größere Wurzel? Verwenden Sie 1,12−4⋅0,28=0,3\sqrt{1{,}1^2-4\cdot0{,}28}=0{,}31,12−4⋅0,28​=0,3 und geben Sie eine Dezimalstelle an.

l=l =l=akzeptiert: ± 0.01
4Numerische Antwort15 Pkt.Nicht beantwortet

Das Weglassen der Verzögerungen von yyy erhöht die Residuenquadratsumme von RSSu=10\mathrm{RSS}_u=10RSSu​=10 auf RSSr=12\mathrm{RSS}_r=12RSSr​=12. Mit q=2q=2q=2 Restriktionen und n−k=40n-k=40n−k=40: Wie groß ist 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=akzeptiert: ± 0.05
5Numerische Antwort15 Pkt.Nicht beantwortet

Ein VAR(2) hat 333 Variablen. Jede seiner 333 Gleichungen hat eine Konstante plus 222 Verzögerungen jeder der 333 Variablen. Wie viele Koeffizienten hat das ganze Modell?

n=n =n=akzeptiert: ± 0.1
6Auswahl10 Pkt.Nicht beantwortet

Der Test, ob der Renditespread die monatlichen SPY-Renditen im Granger-Sinn verursacht, ergibt F = 0,22 mit p = 0,81. Was sagt das?

7In Worten erklären20 Pkt.Nicht beantwortet

Angenommen, yyy verursacht xxx im Granger-Sinn mit sehr kleinem p-Wert. Erklären Sie in ein bis zwei Sätzen, warum das trotzdem nicht beweist, dass yyy xxx bewegt.

Mathe-ModusEingabe per Tastatur oder LaTeXAufgabe 1 von 7 · Strg ⏎ zum Prüfen

Melde dich an, um deinen Code auszuführen.

Alle Lektionen kannst du kostenlos lesen. Erstelle ein Konto oder melde dich an, um Code auszuführen, Tests zu starten und dieses Lab abzugeben.

  • IDE im Browser mit Bash-Terminal — nichts zu installieren
  • Sofortige Bewertung, XP und dein Platz auf der Skill
Konto erstellenAnmelden
ODER
Weiter mit GitHub