https://github.com/ekg/freebayes
Tip revision: b5ccca1e63a9b681ab1f1a6dfd4cae3058422522 authored by Erik Garrison on 01 February 2018, 12:25:19 UTC
allow calling on N-matching bases
allow calling on N-matching bases
Tip revision: b5ccca1
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)