scipy.stats.mstats.

hdquantiles#

scipy.stats.mstats.hdquantiles(data, prob=[0.25, 0.5, 0.75], axis=None, var=False)[source]#

使用 Harrell-Davis 方法计算分位数估计。

分位数估计被计算为顺序统计量的加权线性组合。

参数:

data类数组: 数据数组。
prob序列，可选: 用于计算分位数的概率序列。
axis整数或 None，可选: 计算分位数的轴。如果为 None，则使用扁平化数组。
var布尔值，可选: 是否返回估计值的方差。

返回:

hdquantiles掩码数组: 一个 (p,) 数组的分位数（如果 var 为 False），或者一个 (2,p) 数组的分位数和方差（如果 var 为 True），其中 p 是分位数的数量。

另请参见

hdquantiles_sd

示例

>>> import numpy as np
>>> from scipy.stats.mstats import hdquantiles
>>>
>>> # Sample data
>>> data = np.array([1.2, 2.5, 3.7, 4.0, 5.1, 6.3, 7.0, 8.2, 9.4])
>>>
>>> # Probabilities at which to compute quantiles
>>> probabilities = [0.25, 0.5, 0.75]
>>>
>>> # Compute Harrell-Davis quantile estimates
>>> quantile_estimates = hdquantiles(data, prob=probabilities)
>>>
>>> # Display the quantile estimates
>>> for i, quantile in enumerate(probabilities):
...     print(f"{int(quantile * 100)}th percentile: {quantile_estimates[i]}")
25th percentile: 3.1505820231763066 # may vary
50th percentile: 5.194344084883956
75th percentile: 7.430626414674935