[BioC] hyperGTest(paramsGO)

[James W. MacDonald]

Hi Yu,

Yuan Jian wrote:
> hello there,
>  
> when I run a toy example
>>  paramsGO <- new("GOHyperGParams", geneIds = "100113407",
>            universeGeneIds = c("6310","100113407"), annotation = "org.Hs.eg.db", ontology = "BP",
>            pvalueCutoff = 0.05, conditional = FALSE, testDirection = "over")
>>  hgGO_Over <- hyperGTest(paramsGO)
> I got error below. 
> "Error in getGoToEntrezMap_db(p) : 
>   The genes you are testing do not have any corresponding GO terms for the ontology you are searching."
>  
> In my script, geneIDs are generated at runtime. How to make the script continue running?
> how to catch the error?

?try
?tryCatch

Best,

Jim




>  
> thanks
> YU
>  
> 
>  
> 
> 
>       
> 	[[alternative HTML version deleted]]
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> 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

-- 
James W. MacDonald, M.S.
Biostatistician
Douglas Lab
University of Michigan
Department of Human Genetics
5912 Buhl
1241 E. Catherine St.
Ann Arbor MI 48109-5618
734-615-7826
**********************************************************
Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues