[SciPy-User] scoreatpercentile behaviour

[Andreas Hilboll]

I just had a quick look into scipy.stats.scoreatpercentile, and was
disappointed to see that it's currently not possible to do the
calculation for more than one percentile at a time (``per`` is scalar).
So I had a quick look into the sources, and was surprised to see that
apprently, the function expects ordered input ``a``, which is not noted
in the docstring. (Or maybe it's just my misunderstanding of the word
'percentile'. I had expected the function to work on the input's
**values**, not on the indices.

Is this a bug or a feature? If it's a feature, this should be very
explicitly noted in the docstring, I think. I'm willing to do so if you
can confirm that the current behaviour is actually wanted.

In the sources' TODO, it's stated that a more general percentile
implementation would be welcome. I might be able to contribute something
here; any hints on where to start?

Andreas.