Entering edit mode
cjsifuen
•
0
@cjsifuen-12432
Last seen 7.8 years ago
I am having an issue with using the plotGOgraph() command on the gseaResult object created from the gseGO() command in clusterProfiler. The command works without any issues when using an enrichResult object however.
With enrichGO():
> GO.ora.bp.filt <- enrichGO(gene = gene.filtered, universe = gene, OrgDb = "org.Xl.eg.db", ont = "BP", pAdjustMethod = "BH", pvalueCutoff = 1, qvalueCutoff = 1, readable = TRUE, keytype = "ENTREZID")
> plotGOgraph(GO.ora.bp.filt)
groupGOTerms: GOBPTerm, GOMFTerm, GOCCTerm environments built.
Building most specific GOs .....
( 2151 GO terms found. )
Build GO DAG topology ..........
( 2152 GO terms and 4738 relations. )
Annotating nodes ...............
( 2932 genes annotated to the GO terms. )
$dag
A graphNEL graph with directed edges
Number of Nodes = 56
Number of Edges = 84
$complete.dag
[1] "A graph with 56 nodes."
With gseGO():
> GO.bp <- gseGO(geneList = geneList, ont = "BP", OrgDb = "org.Xl.eg.db", keytype = "ENTREZID", nPerm = 1000, exponent = 1, minGSSize = 10, pvalueCutoff = 1, pAdjustMethod = "BH", verbose = TRUE, seed = FALSE, by = "DOSE") preparing geneSet collections... calculating observed enrichment scores... calculating permutation scores... calculating p values... leading edge analysis... done...
> plotGOgraph(GO.bp)
groupGOTerms: GOBPTerm, GOMFTerm, GOCCTerm environments built. Error in .local(.Object, ...) : no slot of name "ontology" for this object of class "gseaResult"
Any help with this would be greatly appreciated.
Thanks,
C