Stacjonarne szeregi czasowe: ARMA

Modele ARIMA w R

David Stoffer

Professor of Statistics at the University of Pittsburgh

Dekompozycja Wolda

wold.png Wold udowodnił, że każdy stacjonarny szereg czasowy można przedstawić jako liniową kombinację białego szumu:

$$X_t = W_t + a_1 W_{t-1} + a_2 W_{t-2} + ...$$

Dla stałych $ \ a_1,a_2,...$

Każdy model ARMA ma tę postać, co oznacza, że nadaje się do modelowania szeregów czasowych.

Uwaga: Szczególnym przypadkiem jest MA(q), który już ma tę postać, gdzie stałe po q-tym wyrazie wynoszą 0.

Modele ARIMA w R

Generowanie ARMA za pomocą arima.sim()

  • Podstawowa składnia:
arima.sim(model, n, ...)
  • model to lista z rzędem modelu jakoc(p, d, q) i współczynnikami
  • n to długość szeregu
Modele ARIMA w R

Generowanie i wykres MA(1)

ch1_3.013.png

Modele ARIMA w R

Generowanie i wykres MA(1)

ch1_3.014.png

x <- arima.sim(list(order = c(0, 0, 1), ma = 0.9), n = 100)
plot(x)
Modele ARIMA w R

Generowanie i wykres AR(2)

ch1_3.016.png

Modele ARIMA w R

Generowanie i wykres AR(2)

ch1_3.017.png

x <- arima.sim(list(order = c(2, 0, 0), ar = c(0, -0.9)), n = 100)
plot(x)
Modele ARIMA w R

Ćwiczmy!

Modele ARIMA w R

Preparing Video For Download...