[Bioperl-l] problem parsing a newick format

Chrysanthi A. chrysain at gmail.com
Thu Feb 12 02:13:29 UTC 2009


Is the code below correct?? Why it does not print anything???
use strict;

use Bio::TreeIO;


my $input = new Bio::TreeIO(-fh =>"nexusCytochrome7R.newick",
                            -format => "newick");

while(my $tree = $input->next_tree){
    for my $node(grep{!$_->is_Leaf}$tree->get_nodes){
next if !$node->ancestor;
print "Node:", $node->id, "length:", $node->branch_length, " ";
for my $child($node->get_Descendents){
    print "child:", $child->id, "", $child->branch_length, " ";
}
    print "\n";
    }
}

Any ideas? I want to read a tree and mainly get the duplication events.
Could someone help me?

Thanks a lot,

Chrysanthi



More information about the Bioperl-l mailing list