Median including bayesian credible interval.

def median_interval(data, alpha=_alpha):
    """
    Median including bayesian credible interval.
    """
    q = [100*alpha/2., 50, 100*(1-alpha/2.)]
    lo,med,hi = np.percentile(data,q)
    return interval(med,lo,hi)