python - scipy.stats 是否对 iqr 进行了错误的计算？

Question

我正在对数据集 [23,25,28,28,32,33,35] 进行编码

根据 wiki和 scipy doc

IQR = Q3 - Q1 = 33 - 25 = 8

当我在数据集上运行 IQR 时，结果 (6) 与预期的结果 (8) 不同。

我在 https://stackoverflow.com/a/23229224 中尝试了另一种方法, 结果为 6.

这是我的代码

import numpy as np
from scipy.stats import iqr
x = np.array([23,25,28,28,32,33,35])
print(iqr(x, axis=0))

是什么导致了这个问题？

Answer 1

scipy.stats.iqr 似乎不遵循维基百科中记录的递归算法。相反，它只是执行 np.percentile(x, 75) - np.percentile(x, 25) 这不排除中位数，它包含在内，所以你得到 (32 + 33 )/2 - (25 + 28)/2 = 6

如果您想在维基百科中使用该算法，您需要执行以下操作:

def iqr_(m):
    m = np.array(m)
    n = m.size//2
    m_ = np.partition(m.ravel(), n + 1)
    return np.median(m_[n + m.size%2:]) - np.median(m_[:n])

iqr_([23,25,28,28,32,33,35])
8.0

编辑:在 talk page 上of wikipedia 提出所提供的算法不是确定的，实际上 scipy.stats.iqr 的方法也是可以接受的。参见确定四分位数的三种方法Here