Entering edit mode
2.1 years ago
Priya Rao
•
0
This is my code:
library(topGO)
hs_topGOdata <- new("topGOdata", description = "Human proteins", ontology = c("BP"), allGenes = allGenes.hs, geneSel = interestingpValues.hs, nodeSize=5, annot = annFUN.gene2GO, gene2GO = geneID2GO.hs)
resultFisher <- runTest(hs_topGOdata, algorithm = "classic", statistic = "fisher")
resultFisher
resultFisher.elim <- runTest(hs_topGOdata, algorithm = "elim", statistic = "fisher")
resultFisher.elim
resultKS <- runTest(hs_topGOdata, algorithm = "classic", statistic = "ks")
resultKS
resultKS.elim <- runTest(hs_topGOdata, algorithm = "elim", statistic = "ks")
resultKS.elim
All my other tests are running well except `resultKS.elim` which gives me the following error:
```r
Error in ks.test.default(x.a, seq_len(N)[-x.a], alternative = "greater") :
not enough 'y' data
I have searched for the NA
s and there are no NA
s!