Series de tiempo

Concepto

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:
1. Precio de acciones.
2. Niveles de inventario.
3. Rotación de personal.
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 las series de tiempo

#Ejemplo: Los siguientes datos de producción bimestral inician en el primer trimestre del 2020
#Se busca promocionar la producción de los siguientes 5 trimestres.

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

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

3. Crear el modelo ARIMA

ARIMA significa AutoRegressive Integrated Moving Average o Modelo Autoregresivo Integrado de Promedio Movil

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 el pronóstico

pronostico <- forecast(modelo, level=c(95), h=5)
pronostico
##         Point Forecast    Lo 95    Hi 95
## 2022 Q1             61 58.22127 63.77873
## 2022 Q2             63 60.22127 65.77873
## 2022 Q3             61 58.22127 63.77873
## 2022 Q4             66 63.22127 68.77873
## 2023 Q1             67 63.07028 70.92972
plot(pronostico)

LS0tCnRpdGxlOiAiU2VyaWVzIGRlIFRpZW1wbyIKYXV0aG9yOiAiTWFyaWFuYSBHIgpkYXRlOiAiMjAyMy0xMC0wMyIKb3V0cHV0OiAKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiBUUlVFCiAgICB0b2NfZmxvYXQ6IFRSVUUKICAgIGNvZGVfZG93bmxvYWQ6IFRSVUUKICAgIHRoZW1lOiAiYm9vdHN0cmFwIgogICAgaGlnaGxpZ2h0OiAiZXNwcmVzc28iCi0tLQoKIVtdKC9Vc2Vycy9tYXJpYW5hZ3VldmFyYS9EZXNrdG9wL0Fkb2JlU3RvY2tfMzAwMTY3MTYyLXNjYWxlZC5qcGcpCgojIDxzcGFuIHN0eWxlPSJjb2xvcjogYmxhY2s7Ij5TZXJpZXMgZGUgdGllbXBvPC9zcGFuPgoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBibGFjazsiPkNvbmNlcHRvPC9zcGFuPgoKVW5hICoqc2VyaWUgZGUgdGllbXBvKiogZXMgdW5hIGNvbGVjY2nDs24gZGUgb2JzZXJ2YWNpb25lcyBzb2JyZSB1biBkZXRlcm1pbmFkbyBmZW7Ds21lbm8gZWZlY3R1YWRhcyBlbiBtb21lbnRvcyBkZSB0aWVtcG8gc3VjZXNpdm9zLCB1c3VhbG1lbnRlIGVxdWllc3BhY2lhZG9zLgoKRWplbXBsb3MgZGUgc2VyaWVzIGRlIHRpZW1wbzogIAoxLiBQcmVjaW8gZGUgYWNjaW9uZXMuICAKMi4gTml2ZWxlcyBkZSBpbnZlbnRhcmlvLiAgCjMuIFJvdGFjacOzbiBkZSBwZXJzb25hbC4gIAo0LiBWZW50YXMuICAKNS4gUElCIChHRFApLiAgCgoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBibGFjazsiPjEuIEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcsOtYXM8L3NwYW4+CmBgYHtyfQojaW5zdGFsbC5wYWNrYWdlcygiZm9yZWNhc3QiKQpsaWJyYXJ5KGZvcmVjYXN0KQpgYGAKCgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IGJsYWNrOyI+Mi4gQ3JlYXIgbGFzIHNlcmllcyBkZSB0aWVtcG88L3NwYW4+CmBgYHtyfQojRWplbXBsbzogTG9zIHNpZ3VpZW50ZXMgZGF0b3MgZGUgcHJvZHVjY2nDs24gYmltZXN0cmFsIGluaWNpYW4gZW4gZWwgcHJpbWVyIHRyaW1lc3RyZSBkZWwgMjAyMAojU2UgYnVzY2EgcHJvbW9jaW9uYXIgbGEgcHJvZHVjY2nDs24gZGUgbG9zIHNpZ3VpZW50ZXMgNSB0cmltZXN0cmVzLgoKcHJvZHVjY2lvbiA8LSBjKDUwLDUzLDU1LDU3LDU1LDYwKQoKc2VyaWVfZGVfdGllbXBvIDwtIHRzKGRhdGE9cHJvZHVjY2lvbiwgc3RhcnQgPSBjKDIwMjAsMyksIGZyZXF1ZW5jeSA9IDQpCnNlcmllX2RlX3RpZW1wbwpgYGAKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmxhY2s7Ij4zLiBDcmVhciBlbCBtb2RlbG8gQVJJTUE8L3NwYW4+CioqQVJJTUEqKiBzaWduaWZpY2EgKkF1dG9SZWdyZXNzaXZlIEludGVncmF0ZWQgTW92aW5nIEF2ZXJhZ2UqIG8gKk1vZGVsbyBBdXRvcmVncmVzaXZvIEludGVncmFkbyBkZSBQcm9tZWRpbyBNb3ZpbCoKCmBgYHtyfQptb2RlbG8gPC0gYXV0by5hcmltYShzZXJpZV9kZV90aWVtcG8sIEQ9MSkKbW9kZWxvCnN1bW1hcnkobW9kZWxvKQpgYGAKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmxhY2s7Ij40LiBSZWFsaXphciBlbCBwcm9uw7NzdGljbzwvc3Bhbj4KCmBgYHtyfQpwcm9ub3N0aWNvIDwtIGZvcmVjYXN0KG1vZGVsbywgbGV2ZWw9Yyg5NSksIGg9NSkKcHJvbm9zdGljbwpwbG90KHByb25vc3RpY28pCmBgYAoKCgoKCgoKCgoKCgo=