library(ReactomePA)
library(tidyverse)
library(data.table)
library(org.Hs.eg.db)
library(clusterProfiler)
library(biomaRt)
library(enrichplot)
# 基因名转换
gene_map <- select(org.Hs.eg.db, keys=rownames(Myo_DE), keytype="SYMBOL", columns=c("ENTREZID"))
Myo_DE$ENTREZID <- gene_map$ENTREZID
Myo_DE_select <- Myo_DE %>%
arrange(desc(avg_log2FC)) %>%
filter(avg_log2FC > 0.5)
geneList <- Myo_DE_select$avg_log2FC
names(geneList) <- Myo_DE_select$ENTREZID
#GSEA分析——GO
Go_gseresult <- gseGO(geneList, 'org.Hs.eg.db', keyType = "ENTREZID", ont="all", nPerm = 1000,
minGSSize = 10, maxGSSize = 1000, pvalueCutoff=1)
#GSEA分析——KEGG
KEGG_gseresult <- gseKEGG(geneList, nPerm = 1000, minGSSize = 10, maxGSSize = 1000, pvalueCutoff=1)
#GSEA分析——Reactome
Go_Reactomeresult <- gsePathway(geneList, nPerm = 1000, minGSSize = 10, maxGSSize = 1000, pvalueCutoff=1)
#波浪图
ridgeplot(Go_gseresult,10)
ridgeplot(KEGG_gseresult,10)
ridgeplot(Go_Reactomeresult,20)
gseaplot(Go_Reactomeresult,1,pvalue_table = TRUE)
gseaplot2(Go_Reactomeresult,1,pvalue_table = TRUE)
gseaplot2(Go_Reactomeresult, 1:4, pvalue_table = TRUE)
生信笔记16-转录组下游分析之：GSEA
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
- 文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
推荐阅读更多精彩内容
- 对基因counts进行校正定量一般有RPKM、FPKM 、TPM和CPM这几种方法，StatQuest网站中对RP...
- WGCNA是什么 WGCNA全称为weighted gene co-expression network anal...
- 第一步，对所有样本进行中心化 第二步，求特征协方差矩阵 第三步，求协方差矩阵的特征值和特征向量 第四步，将特征值按...
- pearson相关系数：适用于连续性变量，且变量服从正态分布的情况，为参数性的相关系数。 spearman相关系数...
