def hdmedian(data, axis=(-1), var=False):
 result = hdquantiles(data, [0.5], axis=axis, var=var)
 return result.squeeze()

def hdmedian ( data , axis = ( - 1 ) , var = False ) : result = hdquantiles ( data , [ 0.5 ] , axis = axis , var = var ) return result . squeeze ( )
returns the harrell-davis estimate of the median along the given axis .