Gene set enrichment analysis
0
0
Entering edit mode
16 months ago
Farha ▴ 10

I am trying to perform GSEA in R with data from DESeq2(CSV file containing a list of gene names and log2 fold change values)

1- loaded the required packages

2- prepare input (csv file), I extracted the log2 fold change only from the data

3- use gseGO() function from clusterProfiler package to perform GSEA

gse <- gseGO(geneList=gene_list, 
             ont ="ALL", 
             keyType = "ENSEMBL", 
             nPerm = 10000, 
             minGSSize = 3, 
             maxGSSize = 800, 
             pvalueCutoff = 0.05, 
             verbose = TRUE, 
             OrgDb = organism, 
             pAdjustMethod = "none")

4- made Dotplot

require(DOSE)
dotplot(gse, showCategory=10, split=".sign") + facet_grid(.~.sign)

5- and then I am trying to make Encrichment Map

emapplot(gse, showCategory = 10)

but it gives me this error

Error in has_pairsim(x) : 
      Term similarity matrix not available. Please use pairwise_termsim function to deal with the results of enrichment analysis.
NGS GSEA Enrichment • 910 views
ADD COMMENT
0
Entering edit mode

The error message is literally telling you what to do.

ADD REPLY

Login before adding your answer.

Traffic: 2610 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6