python - numpy quantile() 中的线性插值

Question

考虑以下代码：

>>> import numpy as np
>>> l = [21,22,24,24,26,97]
>>> np.quantile(l, 0.25)
22.5

文档说：

线性：i + (j - i) * fraction，其中 fraction 是由 i 和 j 包围的索引的小数部分。

谁能解释一下是什么i，在这个例子中j以及fraction我们是如何得到的22.5？

Answer 1

它的工作原理大致如下：

import numpy as np


def my_quantile(array, q):
    n = len(array)
    index = (n - 1) * q
    if int(index) == index:  # has no fractional part
        return array[index]
    fraction = index - int(index)
    left = int(index)
    right = left + 1
    i, j = array[left], array[right]
    return i + (j - i) * fraction


if __name__ == '__main__':
    arr = [21, 22, 24, 24, 26, 97]
    assert my_quantile(arr, 0.33) == np.quantile(arr, 0.33)
    assert my_quantile(arr, 0.25) == np.quantile(arr, 0.25)
    assert my_quantile(arr, 0.75) == np.quantile(arr, 0.75)