Stationära tidsserier: ARMA

ARIMA-modeller i R

David Stoffer

Professor of Statistics at the University of Pittsburgh

Wolds dekomposition

wold.png Wold visade att varje stationär tidsserie kan representeras som en linjärkombination av vitt brus:

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

För konstanterna $ \ a_1,a_2,...$

Varje ARMA-modell har denna form, vilket innebär att de lämpar sig för modellering av tidsserier.

Notera: Specialfallet MA(q) har redan denna form, där konstanterna är 0 efter den q:te termen.

ARIMA-modeller i R

Generera ARMA med arima.sim()

  • Grundläggande syntax:
arima.sim(model, n, ...)
  • model är en lista med modellens ordning somc(p, d, q) och koefficienterna
  • n är seriens längd
ARIMA-modeller i R

Generera och plotta MA(1)

ch1_3.013.png

ARIMA-modeller i R

Generera och plotta MA(1)

ch1_3.014.png

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

Generera och plotta AR(2)

ch1_3.016.png

ARIMA-modeller i R

Generera och plotta AR(2)

ch1_3.017.png

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

Nu kör vi en övning!

ARIMA-modeller i R

Preparing Video For Download...