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.

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()
library(forecast)
## Warning: package 'forecast' was built under R version 4.3.1
## 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
#inician en el primer trimestre de 2020.
#Se busca pronosticar la producción de los siguiente 5 trimestres.

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

serie_de_tiempo <- ts(data=produccion, start = c(2020,1), frequency = 4)

3. Crear el modelo ARIMA

ARIMA significa AutoReggresive Integrated Moving Average o Modelo AutoRegresivo 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 el pronóstico

pronostico <- forecast(modelo, level=c(95), h=5)
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)

LS0tDQp0aXRsZTogIlNlcmllcyBkZSBUaWVtcG8iDQphdXRob3I6ICJYaW1lbmEgQ2FzdGlsbG8gQTAxMzY5OTQ5Ig0KZGF0ZTogIjIwMjMtMTAtMDMiDQpvdXRwdXQ6DQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiBUUlVFDQogICAgdG9jX2Zsb2F0OiBUUlVFDQogICAgY29kZV9kb3dubG9hZDogVFJVRQ0KICAgIHRoZW1lOiAiYm9vdHN0cmFwIg0KICAgIGhpZ2hsaWdodDogImVzcHJlc3NvIg0KLS0tDQoNCiFbXShDOlxcVXNlcnNcXHhpbWVuXFxEb3dubG9hZHNcXHNlcmllc190aWVtcG8uZ2lmKQ0KDQo8Y2VudGVyPg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjogT3JhbmdlOyI+U2VyaWVzIGRlIFRpZW1wbzwvc3Bhbj4NCg0KDQpVbmEgKnNlcmllIGRlIHRpZW1wbyogZXMgdW5hIGNvbGVjY2nDs24gZGUgb2JzZXJ2YWNpb25lcyBzb2JyZSB1biBkZXRlcm1pbmFkbyBmZW7Ds21lbm8gZWZlY3R1YWRhcyBlbiBtb21lbnRvcyBkZSB0aWVtcG8gc3VjZXNpdm9zLCB1c3VhbG1lbnRlIGVxdWllc3BhY2lhZG9zLg0KDQogIEVqZW1wbG8gZGUgc2VyaWVzIGRlIHRpZW1wbyBzb246DQogIDEuIFByZWNpbyBkZSBhY2Npb25lcw0KICAyLiBOaXZlbCBkZSBpbnZlbnRhcmlvDQogIDMuIFJvdGFjacOzbg0KICA0LiBWZW50YXMuICAgDQogIDUuIFBJQiAoR0RQKS4gICANCg0KICANCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+MS4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hczwvc3Bhbj4NCmBgYHtyfQ0KI2luc3RhbGwucGFja2FnZXMoKQ0KbGlicmFyeShmb3JlY2FzdCkNCmBgYA0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjIuIENyZWFyIGxhIFNlcmllIGRlIFRpZW1wbzwvc3Bhbj4NCmBgYHtyfQ0KI0VqZW1wbG86IExvcyBzaWd1aWVudGVzIGRhdG9zIGRlIHByb2R1Y2Npw7NuIHRyaW1lc3RyYWwNCiNpbmljaWFuIGVuIGVsIHByaW1lciB0cmltZXN0cmUgZGUgMjAyMC4NCiNTZSBidXNjYSBwcm9ub3N0aWNhciBsYSBwcm9kdWNjacOzbiBkZSBsb3Mgc2lndWllbnRlIDUgdHJpbWVzdHJlcy4NCg0KcHJvZHVjY2lvbiA8LSBjKDUwLDUzLDU1LDU3LDU1LDYwKQ0KDQpzZXJpZV9kZV90aWVtcG8gPC0gdHMoZGF0YT1wcm9kdWNjaW9uLCBzdGFydCA9IGMoMjAyMCwxKSwgZnJlcXVlbmN5ID0gNCkNCmBgYA0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjMuIENyZWFyIGVsIG1vZGVsbyBBUklNQTwvc3Bhbj4NCioqQVJJTUEqKiBzaWduaWZpY2EgKkF1dG9SZWdncmVzaXZlIEludGVncmF0ZWQgTW92aW5nIEF2ZXJhZ2UqIG8gKk1vZGVsbyANCkF1dG9SZWdyZXNpdm8gSW50ZWdyYWRvIGRlIFByb21lZGlvIE3Ds3ZpbCouDQoNCmBgYHtyfQ0KbW9kZWxvIDwtIGF1dG8uYXJpbWEoc2VyaWVfZGVfdGllbXBvLCBEPTEpDQptb2RlbG8NCnN1bW1hcnkobW9kZWxvKQ0KYGBgDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+NC4gUmVhbGl6YXIgZWwgcHJvbsOzc3RpY288L3NwYW4+DQpgYGB7cn0NCnByb25vc3RpY28gPC0gZm9yZWNhc3QobW9kZWxvLCBsZXZlbD1jKDk1KSwgaD01KQ0KcHJvbm9zdGljbw0KcGxvdChwcm9ub3N0aWNvKQ0KYGBgDQoNCg==