[Biopython-dev] Bio.Nexus and internal node labels (Bug 2788)

Nick Matzke matzke at berkeley.edu
Mon Jul 20 19:56:18 UTC 2009



Peter wrote:
> On Mon, Jul 20, 2009 at 8:13 PM, Nick Matzke wrote:
>> Hi all, here is my weekly update...
>>
>> 1. Bug fix on Nexus.Tree class is working well so far.  Thanks Brad!!
>>
> 
> Cool. I haven't tried it personally though ;) Frank and/or Cymon - any
> comments regarding Brad checking this in? See Bug 2788 for details.
> 
> I gather you (Nick) are using this on ultrametric Newick trees

yes

  - could
> you supply a sensibly sized example to use as a unit test? Initially
> this can be just to test Brad's patch to Bio.Nexus, but try and pick
> something you can build documentation examples around in future.




This is probable a reasonable size, a subset of my bigger tree.  (Just 
gymnosperms; times are in millions of years.)

(((((Cephalotaxus:125.000000,(Taxus:100.000000,Torreya:100.000000)TT1:25.000000)Taxaceae:90.000000,((((((((Calocedrus:85.000000,Platycladus:85.000000)CP:5.000000,(Cupressus:85.000000,Juniperus:85.000000)CJ:5.000000)CJCP:5.000000,Chamaecyparis:95.000000)CCJCP:5.000000,(Thuja:7.870000,Thujopsis:7.870000)TT2:92.13)CJCPTT:30.000000,((Cryptomeria:120.000000,Taxodium:120.000000)CT:5.000000,Glyptostrobus:125.000000)CTG:5.000000)CupCallTax:5.830000,((Metasequoia:125.000000,Sequoia:125.000000)MS:5.000000,Sequoiadendron:130.000000)Sequoioid:5.830000)STCC:49.060001,Taiwania:184.889999)Taw+others:15.110000,Cunninghamia:200.000000)nonSci:15.000000)Tax+nonSci:10.000000,Sciadopitys:225.000000):25.000000,(((Abies:106.000000,Keteleeria:106.000000)AK:54.000000,(Pseudolarix:156.000000,Tsuga:156.000000)NTP:4.000000)NTPAK:24.000000,((Larix:87.000000,Pseudotsuga:87.000000)LP:81.000000,(Picea:155.000000,Pinus:155.000000)PPC:13.000000)Pinoideae:16.000000)Pinaceae:66.000000)Coniferales:25.000000,Gink
go:275.000000)gymnosperm:75.000000;

I am using this as a test case in my own script, I have put some effort 
into reading up on Unittest but I don't quite get how it all fits 
together yet.


Another important case I will try and come up with is the results of 
pruning a tree, in my experience it is very easy to mess up the tree 
and/or branchlengths when pruning.



> 
> Thanks,
> 
> Peter
> 

-- 
====================================================
Nicholas J. Matzke
Ph.D. Candidate, Graduate Student Researcher
Huelsenbeck Lab
Center for Theoretical Evolutionary Genomics
4151 VLSB (Valley Life Sciences Building)
Department of Integrative Biology
University of California, Berkeley

Lab websites:
http://ib.berkeley.edu/people/lab_detail.php?lab=54
http://fisher.berkeley.edu/cteg/hlab.html
Dept. personal page: 
http://ib.berkeley.edu/people/students/person_detail.php?person=370
Lab personal page: http://fisher.berkeley.edu/cteg/members/matzke.html
Lab phone: 510-643-6299
Dept. fax: 510-643-6264
Cell phone: 510-301-0179
Email: matzke at berkeley.edu

Mailing address:
Department of Integrative Biology
3060 VLSB #3140
Berkeley, CA 94720-3140

-----------------------------------------------------
"[W]hen people thought the earth was flat, they were wrong. When people 
thought the earth was spherical, they were wrong. But if you think that 
thinking the earth is spherical is just as wrong as thinking the earth 
is flat, then your view is wronger than both of them put together."

Isaac Asimov (1989). "The Relativity of Wrong." The Skeptical Inquirer, 
14(1), 35-44. Fall 1989.
http://chem.tufts.edu/AnswersInScience/RelativityofWrong.htm
====================================================



More information about the Biopython-dev mailing list