如何使用Numpy求一个数字在一列数字中的分位数

这个好像没有直接的函数实现,那就自己写吧。因为分位数函数np.percentile()是单调函数,这里使用二分法求解。

import numpy as np


def arg_percentile(series, x):
    # 分位数的启始区间
    a, b = 0, 1
    while True:
        # m是a、b的终点
        m = (a+b)/2
        # 可以打印查看求解过程
        # print(np.percentile(series, 100*m), x)
        if np.percentile(series, 100*m) >= x:
            b = m
        elif np.percentile(series, 100*m) < x:
            a = m
        # 如果区间左右端点足够靠近,则退出循环。
        if np.abs(a-b) <= 0.000001:
            break
    return m

如何使用?

series可以是np.array也可以是pandas.Series,x是待求分位数的数字。举例如下:

s = np.array([1, 3, 4, 4, 6, 7, 8, 9, 10, 11, 13, 15])
x = 12
p = arg_percentile(s, x)
 
print(p)
 
# 结果是 0.8636369705200195

希望对你有帮助。

  • 5
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值