Calculates the credible regions.

def credibleregions(self, probs):
        """ Calculates the credible regions.
        """

        return [brentq(lambda l:  self.pdf[self.pdf > l].sum() - p, 0.0, 1.0) for p in probs]