Series de Tiempo

Concepto

Una Serie de tiempo es una collecion de observaciones sobre un determinado fenomeno efectuadas en momentos de tiempo sucesivos, usualmente equiespaciados.

Ejemplos de series de tiempo: 1. Precio de acciones 2. Niveles de inventario 3. Rotacion de personal 4. Ventas 5. PIB (GDP)

1. Instalar paquetes y llamar librerias

#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 produccion trimestral 
# Inician en el primer trimestre de 2020
# Se busca prinosticar la produccion de los siguientes
# 5 trimestral

produccion <-c(50,53,55,57,55,60)

serie_de_tiempo <- ts(data = produccion, start = c(2020, 1), frequency = 4)
serie_de_tiempo
##      Qtr1 Qtr2 Qtr3 Qtr4
## 2020   50   53   55   57
## 2021   55   60

3. Crear el modelo ARIMA

ARIMA significa AutoRegressive Integrated Moving Average o Modelo Autorregresivo 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
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. 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
plot(pronostico)

LS0tCnRpdGxlOiAiU2VyaWVzIGRlIFRpZW1wbyIKYXV0aG9yOiAiTWVsaXNzYSBHYW5lbSBLaXJjaG5lciIKZGF0ZTogIjIwMjMtMTAtMDMiCm91dHB1dDogCiAgaHRtbF9kb2N1bWVudDoKICAgIHRvYzogVFJVRQogICAgdG9jX2Zsb2F0OiBUUlVFCiAgICBjb2RlX2Rvd25sb2FkOiBUUlVFCiAgICB0aGVtZTogImJvb3RzdHJhcCIKICAgIGhpZ2hsaWdodDogImVzcHJlc3NvIgoKLS0tCgohW10oL1VzZXJzL21lbGlnYW5lbS9Eb2N1bWVudHMvVEVDL1NlbWVzdHJlIDcvQW5hbGl0aWNhIGRlIGRhdG9zL3RpZW1wby5naWYpCgojIDxzcGFuIHN0eWxlPSJjb2xvcjogcHVycGxlOyIgPlNlcmllcyBkZSBUaWVtcG88L3NwYW4+CgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IHB1cnBsZTsiPkNvbmNlcHRvPC9zcGFuPgpVbmEgKlNlcmllIGRlIHRpZW1wbyogZXMgdW5hIGNvbGxlY2lvbiBkZSBvYnNlcnZhY2lvbmVzIHNvYnJlIHVuIApkZXRlcm1pbmFkbyBmZW5vbWVubyBlZmVjdHVhZGFzIGVuIG1vbWVudG9zIGRlIHRpZW1wbyBzdWNlc2l2b3MsIAp1c3VhbG1lbnRlIGVxdWllc3BhY2lhZG9zLgoKRWplbXBsb3MgZGUgc2VyaWVzIGRlIHRpZW1wbzoKMS4gUHJlY2lvIGRlIGFjY2lvbmVzCjIuIE5pdmVsZXMgZGUgaW52ZW50YXJpbwozLiBSb3RhY2lvbiBkZSBwZXJzb25hbAo0LiBWZW50YXMKNS4gUElCIChHRFApCgojIyA8c3BhbiBzdHlsZT0iY29sb3I6IHB1cnBsZTsiPjEuIEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcmlhczwvc3Bhbj4KYGBge3J9CiNpbnN0YWxsLnBhY2thZ2VzKCJmb3JlY2FzdCIpCmxpYnJhcnkoZm9yZWNhc3QpCmBgYAojIyA8c3BhbiBzdHlsZT0iY29sb3I6IHB1cnBsZTsiPjIuIENyZWFyIGxhIHNlcmllIGRlIHRpZW1wbzwvc3Bhbj4KCmBgYHtyfQojIEVqZW1wbG86IExvcyBzaWd1aWVudGVzIGRhdG9zIGRlIHByb2R1Y2Npb24gdHJpbWVzdHJhbCAKIyBJbmljaWFuIGVuIGVsIHByaW1lciB0cmltZXN0cmUgZGUgMjAyMAojIFNlIGJ1c2NhIHByaW5vc3RpY2FyIGxhIHByb2R1Y2Npb24gZGUgbG9zIHNpZ3VpZW50ZXMKIyA1IHRyaW1lc3RyYWwKCnByb2R1Y2Npb24gPC1jKDUwLDUzLDU1LDU3LDU1LDYwKQoKc2VyaWVfZGVfdGllbXBvIDwtIHRzKGRhdGEgPSBwcm9kdWNjaW9uLCBzdGFydCA9IGMoMjAyMCwgMSksIGZyZXF1ZW5jeSA9IDQpCnNlcmllX2RlX3RpZW1wbwpgYGAKCiMjIDxzcGFuIHN0eWxlPSJjb2xvcjogcHVycGxlOyIgPjMuIENyZWFyIGVsIG1vZGVsbyBBUklNQTwvc3Bhbj4KKipBUklNQSoqIHNpZ25pZmljYSAqQXV0b1JlZ3Jlc3NpdmUgSW50ZWdyYXRlZCBNb3ZpbmcgQXZlcmFnZSogbyAKKk1vZGVsbyBBdXRvcnJlZ3Jlc2l2byBJbnRlZ3JhZG8gZGUgUHJvbWVkaW8gTW92aWwqLgoKYGBge3J9Cm1vZGVsbyA8LSBhdXRvLmFyaW1hKHNlcmllX2RlX3RpZW1wbywgRD0xKQptb2RlbG8Kc3VtbWFyeShtb2RlbG8pCmBgYAoKIyMgPHNwYW4gc3R5bGU9ImNvbG9yOiBwdXJwbGU7IiA+NC4gUmVhbGl6YXIgZWwgcHJvbsOzc3RpY288L3NwYW4+CmBgYHtyfQpwcm9ub3N0aWNvIDwtIGZvcmVjYXN0KG1vZGVsbywgbGV2ZWwgPSBjKDk1KSwgaD01KQpwcm9ub3N0aWNvCnBsb3QocHJvbm9zdGljbykKYGBgCgo=