Entering edit mode
I have list of 222 genes with ENSEBL id, gene symbol, log fold Then I executed the following commands
data <- read.table("human_retine.txt", header=TRUE, dec=",")
# Gene Ontology classification
genes <- data$Ensembl_Gene_ID
genes <- as.character(genes)
ggo <- groupGO(gene = genes, Orgdb ="org.Hs.eg.db", keytype = "ENSEMBL", ont = "BP", level = 3, readable = TRUE)
# Enrich GO
genes<-data$Ensembl_Gene_ID[abs(data$log2FC_DeSeq)]
ego <- enrichGO(genes , Orgdb ="org.Hs.eg.db", keytype = "ENSEMBL", , ont= "CC", pAdjustMethod="BH", pvalueCutoff=0.1, qvalueCutoff=0.5, readable= TRUE)
I am getting only 3 genes in ego variable. I need help to get dotplot and enrich barplot for my data.