Lección
La teoría — La media móvil como filtro
Una media móvil de M puntos sustituye cada valor de una serie por la media de M valores consecutivos. Es una convolución con una caja de M pesos de 1/M, así que, por el teorema de convolución, su efecto sobre un ciclo de f ciclos por muestra es una multiplicación por la transformada de Fourier de la caja, la ganancia H(f) = sen(πfM)/(M sen πf). H es cero en f = 1/M, 2/M, …, es negativa entre el primer y el segundo cero, y es real cuando cada media se sitúa en el centro de su ventana; si se sitúan al final de la ventana, los mismos valores conllevan un retraso de (M − 1)/2 muestras.
Qué significa cada símbolo
M- el número de días en cada media
f- la frecuencia de un ciclo, en ciclos por día, de 0 a ½
H(f)- la ganancia: por cuánto multiplica la media a un ciclo de frecuencia f; negativo significa invertido
- Supone
- Un valor por día sin huecos, para que la ventana abarque siempre M días naturales. La cancelación necesita que los propios valores notificados se repitan cada 7 días. Un retraso en la notificación que sea siempre el mismo los mantiene semanales y no cambia nada; un día que falte, un día festivo o un retraso que varíe de una semana a otra la rompen.
Problema resuelto al detalle
-
Promedia la media 4 pasos
Una media de 7 días pone boca abajo una oscilación de 5 días. Un colega sugiere sacar la media de 7 días de la media de 7 días. ¿Qué le pasa ahora a la onda de 5 días, y al ciclo semanal, y cuál es el precio de la segunda pasada? Pulsa primero Una onda de 5 días: la página te da la ganancia de una pasada.
Abre este problema en la herramienta
-
Una pasada, a partir de la página:
-
Una segunda pasada multiplica la onda otra vez por la misma ganancia, y negativo por negativo es positivo:
-
Las dos pasadas juntas son una media sobre 13 días con pesos triangulares, la caja convolucionada consigo misma:
-
Su punto de media potencia se desplaza hacia afuera, de H² = ½ para una pasada a H⁴ = ½ para dos:
Respuesta
La onda de 5 días sale derecha, a 0,0534 de su tamaño, y el ciclo semanal sigue borrado, ya que 0 × 0 = 0. El precio es un suavizador que también atenúa los ciclos más lentos: el periodo cuya potencia se reduce a la mitad pasa de 15,66 días a 21,74, y una versión centrada necesita 6 días a cada lado en lugar de 3, por lo que se detiene 6 días antes que los datos. Hyndman y Athanasopoulos llaman a este tipo de cosas una media móvil de una media móvil, una MM-3×3 en su ejemplo; aquí se trata de una MM-7×7.
-
Ruta de aprendizaje
La transformada de Fourier en la práctica
Referencias (3)
- The source of the page's real series: the daily difference of Germany's cumulative confirmed cases in the JHU CSSE COVID-19 Data Repository, licensed CC BY 4.0, with the dashboard paper its terms of use ask users to cite: COVID-19 Data Repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University, time_series_covid19_confirmed_global.csv (archived 10 March 2023). E. Dong, H. Du and L. Gardner, An interactive web-based dashboard to track COVID-19 in real time. The Lancet Infectious Diseases 20(5), 533–534 (2020). doi:10.1016/S1473-3099(20)30120-1
- The moving average as a convolution with a rectangular pulse, the centred and one-sided forms, and its frequency response H[f] = sin(πfM)/(M sin πf) (equation 15-2), with the verdict the second insight quotes: "an exceptionally good smoothing filter", "an exceptionally bad low-pass filter": S. W. Smith, The Scientist and Engineer's Guide to Digital Signal Processing. California Technical Publishing (online edition, © 1997–2011). Chapter 15, Moving Average Filters.
- The m-MA of a time series, why symmetric moving averages have odd order, the 2×m-MA for an even seasonal period, the statement that "a 7-MA can be used to estimate the trend-cycle of daily data with a weekly seasonality" while other orders leave the estimate contaminated by the seasonality, and moving averages of moving averages (the 3×3-MA), which the worked problem extends to 7×7: R. J. Hyndman and G. Athanasopoulos, Forecasting: Principles and Practice, 3rd edition. OTexts, Melbourne, 2021. Section 3.3, Moving averages.