Stationaire tijdreeksen: ARMA

ARIMA-modellen in R

David Stoffer

Professor of Statistics at the University of Pittsburgh

Wold-decompositie

wold.png Wold bewees dat elke stationaire tijdreeks kan worden weergegeven als een lineaire combinatie van witte ruis:

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

Voor constanten $ \ a_1,a_2,...$

Elk ARMA-model heeft deze vorm, dus ze zijn geschikt voor het modelleren van tijdreeksen.

Opmerking: het bijzondere geval MA(q) heeft al deze vorm, waarbij constanten 0 zijn na de q-de term.

ARIMA-modellen in R

ARMA genereren met arima.sim()

  • Basis-syntax:
arima.sim(model, n, ...)
  • model is een lijst met de orde van het model als c(p, d, q) en de coëfficiënten
  • n is de lengte van de reeks
ARIMA-modellen in R

MA(1) genereren en plotten

ch1_3.013.png

ARIMA-modellen in R

MA(1) genereren en plotten

ch1_3.014.png

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

AR(2) genereren en plotten

ch1_3.016.png

ARIMA-modellen in R

AR(2) genereren en plotten

ch1_3.017.png

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

Laten we oefenen!

ARIMA-modellen in R

Preparing Video For Download...