我试图用BioPython,Phylo模块构建树。
什么我目前所做的就是这种形象:Phylo BioPython构建树
每个名字都有一个四位数字,其次是 - 和一些:这个数字指的倍序列代表的数字。这意味着1578-22,那个节点应该代表22个序列。
对齐序列文件:file
所以,现在我知道如何在节点的每个尺寸变化:file
与构建一棵树的距离的文件。每个节点都有不同的大小,这是很容易做的不同值的数组:
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())
但阵列是任意的,我希望把正确的节点尺寸为右节点,我尝试这样做:
for elem in list_size.keys():
if labels == elem:
Phylo.draw_graphviz(tree_xml, prog="neato", node_size=a)
但是当我使用if语句时什么也没有出现。
反正这么做?
我真的很感激!
谢谢大家
你能否提供你正在使用那棵树的测试文件? – rwilliams 2010-11-02 08:50:07