[BioC] topGO expected value column showing NA

Kishor Tappita kishor.tappita at gmail.com
Mon Mar 11 12:18:26 CET 2013


Dear Bioconductor Users,

I am trying to find enriched GO terms using topGO. I am using Weight01
algorithm for finding the enrichment.
I get significantly enriched GO Terms but I see an NA value in the
expected value column. Below is the code I am
using. Please help me with solving the NA issue.

 sampleGOdata <- new("topGOdata", description = "Simple session",
ontology = "BP",allGenes = geneList, geneSel = topDiffGenes,nodeSize =
10,annot = annFUN.db,affyLib= "hgu133plus2.db")
resultFisher <- runTest(sampleGOdata, algorithm = "classic", statistic
= "fisher")
resultWeight01 <- runTest(sampleGOdata, algorithm ="weight01",
statistic = "fisher")
allRes <- GenTable(sampleGOdata, classicFisher =
resultFisher,weight01=resultWeight01,orderBy = "weight01", ranksOf =
"classicFisher", topNodes = 100,numChar=200)

topDiffGenes <- function(Score) {
  return(Score < 0.05)
}

Thanks,
Kishor

> sessionInfo()
R version 2.15.1 (2012-06-22)
Platform: x86_64-unknown-linux-gnu (64-bit)
locale:
 [1] LC_CTYPE=en_SG       LC_NUMERIC=C         LC_TIME=en_SG
 [4] LC_COLLATE=en_SG     LC_MONETARY=en_SG    LC_MESSAGES=en_SG
 [7] LC_PAPER=C           LC_NAME=C            LC_ADDRESS=C
[10] LC_TELEPHONE=C       LC_MEASUREMENT=en_SG LC_IDENTIFICATION=C
attached base packages:
 [1] grid      stats4    splines   stats     graphics  grDevices utils
 [8] datasets  methods   base
other attached packages:
 [1] hgu133plus2.db_2.8.0 org.Hs.eg.db_2.8.0   Rgraphviz_2.2.1
 [4] BiocInstaller_1.8.3  topGO_2.10.0         SparseM_0.96
 [7] GO.db_2.8.0          RSQLite_0.11.2       DBI_0.2-5
[10] AnnotationDbi_1.20.5 Biobase_2.18.0       BiocGenerics_0.4.0
[13] graph_1.36.2         multicore_0.1-7      coin_1.0-21
[16] modeltools_0.2-19    mvtnorm_0.9-9994     survival_2.36-14
[19] samr_2.0             matrixStats_0.6.2    impute_1.32.0
loaded via a namespace (and not attached):
[1] IRanges_1.16.6    lattice_0.20-6    parallel_2.15.1   R.methodsS3_1.4.2
[5] tools_2.15.1



More information about the Bioconductor mailing list