ComeçarComece de graça

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

Ver curso

Instruções do exercício

  • Plote os dados gerados.
  • Use acf2() do astsa para 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

Editar e executar o código