ComeçarComece de graça

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

Ver curso

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 em x.
  • Plote os dados gerados usando plot().
  • Faça os gráficos do par ACF e PACF amostrais usando o comando acf2() do pacote astsa.
  • Use sarima() do astsa para 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

Editar e executar o código