I have been doing some Gene Ontology enrichment analyses for a non model species using the goseq R package. For visualization I would like to use the fold-enrichment of the respective GO-terms. The output columns are a bit confusing to me though:
category = GO -Identifyer numDEInCat = number of genes in the geneset associated with category. nimInCat = total genes associated with the category (what is called universe/back ground in packages like TOPGO)
Total = total number of genes provided in the bias cdna function ClusterSize =?
I would expect to calculate the fold change based on the number of genes for one go-category to be found in the gene set. Like numInCat/total background genes =expected foldchange would then be numDEInCat/expected.
The total number of background genes is not returned though.
Has anyone else calculated GO-fold enrichment from the goseq output before?
0 answers
No answers yet.
Log in to answer this question.
-
Should gene2cat have all the assayed gene IDs?
written by Kazuna •
-
Functional enrichment analysis of bacteria
written by mickey_95 •
-
topGO on RNA-seq
written by ccstaats •
-
How do I extract the p-value from runTest() in the TopGO package?
written by ktm •
-
Is there a way to get the name or IDs of genes with the goseq table?
written by unawaz •
-
PFAM, KEGG enrichment for non model organism
written by mjp •
-
GOseq returns different number of genes in category for same gene2cat list
written by muhittinozdemir •
-
Why numInCat are different from the number of GOs present in the list it is calculated from?
written by Parham
-
how to extract genes from GO category from a non model organism
written by upendrakumar.devisetty
-
Bioconductor Goseq - Overrepresented P-Values
written by alittleboy
