[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