This is a test version of Biostars. For the public version, visit https://www.biostars.org.
error when using ViSEAGO package

I get the following error when I run the ViSEAGO package

Error in .local(.Object, ...) : allGenes must be a factor with 2 levels

I have a list of differentially expressed genes (from a DE proteomics experiment) and wanted to use the ViSEAGO package to visualise the GO profile. If I understand correctly I need to supply my DE genes as a vector and a list of background genes as a separate vector. Where I am struggling is, what needs to be in these files? Obviously, gene ids need to be present and that's all I supplied. My code is as follows:

str(Downreg_gene_id_v4)

data.frame': 286 obs. of 1 variable:

$ Entrez_Gene_ID: int 2288 5898 7384 3028 11345 5704 7414 1659 84844 10291 ...

#set up gene list to be analysed and background list
write.table(Downreg_gene_id_v4, sep="\t", file="selection.txt")
selection <- read.table("selection.txt")
background=keys(org.Hs.eg.db, keytype ='ENTREZID')
write.table(background, sep="\t", file="background.txt")
background <- read.table("background.txt")
names(background)[1] <- "Entrez_Gene_ID"

from the package vingette

# connect to Entrezgene to get GO annotations
EntrezGene<-ViSEAGO::EntrezGene2GO()

# load GO annotations from EntrezGene
#id = "9606" = homo sapiens
myGENE2GO<-ViSEAGO::annotate(id="9606", EntrezGene)

# create topGOdata for BP
BP<-ViSEAGO::create_topGOdata(
geneSel=selection,
allGenes=background,
gene2GO=myGENE2GO, 
ont="BP",
nodeSize=5
)

It's at this point I get the error message

Error in .local(.Object, ...) : allGenes must be a factor with 2 levels

Any advice is gratefully recieved.

geneontology viseago

0 answers

No answers yet.

Log in to answer this question.