Exemplo de Previsão com o melhor modelo ARIMA para um modelo sazonal

Um analista coletou dados sobre o número de passageiros de companhias aéreas por 108 meses. O analista quer usar um modelo ARIMA para gerar previsões para os dados. O analista analisou anteriormente um gráfico de séries tempoariais dos dados e observou que a variação no ciclo sazonal aumenta ao longo do tempo. O analista concluiu que uma transformação natural dos dados é apropriada. Após a transformação, o analista analisou o gráfico da série temporais dos dados transformados e o gráfico da função de autocorrelação (ACF) dos dados transformados. Ambas as gráficos sugerem que o ponto de partida para o modelo é escolher 1 para a ordem de diferenciamento não sazonal e 1 para a ordem de diferenciação sazonal. O analista pede previsões para os próximos 3 meses.

  1. Abra os dados amostrais PassageirosAereos.mtw.
  2. Escolha Estat > Séries temporais > Previsão com o melhor modelo ARIMA.
  3. Em Séries, digite Número de passageiros.
  4. Em Ordem diferencial d, selecione 1.
  5. Selecione Ajustar modelos sazonais com período e insira 12 para o período.
  6. Em Ordem diferencial sazonal D, selecione 1.
  7. Em Número de previsões, insira 3.
  8. Selecione Opções.
  9. Em Transformação de Box-Cox, selecione λ = 0 (log natural).
  10. Clique em OK em cada caixa de diálogo.

Interprete os resultados

A tabela de seleção de modelos classifica os modelos da pesquisa em ordem pela AICc. O modelo ARIMA(0, 1, 1)(1, 1, 0) tem o menor AICc. Os resultados da ARIMA a seguir são para o modelo ARIMA (0, 1, 1)(1, 1, 0).

Os valores p na tabela de parâmetros mostram que os termos médios móveis são significativos no nível de 0,05. O analista conclui que os coeficientes pertencem ao modelo. Os valores p para as estatísticas modificadas de Box-Pierce (Ljung-Box) são todos insignificantes no nível de 0,05. O ACF dos resíduos e o PACF dos resíduos mostram um pico na defasagem 24. Como um grande pico em um número de atraso alto geralmente é um falso positivo e as estatísticas de teste são todas insignificantes, o analista conclui que o modelo atende à suposição de que os resíduos são independentes. O analista conclui que o exame das previsões é razoável.

* AVISO * Modelos ARIMA (p, d, q) (P, D, Q) que não incluem um termo constante:
(2; 1; 1)(1; 1; 1)

Método

Período sazonal12
Critério para o melhor modeloAICc mínimo
Transformação de Box-Cox 
    λ especificado pelo usuário0
    Série transformada = ln(Número de passageiros) 
Linhas usadas108
Linhas não usadas0

Seleção de modelo

Modelo (d = 1; D = 1)Log-verossimilhançaAICcAICBIC
p = 0; q = 1; P = 1; Q = 0*243,477-480,690-480,954-473,292
p = 2; q = 0; P = 0; Q = 1243,903-479,362-479,806-469,590
p = 1; q = 1; P = 1; Q = 0243,496-478,547-478,992-468,776
p = 0; q = 2; P = 1; Q = 0243,480-478,516-478,961-468,745
p = 2; q = 0; P = 1; Q = 1244,424-478,174-478,848-466,079
p = 0; q = 1; P = 0; Q = 0237,930-471,729-471,859-466,752
p = 1; q = 2; P = 0; Q = 0239,930-471,415-471,859-461,644
p = 1; q = 1; P = 0; Q = 0237,929-469,594-469,858-462,196
p = 0; q = 2; P = 0; Q = 0237,924-469,584-469,848-462,186
p = 1; q = 0; P = 0; Q = 1237,442-468,619-468,883-461,221
p = 1; q = 0; P = 1; Q = 1237,551-466,658-467,102-456,887
p = 2; q = 2; P = 0; Q = 0238,267-465,860-466,534-453,765
p = 2; q = 0; P = 0; Q = 0232,478-458,693-458,957-451,295
p = 0; q = 0; P = 0; Q = 1226,062-447,993-448,124-443,016
p = 0; q = 0; P = 1; Q = 1226,282-446,300-446,563-438,902
p = 2; q = 1; P = 0; Q = 0226,105-443,766-444,211-433,995
p = 1; q = 0; P = 0; Q = 0222,409-440,687-440,818-435,710
p = 2; q = 0; P = 1; Q = 0220,456-432,467-432,911-422,696
p = 0; q = 0; P = 1; Q = 0218,236-432,342-432,472-427,364
p = 1; q = 2; P = 1; Q = 1220,708-428,461-429,416-414,092
p = 0; q = 2; P = 0; Q = 1215,116-421,787-422,232-412,016
p = 0; q = 1; P = 0; Q = 1213,007-419,751-420,015-412,353
p = 2; q = 1; P = 0; Q = 1214,469-418,265-418,939-406,169
p = 1; q = 0; P = 1; Q = 0211,232-416,199-416,463-408,801
p = 2; q = 2; P = 0; Q = 1213,877-414,799-415,754-400,431
p = 2; q = 2; P = 1; Q = 1214,698-414,109-415,397-397,520
p = 1; q = 2; P = 0; Q = 1211,492-412,310-412,984-400,215
p = 1; q = 1; P = 0; Q = 1208,149-407,854-408,299-398,083
p = 0; q = 1; P = 1; Q = 1204,745-401,046-401,490-391,275
p = 0; q = 2; P = 1; Q = 1203,978-397,282-397,956-385,187
p = 1; q = 1; P = 1; Q = 1203,564-396,453-397,127-384,358
p = 1; q = 2; P = 1; Q = 0170,812-330,950-331,624-318,855
p = 2; q = 2; P = 1; Q = 0167,845-322,735-323,690-308,367
p = 2; q = 1; P = 1; Q = 0-202,538415,751415,076427,846
* Melhor modelo com AICc mínimo.  A saída para o melhor modelo está a seguir.

Estimativas Finais de Parâmetros

TipoCoef.EP de CoefValor-TValor-P
ARS  12-0,4030,103-3,920,000
MM   10,87040,051017,080,000
Diferenciação: regular@ 1, sazonal 1 de ordem 12
Número de observações após a diferenciação: 95

Sumário do Modelo

GLSQQMDPMAICcAICBIC
930,03113260,00033480,0003277-480,690-480,954-473,292
MS = variância da série de ruído branco

Estatística Qui-Quadrado de Box-Pierce (Ljung-Box) Modificada

Lag12243648
Qui-Quadrado9,4726,4433,9950,66
GL10223446
Valor-P0,4890,2330,4680,295
* AVISO * Modelos ARIMA (p, d, q) (P, D, Q) que não incluem um termo constante:
(2; 1; 1)(1; 1; 1)

Série original

Período
de tempo

Limites de 95%
PrevisãoInferiorSuperiorReal
109168226641622724217434097 
110208238762008075121587153 
111208267022007744321596450 

Série transformada

Período
de tempo


Limites de 95%
PrevisãoPrevisão SEInferiorSuperiorReal
10916,63810,018296416,602216,6739 
11016,85140,018449516,815316,8876 
11116,85160,018601416,815116,8880