Contexto

La base de datos AirPassengers contiene el número mensual de pasajeros aéreos internacionales entre 1949 y 1960 (144 observaciones mensuales).
Esta base de datos es ampliamente utilizada para el análisis de series de tiempo, reflejando el crecimiento de la industria de la aviación comercial después de la Segunda Guerra Mundial.

Instalar paquetes y llamar librerías

# install.packages("ggplot2") # Para generar gráficas con mejor diseño
library(ggplot2)
# install.packages("forecast") # Para revisar tendencias y generar pronósticos
library(forecast)

Importar base de datos

df <- as.data.frame(AirPassengers)
serie_de_tiempo <- AirPassengers

Identificar patrones o tendencias

Los 4 componentes de las series de tiempo son:

  • Tendencia: Ej. Aumento en las ventas de una tienda en línea debido al crecimiento del comercio electrónico.
  • Ciclo: Ej. Las fluctuaciones en las ventas de automóviles a lo largo de los ciclos económicos (expansión y recesión).
  • Patrón Estacional: Ej. Aumento en las ventas de juguetes durante la navidad.
  • Aleatorio: Ej. Caída de las ventas por un desastre natural.

ggplot(df, aes(x=time(serie_de_tiempo), y = serie_de_tiempo))+
  geom_line(color = "darkgreen") +
  labs(title="Pasajeros Aéreos Internacionales", x="Año", y="Número")
## Don't know how to automatically pick scale for object of type <ts>. Defaulting
## to continuous.
## Don't know how to automatically pick scale for object of type <ts>. Defaulting
## to continuous.

Podemos concluir de la gráfica anterior que la base de datos tiene tendencia creciente y patrón estacional.

Generar pronóstico

# Ajuste del modelo
modelo <- auto.arima(serie_de_tiempo)

# Realizar pronóstico a 12 meses
pronostico <- forecast(modelo, h=12)

# Graficar datos históricos y el pronóstico
autoplot(pronostico) +
  labs(title="Pasajeros Aéreos Internacionales", x="Año", y="Número")

# Generar predicción para un mes en específico
df_prediccion <- data.frame(Fecha = time(pronostico$mean), Prediccion = as.numeric(pronostico$mean))

df_prediccion[df_prediccion$Fecha == 1961.0, ] # Enero de 1961
##   Fecha Prediccion
## 1  1961   445.6349
LS0tDQp0aXRsZTogIkFpclBhc3NhbmdlcnMiDQphdXRob3I6ICJSYXVsIENhbnR1IEEwMTA4NzY4MyINCmRhdGU6ICIyMDI0LTA5LTEyIg0Kb3V0cHV0OiANCiAgaHRtbF9kb2N1bWVudDoNCiAgICB0b2M6IFRSVUUNCiAgICB0b2NfZmxvYXQ6IFRSVUUNCiAgICBjb2RlX2Rvd25sb2FkOiBUUlVFDQogICAgdGhlbWU6IGx1bWVuDQotLS0NCg0KIVtdKEM6XFxVc2Vyc1xccmF1bGNcXE9uZURyaXZlXFxFc2NyaXRvcmlvXFxCb290Y2FtcCBkZSBQcm9ncmFtYWNpw7NuXFxwbGFuZS5naWYpDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOiBncmVlbjsiPkNvbnRleHRvPC9zcGFuPg0KTGEgYmFzZSBkZSBkYXRvcyAqQWlyUGFzc2VuZ2VycyogY29udGllbmUgZWwgbsO6bWVybyBtZW5zdWFsIGRlIHBhc2FqZXJvcyBhw6lyZW9zIGludGVybmFjaW9uYWxlcyBlbnRyZSAxOTQ5IHkgMTk2MCAoMTQ0IG9ic2VydmFjaW9uZXMgbWVuc3VhbGVzKS4gIA0KRXN0YSBiYXNlIGRlIGRhdG9zIGVzIGFtcGxpYW1lbnRlIHV0aWxpemFkYSBwYXJhIGVsIGFuw6FsaXNpcyBkZSBzZXJpZXMgZGUgdGllbXBvLCByZWZsZWphbmRvIGVsIGNyZWNpbWllbnRvIGRlIGxhIGluZHVzdHJpYSBkZSBsYSBhdmlhY2nDs24gY29tZXJjaWFsIGRlc3B1w6lzIGRlIGxhICoqU2VndW5kYSBHdWVycmEgTXVuZGlhbCoqLg0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjogZ3JlZW47Ij5JbnN0YWxhciBwYXF1ZXRlcyB5IGxsYW1hciBsaWJyZXLDrWFzPC9zcGFuPg0KYGBge3IgbWVzc2FnZT1GQUxTRSwgd2FybmluZz1GQUxTRX0NCiMgaW5zdGFsbC5wYWNrYWdlcygiZ2dwbG90MiIpICMgUGFyYSBnZW5lcmFyIGdyw6FmaWNhcyBjb24gbWVqb3IgZGlzZcOxbw0KbGlicmFyeShnZ3Bsb3QyKQ0KIyBpbnN0YWxsLnBhY2thZ2VzKCJmb3JlY2FzdCIpICMgUGFyYSByZXZpc2FyIHRlbmRlbmNpYXMgeSBnZW5lcmFyIHByb27Ds3N0aWNvcw0KbGlicmFyeShmb3JlY2FzdCkNCmBgYA0KDQojIDxzcGFuIHN0eWxlPSJjb2xvcjogZ3JlZW47Ij5JbXBvcnRhciBiYXNlIGRlIGRhdG9zPC9zcGFuPg0KYGBge3J9DQpkZiA8LSBhcy5kYXRhLmZyYW1lKEFpclBhc3NlbmdlcnMpDQpzZXJpZV9kZV90aWVtcG8gPC0gQWlyUGFzc2VuZ2Vycw0KYGBgDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOiBncmVlbjsiPklkZW50aWZpY2FyIHBhdHJvbmVzIG8gdGVuZGVuY2lhczwvc3Bhbj4NCkxvcyA0IGNvbXBvbmVudGVzIGRlIGxhcyBzZXJpZXMgZGUgdGllbXBvIHNvbjogIA0KDQoqICoqVGVuZGVuY2lhOioqIEVqLiBBdW1lbnRvIGVuIGxhcyB2ZW50YXMgZGUgdW5hIHRpZW5kYSBlbiBsw61uZWEgZGViaWRvIGFsIGNyZWNpbWllbnRvIGRlbCBjb21lcmNpbyBlbGVjdHLDs25pY28uDQoqICoqQ2ljbG86KiogRWouIExhcyBmbHVjdHVhY2lvbmVzIGVuIGxhcyB2ZW50YXMgZGUgYXV0b23Ds3ZpbGVzIGEgbG8gbGFyZ28gZGUgbG9zIGNpY2xvcyBlY29uw7NtaWNvcyAoZXhwYW5zacOzbiB5IHJlY2VzacOzbikuDQoqICoqUGF0csOzbiBFc3RhY2lvbmFsOioqIEVqLiBBdW1lbnRvIGVuIGxhcyB2ZW50YXMgZGUganVndWV0ZXMgZHVyYW50ZSBsYSBuYXZpZGFkLg0KKiAqKkFsZWF0b3JpbzoqKiBFai4gQ2HDrWRhIGRlIGxhcyB2ZW50YXMgcG9yIHVuIGRlc2FzdHJlIG5hdHVyYWwuDQoNCiFbXShodHRwczovL3N0YXRpYy52ZWN0ZWV6eS5jb20vc3lzdGVtL3Jlc291cmNlcy90aHVtYm5haWxzLzAzOC85MjYvMTE1L3NtYWxsXzJ4L3RpbWUtc2VyaWVzLWNvbXBvbmVudHMtZm9yLXJhbmRvbS1ncmFwaC10cmVuZC1zZWFzb25hbC1jeWNsZS1ncmFwaC12ZWN0b3IuanBnKQ0KDQpgYGB7cn0NCmdncGxvdChkZiwgYWVzKHg9dGltZShzZXJpZV9kZV90aWVtcG8pLCB5ID0gc2VyaWVfZGVfdGllbXBvKSkrDQogIGdlb21fbGluZShjb2xvciA9ICJkYXJrZ3JlZW4iKSArDQogIGxhYnModGl0bGU9IlBhc2FqZXJvcyBBw6lyZW9zIEludGVybmFjaW9uYWxlcyIsIHg9IkHDsW8iLCB5PSJOw7ptZXJvIikNCmBgYA0KDQpQb2RlbW9zIGNvbmNsdWlyIGRlIGxhIGdyw6FmaWNhIGFudGVyaW9yIHF1ZSBsYSBiYXNlIGRlIGRhdG9zIHRpZW5lICoqdGVuZGVuY2lhIGNyZWNpZW50ZSoqIHkgKipwYXRyw7NuIGVzdGFjaW9uYWwqKi4NCg0KIyA8c3BhbiBzdHlsZT0iY29sb3I6IGdyZWVuOyI+R2VuZXJhciBwcm9uw7NzdGljbzwvc3Bhbj4NCmBgYHtyfQ0KIyBBanVzdGUgZGVsIG1vZGVsbw0KbW9kZWxvIDwtIGF1dG8uYXJpbWEoc2VyaWVfZGVfdGllbXBvKQ0KDQojIFJlYWxpemFyIHByb27Ds3N0aWNvIGEgMTIgbWVzZXMNCnByb25vc3RpY28gPC0gZm9yZWNhc3QobW9kZWxvLCBoPTEyKQ0KDQojIEdyYWZpY2FyIGRhdG9zIGhpc3TDs3JpY29zIHkgZWwgcHJvbsOzc3RpY28NCmF1dG9wbG90KHByb25vc3RpY28pICsNCiAgbGFicyh0aXRsZT0iUGFzYWplcm9zIEHDqXJlb3MgSW50ZXJuYWNpb25hbGVzIiwgeD0iQcOxbyIsIHk9Ik7Dum1lcm8iKQ0KDQojIEdlbmVyYXIgcHJlZGljY2nDs24gcGFyYSB1biBtZXMgZW4gZXNwZWPDrWZpY28NCmRmX3ByZWRpY2Npb24gPC0gZGF0YS5mcmFtZShGZWNoYSA9IHRpbWUocHJvbm9zdGljbyRtZWFuKSwgUHJlZGljY2lvbiA9IGFzLm51bWVyaWMocHJvbm9zdGljbyRtZWFuKSkNCg0KZGZfcHJlZGljY2lvbltkZl9wcmVkaWNjaW9uJEZlY2hhID09IDE5NjEuMCwgXSAjIEVuZXJvIGRlIDE5NjENCmBgYA0KDQo=