Skip to content

Instantly share code, notes, and snippets.

@hyphaltip
Created December 7, 2017 04:52
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save hyphaltip/35130486e2f1036da3af2c713d1b626f to your computer and use it in GitHub Desktop.
Save hyphaltip/35130486e2f1036da3af2c713d1b626f to your computer and use it in GitHub Desktop.
tabulate_gene_tree
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