Finding the proper universe for enrichGO and enrichKEGG - only the genes that are annotated in those ontologies
1
0
Entering edit mode
endrebak85 ▴ 40
@endrebak85-10660
Last seen 5.3 years ago
github.com/endrebak/

I want to find the correct background gene list ("universe") to use with ChIPseeker. Ideally I would like to use all the genes that are annotated in GO as the universe for enrichGO and all the genes that are annotated in KEGG for enrichKEGG. Does anyone know how to find these lists? The current way from the docs is to use the gene list from DOSE, but that list contains all genes that have a disease ontology annotation.

chipseeker ChIPseeker • 4.2k views
ADD COMMENT
2
Entering edit mode
endrebak85 ▴ 40
@endrebak85-10660
Last seen 5.3 years ago
github.com/endrebak/
library(org.Hs.eg.db)

df = as.data.frame(org.Hs.egGO)
go_gene_list = unique(sort(df$gene_id))

dfk = as.data.frame(org.Hs.egPATH)
kegg_gene_list = unique(sort(dfk$gene_id))
ADD COMMENT

Login before adding your answer.

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