Series de Tiempo

Una serie de tiempo es una colección de observaciones sobre un determinado fenómeno efectuadas en momentos de tiempo sucesivos, usualmente equiespaciados.

Ejemplos de series de tiempo son:
- Precio de acciones.
- Niveles de inventario.
- Rotación de personal.
- Ventas.
- PIB (GDP).

1. Llamar librerías

library(forecast)
## Registered S3 method overwritten by 'quantmod':
##   method            from
##   as.zoo.data.frame zoo

2. Crear la salida de tiempo

# Los siguientes datos de producción trimestral incian en el primer trimestre de 2020.  
# Se busca pronosticar la producción de los siguientes 5 trimestres.  

producción <- c(50,53,55,57,55,60)

serie_de_tiempo <- ts(data = producción, start = c(2020,1), frequency = 4)
serie_de_tiempo
##      Qtr1 Qtr2 Qtr3 Qtr4
## 2020   50   53   55   57
## 2021   55   60

3. Crear el modelo ARIMA

ARIMA significa AutoRegressive Integrated Moving Average o Modelo Autoregresivo Integrado de Promedio Móvil.

modelo <- auto.arima(serie_de_tiempo, D = 1) # D=Temporalidad
modelo
## Series: serie_de_tiempo 
## ARIMA(0,0,0)(0,1,0)[4] with drift 
## 
## Coefficients:
##        drift
##       1.5000
## s.e.  0.1768
## 
## sigma^2 = 2.01:  log likelihood = -2.84
## AIC=9.68   AICc=-2.32   BIC=7.06
summary(modelo)
## Series: serie_de_tiempo 
## ARIMA(0,0,0)(0,1,0)[4] with drift 
## 
## Coefficients:
##        drift
##       1.5000
## s.e.  0.1768
## 
## sigma^2 = 2.01:  log likelihood = -2.84
## AIC=9.68   AICc=-2.32   BIC=7.06
## 
## Training set error measures:
##                      ME      RMSE       MAE        MPE      MAPE       MASE
## Training set 0.03333332 0.5787923 0.3666667 0.03685269 0.6429133 0.06111111
##                    ACF1
## Training set -0.5073047

4. Realizar el pronóstico

pronostico <- forecast(modelo, level = c(95), h = 5) # Level=Nivel de confianza h=Número de pronósticos
pronostico
##         Point Forecast    Lo 95    Hi 95
## 2021 Q3             61 58.22127 63.77873
## 2021 Q4             63 60.22127 65.77873
## 2022 Q1             61 58.22127 63.77873
## 2022 Q2             66 63.22127 68.77873
## 2022 Q3             67 63.07028 70.92972
plot(pronostico)

LS0tDQp0aXRsZTogIlNlcmllcyBkZSBUaWVtcG8iDQphdXRob3I6ICJEZWlzeSBMw7NwZXoiDQpkYXRlOiAiMjAyMy0xMC0wMyINCm91dHB1dDogDQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiBUUlVFDQogICAgdG9jX2Zsb2F0OiBUUlVFDQogICAgY29kZV9kb3dubG9hZDogVFJVRQ0KICAgIHRoZW1lOiAiZGFya2x5Ig0KICAgIGhpZ2hsaWdodDogImhhZGRvY2siDQplZGl0b3Jfb3B0aW9uczogDQogIGNodW5rX291dHB1dF90eXBlOiBjb25zb2xlDQotLS0NCiFbXShDOlxcVXNlcnNcXERlaXN5IExvcGV6XFxEb2N1bWVudHNcXDcuIFPDqXB0aW1vIFNlbWVzdHJlIEFnby1EaWMgMjAyM1xcUlN0dWRpb19BTlxcc2VyaWUuZ2lmKQ0KDQojIFNlcmllcyBkZSBUaWVtcG8NCg0KVW5hICoqc2VyaWUgZGUgdGllbXBvKiogZXMgdW5hIGNvbGVjY2nDs24gZGUgb2JzZXJ2YWNpb25lcyBzb2JyZSB1biBkZXRlcm1pbmFkbyBmZW7Ds21lbm8gZWZlY3R1YWRhcyBlbiBtb21lbnRvcyBkZSB0aWVtcG8gc3VjZXNpdm9zLCB1c3VhbG1lbnRlIGVxdWllc3BhY2lhZG9zLiAgDQoNCipFamVtcGxvcyBkZSBzZXJpZXMgZGUgdGllbXBvIHNvbjoqICANCi0gUHJlY2lvIGRlIGFjY2lvbmVzLiAgDQotIE5pdmVsZXMgZGUgaW52ZW50YXJpby4gIA0KLSBSb3RhY2nDs24gZGUgcGVyc29uYWwuICANCi0gVmVudGFzLiAgDQotIFBJQiAoR0RQKS4gIA0KDQojIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBncmF5OyI+MS4gTGxhbWFyIGxpYnJlcsOtYXM8L3NwYW4+DQpgYGB7cn0NCmxpYnJhcnkoZm9yZWNhc3QpDQpgYGANCg0KIyMjIDxzcGFuIHN0eWxlPSJjb2xvcjogZ3JheTsiPjIuIENyZWFyIGxhIHNhbGlkYSBkZSB0aWVtcG88L3NwYW4+DQpgYGB7cn0NCiMgTG9zIHNpZ3VpZW50ZXMgZGF0b3MgZGUgcHJvZHVjY2nDs24gdHJpbWVzdHJhbCBpbmNpYW4gZW4gZWwgcHJpbWVyIHRyaW1lc3RyZSBkZSAyMDIwLiAgDQojIFNlIGJ1c2NhIHByb25vc3RpY2FyIGxhIHByb2R1Y2Npw7NuIGRlIGxvcyBzaWd1aWVudGVzIDUgdHJpbWVzdHJlcy4gIA0KDQpwcm9kdWNjacOzbiA8LSBjKDUwLDUzLDU1LDU3LDU1LDYwKQ0KDQpzZXJpZV9kZV90aWVtcG8gPC0gdHMoZGF0YSA9IHByb2R1Y2Npw7NuLCBzdGFydCA9IGMoMjAyMCwxKSwgZnJlcXVlbmN5ID0gNCkNCnNlcmllX2RlX3RpZW1wbw0KYGBgDQoNCiMjIyA8c3BhbiBzdHlsZT0iY29sb3I6IGdyYXk7Ij4zLiBDcmVhciBlbCBtb2RlbG8gQVJJTUE8L3NwYW4+DQoqKkFSSU1BKiogc2lnbmlmaWNhICpBdXRvUmVncmVzc2l2ZSBJbnRlZ3JhdGVkIE1vdmluZyBBdmVyYWdlKiBvICpNb2RlbG8gQXV0b3JlZ3Jlc2l2byBJbnRlZ3JhZG8gZGUgUHJvbWVkaW8gTcOzdmlsKi4gIA0KDQpgYGB7cn0NCm1vZGVsbyA8LSBhdXRvLmFyaW1hKHNlcmllX2RlX3RpZW1wbywgRCA9IDEpICMgRD1UZW1wb3JhbGlkYWQNCm1vZGVsbw0Kc3VtbWFyeShtb2RlbG8pDQpgYGANCg0KIyMjIDxzcGFuIHN0eWxlPSJjb2xvcjogZ3JheTsiPjQuIFJlYWxpemFyIGVsIHByb27Ds3N0aWNvPC9zcGFuPg0KYGBge3J9DQpwcm9ub3N0aWNvIDwtIGZvcmVjYXN0KG1vZGVsbywgbGV2ZWwgPSBjKDk1KSwgaCA9IDUpICMgTGV2ZWw9Tml2ZWwgZGUgY29uZmlhbnphIGg9TsO6bWVybyBkZSBwcm9uw7NzdGljb3MNCnByb25vc3RpY28NCnBsb3QocHJvbm9zdGljbykNCmBgYA0KDQoNCg0K