https://github.com/ekg/freebayes
Tip revision: 689d6d1336e163bdbda4f4ae431bb23c3e831998 authored by Erik Garrison on 29 May 2019, 10:56:46 UTC
Merge pull request #543 from ekg/coverage-limiting
Merge pull request #543 from ekg/coverage-limiting
Tip revision: 689d6d1
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)