-
Notifications
You must be signed in to change notification settings - Fork 12
Description
Hi,
I am using ClusterProfiler on R to run my GSEA and I am having trouble with the files I am exporting to EnrichmentMap on Cystoscape.
Straight out of R, the file looks like this:
ID setSize enrichmentScore NES pvalue p.adjust qvalue rank leading_edge core_enrichment core_size
GOCC_MITOCHONDRION 1781 0.423194013 2.065931975 3.59326E-46 3.57601E-42 2.69003E-42 3463 tags=34%, list=22%, signal=30% Slc11a2/App/Sdhd/Micu3/... 606
This gives me the cannot parse line 2 error when uploading as GSEA for Analysis type
And then, using the following format as reference: Supplementary_Table8_gsea_report_for_na_pos that I obtained from https://baderlab.github.io/Cytoscape_workflows/EnrichmentMapPipeline/index.html
I changed my headings to:
NAME GS SIZE ES NES pvalue qvalue p.adjust rank leading_edge
GOCC_MITOCHONDRION GOCC_MITOCHONDRION 1781 0.423194013 2.065931975 3.59326E-46 2.69003E-42 3.57601E-42 3463 tags=34%, list=22%, signal=30%
And still got the parsing error
Next, I tried using the generic option for analysis type using the second format mention, and still got the parsing error.
Then, I tried using the headings suggested for generic enrichment results from: https://enrichmentmap.readthedocs.io/en/latest/FileFormats.html, which also gave me the parsing error
ID description p-value FDR gene list
GOCC_MITOCHONDRION GOCC_MITOCHONDRION 3.57601E-42 2.69003E-42 Slc11a2,App,Sdhd,Micu3,...
For all of the above, my gmt file and rnk file stayed the same.
GMT (tab seperated):
GOBERT_CORE_OLIGODENDROCYTE_DIFFERENTIATION NA Abat Abca1
RNK:
gene score
Thy1 6.03204820036826
Slc11a2 5.92889879296482
Overall, I don't know why my enrichment file isn't working.
Thanks,
Nishi