numpy.nanpercentile#
- numpy.nanpercentile(a, q, axis=None, out=None, overwrite_input=False, method='linear', keepdims=<no value>, *, weights=None, interpolation=None)[source]#
计算数据沿指定轴的第 q 个百分位数,同时忽略 NaN 值。
返回数组元素的第 q 个百分位数。
版本 1.9.0 中的新增功能。
- 参数::
- aarray_like
输入数组或可转换为数组的对象,包含要忽略的 NaN 值。
- qfloat 类型的 array_like
要计算的百分位数或百分位数序列,必须介于 0 到 100(含)之间。
- axis{int, int 元组, None}, 可选
计算百分位数的轴或轴。默认情况下,沿数组的扁平化版本计算百分位数。
- outndarray, 可选
放置结果的备用输出数组。它必须与预期输出具有相同的形状和缓冲区长度,但类型(输出的类型)将在必要时被强制转换。
- overwrite_inputbool, 可选
如果为 True,则允许修改输入数组 a,以便在中间计算时节省内存。在这种情况下,此函数完成后输入 a 的内容将是不确定的。
- methodstr, 可选
此参数指定用于估计百分位数的方法。有许多不同的方法,其中一些方法是 NumPy 独有的。有关解释,请参阅注释。按其在 H&F 论文 [1] 中总结的 R 类型排序的选项是
‘inverted_cdf’
‘averaged_inverted_cdf’
‘closest_observation’
‘interpolated_inverted_cdf’
‘hazen’
‘weibull’
‘linear’ (默认)
‘median_unbiased’
‘normal_unbiased’
前三种方法是不连续的。NumPy 进一步定义了默认 ‘linear’ (7.) 选项的以下不连续变体
‘lower’
‘higher’,
‘midpoint’
‘nearest’
更改于版本 1.22.0: 此参数以前称为“interpolation”,并且只提供“linear”默认值和最后四个选项。
- keepdimsbool, 可选
如果将其设置为 True,则缩减的轴将作为大小为一的维度保留在结果中。使用此选项,结果将正确广播到原始数组 a。
如果这与默认值不同,它将被传递给(在空数组的特殊情况下)基础数组的
mean
函数。如果数组是子类并且mean
没有 kwarg keepdims,这将引发 RuntimeError。- weightsarray_like, 可选
与 a 中的值关联的权重数组。 a 中的每个值根据其关联的权重对百分位数做出贡献。权重数组可以是 1-D(在这种情况下,其长度必须是给定轴上 a 的大小)或与 a 形状相同。如果 weights=None,则假定 a 中的所有数据都具有等于 1 的权重。只有 method=”inverted_cdf” 支持权重。
版本 2.0.0 中的新增功能。
- interpolationstr, 可选
method 关键字参数的已弃用名称。
自版本 1.22.0 起已弃用。
- 返回值::
- percentile标量或 ndarray
如果 q 是单个百分位数并且 axis=None,则结果为标量。如果给出多个百分位数,则结果的第一轴对应于百分位数。其他轴是 a 缩减后剩余的轴。如果输入包含整数或小于
float64
的浮点数,则输出数据类型为float64
。否则,输出数据类型与输入相同。如果指定了 out,则返回该数组。
另请参阅
nanmean
nanmedian
相当于
nanpercentile(..., 50)
percentile
,median
,mean
nanquantile
等效于 nanpercentile,除了 q 在 [0, 1] 范围内。
注释
具有百分比 q 的
numpy.nanpercentile
的行为与具有参数q/100
的numpy.quantile
相同(忽略 NaN 值)。有关更多信息,请参阅numpy.quantile
。参考文献
[1]R. J. Hyndman 和 Y. Fan,“统计软件包中的样本分位数”,美国统计学家,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.percentile(a, 50) np.float64(nan) >>> np.nanpercentile(a, 50) 3.0 >>> np.nanpercentile(a, 50, axis=0) array([6.5, 2. , 2.5]) >>> np.nanpercentile(a, 50, axis=1, keepdims=True) array([[7.], [2.]]) >>> m = np.nanpercentile(a, 50, axis=0) >>> out = np.zeros_like(m) >>> np.nanpercentile(a, 50, axis=0, out=out) array([6.5, 2. , 2.5]) >>> m array([6.5, 2. , 2.5])
>>> b = a.copy() >>> np.nanpercentile(b, 50, axis=1, overwrite_input=True) array([7., 2.]) >>> assert not np.all(a==b)