O. Concepto

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

Ejemplo de series de tiempo son:
1. Precio de acciones.
2. Nivel de inventario.
3. Rotación.
4. Ventas. 5. PIB (GDP).

1. Instalar paquetes y llamar librerías

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

2. Crear la serie de tiempo

Ejemplo:
Los siguientes datos de producción trimestral inicial en el primer trimestre del 2020.
Se busca pronosticar la producción de los siguientes 5 trimestres.

produccion <- c(50,53,55,57, 55,60)

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

3. Crear modelo ARIMA

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

modelo <- auto.arima(serie_de_tiempo, D=1)
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 pronóstico

pronostico <- forecast(modelo, level = c(95), h=5)
#h (horizonte de tiempo) es la cantidad de veces que vas a hacer el pronóstico
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)

LS0tCnRpdGxlOiAiU2VyaWVzIGRlIFRpZW1wbyIKYXV0aG9yOiAiQXJhY2VseSBRdWV6YWRhIEEwMTIzNTkzMSIKZGF0ZTogIjIwMjMtMTAtMDMiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogVFJVRQogICAgdG9jX2Zsb2F0OiBUcnVlCiAgICBjb2RlX2Rvd25sb2FkOiBUcnVlCiAgICB0aGVtZTogImJvb3RzdHJhcCIKICAgIGhpZ2hsaWdodDogImVzcHJlc3NvIgotLS0KCiAgIVtdKC9Vc2Vycy9hcmFxdWV6YWRhL0RvY3VtZW50cy9UZWMvN21vIFNlbWVzdHJlL0FuYWxpzIF0aWNhIHBhcmEgTmVnb2Npb3MvTW9kdWxvIDQvdGllbXBvLmdpZikKCgoKCgoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBvcmFuZ2U7Ij5PLiBDb25jZXB0bzwvc3Bhbj4KClVuYSAqKnNlcmllIGRlIHRpZW1wbyoqIGVzIHVuYSBjb2xlY2Npw7NuIGRlIG9ic2VydmFjaW9uZXMgc29icmUgdW4gZGV0ZXJtaW5hZG8gZmVuw7NtZW5vIGVmZWN0dWFkYXMgZW4gbW9tZW50b3MgZGUgdGllbXBvIHN1Y2VzaXZvcywgdXN1YWxtZW50ZSBlcXVpZXNwYWNpYWRvcy4gIAoKRWplbXBsbyBkZSBzZXJpZXMgZGUgdGllbXBvIHNvbjogIAoxLiBQcmVjaW8gZGUgYWNjaW9uZXMuICAgCjIuIE5pdmVsIGRlIGludmVudGFyaW8uICAgCjMuIFJvdGFjacOzbi4gICAKNC4gVmVudGFzLiAKNS4gUElCIChHRFApLiAKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+MS4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hczwvc3Bhbj4KYGBge3J9CiNpbnN0YWxsLnBhY2thZ2VzKCJmb3JlY2FzdCIpCmxpYnJhcnkoZm9yZWNhc3QpCmBgYAoKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+Mi4gQ3JlYXIgbGEgc2VyaWUgZGUgdGllbXBvPC9zcGFuPgoKKipFamVtcGxvOioqICAgCkxvcyBzaWd1aWVudGVzIGRhdG9zIGRlIHByb2R1Y2Npw7NuIHRyaW1lc3RyYWwgaW5pY2lhbCBlbiBlbCBwcmltZXIgdHJpbWVzdHJlIGRlbCAyMDIwLiAgClNlIGJ1c2NhIHByb25vc3RpY2FyIGxhIHByb2R1Y2Npw7NuIGRlIGxvcyBzaWd1aWVudGVzIDUgdHJpbWVzdHJlcy4gCgpgYGB7cn0KcHJvZHVjY2lvbiA8LSBjKDUwLDUzLDU1LDU3LCA1NSw2MCkKCnNlcmllX2RlX3RpZW1wbyA8LSB0cyhkYXRhPXByb2R1Y2Npb24sIHN0YXJ0ID0gMjAyMCwgZnJlcXVlbmN5ID0gNCkKc2VyaWVfZGVfdGllbXBvCmBgYAoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBvcmFuZ2U7Ij4zLiBDcmVhciBtb2RlbG8gQVJJTUE8L3NwYW4+CgoqKkFSSU1BKiogc2lnbmlmaWNhICpBdXRvUmVncmVzc2l2ZSBJbnRlZ3JhdGVkIE1vdmluZyBBdmVyYWdlKiBvICpNb2RlbG8gQXV0b3JyZWdyZXNpdm8gSW50ZWdyYWRvIGRlIFByb21lZGlvIE3Ds3ZpbCoKCmBgYHtyfQptb2RlbG8gPC0gYXV0by5hcmltYShzZXJpZV9kZV90aWVtcG8sIEQ9MSkKbW9kZWxvCnN1bW1hcnkobW9kZWxvKQpgYGAKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBvcmFuZ2U7Ij40LiBSZWFsaXphciBwcm9uw7NzdGljbzwvc3Bhbj4KCmBgYHtyfQpwcm9ub3N0aWNvIDwtIGZvcmVjYXN0KG1vZGVsbywgbGV2ZWwgPSBjKDk1KSwgaD01KQojaCAoaG9yaXpvbnRlIGRlIHRpZW1wbykgZXMgbGEgY2FudGlkYWQgZGUgdmVjZXMgcXVlIHZhcyBhIGhhY2VyIGVsIHByb27Ds3N0aWNvCnByb25vc3RpY28KcGxvdChwcm9ub3N0aWNvKQpgYGAKCgoKCg==