Fundamentos y Descomposición
Definición formal, componentes clásicas (tendencia, estacionalidad, ciclo, error) y modelos de descomposición aditiva y multiplicativa. Análisis descriptivo con R.
¿Por qué estudiar series de tiempo?
Imagina que eres el analista de datos de una aerolínea y tu jefe te pregunta cuántos pasajeros esperar el próximo diciembre. Tienes registros mensuales de los últimos diez años. ¿Cómo respondes?
Los métodos estadísticos clásicos —regresión ordinaria, pruebas t, ANOVA— asumen que las observaciones son independientes. Pero en una serie de tiempo, el valor de diciembre pasado sí contiene información sobre el de diciembre próximo. Ignorar eso equivale a tirar información valiosa a la basura.
Las series de tiempo tienen su propio lenguaje: tendencia, estacionalidad, correlación serial, raíces unitarias. Este módulo te enseña ese lenguaje desde cero. En los módulos siguientes lo usarás para construir modelos que predicen, descomponen y comunican el comportamiento de datos reales.
1.1 Definición de serie de tiempo
Una serie de tiempo es una secuencia de observaciones registradas en momentos equidistantes. La notamos como , donde indexa el tiempo. La clave no es solo que los datos estén ordenados: es que la dependencia entre observaciones cercanas encierra información que un modelo bien construido puede aprovechar.
Ejemplos que encontrarás en el curso
- AirPassengers — pasajeros aéreos mensuales (1949–1960). Tendencia creciente + estacionalidad multiplicativa.
- co2 — concentración de CO₂ atmosférico (Mauna Loa). Tendencia lineal + ciclo anual casi perfecto.
- nottem — temperaturas medias mensuales en Nottingham. Estacionalidad pura, sin tendencia.
- JohnsonJohnson — ventas trimestrales de J&J. Crecimiento exponencial + estacionalidad multiplicativa.
- Ventas de licor (EE.UU.) — series mensuales usadas en los apuntes de la profesora Nelfi González (cap. 3): oscilaciones de diciembre que crecen con el nivel → modelo multiplicativo.
La característica esencial que distingue una serie de tiempo de un conjunto de datos transversal es que las observaciones no son independientes: el valor de ayer contiene información sobre el de hoy, y el de hoy sobre el de mañana. Esta dependencia se mide con la función de autocorrelación (ACF), herramienta que explorarás en los módulos siguientes.
1.2 Componentes clásicas
Antes de ajustar cualquier modelo conviene descomponer la serie en capas interpretables. La idea central es que lo que observas, , es la superposición de fenómenos con naturaleza distinta. Identificarlos por separado permite modelarlos con las herramientas adecuadas.
Movimiento suave y sostenido de largo plazo. Captura si la serie crece, decrece o permanece estable en el tiempo. Se modela con polinomios, logaritmos o suavizamiento.
Patrón periódico que se repite dentro del año. El periodo s = 12 para datos mensuales, s = 4 para trimestrales, s = 52 para semanales. Se modela con variables dummy o series de Fourier.
Fluctuaciones de largo plazo no periódicas, asociadas a ciclos económicos. Su duración es variable (2–10 años). En este curso se incorpora al término de error.
Componente aleatoria residual. En un buen modelo debe comportarse como ruido blanco: Et ~ iid N(0, σ²). Si no lo es, hay estructura que el modelo aún no capturó.
Convención del curso
Nos concentramos en tendencia y estacionalidad. El ciclo se absorbe en . Si los residuos muestran correlación serial —visible en el correlograma— es señal de que hay estructura cíclica que el modelo aún no capturó. Eso lo resolveremos con modelos ARIMA/SARIMA en módulos avanzados.
1.3 Modelos de descomposición
La manera en que se combinan las componentes define el tipo de modelo. Elegir correctamente entre aditivo y multiplicativo no es un detalle menor: afecta la interpretación de cada parámetro y la forma en que se recuperan las predicciones. El primer paso para decidir es siempre graficar la serie y mirar si las oscilaciones estacionales crecen proporcionalmente con el nivel.
Modelo aditivo
Apropiado cuando la variabilidad estacional es constante —las oscilaciones no se amplían aunque la tendencia suba:
Aquí , y están en las mismas unidades que . Si la serie mide toneladas, la tendencia y la estacionalidad también se interpretan en toneladas.
Modelo completamente multiplicativo
Cuando la amplitud estacional crece proporcionalmente con el nivel de la serie (caso típico: ventas, producción industrial, pasajeros aéreos), el modelo correcto es:
El término es el factor de error multiplicativo: siempre positivo, y con valor neutro exactamente 1 cuando (su media). Esto garantiza que y representen el nivel estructural de la serie sin ningún sesgo adicional: cuando no hay perturbación aleatoria, exactamente.
¿Por qué no escribir simplemente Y_t = T_t × S_t × E_t?
Una notación muy frecuente usa un error multiplicativo genérico sin especificar su distribución. El problema: para que represente el nivel esperado de la serie, necesitaríamos . Pero si además asumimos que —como exige el modelo lineal en log—, entonces y su media es . Las dos condiciones no pueden satisfacerse simultáneamente. La formulación con resuelve esta ambigüedad declarando desde el principio que el error gaussiano está en escala logarítmica.
Tomando logaritmo natural de ambos lados se obtiene directamente un modelo lineal aditivo estándar:
Nótese que el error aparece solo —sin transformar— en la ecuación log. No hay ambigüedad: es un error aditivo gaussiano con media cero, exactamente lo que requieren los métodos de regresión lineal. Bajo este modelo, la serie en escala log es , y se puede estimar por mínimos cuadrados con todas las herramientas habituales.
⚠️ Corrección por sesgo al volver a escala original
Una vez ajustado el modelo, la predicción está en escala log. Aplicar de forma directa estima la mediana de , no su media. La razón se obtiene directamente del modelo:
El factor proviene de la propiedad de la distribución log-normal. Para predecir la media, la corrección correcta —conocida como corrección de Duan o smearing estimator— es:
donde es el error cuadrático medio del ajuste en escala log. Si , el factor es : una corrección del 2%. Si , el factor es : un 13% —ya no es despreciable. En el Módulo 4 verás cómo implementar esto automáticamente en R.
Modelo parcialmente multiplicativo
Cuando tendencia y estacionalidad interactúan multiplicativamente pero el error es aditivo:
Este modelo es útil cuando la varianza del error es estable pero la amplitud estacional sí escala con la tendencia. Es menos común pero aparece en series de consumo energético y producción agrícola.
Regla práctica para elegir el tipo de modelo
Paso 1 — Grafique la serie original. ¿Las oscilaciones de verano (o diciembre, o el pico que sea) se ven iguales de altas a lo largo del tiempo, o van creciendo?
Paso 2 — Si crecen proporcionalmente → trabaje con log(Y). Si son constantes → trabaje con Y directamente.
Paso 3 — Después de ajustar el modelo, grafique los residuos. Si la varianza de los residuos cambia con el tiempo (más dispersos al final), la transformación logarítmica probablemente era necesaria.
Este es el ciclo que repetirás en cada módulo: graficar → decidir transformación → ajustar → revisar residuos → iterar. Modelar es esencialmente un proceso visual.
1.4 Análisis descriptivo con R
El primer paso de cualquier análisis de series de tiempo es siempre el mismo: graficar. Antes de correr una prueba de hipótesis, antes de estimar un parámetro, antes de elegir un modelo — grafica la serie y déjala hablar. La función decompose() hace eso de manera automática: descompone la serie en sus cuatro componentes usando medias móviles para la tendencia y promedios por período para la estacionalidad.
¿Qué vas a ver al ejecutar el código?
El siguiente bloque produce tres gráficos:
- Serie original — identifica tendencia y estacionalidad a ojo.
- Descomposición multiplicativa — cuatro paneles: serie observada, tendencia suavizada, índices estacionales (factores que multiplican la tendencia cada mes), y residuos.
- log(serie) — observa cómo la transformación estabiliza la varianza: las oscilaciones se vuelven de amplitud constante.
💡 Cambia 'multiplicative' por 'additive' y observa la diferencia en el panel seasonal: en el aditivo verás valores absolutos (ej. +60 pasajeros en julio), en el multiplicativo verás factores relativos (ej. ×1.35). También prueba con nottem, co2 o JohnsonJohnson.
¿Qué deberías ver en los resultados?
En decompose(AirPassengers, type="multiplicative"):
- Tendencia (trend): curva suavizada que crece de ~125 a ~430. Los primeros y últimos 6 valores serán
NA—esto es normal, es el costo de usar medias móviles. - Estacionalidad (seasonal): un patrón que se repite igual cada año. Los valores en julio serán cerca de 1.35 (35% sobre la tendencia); en noviembre cerca de 0.88 (12% bajo la tendencia).
- Residuos (random): deben oscilar alrededor de 1.0 sin patrón visible. Si ves una curva en los residuos, el modelo de tendencia no capturó algo.
Error frecuente: confundir aditivo y multiplicativo
Si usas type="additive" en una serie multiplicativa como AirPassengers, los residuos mostrarán un patrón en forma de embudo: serán pequeños al inicio y grandes al final. Eso es evidencia de que la descomposición es incorrecta —hay varianza no capturada. La corrección es cambiar a type="multiplicative" o trabajar con log(serie).
Hay una serie que ilustra el caso multiplicativo de manera aún más dramática que AirPassengers: las ventas mensuales de licor en Estados Unidos. En los apuntes de la profesora Nelfi González (Figura 3.11, Capítulo 3) puedes ver cómo las oscilaciones de diciembre —el mes pico— se hacen progresivamente más grandes a medida que crece el nivel de ventas. Aplicando log(), la serie se transforma en un patrón aditivo perfectamente tratable con regresión lineal. Ese es el poder de la transformación logarítmica.
💡 En el ggseasonplot observa si las líneas de distintos años son paralelas (aditivo) o si se van separando (multiplicativo). En el ggsubseriesplot la línea horizontal de cada mes es su media histórica — una 'guía de lo normal' para ese mes.
Cómo interpretar el ggseasonplot
Cada línea de color representa un año. Si las líneas son paralelas —misma forma, distintas alturas— la estacionalidad es aditiva: el efecto de julio es siempre +60 pasajeros, sin importar el año. Si las líneas se van separando hacia los meses pico, la estacionalidad es multiplicativa: el efecto de julio en 1960 es proporcional al nivel de 1960, que es mayor que el de 1950.
Esta es la lectura que harás en cada nueva serie que analices: graficar primero, decidir transformación después.
1.5 Visualización conceptual
El siguiente diagrama ilustra cómo las componentes se combinan en la serie observada bajo el modelo aditivo . La clave para leer una serie de tiempo es aprender a separar mentalmente estas capas antes de modelar.
Figura 1.1 — Descomposición aditiva: la serie observada es la suma de tendencia (azul), estacionalidad (verde) y error (la diferencia entre la curva negra y la suma de azul+verde).
Limitación de decompose(): los NA en los extremos
decompose() estima la tendencia con una media móvil centrada de longitud (el período estacional). Eso significa que necesita observaciones a cada lado de cada punto para calcular la media. El resultado: los primeros y últimos valores de la tendencia son NA. Para datos mensuales (s=12) pierdes 6 observaciones en cada extremo.
Esto no es un error —es una limitación matemática inherente a las medias móviles. Para estimaciones completas hasta el último dato disponible, usa stl() (Módulo 5) o modela directamente con regresión (Módulos 2–4).
1.6 Resumen y conexión con los módulos siguientes
En este módulo estableciste el vocabulario fundamental: serie de tiempo, tendencia, estacionalidad, modelo aditivo vs multiplicativo, y descomposición exploratoria. Estos conceptos son los cimientos sobre los que se construyen todos los modelos del curso.
El camino que seguiremos es siempre el mismo ciclo:
En el Módulo 2 modelarás la tendencia con polinomios y funciones logarítmicas. En el Módulo 3 añadirás la estacionalidad con variables dummy y series de Fourier. En el Módulo 4 ensamblarás el modelo completo, aprenderás a comparar alternativas con AIC/BIC, y aplicarás la corrección de sesgo al predecir en escala original. El diagnóstico de residuos —el paso más importante y el más omitido— lo profundizarás en el Módulo 8.
La pregunta que guía cada módulo
Antes de ajustar un modelo nuevo, hazte siempre esta pregunta: ¿los residuos del modelo anterior tienen estructura? Si tienen tendencia, falta modelar la tendencia. Si tienen patrón estacional, falta modelar la estacionalidad. Si tienen autocorrelación, falta un componente ARIMA. Un buen modelo es aquel cuyos residuos son — finalmente — ruido blanco.