def hdmedian(data, axis=-<NUM_LIT:1>, var=False):
result = hdquantiles(data,[<NUM_LIT:0.5>], axis=axis, var=var)<EOL>return result.squeeze()<EOL>
Returns the Harrell-Davis estimate of the median along the given axis.

Parameters
----------
data : ndarray
    Data array.
axis : int
    Axis along which to compute the quantiles. If None, use a flattened
    array.
var : boolean
    Whether to return the variance of the estimate.