
O. Concepto
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("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 inicial en el primer
trimestre del 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 = 2020, frequency = 4)
serie_de_tiempo
## Qtr1 Qtr2 Qtr3 Qtr4
## 2020 50 53 55 57
## 2021 55 60
3. Crear modelo ARIMA
ARIMA significa AutoRegressive 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 pronóstico
pronostico <- forecast(modelo, level = c(95), h=5)
#h (horizonte de tiempo) es la cantidad de veces que vas a hacer el pronóstico
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

LS0tCnRpdGxlOiAiU2VyaWVzIGRlIFRpZW1wbyIKYXV0aG9yOiAiQXJhY2VseSBRdWV6YWRhIEEwMTIzNTkzMSIKZGF0ZTogIjIwMjMtMTAtMDMiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogVFJVRQogICAgdG9jX2Zsb2F0OiBUcnVlCiAgICBjb2RlX2Rvd25sb2FkOiBUcnVlCiAgICB0aGVtZTogImJvb3RzdHJhcCIKICAgIGhpZ2hsaWdodDogImVzcHJlc3NvIgotLS0KCiAgIVtdKC9Vc2Vycy9hcmFxdWV6YWRhL0RvY3VtZW50cy9UZWMvN21vIFNlbWVzdHJlL0FuYWxpzIF0aWNhIHBhcmEgTmVnb2Npb3MvTW9kdWxvIDQvdGllbXBvLmdpZikKCgoKCgoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBvcmFuZ2U7Ij5PLiBDb25jZXB0bzwvc3Bhbj4KClVuYSAqKnNlcmllIGRlIHRpZW1wbyoqIGVzIHVuYSBjb2xlY2Npw7NuIGRlIG9ic2VydmFjaW9uZXMgc29icmUgdW4gZGV0ZXJtaW5hZG8gZmVuw7NtZW5vIGVmZWN0dWFkYXMgZW4gbW9tZW50b3MgZGUgdGllbXBvIHN1Y2VzaXZvcywgdXN1YWxtZW50ZSBlcXVpZXNwYWNpYWRvcy4gIAoKRWplbXBsbyBkZSBzZXJpZXMgZGUgdGllbXBvIHNvbjogIAoxLiBQcmVjaW8gZGUgYWNjaW9uZXMuICAgCjIuIE5pdmVsIGRlIGludmVudGFyaW8uICAgCjMuIFJvdGFjacOzbi4gICAKNC4gVmVudGFzLiAgIAo1LiBQSUIgKEdEUCkuIAoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBvcmFuZ2U7Ij4xLiBJbnN0YWxhciBwYXF1ZXRlcyB5IGxsYW1hciBsaWJyZXLDrWFzPC9zcGFuPgpgYGB7cn0KI2luc3RhbGwucGFja2FnZXMoImZvcmVjYXN0IikKbGlicmFyeShmb3JlY2FzdCkKYGBgCgoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBvcmFuZ2U7Ij4yLiBDcmVhciBsYSBzZXJpZSBkZSB0aWVtcG88L3NwYW4+CgoqKkVqZW1wbG86KiogICAKTG9zIHNpZ3VpZW50ZXMgZGF0b3MgZGUgcHJvZHVjY2nDs24gdHJpbWVzdHJhbCBpbmljaWFsIGVuIGVsIHByaW1lciB0cmltZXN0cmUgZGVsIDIwMjAuICAKU2UgYnVzY2EgcHJvbm9zdGljYXIgbGEgcHJvZHVjY2nDs24gZGUgbG9zIHNpZ3VpZW50ZXMgNSB0cmltZXN0cmVzLiAKCmBgYHtyfQpwcm9kdWNjaW9uIDwtIGMoNTAsNTMsNTUsNTcsIDU1LDYwKQoKc2VyaWVfZGVfdGllbXBvIDwtIHRzKGRhdGE9cHJvZHVjY2lvbiwgc3RhcnQgPSAyMDIwLCBmcmVxdWVuY3kgPSA0KQpzZXJpZV9kZV90aWVtcG8KYGBgCgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjMuIENyZWFyIG1vZGVsbyBBUklNQTwvc3Bhbj4KCioqQVJJTUEqKiBzaWduaWZpY2EgKkF1dG9SZWdyZXNzaXZlIEludGVncmF0ZWQgTW92aW5nIEF2ZXJhZ2UqIG8gKk1vZGVsbyBBdXRvcnJlZ3Jlc2l2byBJbnRlZ3JhZG8gZGUgUHJvbWVkaW8gTcOzdmlsKgoKYGBge3J9Cm1vZGVsbyA8LSBhdXRvLmFyaW1hKHNlcmllX2RlX3RpZW1wbywgRD0xKQptb2RlbG8Kc3VtbWFyeShtb2RlbG8pCmBgYAojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjQuIFJlYWxpemFyIHByb27Ds3N0aWNvPC9zcGFuPgoKYGBge3J9CnByb25vc3RpY28gPC0gZm9yZWNhc3QobW9kZWxvLCBsZXZlbCA9IGMoOTUpLCBoPTUpCiNoIChob3Jpem9udGUgZGUgdGllbXBvKSBlcyBsYSBjYW50aWRhZCBkZSB2ZWNlcyBxdWUgdmFzIGEgaGFjZXIgZWwgcHJvbsOzc3RpY28KcHJvbm9zdGljbwpwbG90KHByb25vc3RpY28pCmBgYAoKCgoK