[Bioperl-l] phylogeny-trait association methods into BioPerl

Hilmar Lapp hlapp at gmx.net
Wed Sep 10 13:44:27 UTC 2008


Sounds great Heikki!

Just FYI, there is a considerable amount of code and packages for  
comparative analysis in R. You might want to look into the resources  
linked to from the http://r-phylo.org page. There is also a R-SIG- 
Phylo special interest group mailing list (should be linked from the  
aforementioned site).

	-hilmar

On Sep 10, 2008, at 5:32 AM, Heikki Lehvaslaiho wrote:

> FYI,
>
> I've been recently writing code to analyse phylogeny-trait  
> associations. These
> traits are typically geographical location of the sequence but they  
> can be any
> phenotypic characters associated with the sequences.
>
> This involves trees, i.e. Bio::Tree::Tree and Bio::Tree::Node  
> objects and
> strings describing the traits. I've been using tags to store trait  
> values
> within nodes. The tag methods are:
>
> Bio::Tree::Node::add_tag_value
> Bio::Tree::Node::get_all_tags
> Bio::Tree::Node::get_tag_values
> Bio::Tree::Node::has_tag
> Bio::Tree::Node::remove_all_tags
> Bio::Tree::Node::remove_tag
>
> Question: Is  there any particular reason why there is no
> set_tag_value(scalar|@array) method?
>
> I am getting tired of writing:
>  $node->remove_tag($key);
>  map {$node->add_tag_value($key)} @values ;
> so I am going to implement that unless there is are strong objections.
>
> Otherwise it has been smooth sailing. I am going to add
> Bio::Tree::TreeFunctions::is_binary() and start populating
> Bio::Tree::Statistics soon with these methods:
>
> ps() - Parsimony Score (PS) from Fitch 1971
> ai() - Association index (AI) of Whang et al. 2001
> mc() - Monophyletic Clade (MC) size statistics by Salemi at al. 2005
> cherries() - number of leaf node pairs
>
> If you have any comments, please feel free to post them here.
>
>    -Heikki
>
> -- 
> ______ _/      _/_____________________________________________________
>      _/      _/
>     _/  _/  _/  Heikki Lehvaslaiho    heikki at_sanbi _ac _za
>    _/_/_/_/_/  Senior Scientist    skype: heikki_lehvaslaiho
>   _/  _/  _/  SANBI, South African National Bioinformatics Institute
>  _/  _/  _/  University of Western Cape, South Africa
>     _/      Phone: +27 21 959 2096   FAX: +27 21 959 2512
> ___ _/_/_/_/_/________________________________________________________
>
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l

-- 
===========================================================
: Hilmar Lapp  -:-  Durham, NC  -:-  hlapp at gmx dot net :
===========================================================






More information about the Bioperl-l mailing list