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
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 mitdiff()undlog()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 einxts-Objekt. Anschließend (2) plotte die approximativen DJIA-Renditen mitdiff()undlog(). 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)