Ich versuche, einen Baum mit BioPython, Phylo-Modul zu bauen.
Was ich bisher getan habe, ist dieses Bild: Phylo BioPython Gebäude Bäume
jeder Name hat eine vierstellige Zahl, gefolgt von - und eine Zahl: diese Nummer bezieht sich auf die Anzahl der Male, die Sequenz dargestellt wird. Das bedeutet 1578 - 22, dieser Knoten sollte 22 Folgen darstellen.
die Datei mit den ausgerichteten Sequenzen: file
die Datei mit dem Abstand, einen Baum zu erstellen: file
I So jetzt bekannt, wie jede Größe des Knotens zu ändern. Jeder Knoten eine andere Größe hat, das einfach ein Array von den verschiedenen Werten tut:
fh = open(MEDIA_ROOT + "groupsnp.txt")
list_size = {}
for line in fh:
if '>' in line:
labels = line.split('>')
label = labels[-1]
label = label.split()
num = line.split('-')
size = num[-1]
size = size.split()
for lab in label:
for number in size:
list_size[lab] = int(number)
a = array(list_size.values())
Aber das Array ist willkürlich, ich möchte den richtigen Knotengröße in den rechten Knoten setzen, habe ich versucht, dies:
for elem in list_size.keys():
if labels == elem:
Phylo.draw_graphviz(tree_xml, prog="neato", node_size=a)
aber nichts erscheint, wenn ich die if-Anweisung verwende.
Wie auch immer dabei?
würde ich wirklich schätzen!
Danke allen
Können Sie die Testdateien zur Verfügung stellen Sie für diesen Baum mit? – rwilliams