[BioC] GOstats : Error in getUniverseHelper - No Entrez Gene ids left in universe

Robert Gentleman rgentlem at fhcrc.org
Wed Aug 5 23:37:51 CEST 2009


Hi Ben,
   First you should update your R (2.9.1 is the current release).
Next, it looks like the IDs in your universe might not be correct (or 
something of that nature). Are you sure they are EG IDs - they don't 
look like they are.

   best wishes
    Robert


Benoit Ballester wrote:
> Hi,
> 
> I am getting this error when I run hyperGTest from GOstats. Any idea on 
> what is going wrong would be appreciated ?
> 
> R version 2.8.1 (2008-12-22)
> GOstats  1.7.4
> I am also attaching the files genefile/universe.txt
> 
> write.table(ori$Accession, "universe.txt", sep = "\t", 
> quote=FALSE,row.names=FALSE,col.names=FALSE)
> 
> geneFile  <- "genefile.txt"
> uniFile   <- "universe.txt"
> 
> universe  <- scan(uniFile, what="character")
> genelist  <- scan(geneFile, what="character")
> 
> pcut    <- 0.05
> catSize <- 5
> 
> 
> params    <- new("GOHyperGParams", geneIds=genelist, 
> universeGeneIds=universe,
>               annotation="org.Mm.eg", ontology="BP", pvalueCutoff=pcut,
>               conditional=FALSE, testDirection="over")
> 
> 
> results <- hyperGTest(params)
> Error in getUniverseHelper(probes, datPkg, entrezIds) :
>   No Entrez Gene ids left in universe
> 
> 
> 
> 
> Thanks,
> 
> Ben
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor

-- 
Robert Gentleman, PhD
Program in Computational Biology
Division of Public Health Sciences
Fred Hutchinson Cancer Research Center
1100 Fairview Ave. N, M1-B514
PO Box 19024
Seattle, Washington 98109-1024
206-667-7700
rgentlem at fhcrc.org



More information about the Bioconductor mailing list