LESMTH - Suavizado exponencial lineal de Brown

Devuelve el estimado de pronóstico fuera de la muestra del suavizado exponencial lineal (aka Doble Brown).

 

Sintaxis

LESMTH(X, Orden, Alfa, Optimizar, T, Tipo de retorno)

X es la serie univariada de tiempo de datos (un despliegue de celdas unidimensional (e.g. fila o columna).

Orden es el orden del tiempo en la serie de datos (i.e. la primera fecha correspondiente del punto de datos (la fecha más reciente=1 (aleatoria), la fecha más tardía=0)).

Orden Descripción
1 ascendente (el primer punto de datos corresponde a la fecha más temprana) (defecto)
0 descendente (el primer punto de datos corresponde a la ultima fecha)

Alfa es el factor del nivel de suavizado (alfa debe estar entre cero(0) y uno(1) (exclusivamente)). Si este valor faltara o se omitiera, se usaría 0.333.

Optimizar es una bandera de (Veradero/Falso) para buscar y usar el valor óptimo del factor de suavizado. Si el valor faltara o se omitiera, la optimización se asumiría como falsa.

T es el estimado tiempo-horizonte más allá del final de X. Si faltara, un valor estándar de 0 (más reciente o al final de X) se asumirá.

Tipo de retornoes un número que determina el tipo de valor de retorno: 0 (o faltante) = Pronóstico, 1=Alfa, 2=nivel de componente (series), 3=componente de tendencia (series), 4=pronósticos de un paso (series).

TIPO DE RETORNO NUMERO RETORNADO O DEVUELTO
0 or omitted Pronóstico
1 el factor del nivel de suavizado (Alfa)
2 Nivel de componente (series)
3 Componente de tendencia (series)
4 Pronósticos de un paso (series)
 

Observaciones

  1. La serie de tiempo es homogenea o igualmente espaciada.
  2. La serie de tiempo puede incluir los valores faltantes (e.g. #N/A) sólo en cada extremo.
  3. El suavizado exponencial lineal se aplica mejor a las series de tiempo que exhiben una tendencia prevalente pero no muestran estaciones.
  4. La forma recursiva de la ecuación de suavizado exponencial de Brown se expresa así:

    $$ \begin{array}{l} S_{t\succ 1}^{'} = \alpha\times X_t+(1-\alpha)\times S_{t-1}^{'} \\ S_{t\succ 1}^{''} = \alpha\times S_t^{'}+(1-\alpha)\times S_{t-1}^{''} \\ S_t = 2\times S_{t}^{'} - S_{t}^{''} \\ b_t = \frac{\alpha}{1-\alpha} (S_{t}^{'} - S_{t}^{''})\\ \\ \hat{F}_t(m)=S_t + m\times b_t \end{array} $$
    Donde:
    • $X_t$ es el valor de las series de tiempo en el tiempo t.
    • $S_{t}^{'}$ es el valor de suavizado exponencial simple de X.
    • $S_{t}^{''}$ es el doble valor del suavizado exponencial.
    • $S_{t}$ es el estimado del nivel.
    • $b_{t}$ es el estimado de la pendiente.
    • $\alpha$ es el coeficiente de suavizado por nivel.
    • $\hat{F}_t(m)$ es el valor pronóstico de suavizado en el paso adelantado m para X en el tiempo t.
  5. 5. Para LESMTH estamos aplicando el suavizado exponencial dos veces y calculando el nivel y la pendiente para esas dos series:
    • Usar $X_t$ como entrada y $\alpha$ como coeficiente de suavizado, genera $S_{t}^{'}$.
    • Usar $S_{t}^{'}$ como entrada y $\alpha$ como coeficiente de suavizado, genera $S_{t}^{''}$.
  6. Para $\alpha = 1$, el suavizado exponencial doble es equivalente al método de extrapolación inmodificable primitiva (NCE), o simplemente una caminata aleatoria. En este caso, la división por cero ocurre cuando calculamos la pendiente, y el pronóstico es indefinido.
  7. Para $\alpha = 0$, el pronóstico será una constante tomando sus valores desde el valor inicial para $S_{t}^{'}$ y $S_{t}^{''}$:
    $$ \hat{F}_t(m) = 2\times S_{1}^{'} - S_{1}^{''} $$
  8. LESMTH calcula el pronóstico por puntos. No hay un modelo probalístico asumido para el suavizado exponencial simple, entonces no podemos derivar un intervalo estadístico confiable para los valores calculados.
  9. En la práctica, el error cuadrático medio (MSE, por sus siglas en inglés) para valores anteriores de pronóstico fuera de la muestra se usa, con mayor frecuencia, como indicador de incertidumbre (i.e. varianza) en el valor de pronóstico más reciente.
  10. Este método requiere dos valores de inicio ($S_{1}^{'},S_{1}^{''}$) para arrancar con la actualización recursiva de la ecuación. En NumXL, hemos agrupado estos valores tal como se muestra a continuación:
    • $S_{1}^{'}$ está establecido con la media de las primeras cuatro oservaciones, y para una serie de tiempo corta se ha establecido como el valor de la primera observación.

      $$ S_{1}^{'}=\left\{\begin{array}{l} X_1\\ \frac{\sum_{t=1}^4 X_t}{4} \end{array}\right. \begin{array}{r} N \leq 4\\ N \gt 4 \end{array} $$
    • $S_{1}^{''}$ está establecido como la media de las primeras cuatro observaciones de $S_{t}^{'}$, y para una serie de tiempo corta se ha establecido como el valor de la primera observación.

      $$ S_{1}^{''}=\left\{\begin{array}{l} S_{1}^{'}\\ \frac{\sum_{t=1}^4 S_{t}^{'}}{4} \end{array}\right. \begin{array}{r} N \leq 4\\ N \gt 4 \end{array} $$
  11. Empezando desde la version 1.63 de NumXL, SESMTH tiene un optimizador incorporado para encontrar el mejor valor de $\alpha$ que minimice el SSE (función de pérdida ($U(.)$) para el pronóstico de un paso calculado dentro de la muestra.

    $$ \begin{array}{l} U(\alpha)=\mathrm{SSE}=\sum_{t=1}^{N-1}(X_{t+1}-\hat{F}_t(1))^2\\ \min_{\alpha \in (0,1)} U(\alpha) \end{array} $$
  12. Para los valores iniciales, el optimizador de NumXL utilizará el valor de registro de alfa (si está disponible) en el problema de minimización, y los valores iniciales para las dos series de suavizado ($S_{1}^{'},S_{1}^{''}$) se calculan a partir de los datos de registro.
  13. Empezando con la version 1.65 de NumXL, la función SESMTH devuelve el valor óptimo encontrado a alfa, y a las series de suavizado de un paso correspondientes de nivel, de tendencia y al pronóstico calculado dentro de la muestra.
  14. Las series de tiempo deben tener al menos tres (3) observaciones con valores no faltantes para usar el optimizador integrado.
  15. NumXL implementa el método de pendiente espectral proyectada (SPG) para encontrar el mínimo con una frontera dentro del cuadro.
    • El SPG requiere un valor de la función de pérdida y el primer derivativo. NumXL implementa la fórmula exacta derivativa (vs. la aproximación numérica) para propósitos de funcionamiento.

      $$ \begin{array}{l} \frac{\partial U}{\partial \alpha}=-2\times\sum_{t=1}^{N-1}(X_{t+1}-\hat{F}_t(1))\times \frac{\partial \hat{F}_t}{\partial \alpha} \\ \\ \frac{\partial \hat{F}_t}{\partial \alpha}=\frac{\partial S_t}{\partial \alpha}+\frac{\partial b_t}{\partial \alpha}\\ \\ \frac{\partial S_t}{\partial \alpha}= 2\times \frac{\partial S_t^{'}}{\partial \alpha}-\frac{\partial S_t^{''}}{\partial \alpha}\\ \frac{\partial b_t}{\partial \alpha} = \frac{S_t^{'}-S_t^{''}}{(1-\alpha)^2}+\frac{\alpha}{1-\alpha}\times(\frac{\partial S_t^{'}}{\partial \alpha}-\frac{\partial S_t^{''}}{\partial \alpha}) \\ \frac{\partial S_t^{'}}{\partial \alpha}=X_{t+1}+(1-\alpha)\times \frac{\partial S_{t-1}^{'}}{\partial \alpha}-S_{t-1}^{'}\\ \frac{\partial S_t^{''}}{\partial \alpha}=S_t^{'}+\alpha\frac{\partial S_t^{'}}{\partial \alpha} +(1-\alpha)\times \frac{\partial S_{t-1}^{''}}{\partial \alpha}-S_{t-1}^{''}\\ \\ \end{array} $$
    • Internamente, durante la optimización, NumXL calcula de manera recursiva tanto el suavizado de series de tiempo, de nivel y de tendencia como los derivados dentro de la muestra, los cuales se usan para la función de pérdida y sus derivadas.
    • El SPG es un método reiterativo (recursivo), y es posible que el mínimo no se pueda encontrar en el número de reiteraciones permitidas y/o toleradas. En este caso, NumXL no fracasará, en vez, NumXL usará el mejor valor de alfa encontrado hasta el momento.
    • El SPG no tiene provisiones para detectar o evadir la trampa minima local. No existe garantía del mínimo global.
  16. En general, la función SSE en LESMTH cede el paso a una curva convexa monótona de suavizado continuo, ese minimizador de SPG casi siempre encuentra una solución óptima en muy pocas reiteraciones.

Ejemplos

Ejemplo 1:

 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
A B C
Date Data LESMTH
Enero 10, 2008 -0.30 #N/A
January 11, 2008 -1.28 -0.3
January 12, 2008 0.24 -0.89
January 13, 2008 1.28 -0.30
January 14, 2008 1.20 0.66
January 15, 2008 1.73 1.14
January 16, 2008 -2.18 1.70
January 17, 2008 -0.23 -0.37
January 18, 2008 1.10 -0.38
January 19, 2008 -1.09 0.43
January 20, 2008 -0.69 -0.43
January 21, 2008 -1.69 -0.67
January 22, 2008 -1.85 -1.39
January 23, 2008 -0.98 -1.86
January 24, 2008 -0.77 -1.57
January 25, 2008 -0.30 -1.25
January 26, 2008 -1.28 -0.77
January 27, 2008 0.24 -1.08
January 28, 2008 1.28 -0.34
January 29, 2008 1.20 0.70
January 30, 2008 1.73 1.22
January 31, 2008 -2.18 1.79
Febrero 1, 2008 -0.23 -0.29
Febrero 2, 2008 1.10 -0.30
Febrero 3, 2008 -1.09 0.49
Febrero 4, 2008 -0.69 -0.38
Febrero 5, 2008 -1.69 -0.63
Febrero 6, 2008 -1.85 -1.35
Febrero 7, 2008 -0.98 -1.84

Ejemplos de archivos

Referencias

¿Tiene más preguntas? Enviar una solicitud

0 Comentarios