ComeçarComece de graça

Ajustando um modelo ARMA

Agora você está pronto para unir o modelo AR e o modelo MA no modelo ARMA. Geramos dados do modelo 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). Observe os dados simulados e o par ACF e PACF amostrais para identificar um possível modelo.

Lembre-se de que, para modelos ARMA(\(p, q\)), tanto a ACF teórica quanto a PACF teórica apresentam cauda (tail off). Nesse caso, as ordens são difíceis de discernir a partir dos dados e pode não estar claro se a ACF ou a PACF amostrais estão se truncando ou apresentando cauda. Aqui, você conhece as ordens reais do modelo, então ajuste um ARMA(2,1) aos dados gerados. Estratégias gerais de modelagem serão discutidas mais adiante no curso.

Este exercício faz parte do curso

Modelos ARIMA em R

Ver curso

Instruções do exercício

  • O pacote astsa já está carregado. Há 250 observações ARMA(2,1) em x.
  • Como nos exercícios anteriores, use plot() para visualizar os dados gerados em x e use acf2() para ver os pares ACF e PACF amostrais.
  • Use sarima() para ajustar um ARMA(2,1) aos dados gerados. Examine a t-table e compare as estimativas com os valores verdadeiros.

Exercício interativo prático

Experimente este exercício completando este código de exemplo.

# 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

Editar e executar o código