I am having trouble getting topGO to with with my data. It has something to due with my GO terms being in an inverse order, but from what I can tell that really shouldn't be an issue.

Anyone encounter this problem before? Below is the code and error associated with it.

geneID2GO <- readMappings(file = "aust_GO.txt") GO2geneID <- inverseList(geneID2GO) str(head(GO2geneID)) List of 6 $ GO:0000002: chr [1:2] "Caust.v2_007320" "Caust.v2_019771" $ GO:0000003: chr "Caust.v2_016664" $ GO:0000012: chr "Caust.v2_020298" $ GO:0000014: chr "Caust.v2_001917" $ GO:0000015: chr [1:2] "Caust.v2_001229" "Caust.v2_021398" $ GO:0000022: chr [1:5] "Caust.v2_000558" "Caust.v2_002098" "Caust.v2_002701" "Caust.v2_010710" ... geneUniverse <- names(GO2geneID) head(GO2geneID) $GO:0000002 [1] "Caust.v2_007320" "Caust.v2_019771"

$GO:0000003 [1] "Caust.v2_016664"

$GO:0000012 [1] "Caust.v2_020298"

$GO:0000014 [1] "Caust.v2_001917"

$GO:0000015 [1] "Caust.v2_001229" "Caust.v2_021398"

$GO:0000022 [1] "Caust.v2_000558" "Caust.v2_002098" "Caust.v2_002701" "Caust.v2_010710" "Caust.v2_019642"

genesOfInterest <- read.table("GO_MB1.txt",header=FALSE) genesOfInterest <- as.vector(genesOfInterest$V2) geneList <- factor(as.integer(geneUniverse %in% genesOfInterest)) names(geneList) <- geneUniverse str(geneList) Factor w/ 2 levels "0","1": 1 1 1 1 1 1 1 1 1 1 ... - attr(*, "names")= chr [1:3327] "GO:0000002" "GO:0000003" "GO:0000012" "GO:0000014" ... myGOdata <- new("topGOdata", description="My project", ontology="BP", allGenes=geneList, annot = annFUN.GO2genes, GO2gene = GO2geneID)

Building most specific GOs ..... ( 0 GO terms found. )

Build GO DAG topology .......... ( 0 GO terms and 0 relations. ) Error in if is.na(index) || index < 0 || index > length(nd)) stop("vertex is not in graph: ", : missing value where TRUE/FALSE needed

Thanks, Mike