
0. 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 son:
Precio de Acciones
Niveles de inventario
Rotación de personal.
Ventas.
PIB (GDP)
1. Cargar la librería
#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.
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 “Auroregressive 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
## 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

LS0tCnRpdGxlOiAiU2VyaWVzIGRlIFRpZW1wbyIKYXV0aG9yOiAiQW5kcmVhIFPDoW5jaGV6IgpkYXRlOiAiMjAyMy0xMC0wMyIKb3V0cHV0OiAKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiBUUlVFCiAgICB0b2NfZmxvYXQ6IFRSVUUKICAgIGNvZGVfZG93bmxvYWQ6IFRSVUUKICAgIHRoZW1lOiAicmVhZGFibGUiCiAgICBoaWdobGlnaHQ6ICJweWdtZW50cyIgCi0tLQoKIVtdKC9Vc2Vycy9hbmRyZWFzYW5jaGV6Z29tZXovRG9jdW1lbnRzL1Ivc2VyaWVzIGRlIHRpZW1wby5naWYpCgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IGdvbGRlbnJvZDsiPiAwLiBDb25jZXB0byA8L3NwYW4+CgpVbmEgKipzZXJpZSBkZSB0aWVtcG8qKiBlcyB1bmEgY29sZWNjacOzbiBkZSBvYnNlcnZhY2lvbmVzIHNvYnJlIHVuIGRldGVybWluYWRvIGZlbsOzbWVubyBlZmVjdHVhZGFzIGVuIG1vbWVudG9zIGRlIHRpZW1wbyBzdWNlc2l2b3MsIHVzdWFsbWVudGUgZXF1aWVzcGFjaWFkb3MuIAoKRWplbXBsb3MgZGUgc2VyaWVzIGRlIHRpZW1wbyBzb246CgoxLiBQcmVjaW8gZGUgQWNjaW9uZXMgCgoyLiBOaXZlbGVzIGRlIGludmVudGFyaW8KCjMuIFJvdGFjacOzbiBkZSBwZXJzb25hbC4gCgo0LiBWZW50YXMuIAoKNS4gUElCIChHRFApCgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IGdvbGRlbnJvZDsiPiAxLiBDYXJnYXIgbGEgbGlicmVyw61hIDwvc3Bhbj4KYGBge3J9CiNpbnN0YWxsLnBhY2thZ2VzKCJmb3JlY2FzdCIpCmxpYnJhcnkoZm9yZWNhc3QpCmBgYAoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBnb2xkZW5yb2Q7Ij4gMi4gQ3JlYXIgbGEgc2VyaWUgZGUgdGllbXBvIDwvc3Bhbj4KYGBge3J9CiMgRWplbXBsbzogTG9zIHNpZ3VpZW50ZXMgZGF0b3MgZGUgcHJvZHVjY2nDs24gdHJpbWVzdHJhbCBpbmljaWFuIGVuIGVsIHByaW1lciB0cmltZXN0cmUgZGUgMjAyMC4gCiMgU2UgYnVzY2EgcHJvbm9zdGljYXIgbGEgcHJvZHVjY2nDs24gZGUgbG9zIHNpZ3VpZW50ZXMgNSB0cmltZXN0cmVzLiAKCnByb2R1Y2Npb24gPC0gYyg1MCwgNTMsIDU1LCA1NywgNTUsIDYwKQpzZXJpZV9kZV90aWVtcG8gPC0gdHMoZGF0YSA9IHByb2R1Y2Npb24sIHN0YXJ0ID0gYygyMDIwLCAzKSwgZnJlcXVlbmN5ID0gNCkKc2VyaWVfZGVfdGllbXBvCmBgYAoKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogZ29sZGVucm9kOyI+IDMuIENyZWFyIGVsIG1vZGVsbyBBUklNQSA8L3NwYW4+CgoqKiJBUklNQSLCqCoqIHNpZ25pZmljYSAiQXVyb3JlZ3Jlc3NpdmUgSW50ZWdyYXRlZCBNb3ZpbmcgQXZlcmFnZSIgbyAiTW9kZWxvIEF1dG9ycmVncmVzaXZvIEludGVncmFkbyBkZSBQcm9tZWRpbyBNw7N2aWwiCgpgYGB7cn0KbW9kZWxvIDwtIGF1dG8uYXJpbWEoc2VyaWVfZGVfdGllbXBvLCBEPTEpCm1vZGVsbwpzdW1tYXJ5KG1vZGVsbykKYGBgCgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IGdvbGRlbnJvZDsiPiA0LiBSZWFsaXphciBlbCBwcm9uw7NzdGljbyA8L3NwYW4+CmBgYHtyfQpwcm9ub3N0aWNvIDwtIGZvcmVjYXN0KG1vZGVsbywgbGV2ZWwgPSBjKDk1KSwgaD01KQpwcm9ub3N0aWNvCnBsb3QocHJvbm9zdGljbykKYGBgCgo=