Trend i heteroskedastyczność – jak sobie z nimi radzić
W tym ćwiczeniu przekształcimy niestacjonarne dane do postaci stacjonarnej, obliczając stopę zwrotu lub stopę wzrostu w następujący sposób.
Często szeregi czasowe są generowane zgodnie ze wzorem $$X_t = (1 + p_t) X_{t-1}$$, co oznacza, że wartość szeregu zaobserwowana w chwili \(t\) jest równa wartości z chwili \(t-1\) powiększonej o niewielką procentową zmianę \(p_t\) w chwili \(t\).
Prostym, deterministycznym przykładem jest trzymanie pieniędzy na koncie bankowym ze stałym oprocentowaniem \(p\). W takim przypadku \(X_t\) oznacza stan konta w okresie \(t\) przy początkowym depozycie \(X_0\).
Wielkość \(p_t\) jest zazwyczaj określana jako stopa zwrotu lub stopa wzrostu szeregu czasowego – i ten proces jest często stabilny.
Z powodów wykraczających poza zakres tego kursu można pokazać, że stopę wzrostu \(p_t\) da się przybliżyć wzorem $$Y_t = \log X_t - \log X_{t-1} \approx p_t.$$
W R wartość \(p_t\) oblicza się najczęściej jako diff(log(x)), a jej wykres można wygenerować w jednej linii: plot(diff(log(x))).
To ćwiczenie jest częścią kursu
Modele ARIMA w R
Instrukcje do ćwiczenia
- Podobnie jak wcześniej, pakiety astsa i xts są już wczytane.
- Utwórz wykres wielopanelowy, aby: (1) wykreślić kwartalne dane PKB USA (
gnp) i zauważyć, że nie są stacjonarne, oraz (2) wykreślić przybliżoną stopę wzrostu PKB USA przy użyciudiff()ilog(). - Użyj wykresu wielopanelowego, aby: (1) wykreślić dzienne ceny zamknięcia indeksu DJIA (
djia$Close) i zauważyć, że nie są stacjonarne – dane są obiektemxts. Następnie (2) wykreśl przybliżone stopy zwrotu indeksu DJIA przy użyciudiff()ilog(). Jak wyniki wypadają w porównaniu ze stopą wzrostu PKB?
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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)