def log_counts_to_cols_ref(counts, A=4):
"""return number of cols associated given counts"""
N = sum(counts)
# all_cols = 4**N
metacounts = Counter(counts)
log_counts_to_bases = log_fac(A) - sum(log_fac(multiplicity) for multiplicity in metacounts.values())
log_bases_to_pos = (log_fac(N) - sum(log_fac(count) for count in counts))
return log_counts_to_bases + log_bases_to_pos
def log_counts_to_cols_ref(counts, A=4):
"""return number of cols associated given counts"""
N = sum(counts)
# all_cols = 4**N
metacounts = Counter(counts)
log_counts_to_bases = log_fac(A) - sum(
log_fac(multiplicity) for multiplicity in metacounts.values())
log_bases_to_pos = (log_fac(N) - sum(log_fac(count) for count in counts))
return log_counts_to_bases + log_bases_to_pos
def log_counts_to_cols(counts, A=4):
N = sum(counts)
# all_cols = 4**N
metacounts = defaultdict(int)
for c in counts:
metacounts[c] += 1
log_counts_to_bases = log_fac(A) - sum(log_fac(multiplicity) for multiplicity in metacounts.values())
log_bases_to_pos = (log_fac(N) - sum(log_fac(count) for count in counts))
return log_counts_to_bases + log_bases_to_pos
def log_counts_to_cols(counts, A=4):
N = sum(counts)
# all_cols = 4**N
metacounts = defaultdict(int)
for c in counts:
metacounts[c] += 1
log_counts_to_bases = log_fac(A) - sum(
log_fac(multiplicity) for multiplicity in metacounts.values())
log_bases_to_pos = (log_fac(N) - sum(log_fac(count) for count in counts))
return log_counts_to_bases + log_bases_to_pos
