Commit a878dcd6 authored by Yixiong Sun's avatar Yixiong Sun
Browse files

Removed 0.25 added to total sum for incorrect probabilities

parent 650146c8
......@@ -21,7 +21,7 @@ def get_consensus(module_seqs):
#print(scores.index(max(scores)))
#print(np.argmax(scores))
tot = max(sum(scores),4)
consensus.append([x/(tot+0.25) for x in scores])
consensus.append([x/tot for x in scores])
return consensus
def rfam_to_module(positions=[], family_file="",output_name="module_seqs.fasta", as_list=False, seqs=[]):
......@@ -68,7 +68,7 @@ def rfam_to_module(positions=[], family_file="",output_name="module_seqs.fasta",
this_mod = this_mod + str(fullseq[col-1])
module_seqs.append(this_mod)
cons = get_consensus(module_seqs)
fmodule_seqs = []
for s in module_seqs:
this_mod = ""
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment