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

[Chris]

Hi Mark,

I might be missing something, but I cant see how that will work. ps()  takes
as input a tree with the tags already in the leaf nodes and the key for
those tags.  I think your suggested solution places the tags from the @seq
array into ps() which throws an exception error.

Maybe another explanaition;

I would like to calculate the parsimony score at each position of an
alignment. Consider 1 block with only 1 character in for each species;

species1    A
species 2   G
species 3   G
species 4   G

I can do this for 1 position, by adding the tag (which Ive named 'seq')
value pair  to each leaf node of the phylogeny tree for the 4 species, where
the tag is a single character.  Then calculating the parsimony for that
column is simply

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

Now consider an alignment block with more than 1 character
species1    ACG
species 2   GCT
species 3   GGG
species 4   GGT

 I see it is possible to add an array of values as a tag, now, instead of a
single character as a tag on each leaf node, I have an array of nucleotides
as tags at each leaf node. I was hoping for a way to calculate the parsimony
given a tree with such tag value pairs.

Regards, Chris.




On Fri, Oct 23, 2009 at 11:12 PM, Mark A. Jensen <maj at fortinbras.us> wrote:

> up to list...
> ----- Original Message ----- From: "Mark A. Jensen" <maj at fortinbras.us>
> To: "Chris" <coldmeadow at gmail.com>
> Sent: Friday, October 23, 2009 8:14 AM
> Subject: Re: [Bioperl-l] Bio::Tree:Statistics parsimony question
>
>
>
>  Hi Chris,
>> looks like you could simply do
>>
>> @ps = map { $obj->ps($tree, $_) } @seq;
>>
>> or into a hash %ps with tags as keys:
>>
>> @ps{@seq} = map { $obj->ps($tree,$_) } @seq;
>>
>> cheers MAJ
>> ----- Original Message ----- From: "Chris" <coldmeadow at gmail.com>
>> To: <bioperl-l at bioperl.org>
>> Sent: Friday, October 23, 2009 1:10 AM
>> Subject: [Bioperl-l] Bio::Tree:Statistics parsimony question
>>
>>
>>  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.
>>> _______________________________________________
>>> Bioperl-l mailing list
>>> Bioperl-l at lists.open-bio.org
>>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>>>
>>>
>>>  _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>