[BioC] using GOSim with yeast genes
prsmra01 at uniroma2.it
prsmra01 at uniroma2.it
Wed Mar 4 18:16:04 CET 2009
Dear list,
I cannot use getGeneSim() with yeast genes (the code in vignette works
properly for me but the entrez ids in the provided examples are from
human); actually the problem seems linked to the mapping "yeast
entrezIDs-GOterms" as you can see by trying the chunk of code below:
library(org.Sc.sgd.db)
library("GOSim")
x <- org.Sc.sgdENTREZID
setEvidenceLevel(evidences = c("IMP","IGI","ISS","IDA","IEP","IEA"),
organism="yeast", gomap=NULL)
calcICs()
#now I move the ICs file in GOSim/data
setOntology("BP", loadIC=FALSE)
filterGO(Rkeys(x))
filtering out genes not mapping to the currently set GO category ...
===> list of 6199 genes reduced to 0
list()
> sessionInfo()
R version 2.9.0 Under development (unstable) (2009-03-03 r48046)
x86_64-unknown-linux-gnu
locale:
LC_CTYPE=it_IT.UTF-8;LC_NUMERIC=C;LC_TIME=it_IT.UTF-8;LC_COLLATE=it_IT.UTF-8;LC_MONETARY=C;LC_MESSAGES=it_IT.UTF-8;LC_PAPER=it_IT.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=it_IT.UTF-8;LC_IDENTIFICATION=C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] GOSim_1.1.5.4 corpcor_1.5.2 Matrix_0.999375-21
[4] lattice_0.17-20 RBGL_1.19.2 mclust_3.1-10
[7] cluster_1.11.12 topGO_1.11.1 SparseM_0.79
[10] graph_1.21.4 org.Rn.eg.db_2.2.6 org.Pf.plasmo.db_2.2.6
[13] org.Mm.eg.db_2.2.6 org.Hs.eg.db_2.2.6 org.Dm.eg.db_2.2.6
[16] annotate_1.21.4 GO.db_2.2.5 org.Sc.sgd.db_2.2.8
[19] RSQLite_0.7-1 DBI_0.2-4 AnnotationDbi_1.5.17
[22] Biobase_2.3.10
loaded via a namespace (and not attached):
[1] grid_2.9.0 tools_2.9.0 xtable_1.5-4
May you give me some suggestions?
Thanks for your attention,
Best,
Maria
--
Maria Persico, PhD. student
http://cbm.bio.uniroma2.it/~maria/
MINT database group
Universita' di Tor Vergata, via della Ricerca scientifica 11
00133 Roma, Italy
Tel +39 0672594315 (Supervisor's room)
Fax +39 0672594766
Mobile phone: +393479715662
e-mail maria.persico at uniroma2.it
More information about the Bioconductor
mailing list