def _false_positive_probability(threshold, b, r, xq): ''' Compute the false positive probability given the containment threshold. xq is the ratio of x/q. ''' _probability = lambda t: 1 - (1 - (t / (1 + xq - t))**float(r))**float(b) if xq >= threshold: a, err = integrate(_probability, 0.0, threshold) return a a, err = integrate(_probability, 0.0, xq) return a
def _false_negative_probability(threshold, b, r, xq): ''' Compute the false negative probability given the containment threshold ''' _probability = lambda t : 1 - (1 - (1 - (t/(1 + xq - t))**float(r))**float(b)) if xq >= 1.0: a, err = integrate(_probability, threshold, 1.0) return a if xq >= threshold: a, err = integrate(_probability, threshold, xq) return a return 0.0