LoslegenKostenlos loslegen

Umgang mit Trend und Heteroskedastizität

Hier zwingen wir nichtstationäre Daten zur Stationarität, indem wir die Rendite bzw. Wachstumsrate wie folgt berechnen.

Häufig werden Zeitreihen durch $$X_t = (1 + p_t) X_{t-1}$$ erzeugt. Das bedeutet: Der zum Zeitpunkt \(t\) beobachtete Wert der Zeitreihe entspricht dem Wert zum Zeitpunkt \(t-1\) plus einer kleinen prozentualen Änderung \(p_t\) zum Zeitpunkt \(t\).

Ein einfaches deterministisches Beispiel ist das Einzahlen von Geld auf ein Bankkonto mit einem festen Zinssatz \(p\). In diesem Fall ist \(X_t\) der Kontostand im Zeitraum \(t\) bei einer anfänglichen Einzahlung von \(X_0\).

Typischerweise wird \(p_t\) als Rendite oder Wachstumsrate einer Zeitreihe bezeichnet, und dieser Prozess ist oft stabil.

Aus Gründen, die außerhalb des Umfangs dieses Kurses liegen, lässt sich zeigen, dass die Wachstumsrate \(p_t\) näherungsweise durch $$Y_t = \log X_t - \log X_{t-1} \approx p_t$$ beschrieben werden kann.

In R wird \(p_t\) häufig mit diff(log(x)) berechnet, und die Visualisierung gelingt in einer Zeile mit plot(diff(log(x))).

Diese Übung ist Teil des Kurses

ARIMA-Modelle in R

Kurs anzeigen

Anleitung zur Übung

  • Wie zuvor sind die Pakete astsa und xts vorab geladen.
  • Erzeuge ein Mehrfach-Plot, um (1) die vierteljährlichen US-GNP-Daten (gnp) zu plotten und festzustellen, dass sie nicht stationär sind, und (2) die approximative Wachstumsrate des US-GNP mit diff() und log() zu plotten.
  • Verwende ein Mehrfach-Plot, um (1) die täglichen DJIA-Schlusskurse (djia$Close) zu plotten und festzustellen, dass sie nicht stationär sind. Die Daten sind ein xts-Objekt. Anschließend (2) plotte die approximativen DJIA-Renditen mit diff() und log(). Wie verhält sich das im Vergleich zur Wachstumsrate des GNP?

Interaktive Übung

Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.

# astsa and xts are preloaded 

# Plot GNP series (gnp) and its growth rate
par(mfrow = c(2,1))
plot(gnp)


# Plot DJIA closings (djia$Close) and its returns
par(mfrow = c(2,1))
plot(djia$Close)


 
Code bearbeiten und ausführen