SciPy

scipy.stats.mstats.hdmedian

scipy.stats.mstats.hdmedian(data, axis=-1, var=False)[source]

Returns the Harrell-Davis estimate of the median along the given axis.

Parameters:
data : ndarray

Data array.

axis : int, optional

Axis along which to compute the quantiles. If None, use a flattened array.

var : bool, optional

Whether to return the variance of the estimate.

Returns:
hdmedian : MaskedArray

The median values. If var=True, the variance is returned inside the masked array. E.g. for a 1-D array the shape change from (1,) to (2,).