Ajustando um modelo AR(1)
Lembre que você usa o par ACF e PACF para ajudar a identificar as ordens \(p\) e \(q\) de um modelo ARMA. A tabela a seguir resume os resultados:
| AR(\(p\)) | MA(\(q\)) | ARMA(\(p,q\)) | |
|---|---|---|---|
| ACF | Decai gradualmente | Corta após defasagem \(q\) |
Decai gradualmente |
| PACF | Corta após defasagem \(p\) |
Decai gradualmente | Decai gradualmente |
Neste exercício, você vai gerar dados do modelo AR(1), $$X_t = .9 X_{t-1} + W_t,$$ observar os dados simulados e o par ACF e PACF amostrais para determinar a ordem. Depois, você vai ajustar o modelo e comparar os parâmetros estimados com os parâmetros verdadeiros.
Ao longo deste curso, você usará sarima() do pacote astsa para ajustar modelos aos dados de forma simples. O comando gera um gráfico de diagnóstico dos resíduos que pode ser ignorado até discutirmos diagnósticos mais adiante neste capítulo.
Este exercício faz parte do curso
Modelos ARIMA em R
Instruções do exercício
- O pacote astsa já está pré-carregado.
- Use o comando
arima.sim()já escrito para gerar 100 observações de um modelo AR(1) com parâmetro AR igual a 0,9. Salve emx. - Plote os dados gerados usando
plot(). - Faça os gráficos do par ACF e PACF amostrais usando o comando
acf2()do pacoteastsa. - Use
sarima()doastsapara ajustar um AR(1)aos dados gerados anteriormente. Examine a t-table e compare as estimativas com os valores verdadeiros. Por exemplo, se a série temporal estiver emx, para ajustar um AR(1) aos dados, usesarima(x, p = 1, d = 0, q = 0)ou simplesmentesarima(x, 1, 0, 0)`.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Generate 100 observations from the AR(1) model
x <- arima.sim(model = list(order = c(1, 0, 0), ar = .9), n = 100)
# Plot the generated data
# Plot the sample P/ACF pair
# Fit an AR(1) to the data and examine the t-table