Building Enrichment map from a gprofiler R package

66 views
Skip to first unread message

Moutaz Helal

unread,
Jun 10, 2020, 7:02:45 AM6/10/20
to cytoscape-helpdesk
Hello,
Thanks for helping me, I am using Cytoscape Version: 3.8.0 to build an Enrichment map after I used  a gprofiler2  R package to generate GO terms. I used the following chunk of code 
gostres0 <- gost(query = rownames(sc_C0_markers),
                organism
= "mmusculus", ordered_query = TRUE,
                multi_query
= FALSE, significant = TRUE, exclude_iea = TRUE,
                measure_underrepresentation
= FALSE, evcodes = TRUE,
                user_threshold
= 0.05, correction_method = "g_SCS",
                domain_scope
= "annotated", custom_bg = NULL,
                numeric_ns
= "", sources = "GO", as_short_link = FALSE)
sc_C0gprofiler_results
<- gostres0$result


## Order the results by p-adjusted value
sc_C0gprofiler_results
<- sc_C0gprofiler_results[order(sc_C0gprofiler_results$p_value), ]
gostplot
(gostres0, capped = TRUE, interactive = TRUE)
p0
<- gostplot(gostres0, capped = TRUE, interactive = FALSE)
pdf
("sc_C0__gplot.pdf", width = 11, height =11 )


publish_gostplot
(p0, highlight_terms = sc_C0gprofiler_results$term_id[1:19],
                       width
= 20, height = 20, filename = NULL)
 dev
.off()
 
# get the GEM file for cytoscape
 gem0
<- gostres0$result[,c("term_id", "term_name", "p_value", "intersection")]
 colnames
(gem0) <- c("GO.ID", "Description", "p.Val", "Genes")
 gem0$FDR
<- gem0$p.Val
 write
.table(gem0, file = "c0gProfiler_gem.txt", sep = "\t", quote = F, row.names = F)

then I downloaded GMT file from gprofiler website,  when I attempted to build the map it gave me the following Error 



please find files attached if this would help,
gprofiler output written from R (c0gProfiler_gem.txt)
gprofiler GMT file ( gprofiler_full_mmusculus.name.gmt)
Thanks very much

c0gProfiler_gem.txt
Reply all
Reply to author
Forward
0 new messages