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

Chris coldmeadow at gmail.com
Fri Oct 23 01:10:12 EDT 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