基于中值绝对偏差 (MAD) 的异常值检测

Question

我想使用@Joe Kington 的答案应用基于中值绝对偏差 (MAD) 的异常值检测，如下所示：

一维观测数据中检测异常值的 Pythonic 方法

但是，我的代码出了什么问题，我不知道如何将异常值指定为 MY DATA 的 nan 值：

import numpy as np
data = np.array([55,32,4,5,6,7,8,9,11,0,2,1,3,4,5,6,7,8,25,25,25,25,10,11,12,25,26,27,28],dtype=float)
median = np.median(data, axis=0)
diff = np.sum((data - median)**2, axis=-1)
diff = np.sqrt(diff)
med_abs_deviation = np.median(diff)
modified_z_score = 0.6745 * diff / med_abs_deviation
data_without_outliers = data[modified_z_score < 3.5]
?????
print data_without_outliers

Answer 1

使用有什么问题：

data[modified_z_score > 3.5] = np.nan

data请注意，这仅在是浮点数组时才有效（如果您正在计算 MAD，则应该是浮点数组）。