ARIMA - plug and play
Como você viu no vídeo, uma série temporal é chamada de ARIMA(\(p,d,q\)) se a série diferenciada (de ordem \(d\)) for ARMA(\(p,q\)).
Para entender como o modelo funciona, você vai analisar dados simulados do modelo integrado $$ Y_t = .9 Y_{t-1} + W_t\, $$ em que \(Y_t = \nabla X_t = X_t - X_{t-1}\). Nesse caso, o modelo é um ARIMA(1,1,0) porque os dados diferenciados seguem uma autoregressão de ordem 1.
A série temporal simulada está em x e foi gerada no R como
x <- arima.sim(model = list(order = c(1, 1, 0), ar = .9), n = 200).
Você vai traçar os dados gerados e a ACF e PACF amostrais dos dados para ver como dados integrados se comportam. Em seguida, vai diferenciar os dados para torná-los estacionários. Você vai traçar os dados diferenciados e as respectivas ACF e PACF amostrais para ver como a diferenciação faz diferença.
Como antes, o pacote astsa já foi carregado para você. Dados de um ARIMA(1,1,0) com parâmetro AR igual a .9 estão salvos no objeto x.
Este exercício faz parte do curso
Modelos ARIMA em R
Instruções do exercício
- Plote os dados gerados.
- Use
acf2()doastsapara traçar o par P/ACF amostral dos dados gerados. - Plote os dados diferenciados.
- Use outra chamada a
acf2()para ver o par P/ACF amostral dos dados diferenciados. Observe como eles sugerem um modelo AR(1) para os dados diferenciados.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Plot x
# Plot the P/ACF pair of x
# Plot the differenced data
# Plot the P/ACF pair of the differenced data