Ajustarea unui model ARMA
Acum ești pregătit(ă) să combini modelul AR cu modelul MA într-un model ARMA. Am generat date din modelul ARMA(2,1), $$X_t = X_{t-1} - .9 X_{t-2} + W_t + .8 W_{t-1}, $$ x <- arima.sim(model = list(order = c(2, 0, 1), ar = c(1, -.9), ma = .8), n = 250). Analizează datele simulate și perechea de ACF și PACF eșantion pentru a determina un model posibil.
Amintește-ți că, în cazul modelelor ARMA(\(p, q\)), atât ACF teoretic, cât și PACF teoretic descresc treptat. În această situație, ordinele sunt greu de identificat din date și s-ar putea să nu fie clar dacă ACF eșantion sau PACF eșantion se taie brusc sau descrește treptat. Deoarece cunoști ordinele reale ale modelului, ajustează un ARMA(2,1) pe datele generate. Strategiile generale de modelare vor fi discutate în continuare în cadrul cursului.
Acest exercițiu face parte din cursul
Modele ARIMA în R
Instrucțiuni pentru exercițiu
- Pachetul astsa este preîncărcat. 250 de observații ARMA(2,1) se află în
x. - Ca în exercițiile anterioare, folosește
plot()pentru a reprezenta grafic datele generate dinxși foloseșteacf2()pentru a vizualiza perechile de ACF și PACF eșantion. - Folosește
sarima()pentru a ajusta un model ARMA(2,1) pe datele generate. Examinează tabelul t și compară estimările cu valorile reale.
Exercițiu interactiv practic
Încearcă acest exercițiu completând acest cod de exemplu.
# astsa is preloaded
# Plot x
# Plot the sample P/ACF of x
# Fit an ARMA(2,1) to the data and examine the t-table