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
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 emxe useacf2()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