Model ARIMA musiman

Peramalan di R

Rob J. Hyndman

Professor of Statistics at Monash University

Model ARIMA

ch4_vid3_slides.003.png

d = Jumlah beda lag-1
p = Jumlah lag AR biasa:
q = Jumlah lag MA biasa:

Model ARIMA

ch4_vid3_slides.007.png

d = Jumlah beda lag-1
p = Jumlah lag AR biasa:
q = Jumlah lag MA biasa:

Model ARIMA

ch4_vid3_slides.008.png

d = Jumlah beda lag-1
p = Jumlah lag AR biasa: $\ y_{t-1}, y_{t-2},...,y_{t-p}$
q = Jumlah lag MA biasa: $\ \epsilon_{t-1}, \epsilon_{t-2},...,\epsilon_{t-q}$
D = Jumlah beda musiman
P = Jumlah lag AR musiman: $\ y_{t-m}, y_{t-2m},...,y_{t-Pm}$
Q = Jumlah lag MA musiman:$\ \epsilon_{t-m}, \epsilon_{t-2m},...,\epsilon_{t-Qm}$
m = Jumlah observasi per tahun

Contoh: Penggunaan kartu debit ritel bulanan di Islandia

autoplot(debitcards) +
  xlab("Year") + ylab("million ISK") +
  ggtitle("Retail debit card usage in Iceland")

Contoh: Penggunaan kartu debit ritel bulanan di Islandia

fit <- auto.arima(debitcards, lambda = 0)
fit

Series: debitcards
ARIMA(0,1,4)(0,1,1)[12]
Box Cox transformation: lambda= 0
Coefficients:
         ma1    ma2    ma3     ma4    sma1
      -0.796  0.086  0.263  -0.175  -0.814
s.e.   0.082  0.099  0.100   0.080   0.112
sigma^2 estimated as 0.00232:  log likelihood=239.3
AIC=-466.7   AICc=-466.1   BIC=-448.6

Contoh: Penggunaan kartu debit ritel bulanan di Islandia

fit %>%
  forecast(h = 36) %>%
  autoplot() + xlab("Year")

Ayo berlatih!

Peramalan di R