Ajuste um modelo puramente sazonal
Assim como outros modelos, você pode ajustar modelos sazonais no R usando o comando sarima() do pacote astsa.
Para entender como funcionam os modelos puramente sazonais, é melhor considerar dados simulados. Geramos 250 observações de um modelo puramente sazonal dado por $$X_t = .9 X_{t-12} + W_t + .5 W_{t-12}\,,$$ que denotaríamos como um SARMA(P = 1, Q = 1)S = 12. Três anos de dados e a ACF e PACF do modelo estão plotados para você.
Você vai comparar os valores amostrais de ACF e PACF dos dados gerados com os valores verdadeiros exibidos.
O pacote astsa já está carregado e os dados gerados estão em x.
Este exercício faz parte do curso
Modelos ARIMA em R
Instruções do exercício
- Use
acf2()para plotar a ACF e a PACF amostrais dos dados gerados até a defasagem 60 e compare com os valores reais. Para estimar até a defasagem 60, defina o argumentomax.lagigual a60. - Ajuste o modelo aos dados gerados usando
sarima(). Além dos argumentosp,deqno seu comandosarima(), especifiqueP,D,QeS(note que R diferencia maiúsculas de minúsculas).
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Plot sample P/ACF to lag 60 and compare to the true values
acf2(___, max.lag = ___)
# Fit the seasonal model to x
sarima(x, p = 0, d = 0, q = 0, P = ___, D = 0, Q = ___, S = ___)