
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
## 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

LS0tCnRpdGxlOiAiU2VyaWVzIGRlIFRpZW1wbyIKYXV0aG9yOiAiTWFyaWFuYSBHIgpkYXRlOiAiMjAyMy0xMC0wMyIKb3V0cHV0OiAKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiBUUlVFCiAgICB0b2NfZmxvYXQ6IFRSVUUKICAgIGNvZGVfZG93bmxvYWQ6IFRSVUUKICAgIHRoZW1lOiAiYm9vdHN0cmFwIgogICAgaGlnaGxpZ2h0OiAiZXNwcmVzc28iCi0tLQoKIVtdKC9Vc2Vycy9tYXJpYW5hZ3VldmFyYS9EZXNrdG9wL0Fkb2JlU3RvY2tfMzAwMTY3MTYyLXNjYWxlZC5qcGcpCgojIDxzcGFuIHN0eWxlPSJjb2xvcjogYmxhY2s7Ij5TZXJpZXMgZGUgdGllbXBvPC9zcGFuPgoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBibGFjazsiPkNvbmNlcHRvPC9zcGFuPgoKVW5hICoqc2VyaWUgZGUgdGllbXBvKiogZXMgdW5hIGNvbGVjY2nDs24gZGUgb2JzZXJ2YWNpb25lcyBzb2JyZSB1biBkZXRlcm1pbmFkbyBmZW7Ds21lbm8gZWZlY3R1YWRhcyBlbiBtb21lbnRvcyBkZSB0aWVtcG8gc3VjZXNpdm9zLCB1c3VhbG1lbnRlIGVxdWllc3BhY2lhZG9zLgoKRWplbXBsb3MgZGUgc2VyaWVzIGRlIHRpZW1wbzogIAoxLiBQcmVjaW8gZGUgYWNjaW9uZXMuICAKMi4gTml2ZWxlcyBkZSBpbnZlbnRhcmlvLiAgCjMuIFJvdGFjacOzbiBkZSBwZXJzb25hbC4gIAo0LiBWZW50YXMuICAKNS4gUElCIChHRFApLiAgCgoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBibGFjazsiPjEuIEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcsOtYXM8L3NwYW4+CmBgYHtyfQojaW5zdGFsbC5wYWNrYWdlcygiZm9yZWNhc3QiKQpsaWJyYXJ5KGZvcmVjYXN0KQpgYGAKCgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IGJsYWNrOyI+Mi4gQ3JlYXIgbGFzIHNlcmllcyBkZSB0aWVtcG88L3NwYW4+CmBgYHtyfQojRWplbXBsbzogTG9zIHNpZ3VpZW50ZXMgZGF0b3MgZGUgcHJvZHVjY2nDs24gYmltZXN0cmFsIGluaWNpYW4gZW4gZWwgcHJpbWVyIHRyaW1lc3RyZSBkZWwgMjAyMAojU2UgYnVzY2EgcHJvbW9jaW9uYXIgbGEgcHJvZHVjY2nDs24gZGUgbG9zIHNpZ3VpZW50ZXMgNSB0cmltZXN0cmVzLgoKcHJvZHVjY2lvbiA8LSBjKDUwLDUzLDU1LDU3LDU1LDYwKQoKc2VyaWVfZGVfdGllbXBvIDwtIHRzKGRhdGE9cHJvZHVjY2lvbiwgc3RhcnQgPSBjKDIwMjAsMyksIGZyZXF1ZW5jeSA9IDQpCnNlcmllX2RlX3RpZW1wbwpgYGAKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmxhY2s7Ij4zLiBDcmVhciBlbCBtb2RlbG8gQVJJTUE8L3NwYW4+CioqQVJJTUEqKiBzaWduaWZpY2EgKkF1dG9SZWdyZXNzaXZlIEludGVncmF0ZWQgTW92aW5nIEF2ZXJhZ2UqIG8gKk1vZGVsbyBBdXRvcmVncmVzaXZvIEludGVncmFkbyBkZSBQcm9tZWRpbyBNb3ZpbCoKCmBgYHtyfQptb2RlbG8gPC0gYXV0by5hcmltYShzZXJpZV9kZV90aWVtcG8sIEQ9MSkKbW9kZWxvCnN1bW1hcnkobW9kZWxvKQpgYGAKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogYmxhY2s7Ij40LiBSZWFsaXphciBlbCBwcm9uw7NzdGljbzwvc3Bhbj4KCmBgYHtyfQpwcm9ub3N0aWNvIDwtIGZvcmVjYXN0KG1vZGVsbywgbGV2ZWw9Yyg5NSksIGg9NSkKcHJvbm9zdGljbwpwbG90KHByb25vc3RpY28pCmBgYAoKCgoKCgoKCgoKCgo=