Tutorial:Easy way to do GO enrichment in R
1
11
Entering edit mode
5.0 years ago
ahmad mousavi ▴ 730

Hi friends

After all I found a way to do GO enrichment in a programming and very easy way :

You should install enrichR package and then use following code:

library(enrichR)
dbs <- listEnrichrDbs()
dbs <- c("GO_Molecular_Function_2015", "GO_Cellular_Component_2015", "GO_Biological_Process_2015" , "ChEA_2016" ,"KEGG_2016")
enriched <- enrichr(c("Runx1", "Gfi1", "Gfi1b", "Spi1", "Gata1", "Kdr"), dbs)
printEnrich(enriched, "output.txt" , sep = "\t", columns = c(1:9))

bp <- enriched[["GO_Biological_Process_2015"]]

Thanks guys who developed this package !!

Cheers :D

enrichr gene RNA-Seq R Tutorial • 10k views
ADD COMMENT
1
Entering edit mode

Enrichr R is very useful because of multiple enrichment analyses in one place....however, when you go through the paper you will notice that part of their algorithm for producing the combined score is "intuition".

Because I don;t like black boxes...I stick to g:Profiler

ADD REPLY
0
Entering edit mode

enrichr was published in a good journal (NAR) and is not exactly a black box. While they provide a correction to the Fisher exact test, you don't have to use it. They also give you the uncorrected p-values. Of course, you should use whatever tool you prefer, but I don't want to give other people the false impression that enrichr is statistically unprincipled.

ADD REPLY
0
Entering edit mode

How to add or input nearly 8000 genes for annotation in enrichR.

ADD REPLY
0
Entering edit mode

just paste your gene in enrichR website; in R you could input your genes as a vector in below code :

enriched <- enrichr(genes8000, dbs)
ADD REPLY
2
Entering edit mode
6 months ago

topGo is one program - see my answers, here:

Kevin

ADD COMMENT

Login before adding your answer.

Traffic: 1876 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