Entering edit mode
7.8 years ago
miyuki_iii
•
0
Hi, I have a SummarizeExperiment data of Mus Musculus. I'm trying to do a hyperGTest , buit I obtain an error.
library("GO.db") library("Category") library("GOstats") library("org.Mm.eg.db")
1) The gene universe:
G = unlist(rownames(se), org.Mm.egSYMBOL2EG)
anyDuplicated(G)
2)
genes.c= c(genesUp.c, genesDown.c)
genes.t= c(genesUp.t, genesDown.t)
params.c= new("GOHyperGParams", geneIds= genes.c, universeGeneIds= G,
annotation= "org.Mm.eg", ontology = "BP", pvalueCutoff= 0.001,
conditional= FALSE, testDirection= "over")
overRepresented.c= hyperGTest(params.c)
ERROR -> Error in getUniverseHelper(probes, datPkg, entrezIds) : After filtering, there are no valid IDs that can be used as the Gene universe. Check input values to confirm they are the same type as the central ID used by your annotation package. For chip packages, this will still mean the central GENE identifier used by the package (NOT the probe IDs).