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

LS0tDQp0aXRsZTogIlNlcmllcyBkZSBUaWVtcG8iDQphdXRob3I6ICJYaW1lbmEgQ2FzdGlsbG8gQTAxMzY5OTQ5Ig0KZGF0ZTogIjIwMjMtMTAtMDMiDQpvdXRwdXQ6DQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiBUUlVFDQogICAgdG9jX2Zsb2F0OiBUUlVFDQogICAgY29kZV9kb3dubG9hZDogVFJVRQ0KICAgIHRoZW1lOiAiYm9vdHN0cmFwIg0KICAgIGhpZ2hsaWdodDogImVzcHJlc3NvIg0KLS0tDQoNCiFbXShDOlxcVXNlcnNcXHhpbWVuXFxEb3dubG9hZHNcXHNlcmllc190aWVtcG8uZ2lmKQ0KDQo8Y2VudGVyPg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjogT3JhbmdlOyI+U2VyaWVzIGRlIFRpZW1wbzwvc3Bhbj4NCg0KDQpVbmEgKnNlcmllIGRlIHRpZW1wbyogZXMgdW5hIGNvbGVjY2nDs24gZGUgb2JzZXJ2YWNpb25lcyBzb2JyZSB1biBkZXRlcm1pbmFkbyBmZW7Ds21lbm8gZWZlY3R1YWRhcyBlbiBtb21lbnRvcyBkZSB0aWVtcG8gc3VjZXNpdm9zLCB1c3VhbG1lbnRlIGVxdWllc3BhY2lhZG9zLg0KDQogIEVqZW1wbG8gZGUgc2VyaWVzIGRlIHRpZW1wbyBzb246DQogIDEuIFByZWNpbyBkZSBhY2Npb25lcw0KICAyLiBOaXZlbCBkZSBpbnZlbnRhcmlvDQogIDMuIFJvdGFjacOzbg0KICA0LiBWZW50YXMuICAgDQogIDUuIFBJQiAoR0RQKS4gICANCg0KICANCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+MS4gSW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hczwvc3Bhbj4NCmBgYHtyfQ0KI2luc3RhbGwucGFja2FnZXMoKQ0KbGlicmFyeShmb3JlY2FzdCkNCmBgYA0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjIuIENyZWFyIGxhIFNlcmllIGRlIFRpZW1wbzwvc3Bhbj4NCmBgYHtyfQ0KI0VqZW1wbG86IExvcyBzaWd1aWVudGVzIGRhdG9zIGRlIHByb2R1Y2Npw7NuIHRyaW1lc3RyYWwNCiNpbmljaWFuIGVuIGVsIHByaW1lciB0cmltZXN0cmUgZGUgMjAyMC4NCiNTZSBidXNjYSBwcm9ub3N0aWNhciBsYSBwcm9kdWNjacOzbiBkZSBsb3Mgc2lndWllbnRlIDUgdHJpbWVzdHJlcy4NCg0KcHJvZHVjY2lvbiA8LSBjKDUwLDUzLDU1LDU3LDU1LDYwKQ0KDQpzZXJpZV9kZV90aWVtcG8gPC0gdHMoZGF0YT1wcm9kdWNjaW9uLCBzdGFydCA9IGMoMjAyMCwxKSwgZnJlcXVlbmN5ID0gNCkNCmBgYA0KDQojIyA8c3BhbiBzdHlsZT0iY29sb3I6IG9yYW5nZTsiPjMuIENyZWFyIGVsIG1vZGVsbyBBUklNQTwvc3Bhbj4NCioqQVJJTUEqKiBzaWduaWZpY2EgKkF1dG9SZWdncmVzaXZlIEludGVncmF0ZWQgTW92aW5nIEF2ZXJhZ2UqIG8gKk1vZGVsbyANCkF1dG9SZWdyZXNpdm8gSW50ZWdyYWRvIGRlIFByb21lZGlvIE3Ds3ZpbCouDQoNCmBgYHtyfQ0KbW9kZWxvIDwtIGF1dG8uYXJpbWEoc2VyaWVfZGVfdGllbXBvLCBEPTEpDQptb2RlbG8NCnN1bW1hcnkobW9kZWxvKQ0KYGBgDQoNCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogb3JhbmdlOyI+NC4gUmVhbGl6YXIgZWwgcHJvbsOzc3RpY288L3NwYW4+DQpgYGB7cn0NCnByb25vc3RpY28gPC0gZm9yZWNhc3QobW9kZWxvLCBsZXZlbD1jKDk1KSwgaD01KQ0KcHJvbm9zdGljbw0KcGxvdChwcm9ub3N0aWNvKQ0KYGBgDQoNCg==