如何利用Python NumPy快速计算数据的分位数?

云静大大_5786

云静大大_5786

2026-07-02

746人浏览

原创

优先用 np.quantile 而非 np.percentile,因其参数 q 为 0.0–1.0 分位数,语义清晰、贴合统计直觉且与 sklearn/scipy 一致;需显式指定 interpolation 控制插值,用 keepdims 避免维度压缩,含 nan 时须改用 np.nanquantile。

如何利用python numpy快速计算数据的分位数?

用 np.quantile 直接算,别用 np.percentile 混淆语义

两者功能几乎一样,但参数单位不同:np.quantile 接收 0–1 区间的分位数(比如中位数是 0.5),np.percentile 接收 0–100 的百分位(中位数是 50)。实际项目里混用容易写错,尤其当逻辑要和统计教材或下游系统对齐时。建议统一用 np.quantile,语义更清晰。

常见错误:传入 q=95 给 np.quantile,结果远超数据范围——其实是把百分位当成分位数用了。

  • 正确写法:np.quantile(data, 0.95)(95% 分位)
  • 错误写法:np.quantile(data, 95)(会当作索引越界处理,可能报 IndexError 或返回意外值)
  • 若必须用 np.percentile,记得除以 100:np.percentile(data, 95)

多维数组分位数计算必须指定 axis,否则默认展平

这是最常踩的坑:对二维数组直接调用 np.quantile(data, 0.5),你以为是在每列上算中位数,结果它把整个矩阵拉成一维再算——得到的是全局中位数,不是你想要的列/行统计。

典型场景:你想知道每只股票每日收益率的 25% 分位数(即每列一个值),就得明确告诉 NumPy 按哪一维操作。

huaweicloud-python-sdk-helper
huaweicloud-python-sdk-helper

华为云 Python SDK 查询工具。自动发现已安装 SDK 的服务包、Client类、API 方法、入参详情、类信息(__init__ 参数、属性、类型、描述),支持关键字搜索,快速定位目标 API 使用方法。

下载
  • 按列(每列一个分位数):np.quantile(data, 0.25, axis=0)
  • 按行(每行一个分位数):np.quantile(data, 0.25, axis=1)
  • 不指定 axis → 等价于 axis=None → 全局计算
  • 如果想保留原维度形状(比如输出和输入同 shape),加 keepdims=True

interpolation 参数影响结果精度,尤其小样本时不能忽略

当分位点落在两个数据点之间时,NumPy 需插值。默认插值方式是 'linear',但还有 'lower'、'higher'、'midpoint'、'nearest' 可选。不同业务对“保守估计”或“向上取整”有硬性要求(比如风控阈值必须不低于某分位),这时必须显式指定。

例如:某监控指标要求告警阈值取“不低于 99% 分位的最小值”,就该用 interpolation='higher';若做学术报告需严格匹配 R 或 SPSS 的默认行为,得查清它们用的是哪种插值。

  • 小样本(
  • 'linear' 是线性插值(最常用),'nearest' 返回最近的原始观测值
  • 注意:NumPy 1.22+ 已弃用 interpolation 参数名,改用 method,但旧版本仍认 interpolation;建议检查 np.__version__ 后适配

大数据量时优先用 np.quantile 而非排序后索引取值

有人图省事手写 sorted_data[int(0.95 * len(data))],这在小数据上可行,但问题不少:没处理插值、不支持多维、无法指定 axis、且对大数组排序开销远高于 NumPy 内置算法。NumPy 的 quantile 实际用的是 introselect 类算法,平均时间复杂度 O(n),比全排序 O(n log n) 更快。

  • 100 万点数组,np.quantile 通常比 np.sort + 索引快 3–5 倍
  • 内存友好:不需要额外存完整排序数组
  • 若只需单个分位数,别用 np.percentile 传数组如 [25, 50, 75]——它内部会重复扫描,不如循环调用单次 np.quantile

分位数计算看着简单,但 axis、插值方式、参数单位这三个点,任一个设错都会让结果偏离预期。特别是生产环境跑批任务时,最好在测试数据上打印中间结果,确认维度和数值都符合业务定义。

Python免费学习笔记(深入):立即使用
在学习笔记中,你将探索 Python 的核心概念和高级技巧!

相关文章

PHP速学视频免费教程(入门到精通)
PHP速学视频免费教程(入门到精通)

PHP怎么学习?PHP怎么入门?PHP在哪学?PHP怎么学才快?不用担心,这里为大家提供了PHP速学教程(入门到精通),有需要的小伙伴保存下载就能学习啦!

下载

相关标签:

python

本站声明:本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系admin@php.cn

相关专题

更多
python打包成可执行文件
python打包成可执行文件

本专题为大家带来python打包成可执行文件相关的文章,大家可以免费的下载体验。

2023.07.20

1771

4

python能做什么
python能做什么

python能做的有:可用于开发基于控制台的应用程序、多媒体部分开发、用于开发基于Web的应用程序、使用python处理数据、系统编程等等。本专题为大家提供python相关的各种文章、以及下载和课程。

2023.07.25

4604

7

format在python中的用法
format在python中的用法

Python中的format是一种字符串格式化方法,用于将变量或值插入到字符串中的占位符位置。通过format方法,我们可以动态地构建字符串,使其包含不同值。php中文网给大家带来了相关的教程以及文章,欢迎大家前来阅读学习。

2023.07.31

1729

3

python教程
python教程

Python已成为一门网红语言,即使是在非编程开发者当中,也掀起了一股学习的热潮。本专题为大家带来python教程的相关文章,大家可以免费体验学习。

2023.08.03

27237

23

python环境变量的配置
python环境变量的配置

Python是一种流行的编程语言,被广泛用于软件开发、数据分析和科学计算等领域。在安装Python之后,我们需要配置环境变量,以便在任何位置都能够访问Python的可执行文件。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3287

5

python eval
python eval

eval函数是Python中一个非常强大的函数,它可以将字符串作为Python代码进行执行,实现动态编程的效果。然而,由于其潜在的安全风险和性能问题,需要谨慎使用。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.04

3307

5

scratch和python区别
scratch和python区别

scratch和python的区别:1、scratch是一种专为初学者设计的图形化编程语言,python是一种文本编程语言;2、scratch使用的是基于积木的编程语法,python采用更加传统的文本编程语法等等。本专题为大家提供scratch和python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

1203

5

python合并两个列表
python合并两个列表

Python是一种强大的编程语言,具有许多方便的功能和工具。在Python中,有多种方法可以合并两个列表。php中文网给大家带来了相关的教程以及文章,欢迎大家前来学习阅读。

2023.08.10

596

4

python是前端还是后端
python是前端还是后端

Python属于前端也属于后端,其灵活性和丰富的生态系统使得开发人员能够在不同的领域中灵活运用。本专题为大家提供python相关的文章、下载、课程内容,供大家免费下载体验。

2023.08.11

2563

5

热门下载

更多
网站特效
/
网站源码
/
网站素材
/
前端模板

精品课程

更多
相关推荐
/
热门推荐
/
最新课程
Uvicorn官方安装与入门手册
Uvicorn官方安装与入门手册

共0课时 | 0人学习

Python uv官方项目管理指南
Python uv官方项目管理指南

共0课时 | 0人学习

Python uv官方安装与入门文档
Python uv官方安装与入门文档

共0课时 | 0人学习