ADTK básico

Actualizado 2020-01-29 Lectura 2 min

Resumen

Organizar objetivo, supuestos y condiciones de uso del metodo.
Revisar como reglas de actualizacion o criterios de division afectan el comportamiento.
Usar ejemplos de implementacion para concretar decisiones de ajuste de parametros.

Intuicion #

Vamos a realizar detección de anomalías utilizando Anomaly Detection Toolkit (ADTK). Los datos originales provienen de Numenta Anomaly Benchmark.

Explicacion Detallada #

Vamos a realizar detección de anomalías utilizando Anomaly Detection Toolkit (ADTK). Los datos originales provienen de Numenta Anomaly Benchmark.

1
2
3
4
5
6
7
8
9
import pandas as pd

s_train = pd.read_csv(
    "./training.csv", index_col="timestamp", parse_dates=True, squeeze=True
)
from adtk.data import validate_series

s_train = validate_series(s_train)
print(s_train)

timestamp
2014-04-01 00:00:00    18.090486
2014-04-01 00:05:00    20.359843
2014-04-01 00:10:00    21.105470
2014-04-01 00:15:00    21.151585
2014-04-01 00:20:00    18.137141
                         ...    
2014-04-14 23:35:00    18.269290
2014-04-14 23:40:00    19.087351
2014-04-14 23:45:00    19.594689
2014-04-14 23:50:00    19.767817
2014-04-14 23:55:00    20.479156
Freq: 5T, Name: value, Length: 4032, dtype: float64


  s_train = pd.read_csv("./training.csv", index_col="timestamp", parse_dates=True, squeeze=True)

1
2
3
from adtk.visualization import plot

plot(s_train)

csv", index_col=“timestamp”, parse_dates=True, squeeze=True) (figura)

Comparación de métodos de detección de anomalías #

Realizaremos la detección de anomalías utilizando SeasonalAD. Para otros métodos, consulte Detector.

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
import matplotlib.pyplot as plt
from adtk.detector import (
    AutoregressionAD,
    InterQuartileRangeAD,
    LevelShiftAD,
    PersistAD,
    SeasonalAD,
)

model_dict = {
    "LevelShiftAD": LevelShiftAD(window=5),
    "SeasonalAD": SeasonalAD(),
    "PersistAD": PersistAD(c=3.0, side="positive"),
    "InterQuartileRangeAD": InterQuartileRangeAD(c=1.5),
    "AutoregressionAD": AutoregressionAD(n_steps=7 * 2, step_size=24, c=3.0),
}

for model_name, model in model_dict.items():
    anomalies = model.fit_detect(s_train)
    plot(s_train, anomaly=anomalies, anomaly_color="red", anomaly_tag="marker")
    plt.title(model_name)
    plt.show()

Para otros métodos, consulte Detector (figura)