Předpovědi s modely ARIMA
Automatická metoda z předchozího cvičení zvolila pro data austa model ARIMA(0,1,1) s driftem, tedy:
\(y_t = c + y_{t-1} + \theta e_{t-1} + e_t.\)
Teď si vyzkoušíš různé další modely ARIMA na stejných datech a uvidíš, jak ovlivňují výsledné předpovědi.
Funkce Arima() umožňuje zvolit konkrétní model ARIMA. Její první argument order se nastaví na vektor určující hodnoty \(p\), \(d\) a \(q\). Druhý argument include.constant je logická hodnota, která určuje, zda má být do modelu zahrnuta konstanta \(c\), neboli drift. Níže je příklad řetězce funkcí, který vykreslí předpovědi pro usnetelec z modelu ARIMA(2,1,2) s driftem:
> usnetelec %>%
Arima(order = c(2,1,2), include.constant = TRUE) %>%
forecast() %>%
autoplot()
Sleduj, jak různé modely ovlivňují předpovědi i intervaly spolehlivosti. Data austa máš připravená v pracovním prostředí.
Toto cvičení je součástí kurzu
Forecasting v R
Pokyny k cvičení
- Vykresli předpovědi z modelu ARIMA(0,1,1) bez driftu.
- Vykresli předpovědi z modelu ARIMA(2,1,3) s driftem.
- Vykresli předpovědi z modelu ARIMA(0,0,1) s konstantou.
- Vykresli předpovědi z modelu ARIMA(0,2,1) bez konstanty.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# Plot forecasts from an ARIMA(0,1,1) model with no drift
austa %>% Arima(order = c(___, ___, ___), include.constant = ___) %>% ___ %>% ___
# Plot forecasts from an ARIMA(2,1,3) model with drift
austa %>% Arima(___, ___) %>% ___ %>% ___
# Plot forecasts from an ARIMA(0,0,1) model with a constant
austa %>% Arima(___, ___) %>% ___ %>% ___
# Plot forecasts from an ARIMA(0,2,1) model with no constant
austa %>% Arima(___, ___) %>% ___ %>% ___