# Licensed under a 3-clause BSD style license - see LICENSE.rst
# -*- coding: utf-8 -*-
"""
==============
desiutil.stats
==============
Just contains a dead-simple wrapper on :func:`numpy.percentile`.
"""
[docs]def perc(x, per=68.2):
"""Calculate the percentile bounds of a distribution,
*i.e.* for per=68, the code returns the upper and lower bounds
that encompass 68 percent of the distribution.
Uses simple interpolation.
Parameters
----------
x : :class:`numpy.ndarray`
numpy array of values
per : :class:`float`, optional
Percentile for the calculation [0-100].
Returns
-------
:class:`numpy.ndarray`
Value at lower, value at upper.
"""
from numpy import percentile
return percentile(x, [50-per/2.0, 50+per/2.0])