[Bioperl-l] Bio::Tree:Statistics parsimony question

Chris coldmeadow at gmail.com
Fri Oct 23 05:10:12 UTC 2009


Hello,
I have been able to use the ps() method in Bio::Tree::Statistics to calulate
the parsimony for tag values from a tree if they are strings.

This is the way I have been implementing the method:

my $obj = Bio::Tree::Statistics->new();
my $pars = $obj->ps($tree,'seq');
print "Parsimony = $pars \n";

Is it possible to do so if they are an array?

say I have an array @seq = qw(A,A,C,G), I can add these as values to a leaf
node of a tree using;

foreach my $char (@seq){
   $node->add_tag_value('seq',$char);
}

Does anybody have a suggestion to get the parsimony for each label of the
array of labels for such a tree ?

Thanks in advance,
Chris.



More information about the Bioperl-l mailing list