https://github.com/ekg/freebayes
Tip revision: a676799f5a64c286b3508997bd3cfe92ed4f759b authored by Erik Garrison on 07 January 2020, 17:28:13 UTC
be more sensitive to variant input
be more sensitive to variant input
Tip revision: a676799
factorialln.py
from scipy.special import gamma, gammaln
def factorialln(n):
if n == 1:
return 0
elif n == 0:
return 0
elif n < 0:
raise Exception("factorial is not defined for n < 0")
else:
return gammaln(n + 1)