Modelos ARIMA automáticos para series temporales no estacionales
En el vídeo, viste que la función auto.arima() selecciona un modelo autoregressive integrated moving average (ARIMA) adecuado dada una serie temporal, del mismo modo que ets() lo hace para los modelos ETS. La función summary() puede aportar información adicional:
> # p = 2, d = 1, p = 2
> summary(fit)
Series: usnetelec
ARIMA(2,1,2) with drift
...
En este ejercicio, elegirás automáticamente un modelo ARIMA para la serie austa (ya cargada), que contiene el número anual de visitantes internacionales a Australia de 1980 a 2015. Después comprobarás los residuos (recuerda que un p-valor mayor que 0.05 indica que los datos se parecen a ruido blanco) y generarás algunos pronósticos. Aparte de la función de modelado, esto es idéntico a lo que hiciste con el pronóstico ETS.
Este ejercicio forma parte del curso
Forecasting en R
Instrucciones del ejercicio
- Ajusta un modelo ARIMA automático a la serie
austausando la función recién presentada. Guarda el resultado enfit. - Usa la función adecuada para comprobar que los residuos del modelo resultante se comportan como ruido blanco. Asigna
TRUE(si los residuos parecen ruido blanco) oFALSE(si no) aresidualsok. - Aplica
summary()al modelo para ver los coeficientes ajustados. - Según los resultados de
summary(), ¿cuál es el valor de AICc con dos decimales? ¿Cuántas diferenciaciones se usaron? Asigna estos valores aAICcyd, respectivamente. - Por último, usando el operador pipe, traza los pronósticos de los próximos 10 periodos a partir del modelo elegido.
ejercicio interactivo práctico
Prueba este ejercicio completando este código de ejemplo.
# Fit an automatic ARIMA model to the austa series
fit <- ___
# Check that the residuals look like white noise
___
residualsok <- ___
# Summarize the model
___
# Find the AICc value and the number of differences used
AICc <- ___
d <- ___
# Plot forecasts of fit
fit %>% forecast(h = ___) %>% ___()