[Bioperl-l] perl for GO

Ross KK Leung ross at cuhk.edu.hk
Sat Sep 25 08:30:39 UTC 2010


Given a set of GO IDs, e.g. 

GO:0008150
GO:0005750
GO:0006122
GO:0008121
GO:0003674
GO:0005575
GO:0008150
GO:0009507
GO:0009535
GO:0009567
GO:0009977
GO:0010027
GO:0031361

from
http://www.geneontology.org/ontology/obo_format_1_2/gene_ontology_ext.obo

one can manually examine the hierarchy. Although there is go-perl
(http://search.cpan.org/~cmungall/go-perl/) and go-db-perl
(http://search.cpan.org/~cmungall/go-db-perl/), as a life science student
who just learns Perl, I find it difficult to draw a hierarchy tree (or
simply make it a table to count the occurrence) to produce something like:


biological_process (4)
*** cellular process (4)
****** cell adhesion (1)
****** cell differention (3)
Molecular function (4)
<snip>
Cellular component (4)
<snip>

Can anybody advise? I don't need any fancy figures at all...




More information about the Bioperl-l mailing list