Instalar paquetes y llamar librerias

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

Crear serie de tiempo

produccion <- c(100, 200, 300, 400, 500)

produccion_st <- ts(data = produccion, start = c(2018, 1), frequency = 1)
produccion_st
## Time Series:
## Start = 2018 
## End = 2022 
## Frequency = 1 
## [1] 100 200 300 400 500

ARIMA: Modelo Autoregresivo Integrado de Media Movil (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 se usa?

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

Ejemplo: Tipo de cambio

modelo <- auto.arima(produccion_st)
modelo
## Series: produccion_st 
## ARIMA(0,1,0) with drift 
## 
## Coefficients:
## drift  
##   100  
## 
## sigma^2 = 0:  log likelihood = Inf
## AIC=-Inf   AICc=-Inf   BIC=-Inf
pronostico <- forecast(modelo, level = c(95), h=5)
pronostico
##      Point Forecast Lo 95 Hi 95
## 2023            600   600   600
## 2024            700   700   700
## 2025            800   800   800
## 2026            900   900   900
## 2027           1000  1000  1000
plot(pronostico)

Conclusión:

Con esta herramienta pudimos predecir las piezas de producción en los próximos años, teniendo como información la producción actual desde 2018. En la gráfica realizada podemos ver que la producción tiene un crecimiento lineal positivo, creciendo 100 piezas por año. Por lo que para 2027, se espera una producción de 1000 piezas.

LS0tDQp0aXRsZTogIDxzcGFuIHN0eWxlPSJDb2xvcjpCbHVlIiA+IFByb27Ds3N0aWNvcw0KYXV0aG9yOiAiTWFyaWFuYSBVbGxvYSAtIEEwMTI1MzA3NiINCmRhdGU6ICIyMDIyLTEwLTA0Ig0Kb3V0cHV0OiANCiAgaHRtbF9kb2N1bWVudDoNCiAgICB0b2M6IHRydWUNCiAgICB0b2NfZmxvYXQ6IHRydWUNCiAgICBjb2RlX2Rvd25sb2FkOiB0cnVlDQotLS0NCjxpbWcgc3JjPSAiQzpcXFVzZXJzXFxtYXJpYVxcRG9jdW1lbnRzXFxJVEVTTSBMQUVUXFxTZW1lc3RyZSA3XFxNM1xcaW1hZ2VuZXNcXHByb25vc3RpY28ucG5nIj4NCg0KIyMjIEluc3RhbGFyIHBhcXVldGVzIHkgbGxhbWFyIGxpYnJlcmlhcw0KYGBge3J9DQojIGluc3RhbGwucGFja2FnZXMoImZvcmVjYXN0IikNCmxpYnJhcnkoZm9yZWNhc3QpDQpgYGANCg0KIyMjIENyZWFyIHNlcmllIGRlIHRpZW1wbw0KYGBge3J9DQpwcm9kdWNjaW9uIDwtIGMoMTAwLCAyMDAsIDMwMCwgNDAwLCA1MDApDQoNCnByb2R1Y2Npb25fc3QgPC0gdHMoZGF0YSA9IHByb2R1Y2Npb24sIHN0YXJ0ID0gYygyMDE4LCAxKSwgZnJlcXVlbmN5ID0gMSkNCnByb2R1Y2Npb25fc3QNCmBgYA0KDQojIyMgQVJJTUE6IE1vZGVsbyBBdXRvcmVncmVzaXZvIEludGVncmFkbyBkZSBNZWRpYSBNb3ZpbCAoQXV0b3JlZ3Jlc3NpdmUgSW50ZWdyYXRlZCBNb3ZpbmcgQXZlcmFnZSkNCioqQVJJTUEgKHAsZCxxKSoqDQoNCnAgPSBvcmRlbiBkZSBhdXRvLXJlZ3Jlc2nDs24NCg0KZCA9IG9yZGVuIGRlIGludGVncmFjacOzbiAoZGlmZXJlbmNpYWNpw7NuKQ0KDQpxID0gb3JkZW4gZGVsIHByb21lZGlvIG3Ds3ZpbA0KDQoqKsK/Q3XDoW5kbyBzZSB1c2E/KioNCg0KQ3VhbmRvIGxhcyBlc3RpbWFjaW9uZXMgZnV0dXJhcyBzZSBleHBsaWNhbiBwb3IgbG9zIGRhdG9zIGRlbCBwYXNhZG8geSBubyBwb3IgdmFyaWFibGVzIGluZGVwZW5kaWVudGVzLg0KDQpFamVtcGxvOiAqKlRpcG8gZGUgY2FtYmlvKioNCg0KYGBge3J9DQptb2RlbG8gPC0gYXV0by5hcmltYShwcm9kdWNjaW9uX3N0KQ0KbW9kZWxvDQoNCnByb25vc3RpY28gPC0gZm9yZWNhc3QobW9kZWxvLCBsZXZlbCA9IGMoOTUpLCBoPTUpDQpwcm9ub3N0aWNvDQpwbG90KHByb25vc3RpY28pDQpgYGANCg0KIyMjIENvbmNsdXNpw7NuOg0KQ29uIGVzdGEgaGVycmFtaWVudGEgcHVkaW1vcyBwcmVkZWNpciBsYXMgcGllemFzIGRlIHByb2R1Y2Npw7NuIGVuIGxvcyBwcsOzeGltb3MgYcOxb3MsIHRlbmllbmRvIGNvbW8gaW5mb3JtYWNpw7NuIGxhIHByb2R1Y2Npw7NuIGFjdHVhbCBkZXNkZSAyMDE4LiBFbiBsYSBncsOhZmljYSByZWFsaXphZGEgcG9kZW1vcyB2ZXIgcXVlIGxhICoqcHJvZHVjY2nDs24gdGllbmUgdW4gY3JlY2ltaWVudG8gbGluZWFsIHBvc2l0aXZvKiosICpjcmVjaWVuZG8gMTAwIHBpZXphcyBwb3IgYcOxbyouIFBvciBsbyBxdWUgcGFyYSAyMDI3LCBzZSBlc3BlcmEgdW5hIHByb2R1Y2Npw7NuIGRlIDEwMDAgcGllemFzLg==