Created
December 7, 2017 04:52
-
-
Save hyphaltip/35130486e2f1036da3af2c713d1b626f to your computer and use it in GitHub Desktop.
tabulate_gene_tree
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
def tabulate_names(tree): | |
cladenames={} | |
for idx, clade in enumerate(tree.find_clades()): | |
if clade.name: | |
clade.name = '%d_%s' % (idx, clade.name) | |
else: | |
clade.name = str(idx) | |
cladenames[clade.name] = clade | |
return cladenames | |
genefamily = {} | |
names = lookup_by_names(tree) | |
for gene in tree.get_terminals(): | |
clade = names[gene] | |
if clade in genefamily: | |
genefamily[clade] +=1 | |
else: | |
names[clade] =1 |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment