bd9bdcb334be151fe0195d8685d575355df5a551,third_party/nucleus/util/genomics_math.py,,log10_binomial,#Any#Any#Any#,161

Before Change

  Returns:
    log10 probability of seeing k successes in n trials with p.
  
  return scipy.stats.binom.logpmf(k, n, p) / LOG_E_OF_10


def log10sumexp(log10_probs):
  Returns log10(sum(10^log10_probs)) computed in a numerically-stable way.

After Change

  Returns:
    log10 probability of seeing k successes in n trials with p.
  
  r = math.lgamma(n + 1) - (math.lgamma(k + 1) + math.lgamma(n - k + 1))
  if k > 0:
    r += k * math.log(p)
  if n > k:
    r += (n-k) * math.log1p(-p)
  return r / LOG_E_OF_10


def log10sumexp(log10_probs):
  Returns log10(sum(10^log10_probs)) computed in a numerically-stable way.

Instances