BioPython、Phyloモジュールでツリーを構築しようとしています。
私がこれまでに行ったことは、この画像です。
各名前には、4桁の数字の後に-と数字が続きます。この数字は、そのシーケンスが表される回数を示します。つまり、1578-22、そのノードは22シーケンスを表す必要があります。
シーケンスが整列され
たファイル:ツリーを構築する距離のファイルをファイルします: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ステートメントを使用しても何も表示されません。
とにかくこれを行う?
本当にありがたいです!
みんなありがとう