Instalar paquetes y llamar librerías

# install.packages("forecast")
library(forecast)
## Registered S3 method overwritten by 'quantmod':
##   method            from
##   as.zoo.data.frame zoo

Crear serie de Tiempo

produccion <- c(19,18,15,20,18,22,20)

produccion_st <- ts(data = produccion, start = c(2022, 2), frequency = 12)
produccion_st
##      Feb Mar Apr May Jun Jul Aug
## 2022  19  18  15  20  18  22  20

ARIMA: Modelo Autorregresivo Integrado de Media Móvil (AutoRegressive Integrated Moving Average)

ARIMA (p, d, q)

p = orden de auto-regresión

d = orden de integración (diferenciación)

q = orden del promedio móvil

¿Cuándo de usa?

Cuando las estimaciones futuras se explican por los datos del pasado y no por variables independientes.

Ejemplo: Tipo de Cambio

Modelo

modelo <- auto.arima(produccion_st)
modelo
## Series: produccion_st 
## ARIMA(0,0,0) with non-zero mean 
## 
## Coefficients:
##          mean
##       18.8571
## s.e.   0.7674
## 
## sigma^2 = 4.81:  log likelihood = -14.89
## AIC=33.78   AICc=36.78   BIC=33.67

Pronóstico

pronostico <- forecast(modelo, level=c(95), h=5)
pronostico
##          Point Forecast    Lo 95    Hi 95
## Sep 2022       18.85714 14.55882 23.15547
## Oct 2022       18.85714 14.55882 23.15547
## Nov 2022       18.85714 14.55882 23.15547
## Dec 2022       18.85714 14.55882 23.15547
## Jan 2023       18.85714 14.55882 23.15547
plot(pronostico)

Conclusiones

En este código trabajamos con Pronósticos que nos sirve para poder predecir las variables que nos gustaria conocer, por ejemplo, en el caso de tema de negocios nos encontramos con pronosticos que nos ayudan a conocer futuras ventas, producción, exportaciones, etc. A pesar de que no son perfectos podemos tener estimaciones que pueden acercarse a la realidad.

En esta caso creamos un modelo Autorregresivo para poder obtener la predicción de producción de los siguientes meses dando como resultado una tendencia constante para los siguientes meses, imaginando que son 18 piezas de producción y que nuestro intervalo de confianza nos dice que se podría mantener entre 14.55882-23.15547.

LS0tCnRpdGxlOiA8c3BhbiBzdHlsZT0iY29sb3I6Z29sZCI+UHJvbm9zdGljb3M8L3NwYW4+CmF1dGhvcjogIkthcmxhIEdlb3JnaW5hIE1hcnTDrW5leiBHb256w6FsZXogQTAwODI3NTAwIgpkYXRlOiAiMjAyMi0xMC0wNSIKb3V0cHV0OiAKICBodG1sX2RvY3VtZW50OgogICAgdG9jOiB0cnVlCiAgICB0b2NfZmxvYXQ6IHRydWUKICAgIGNvZGVfZG93bmxvYWQ6IHRydWUKLS0tCgo8aW1nIHNyYz0gIi9Vc2Vycy9nZW9yZ2luYW1hcnRpbmV6L0RvY3VtZW50cy9UZWMvU2XMgXB0aW1vIFNlbWVzdHJlL0FuYWxpzIF0aWNhIHBhcmEgbmVnb2Npb3MsIGRlIGxvcyBkYXRvcyBhIGRlY2lzaW9uZXMvcHJvbm9zdGljb3MuanBlZyI+CgojIDxzcGFuIHN0eWxlPSJjb2xvcjpkYXJrYmx1ZSI+SW5zdGFsYXIgcGFxdWV0ZXMgeSBsbGFtYXIgbGlicmVyw61hczwvc3Bhbj4KYGBge3J9CiMgaW5zdGFsbC5wYWNrYWdlcygiZm9yZWNhc3QiKQpsaWJyYXJ5KGZvcmVjYXN0KQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmRhcmtibHVlIj5DcmVhciBzZXJpZSBkZSBUaWVtcG88L3NwYW4+CmBgYHtyfQpwcm9kdWNjaW9uIDwtIGMoMTksMTgsMTUsMjAsMTgsMjIsMjApCgpwcm9kdWNjaW9uX3N0IDwtIHRzKGRhdGEgPSBwcm9kdWNjaW9uLCBzdGFydCA9IGMoMjAyMiwgMiksIGZyZXF1ZW5jeSA9IDEyKQpwcm9kdWNjaW9uX3N0CmBgYAoKIyA8c3BhbiBzdHlsZT0iY29sb3I6ZGFya2JsdWUiPkFSSU1BOiBNb2RlbG8gQXV0b3JyZWdyZXNpdm8gSW50ZWdyYWRvIGRlIE1lZGlhIE3Ds3ZpbCAoQXV0b1JlZ3Jlc3NpdmUgSW50ZWdyYXRlZCBNb3ZpbmcgQXZlcmFnZSk8L3NwYW4+CiMjIEFSSU1BIChwLCBkLCBxKSAgCiMjIyBwID0gb3JkZW4gZGUgYXV0by1yZWdyZXNpw7NuICAKIyMjIGQgPSBvcmRlbiBkZSBpbnRlZ3JhY2nDs24gKGRpZmVyZW5jaWFjacOzbikgIAojIyMgcSA9IG9yZGVuIGRlbCBwcm9tZWRpbyBtw7N2aWwgIAoKIyMgwr9DdcOhbmRvIGRlIHVzYT8gIAojIyMgQ3VhbmRvIGxhcyBlc3RpbWFjaW9uZXMgZnV0dXJhcyBzZSBleHBsaWNhbiBwb3IgbG9zIGRhdG9zIGRlbCBwYXNhZG8geSBubyBwb3IgdmFyaWFibGVzIGluZGVwZW5kaWVudGVzLiAgCgojIyMgRWplbXBsbzogVGlwbyBkZSBDYW1iaW8gIAoKIyMgTW9kZWxvCmBgYHtyfQptb2RlbG8gPC0gYXV0by5hcmltYShwcm9kdWNjaW9uX3N0KQptb2RlbG8KYGBgCgojIyBQcm9uw7NzdGljbwpgYGB7cn0KcHJvbm9zdGljbyA8LSBmb3JlY2FzdChtb2RlbG8sIGxldmVsPWMoOTUpLCBoPTUpCnByb25vc3RpY28KcGxvdChwcm9ub3N0aWNvKQpgYGAKCiMgPHNwYW4gc3R5bGU9ImNvbG9yOmRhcmtibHVlIj5Db25jbHVzaW9uZXM8L3NwYW4+CkVuIGVzdGUgY8OzZGlnbyB0cmFiYWphbW9zIGNvbiA8c3BhbiBzdHlsZT0iY29sb3I6Ymx1ZSI+KipQcm9uw7NzdGljb3MqKjwvc3Bhbj4gcXVlIG5vcyBzaXJ2ZSBwYXJhIHBvZGVyIHByZWRlY2lyIGxhcyB2YXJpYWJsZXMgcXVlIG5vcyBndXN0YXJpYSBjb25vY2VyLCBwb3IgZWplbXBsbywgZW4gZWwgY2FzbyBkZSB0ZW1hIGRlIG5lZ29jaW9zIG5vcyBlbmNvbnRyYW1vcyBjb24gcHJvbm9zdGljb3MgcXVlIG5vcyBheXVkYW4gYSBjb25vY2VyIGZ1dHVyYXMgdmVudGFzLCBwcm9kdWNjacOzbiwgZXhwb3J0YWNpb25lcywgZXRjLiBBIHBlc2FyIGRlIHF1ZSBubyBzb24gcGVyZmVjdG9zIHBvZGVtb3MgdGVuZXIgPHNwYW4gc3R5bGU9ImNvbG9yOnJlZCI+KmVzdGltYWNpb25lcyo8L3NwYW4+IHF1ZSBwdWVkZW4gYWNlcmNhcnNlIGEgbGEgcmVhbGlkYWQuICAKCkVuIGVzdGEgY2FzbyBjcmVhbW9zIHVuIG1vZGVsbyBBdXRvcnJlZ3Jlc2l2byBwYXJhIHBvZGVyIG9idGVuZXIgbGEgcHJlZGljY2nDs24gZGUgcHJvZHVjY2nDs24gZGUgbG9zIHNpZ3VpZW50ZXMgbWVzZXMgZGFuZG8gY29tbyByZXN1bHRhZG8gdW5hIHRlbmRlbmNpYSBjb25zdGFudGUgcGFyYSBsb3Mgc2lndWllbnRlcyBtZXNlcywgaW1hZ2luYW5kbyBxdWUgc29uIDE4IHBpZXphcyBkZSBwcm9kdWNjacOzbiB5IHF1ZSBudWVzdHJvIDxzcGFuIHN0eWxlPSJjb2xvcjpyZWQiPmludGVydmFsbyBkZSBjb25maWFuemE8L3NwYW4+IG5vcyBkaWNlIHF1ZSBzZSBwb2Ryw61hIG1hbnRlbmVyIGVudHJlIDxzcGFuIHN0eWxlPSJjb2xvcjpyZWQiPjE0LjU1ODgyLTIzLjE1NTQ3Ljwvc3Bhbj4KCg==