numpy.nanquantile#

numpy.nanquantile(a, q, axis=None, out=None, overwrite_input=False, method='linear', keepdims=<no value>, *, weights=None, interpolation=None)[源代码]#

计算沿指定轴的数据的第 q 个分位数,同时忽略 nan 值.返回数组元素的第 q 个分位数.

在 1.15.0 版本加入.

参数:
aarray_like

可以转换为数组的输入数组或对象,包含要忽略的 nan 值

q类似数组的浮点数

计算分位数所需的概率或概率序列.值必须在0到1之间（包括0和1）.

axis{int, tuple of int, None}, 可选

计算分位数所沿的轴或轴.默认是计算数组展平版本上的分位数.

outndarray, 可选

要在其中放置结果的替代输出数组.它必须具有与预期输出相同的形状和缓冲区长度,但如果需要,输出类型将被强制转换.

overwrite_inputbool, 可选

如果为真,则允许输入数组 a 被中间计算修改,以节省内存.在这种情况下,此函数完成后输入 a 的内容是未定义的.

methodstr, 可选

此参数指定用于估计分位数的方法.有许多不同的方法,其中一些是NumPy特有的.请参阅注释以获取解释.按其在H&F论文[R02de30f409d2-1]_中总结的R类型排序的选项有:

  1. ‘inverted_cdf’

  2. ‘averaged_inverted_cdf’

  3. ‘closest_observation’

  4. ‘interpolated_inverted_cdf’

  5. ‘hazen’

  6. ‘weibull’

  7. ‘linear’ （默认）

  8. ‘中位数无偏’

  9. ‘normal_unbiased’

前三种方法是间断的.NumPy 进一步定义了以下默认 ‘linear’ (7.) 选项的间断变体:

  • ‘lower’

  • ‘更高’,

  • ‘中点’

  • ‘nearest’

在 1.22.0 版本发生变更: 这个参数之前被称为”插值”,并且只提供”线性”默认和最后四个选项.

keepdimsbool, 可选

如果设置为 True,被减少的轴将作为尺寸为1的维度保留在结果中.使用此选项,结果将正确地与原始数组 a 进行广播.

如果这不是默认值,它将被传递（在空数组的特殊情况下）到基础数组的 mean 函数.如果数组是子类,并且 mean 没有 keepdims kwarg,这将引发 RuntimeError.

weights类似数组, 可选

a 中值相关的一组权重.`a` 中的每个值根据其关联的权重对分位数做出贡献.权重数组可以是 1-D（在这种情况下,其长度必须与 a 沿给定轴的大小相同）或与 a 的形状相同.如果 weights=None,则假定 a 中的所有数据具有等于一的权重.只有 method=”inverted_cdf” 支持权重.

在 2.0.0 版本加入.

interpolationstr, 可选

该方法的关键字参数的弃用名称.

自 1.22.0 版本弃用.

返回:
quantile标量或ndarray

如果 q 是一个单一的概率并且 axis=None,那么结果是一个标量.如果给出了多个概率水平,结果的第一个轴对应于分位数.其他的轴是 a 减少后剩余的轴.如果输入包含小于 float64 的整数或浮点数,则输出数据类型是 float64.否则,输出数据类型与输入相同.如果指定了 out,则返回该数组.

参见

quantile
nanmean, nanmedian
nanmedian

等同于 nanquantile(..., 0.5)

nanpercentile

与 nanquantile 相同,但 q 的范围在 [0, 100] 之间.

备注

numpy.nanquantile 的行为与 numpy.quantile 相同（忽略 nan 值）.更多信息请参见 numpy.quantile.

参考文献

[1]

R. J. Hyndman and Y. Fan, “Sample quantiles in statistical packages,” The American Statistician, 50(4), pp. 361-365, 1996

示例

>>> import numpy as np
>>> a = np.array([[10., 7., 4.], [3., 2., 1.]])
>>> a[0][1] = np.nan
>>> a
array([[10.,  nan,   4.],
      [ 3.,   2.,   1.]])
>>> np.quantile(a, 0.5)
np.float64(nan)
>>> np.nanquantile(a, 0.5)
3.0
>>> np.nanquantile(a, 0.5, axis=0)
array([6.5, 2. , 2.5])
>>> np.nanquantile(a, 0.5, axis=1, keepdims=True)
array([[7.],
       [2.]])
>>> m = np.nanquantile(a, 0.5, axis=0)
>>> out = np.zeros_like(m)
>>> np.nanquantile(a, 0.5, axis=0, out=out)
array([6.5, 2. , 2.5])
>>> m
array([6.5,  2. ,  2.5])
>>> b = a.copy()
>>> np.nanquantile(b, 0.5, axis=1, overwrite_input=True)
array([7., 2.])
>>> assert not np.all(a==b)
On this page