Dotplot from enricher outputting very few genes
0
0
Entering edit mode
Lucía • 0
@16997962
Last seen 3 days ago

Hi,

I am running the generic version of enricher supplying my own GO mappings. When I create a dotplot from this object, I get a dotplot with less genes than the total amount of significant genes I supplied in the first place. Here's my code:

library(clusterProfiler)
library(GO.db)

aga_genes <- aga_genes\$X

aga_genes <- as.character(aga_genes)

go_terms <- go_mappings

goterms <- Term(GOTERM)

term2name <- data.frame("GOID"=names(goterms),"term"=goterms )

ego_bp_aga <- enricher(gene = aga_genes,
TERM2GENE = go_terms,
TERM2NAME = term2name)

dotplot(ego_bp_aga, showCategory=30) + ggtitle("dotplot for ORA")


where treat_vs_control_aga.csv contains my DEGs from DESeq2 and go_mappings are the gene and respective GO term that I'm supplying.

My total number of DEGs is about 500, and the dotplot shows about 100 if I add up all the categories. (Is there a way to include my image?)

Why is this? Thanks

clusterProfiler enrichplot • 62 views
0
Entering edit mode

It is difficult to guess without exactly knowing your data objects, but realize that by default GO categories containing less than 10 genes are excluded from the analyses. See the argument minGSSize = 10. Related to this, the result of sort( table(go_terms[,1]) ) will show you the number of genes per category. Also realize that only significant categories are plotted in the dotplot. See arguments pvalueCutoff, pAdjustMethod and qvalueCutoff. Set the 1st and 3rd argument to 1, and minGSSize = 2, to be most permissive.