πŸ“Š ARIMA y Box-Jenkins

Modelos autorregresivos, medias móviles y la metodología sistemática de Box-Jenkins

AR(p) MA(q) ARIMA(p,d,q) ACF/PACF
Módulo 2 de 12
πŸ’‘ PiΓ©nsalo asΓ­: ARIMA es como un chef que predice el siguiente ingrediente: los modelos AR usan ingredientes anteriores (pasado), MA corrige errores previos (como ajustar la sal), y la diferenciación elimina tendencias (como nivelar una mesa coja). Box-Jenkins es la receta paso a paso para encontrar el mejor modelo.
  • Comprender los modelos AR(p), MA(q) y ARIMA(p,d,q) desde cero.
  • Interpretar los gráficos ACF y PACF para identificar el orden del modelo.
  • Aplicar la metodología Box-Jenkins en sus 4 pasos: identificación, estimación, validación y pronóstico.
  • Diferenciar entre series estacionarias y no estacionarias, y saber cuándo aplicar diferenciación.
  • Modelos ARIMA

    πŸ“ˆ AR(p) β€” Autorregresivo

    El modelo AR(p) asume que el valor actual de una serie depende linealmente de sus propios valores pasados. Piensa en la inercia: si hoy hace calor, es probable que mañana también. La "p" indica cuántos valores pasados usamos.

    La fórmula es: yβ‚œ = φ₁yβ‚œβ‚‹β‚ + Ο†β‚‚yβ‚œβ‚‹β‚‚ + ... + Ο†β‚šyβ‚œβ‚‹β‚š + Ξ΅β‚œ. El término Ξ΅β‚œ es el error aleatorio (ruido blanco) que no podemos predecir. Cada coeficiente Ο†α΅’ nos dice qué tan importante es ese valor pasado.

    Ejemplo real: Para modelizar el PIB mensual, un AR(2) usaría el PIB de los últimos 2 meses para predecir el siguiente. Si φ₁ = 0.8 y Ο†β‚‚ = -0.3, el PIB de hace 1 mes influye positivamente y el de hace 2 meses tiene un efecto de corrección.

    Peligro común: Si la suma de los coeficientes AR es cercana a 1, la serie puede ser no estacionaria (raíz unitaria). En AR(1), si Ο† = 1, la serie explota. Siempre verifica que |Ο†| < 1 para AR(1).

    yβ‚œ = φ₁yβ‚œβ‚‹β‚ + ... + Ο†β‚šyβ‚œβ‚‹β‚š + Ξ΅β‚œ

    🌊 MA(q) β€” Media Móvil

    El modelo MA(q) usa los errores pasados para predecir el presente. Cada vez que el modelo se equivoca (el error Ξ΅β‚œ), esa información se usa para corregir la predicción de los próximos períodos. Es como aprender de los errores.

    Fórmula: yβ‚œ = Ξ΅β‚œ + ΞΈβ‚Ξ΅β‚œβ‚‹β‚ + ΞΈβ‚‚Ξ΅β‚œβ‚‹β‚‚ + ... + ΞΈβ‚™Ξ΅β‚œβ‚‹q. Los errores Ξ΅ son ruido blanco con media cero y varianza constante. Los coeficientes ΞΈ miden el impacto de cada error pasado.

    Ejemplo real: En finanzas, un MA(1) captura el efecto de "noticia de último momento": si hoy hubo un shock (buena/mala noticia), su efecto se disipa rápidamente en el siguiente período. ΞΈ negativo significa que el shock se revierte parcialmente.

    Peligro común: Los modelos MA requieren que los coeficientes cumplan la "invertibilidad" (|ΞΈ| < 1 para MA(1)). Si no, el modelo no puede representarse como un AR infinito y es inestable.

    yβ‚œ = Ξ΅β‚œ + θ₁εβ‚œβ‚‹β‚ + ... + θβ‚™εβ‚œβ‚‹q

    πŸ”§ ARIMA(p,d,q) β€” El modelo completo

    ARIMA integra tres componentes: AR (autorregresivo), I (diferenciación, del inglés "Integrated") y MA (media móvil). La "d" indica cuántas diferencias aplicamos para hacer la serie estacionaria.

    Intuición de la diferenciación: Si una serie sube constantemente (tendencia), sus diferencias (cambios) son estables. La primera diferencia es Ξ”yβ‚œ = yβ‚œ - yβ‚œβ‚‹β‚. Si aún hay tendencia, aplicamos segunda diferencia: Δ²yβ‚œ = Ξ”yβ‚œ - Ξ”yβ‚œβ‚‹β‚.

    Ejemplo: Para el PIB argentino (no estacionario en niveles pero sí en diferencias), usar d=1. Luego modelamos las diferencias con ARMA(p,q). El modelo final es ARIMA(p,1,q).

    Peligro común: Sobrediferenciar (d muy grande) introduce correlación espuria y empeora el pronóstico. Usa solo la d mínima necesaria para lograr estacionariedad, verificada con tests ADF o KPSS.

    (1-B)ᡈyβ‚œ = ARMA(p,q)

    πŸ“Š ACF y PACF

    La función de autocorrelación (ACF) mide la correlación entre yβ‚œ e yβ‚œβ‚‹β‚–. Si el ACF decae lentamente, la serie probablemente es no estacionaria o tiene un componente AR fuerte.

    La función de autocorrelación parcial (PACF) mide la correlación entre yβ‚œ e yβ‚œβ‚‹β‚– después de eliminar el efecto de los valores intermedios. ¡Es como la correlación "pura" entre dos puntos!

    Regla de identificación: Si el PACF se corta abruptamente en el rezago p y el ACF decae β†’ AR(p). Si el ACF se corta en q y el PACF decae β†’ MA(q). Si ambos decaen β†’ ARMA(p,q).

    Peligro común: No confíes ciegamente en ACF/PACF con muestras pequeñas (T<50). Las bandas de confianza (Β±1.96/√T) son aproximadas; úsalas como guía, no como verdad absoluta.

    ACF(k) = Corr(yβ‚œ, yβ‚œβ‚‹β‚–)
    Metodología Box-Jenkins

    1. Identificación

    Usar ACF/PACF para identificar p, d, q tentativos. Examinar estacionariedad con tests ADF o KPSS. Si la serie es no estacionaria, diferenciar hasta que lo sea.

    2. Estimación

    Estimar parámetros por máxima verosimilitud o mínimos cuadrados condicionales. Los coeficientes deben ser estadísticamente significativos (p<0.05).

    3. Validación

    Verificar residuos tipo ruido blanco. Tests de Ljung-Box (Hβ‚€: no autocorrelación), normalidad y heterocedasticidad. Si los residuos tienen patrón, el modelo es inadecuado.

    4. Pronóstico

    Generar predicciones puntuales e intervalos de confianza que se ensanchan con el horizonte. Comparar con modelos alternativos usando AIC o BIC.

    Simulación: Modelos ARMA

    Genera datos de AR(1), MA(1) o ARMA(1,1) y observa ACF/PACF. Identifica el modelo correcto.

    Presiona "Generar" para simular.
    Ejemplo Paso a Paso: AR(1) Numérico
    Paso 1: Definir el modelo

    Supón un AR(1): yβ‚œ = 0.7Β·yβ‚œβ‚‹β‚ + Ξ΅β‚œ, con Ξ΅β‚œ ~ N(0,1). El coeficiente Ο† = 0.7 indica que el valor actual retiene el 70% del valor anterior. Como |0.7| < 1, el proceso es estacionario.

    Paso 2: Generar datos iniciales

    Si yβ‚€ = 0 y Ρ₁ = 0.5, entonces y₁ = 0.7(0) + 0.5 = 0.5. Si Ξ΅β‚‚ = -0.3, entonces yβ‚‚ = 0.7(0.5) + (-0.3) = 0.35 - 0.3 = 0.05. ¡Así se construye la serie paso a paso!

    Paso 3: Calcular el ACF teórico

    Para un AR(1), el ACF decae geométricamente: ρ(k) = φᡏ. Con Ο†=0.7: ρ(1)=0.7, ρ(2)=0.49, ρ(3)=0.343, ρ(4)=0.24. El PACF tiene un solo pico en k=1 (valor 0.7) y ceros después. Esta es la "firma" de un AR(1).

    Paso 4: Pronosticar

    Para pronosticar yβ‚œβ‚Šβ‚ desde el último valor observado yβ‚œ: Ε·β‚œβ‚Šβ‚ = 0.7Β·yβ‚œ. El intervalo de confianza al 95% es Ε·β‚œβ‚Šβ‚ Β± 1.96Β·Οƒβ‚‘, donde Οƒβ‚‘ = 1 (la desviación del error). Para horizontes más lejanos, el pronóstico converge a 0 (la media incondicional).

    Mini Quiz: ARIMA y Box-Jenkins
    πŸ’‘ Responde las preguntas sobre modelos ARIMA y Box-Jenkins.

    πŸ“Œ Resumen del módulo

  • AR(p) modela dependencia con valores pasados; MA(q) modela dependencia con errores pasados.
  • ARIMA(p,d,q) combina AR + diferenciación + MA para modelar series no estacionarias.
  • El ACF identifica el orden MA(q); el PACF identifica el orden AR(p).
  • Box-Jenkins tiene 4 pasos: identificar, estimar, validar, pronosticar. Nunca saltarse la validación.