¿Cómo manejar el pronóstico de series de tiempo en línea?

He estado lidiando con el siguiente problema. Tengo una especie de sistema de tiempo real y cada período de tiempo leo su valor actual, creando una serie de tiempo (como 1, 12, 2, 3, 5, 9, 1, ...). Me gustaría conocer métodos (estadística y aprendizaje automático) para pronosticar el siguiente valor en línea (es decir, cada vez que se lee un nuevo valor). Probé ARIMA de R y SMOreg de Weka, dan buenas predicciones, pero son algo lentas ya que cada vez que aparece un nuevo valor tengo que rehacer las matemáticas.

PD: Sería genial también si el método tuviera un intervalo de confianza.

En su sistema de tiempo real, ¿son los tiempos de observación no homogéneos y los datos no estacionarios? Si desea algo simple y rápido, sugiero utilizar los operadores de tipo EMA no homogéneos:

Operadores en series de tiempo no homogéneas

Actualizan la EMA ( ) con cada nueva observación de acuerdo con,$\text{value}$

con un parámetro de suavizado / ajuste. Es una forma simple de estimar una expectativa.$\tau$

También se puede crear una estimación mediana simple en línea a través de la actualización

En la práctica, desea pequeño (o en descomposición con más observaciones). Idealmente, debería depender de qué tan asimétricas sean las actualizaciones; es decir, si realidad es igual a la mediana, entonces$\epsilon$$\epsilon$$\text{med}$$\text{sg}$ debe ser uniforme en $\{-1,1\}$. Luego puede extender esto a una profundidad$d$ estructura de tipo de árbol binario equilibrado para obtener $2^{d+1}-1$ cuantiles uniformemente espaciados.

La combinación de lo anterior debería proporcionarle una distribución en línea decente de sus datos. El árbol es difícil de entender, tengo implementaciones de ambos en C ++ si estás interesado. Utilizo mucho ambos en la práctica (datos de tick financieros en tiempo real) y han funcionado bien.

Primero debe hacer una incrustación de sus datos en el tiempo. Por ejemplo, tome como primera entrada [1, 12, 2, 3] y la salida correspondiente [5], y como segunda entrada tome [12, 2, 3, 5] y la salida correspondiente [9]. (Esto está incrustado con el retraso 4, pero puede elegir otro valor que se ajuste mejor).

Ahora tiene un problema de predicción válido. A estos datos puede aplicar los procesos gaussianos en línea . Este es un método de aprendizaje automático que hace exactamente lo que usted describe y proporciona intervalos de confianza.

Si su modelo no es estacionario, puede probar la extensión no estacionaria, el rastreador recursivo de mínimos cuadrados del núcleo . Por cierto, ese documento incluye el código de Matlab para casos estacionarios y no estacionarios.

Estos métodos son razonablemente rápidos: su complejidad computacional es cuadrática en términos de la cantidad de datos que almacena en la memoria (que generalmente es una parte pequeña y representativa de todos los datos procesados). Para métodos más rápidos, recomiendo, por ejemplo, el método de mínimos cuadrados del núcleo, pero su precisión es menor.

El filtro de Kalman es un algoritmo recursivo. Toma la nueva observación y la combina con la predicción previa. Sería bueno usarlo, pero solo si es un modelo apropiado para sus datos. No estoy seguro de lo fácil que es actualizar el intervalo de predicción.

No sé si intentaste esto, pero en R cuando usas la función Arima puedes especificar el modelo como entrada. Entonces, si inicialmente encontró un modelo de arima, digamos Arima (1,2,1) con sus respectivos componentes de suavizado, puede arreglar el modelo en iteraciones posteriores para que no intente reajustar un modelo. Si sus datos son estacionarios en ese caso, las predicciones pueden ser lo suficientemente buenas para usted y mucho más rápidas.

Espero que esto ayude..