σμ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
  • 日本語
  • 中文
Particle filtering and non-stationary time series-Skill›Modul 2 · Partikelfilter›Thema 2.5·~1 Std.Fortgeschritten

Die Parameter lernen

Position in der SkillModul 2 von 3 · Thema 5 von 5
2.5 Die Parameter lernen
—Lab-Punkte— / 100XP gesamt—Serie—

Die Zahlen des Modells selbst lernen

Bisher waren die Zahlen des Modells, etwa qqq, rrr, μ\muμ, φ\varphiφ und σ\sigmaσ, gegeben. In der Praxis muss man sie aus denselben Daten lernen. Ein Partikelfilter hilft auf zwei Wegen. Erstens schätzt er die Likelihood der Daten für eine Wahl θ\thetaθ der Zahlen: An jedem Tag schätzt der Durchschnitt der Partikelgewichte, wie wahrscheinlich die Rendite des Tages war, und die Summe der Logarithmen dieser Durchschnitte ist die Log-Likelihood. Man lässt den Filter für mehrere Wahlen von θ\thetaθ auf einem Gitter laufen, stets mit denselben Zufallszahlen, damit der Vergleich glatt ist, und behält die beste. Die Schätzung ist verrauscht, mit einer Streuung von etwa 333 bei einer Log-Likelihood nahe 17 75017\,75017750, nur größere Abstände bedeuten etwas.

log⁡L^(θ)=∑tlog⁡(1N∑iwt(i)),h2=1−(3δ−12δ)2,a=1−h2,θ′=a θ+(1−a) θˉ+hV z\log\hat L(\theta)=\sum_t\log\left(\frac1N\sum_i w_t^{(i)}\right),\qquad h^{2}=1-\left(\frac{3\delta-1}{2\delta}\right)^{2},\quad a=\sqrt{1-h^{2}},\quad \theta'=a\,\theta+(1-a)\,\bar\theta+h\sqrt V\,zlogL^(θ)=t∑​log(N1​i∑​wt(i)​),h2=1−(2δ3δ−1​)2,a=1−h2​,θ′=aθ+(1−a)θˉ+hV​z(1)

Der zweite Weg legt die unbekannte Zahl in den Zustand: Jedes Partikel trägt seinen eigenen Wert von θ\thetaθ. Das scheitert, wenn θ\thetaθ nie bewegt wird, denn das Neuziehen kopiert dann wenige Werte, und die θ\thetaθ-Wolke schrumpft auf einen Punkt. Die Liu-West-Lösung bewegt das θ\thetaθ jedes Partikels nach jedem Neuziehen ein wenig, aber so, dass Mittelwert und Streuung der Wolke erhalten bleiben. Man wählt einen Diskontfaktor δ\deltaδ nahe 111 (etwa 0,950{,}950,95). Dann ist h2=1−((3δ−1)/(2δ))2h^{2}=1-\big((3\delta-1)/(2\delta)\big)^{2}h2=1−((3δ−1)/(2δ))2 und a=1−h2a=\sqrt{1-h^{2}}a=1−h2​. Jeder Wert wird zuerst mit dem Faktor aaa zum Wolkenmittel θˉ\bar\thetaθˉ gezogen und dann mit Rauschen der Varianz h2Vh^{2}Vh2V gestört, wobei VVV die Wolkenvarianz ist. Ein kleineres δ\deltaδ bedeutet mehr Störung.

Zahlen für Liu-West und den SPYMit δ=0,95\delta=0{,}95δ=0,95: (3δ−1)/(2δ)=0,9737(3\delta-1)/(2\delta)=0{,}9737(3δ−1)/(2δ)=0,9737, also h2=1−0,97372=0,0519h^{2}=1-0{,}9737^{2}=0{,}0519h2=1−0,97372=0,0519 und a=0,9737a=0{,}9737a=0,9737. Ein Partikel mit θ=0,90\theta=0{,}90θ=0,90 in einer Wolke mit Mittel 0,970{,}970,97 wird auf 0,9737×0,90+0,0263×0,97=0,9020{,}9737\times0{,}90+0{,}0263\times0{,}97=0{,}9020,9737×0,90+0,0263×0,97=0,902 gezogen. Ist die Wolkenvarianz V=0,0004V=0{,}0004V=0,0004, hat die Störung die Varianz 0,0519×0,00040{,}0519\times0{,}00040,0519×0,0004 und die Standardabweichung 0,00460{,}00460,0046. Beim SPY erreicht die Log-Likelihood der stochastischen Volatilität auf einem Gitter ihr Maximum bei φ=0,98\varphi=0{,}98φ=0,98, σ=0,25\sigma=0{,}25σ=0,25 mit 17 79317\,79317793, gegenüber 16 25216\,25216252 bei konstanter Volatilität. Ihr BIC mit 333 Zahlen, −2×17 793,1+3ln⁡5404=−35 560-2\times17\,793{,}1+3\ln5404=-35\,560−2×17793,1+3ln5404=−35560, schlägt das Hidden-Markov-Modell mit 3 Zuständen (−35 485-35\,485−35485) und sogar das mit 4 (−35 557-35\,557−35557).
  • Likelihood durch Filtern: die Summe der Logarithmen der mittleren Gewichte, mit denselben Zufallszahlen für jede Wahl.
  • Parameter im Zustand: braucht nach dem Neuziehen eine kleine Störung, sonst kollabiert die Parameterwolke.
  • Modelle fair vergleichen: Eine Log-Likelihood braucht eine Strafe für die Parameterzahl, etwa das BIC.

Lab: Was bewertet wird

  • Der Log-Likelihood-Gewinn10 Pkt.
  • Die Liu-West-Konstante h zum Quadrat15 Pkt.
  • Der Schrumpfungsfaktor a10 Pkt.
  • Ein geschrumpftes Partikel15 Pkt.
  • Die Größe der Störung10 Pkt.
  • Ein BIC15 Pkt.
  • Parameter ohne Störung10 Pkt.
  • Warum man mit einem Abstand von etwa 3 vergleicht15 Pkt.
0/ 8 beantwortet
Aufgaben
1Numerische Antwort10 Pkt.Nicht beantwortet

Beim SPY beträgt die beste Log-Likelihood der stochastischen Volatilität 17 793,117\,793{,}117793,1, konstante Volatilität ergibt 16 251,816\,251{,}816251,8. Wie groß ist der Gewinn, auf die nächste ganze Zahl?

g=g =g=akzeptiert: ± 2
2Numerische Antwort15 Pkt.Nicht beantwortet

Berechnen Sie mit δ=0,95\delta=0{,}95δ=0,95 den Wert h2=1−((3δ−1)/(2δ))2h^{2}=1-\big((3\delta-1)/(2\delta)\big)^{2}h2=1−((3δ−1)/(2δ))2. Geben Sie vier Dezimalstellen an.

h2=h2 =h2=akzeptiert: ± 0.0003
3Numerische Antwort10 Pkt.Nicht beantwortet

Mit h2=0,0519h^{2}=0{,}0519h2=0,0519: Wie groß ist a=1−h2a=\sqrt{1-h^{2}}a=1−h2​? Geben Sie vier Dezimalstellen an.

a=a =a=akzeptiert: ± 0.0003
4Numerische Antwort15 Pkt.Nicht beantwortet

Ein Partikel hat θ=0,90\theta=0{,}90θ=0,90, das Wolkenmittel ist θˉ=0,97\bar\theta=0{,}97θˉ=0,97 und a=0,9737a=0{,}9737a=0,9737. Wie groß ist der geschrumpfte Wert a θ+(1−a) θˉa\,\theta+(1-a)\,\bar\thetaaθ+(1−a)θˉ vor jeder Störung? Geben Sie drei Dezimalstellen an.

t=t =t=akzeptiert: ± 0.001
5Numerische Antwort10 Pkt.Nicht beantwortet

Die Störung hat die Varianz h2V=0,0519×0,0004=0,0000208h^{2}V=0{,}0519\times0{,}0004=0{,}0000208h2V=0,0519×0,0004=0,0000208. Wie groß ist ihre Standardabweichung? Geben Sie vier Dezimalstellen an.

s=s =s=akzeptiert: ± 0.0002
6Numerische Antwort15 Pkt.Nicht beantwortet

Das Modell der stochastischen Volatilität hat k=3k=3k=3 Zahlen, eine Log-Likelihood von 17 793,117\,793{,}117793,1 und n=5404n=5404n=5404 Renditen. Berechnen Sie −2log⁡L+kln⁡n-2\log L+k\ln n−2logL+klnn mit ln⁡5404=8,595\ln5404=8{,}595ln5404=8,595. Geben Sie eine Dezimalstelle an.

bic=bic =bic=akzeptiert: ± 1
7Auswahl10 Pkt.Nicht beantwortet

Jedes Partikel trägt seinen eigenen Parameterwert, aber die Werte werden nach dem Neuziehen nie bewegt. Was geschieht im Lauf der Zeit?

8In Worten erklären15 Pkt.Nicht beantwortet

Zwei Einstellungen liefern Partikelfilter-Log-Likelihoods von 17 793,117\,793{,}117793,1 und 17 789,317\,789{,}317789,3 für φ=0,98\varphi=0{,}98φ=0,98 und φ=0,99\varphi=0{,}99φ=0,99. Erklären Sie in ein bis zwei Sätzen, warum man nicht schließen sollte, dass 0,980{,}980,98 wirklich besser ist.

Mathe-ModusEingabe per Tastatur oder LaTeXAufgabe 1 von 8 · 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