我在数据集上编码[23,25,28,32,33,35]
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,35])
print(iqr(x,axis=0))
是什么导致了这个问题？
最佳答案
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,35])
8.0
编辑：在维基百科的talk page上提出,所提出的算法不是确定的,实际上scipy.stats.iqr的方法也是可以接受的.参见确定四分位数Here的三种方法