R语言——GSEV可视化

最新推荐文章于 2024-08-19 19:05:44 发布

bobo JLU

最新推荐文章于 2024-08-19 19:05:44 发布

阅读量920

点赞数

分类专栏： R语言学习及分析文章标签： r语言开发语言

本文链接：https://blog.csdn.net/codyliu_bo/article/details/127623474

版权

本文主要介绍了如何使用R语言进行GSEV的可视化，包括气泡图、网络图、热图和富集图的绘制，并深入探讨了高级venn图和GSEA图的使用，特别是通过gseaplot、gseaplot2和gsearank函数展示GSEA结果。

摘要由CSDN通过智能技术生成

环境加载

注意这里面一些参数现在我还没有完全运行好，因此需要注意。
# 清空变量
remove(list = ls())

# 加载所需R包
library(clusterProfiler)
library(ggplot2)

加载GSEA结果，以便用于后续可视化

load(file = "data/gsea_result.Rda")

气泡图

########## 气泡图：dotplot ##########
tiff(file="figures/GSEA_GO_dotplot.tiff",
     width = 35,height = 22,units ="cm",
     compression="lzw",bg="white",res=300)
dotplot(gsego, showCategory=30) + ggtitle("dotplot for GSEA")
dev.off()

dotplot(gsekk, showCategory=30)
dotplot(gsegmt, showCategory=30)

网络图

########## 网络图：cnetplot ##########
tiff(file="figures/GSEA_KEGG_cnetplot.tiff",
     width = 35,height = 22,units ="cm",
     compression="lzw",bg="white",res=300)
cnetplot(gsekk, categorySize="pvalue", foldChange=geneList)
dev.off()

cnetplot(gsekk, foldChange=geneList)
cnetplot(gsekk, foldChange=geneList, circular = TRUE