I've performed a GO enrichment analysis using goseq. The output of
goseq() tells you the enriched categories with how many significant genes are in the categories etc, but it doesn't give you the ID of genes that are in any of those categories.
Obviously you can use
getgo() to retrieve GOs for genes of interest, but is there a way to get a column of the genes in each of the enriched GO terms using the
So my desired output would be:
category over_represented_pvalue under_represented_pvalue numDEInCat numInCat term ontology Ensembl ID GO:0000786 2.143408e-16 1 12 43 nucleosome CC. ENSG00000112655, ENSG00000158483 etc..
With the "Ensembl ID" column is the one I'm looking to add
Any help would be greatly appreciated!