ENH: percentile function with additional parameters and vecorization

The percentile function was enhanced by adding limit and interpolation
parameters to give it similar functionality to SciPy's stats.scoreatpercentile
function.  In addition the function was vecorized along q and rewritten to
use the partition method for better performance.

This pull request is in replacement of PR #2970