EnrichGO
0
0
Entering edit mode
Neha • 0
@3fa075d2
Last seen 2.8 years ago
Finland

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.

R enrichplot • 1.5k views
ADD COMMENT

Login before adding your answer.

Traffic: 810 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6