<p>I'd probably just use <a href="https://bioconductor.org/packages/release/bioc/html/fgsea.html" rel="nofollow noreferrer">fgsea</a>. Your code won't have to change much, it also just uses a ranked list.</p>

CC BY-SA 4.0
2021-07-06T15:57:44.223Z
2021-07-06T15:57:44.223Z
Jared Andrews
<p>I want to draw the standard plot of GSEA on Rstudio. I have a data frame that consists of a list of DEGs as below.
<a href="https://i.sstatic.net/D4kdb.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/D4kdb.png" alt="Structure" /></a></p>
<p>Actually, I put the DEG list directly to PANTHER that resulted in the GO or pathway analysis data. However, I want to use GSEA.</p>
<p>Meanwhile, I found nice online site and follow their codes for the GSEA analysis(<a href="https://bioinformaticsbreakdown.com/how-to-gsea/" rel="nofollow noreferrer">https://bioinformaticsbreakdown.com/how-to-gsea/</a>).</p>
<p>And, here is the my code.</p>
<pre><code>tmp = read.csv(&quot;/data/directory/DEGs_HCvsJSS_DC.csv&quot;, header = TRUE) %&gt;% filter(Gene.Symbol !=&quot;&quot;)

GO_file = &quot;/blue/zhangw/my.kim/HCvsJSS/c5.all.v7.4.symbols.gmt&quot; # got this from http://www.gsea-msigdb.org/gsea/msigdb/collections.jsp

gene_list = tmp$avg_log2FC

names(gene_list) = tmp$Gene.Symbol

gene_list = sort(gene_list, decreasing = TRUE)

gene_list = gene_list[!duplicated(names(gene_list))]

head(gene_list)

res = GSEA(gene_list, GO_file, pval = 0.05) # this caused the error as follows. Error in build_Anno(TERM2GENE, TERM2NAME) :  argument &quot;TERM2GENE&quot; is missing, with no default
</code></pre>
<p>It's really dissapointing because I found this relatively easy protocol for GSEA, but found the error, lol.</p>
<p>Do you guys have any suggests for this troubles? or any recommendations for GSEA analysis that can be done with DEG list generated from seruat scRNAseq data?</p>

<r><scrnaseq><gsea>
How to run GSEA analysis on R studio using DEG file list generated from scRNAseq analysis
# How to run GSEA analysis on R studio using DEG file list generated from scRNAseq analysis

(Asked by: username_0 on 2021-07-02)

I want to draw the standard plot of GSEA on Rstudio. I have a data frame that consists of a list of DEGs as below.

Actually, I put the DEG list directly to PANTHER that resulted in the GO or pathway analysis data. However, I want to use GSEA.
Meanwhile, I found nice online site and follow their codes for the GSEA analysis(https://bioinformaticsbreakdown.com/how-to-gsea/).
And, here is the my code.
tmp = read.csv("/data/directory/DEGs_HCvsJSS_DC.csv", header = TRUE) %>% filter(Gene.Symbol !="")

GO_file = "/blue/zhangw/my.kim/HCvsJSS/c5.all.v7.4.symbols.gmt" # got this from http://www.gsea-msigdb.org/gsea/msigdb/collections.jsp

gene_list = tmp$avg_log2FC

names(gene_list) = tmp$Gene.Symbol

gene_list = sort(gene_list, decreasing = TRUE)

gene_list = gene_list[!duplicated(names(gene_list))]

head(gene_list)

res = GSEA(gene_list, GO_file, pval = 0.05) # this caused the error as follows. Error in build_Anno(TERM2GENE, TERM2NAME) :  argument "TERM2GENE" is missing, with no default

It's really dissapointing because I found this relatively easy protocol for GSEA, but found the error, lol.
Do you guys have any suggests for this troubles? or any recommendations for GSEA analysis that can be done with DEG list generated from seruat scRNAseq data?

**Tags:** r, scrnaseq, gsea

---

## 1 answers

### Answer by username_1 on 2021-07-06. Score: 1

I'd probably just use [fgsea](https://bioconductor.org/packages/release/bioc/html/fgsea.html). Your code won't have to change much, it also just uses a ranked list.


---

