Series de Tiempo

Concepto

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

Ejemplos de series de tiempo son:
1. Precio de acciones.
2. Niveles de inventario.
3. Rotación de personal.
4. Ventas.
5. PIB.

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 inician 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,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 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. Realiar el pronóstico

pronóstico <- forecast(modelo, level=c(95), h=5)
pronóstico
##         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(pronóstico)

LS0tDQp0aXRsZTogIlNlcmllcyBkZSB0aWVtcG8iDQphdXRob3I6ICJMZWlzbHkgQ3J1eiINCmRhdGU6ICIyMDIzLTEwLTAzIg0Kb3V0cHV0Og0KICBodG1sX2RvY3VtZW50Og0KICAgIHRvYzogVFJVRQ0KICAgIHRvY19mbG9hdDogVFJVRQ0KICAgIGNvZGVfZG93bmxvYWQ6IFRSVUUNCiAgICB0aGVtZTogImJvb3RzdHJhcCINCiAgICBoaWdobGlnaHQ6ICJlc3ByZXNzbyINCi0tLQ0KDQohW10oQzpcXFVzZXJzXFxrYXJpbVxcT25lRHJpdmVcXFBpY3R1cmVzXFxzZXJpZXN0aWVtcG8uanBnKQ0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+U2VyaWVzIGRlIFRpZW1wbzwvc3Bhbj4NCg0KIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBvcmFuZ2U7Ij5Db25jZXB0bzwvc3Bhbj4NClVuYSAqc2VyaWUgZGUgdGllbXBvKiBlcyB1bmEgY29sZWNjacOzbiBkZSBvYnNlcnZhY2lvbmVzIHNvYnJlIHVuIGRldGVybWluYWRvIGZlbsOzbWVubyBlZmVjdHVhZGFzIGVuIG1vbWVudG9zIGRlIHRpZW1wb3Mgc3VjZXNpdm9zLCB1c3VhbG1lbnRlIGVxdWllc3BhY2lhZG9zLiAgDQoNCkVqZW1wbG9zIGRlIHNlcmllcyBkZSB0aWVtcG8gc29uOiAgDQoxLiBQcmVjaW8gZGUgYWNjaW9uZXMuICAgDQoyLiBOaXZlbGVzIGRlIGludmVudGFyaW8uICANCjMuIFJvdGFjacOzbiBkZSBwZXJzb25hbC4gIA0KNC4gVmVudGFzLiAgDQo1LiBQSUIuIA0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjEuIEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcsOtYXM8L3NwYW4+DQoNCmBgYHtyfQ0KI2luc3RhbGwucGFja2FnZXMoImZvcmVjYXN0IikNCmxpYnJhcnkoZm9yZWNhc3QpDQpgYGANCg0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjIuIENyZWFyIGxhIHNlcmllIGRlIHRpZW1wbzwvc3Bhbj4NCmBgYHtyfQ0KIyBFamVtcGxvOiBMb3Mgc2lndWllbnRlcyBkYXRvcyBkZSBwcm9kdWNjacOzbiB0cmltZXN0cmFsIGluaWNpYW4gZW4gZWwgcHJpbWVyIHRyaW1lc3RyZSBkZSAyMDIwLg0KIyBTZSBidXNjYSBwcm9ub3N0aWNhciBsYSBwcm9kdWNjacOzbiBkZSBsb3Mgc2lndWllbnRlcyA1IHRyaW1lc3RyZXMuDQoNCnByb2R1Y2Npw7NuIDwtIGMoNTAsNTMsNTUsNTcsNTUsNjApDQpzZXJpZV9kZV90aWVtcG8gPC0gdHMoZGF0YSA9IHByb2R1Y2Npw7NuLCBzdGFydCA9IGMoMjAyMCwzKSwgZnJlcXVlbmN5ID0gNCkNCnNlcmllX2RlX3RpZW1wbw0KYGBgDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+My4gQ3JlYXIgZWwgbW9kZWxvIEFSSU1BPC9zcGFuPg0KDQoqKkFSSU1BKiogc2lnbmlmaWNhICpBdXRvUmVncmVzc2l2ZSBJbnRlZ3JhdGVkIE1vdmluZyBBdmVyYWdlKiBvICpNb2RlbG8gQXV0b3JlZ3Jlc2l2byBJbnRlZ3JhZG8gZGUgUHJvbWVkaW8gTcOzdmlsKi4gIA0KDQoNCmBgYHtyfQ0KbW9kZWxvIDwtIGF1dG8uYXJpbWEoc2VyaWVfZGVfdGllbXBvLCBEPTEpDQptb2RlbG8NCnN1bW1hcnkobW9kZWxvKQ0KYGBgDQoNCg0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjQuIFJlYWxpYXIgZWwgcHJvbsOzc3RpY288L3NwYW4+DQpgYGB7cn0NCnByb27Ds3N0aWNvIDwtIGZvcmVjYXN0KG1vZGVsbywgbGV2ZWw9Yyg5NSksIGg9NSkNCnByb27Ds3N0aWNvDQpwbG90KHByb27Ds3N0aWNvKQ0KYGBgDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQo=