时间序列的离群点检测

2024-04-19 02:09:53 发布

您现在位置:Python中文网/ 问答频道 /正文

我想生成一个非常简单的时间序列异常检测示例。所以我创建了一个非常明显的异常值的样本数据,但是到目前为止,我还没有得到任何方法来可靠地检测异常值。我尝试了局部离群因子,隔离森林和k最近的邻居。从我所读的内容来看,至少其中一种方法应该是合适的。我也尝试了调整参数,但这并没有真正的帮助。

我在这里犯了什么错误?这些方法不合适吗?

下面是一个代码示例。

提前谢谢!

import numpy as np
import matplotlib.pyplot as plt

np.random.seed(1)

t=np.linspace(0,10,101).reshape(-1,1)
y_test=0.5+t+t**2+2*np.random.randn(len(t),1)

y_test[10]=y_test[10]*7

plt.figure(1)
plt.plot(t,y_test)
plt.show;

from sklearn.neighbors import LocalOutlierFactor

clf=LocalOutlierFactor(contamination='auto')
pred=clf.fit_predict(y_test)

plt.figure(2)
plt.plot(t[pred==1],y_test[pred==1],'bx')
plt.plot(t[pred==-1],y_test[pred==-1],'ro')
plt.show

from sklearn.ensemble import IsolationForest

clf=IsolationForest(behaviour='new',contamination='auto')
pred=clf.fit_predict(y_test)

plt.figure(3)
plt.plot(t[pred==1],y_test[pred==1],'bx')
plt.plot(t[pred==-1],y_test[pred==-1],'ro')
plt.show

from pyod.models.knn import KNN

clf = KNN()
clf.fit(y_test)
pred=clf.predict(y_test)

plt.figure(4)
plt.plot(t[pred==0],y_test[pred==0],'bx')
plt.plot(t[pred==1],y_test[pred==1],'ro')
plt.show

Tags: 方法fromtestimportroplotshownp