(R语言)-15-KEGG富集分析气泡图

 重点

geom_point()绘制点图

  • 用aes(size=)将点的大小与某一数值的大小连接起来
  • 点的颜色用aes(color=)

修改默认颜色scale_collor_gradient(low="",hign="")

labs(title="",x="",y="")修改标题,不能用title(main,xlab,ylab)

#KEGG功能注释图的绘制
pathway <- read.csv(file.choose(),header = T)#kegg.csv包
View(pathway)
colnames(pathway)
pp <- ggplot(data=pathway,mapping = aes(x=richFactor,y=Pathway))#richFactor做横轴,Pathway做纵轴
pp
#geom_point()绘制点图,AvsB表示点的大小用aes(size=)表示
pp+geom_point(mapping = aes(size=AvsB))
#点的颜色用aes(color=)
pp+geom_point(mapping = aes(size=AvsB,color=Qvalue))
#修改默认颜色scale_collor_gradient(low="",high="")
pp+geom_point(mapping = aes(size=AvsB,color=Qvalue))+scale_color_gradient(low="green",high="red")
#labs(title="",x="",y="")修改标题
pp+geom_point(mapping = aes(size=AvsB,color=Qvalue))+scale_color_gradient(low="green",high="red")+labs(title = "Top20 of pathway enrichment",x="Rich factor",y="Pathway",color="-log10(PValue)",size="Gene Numbers")+theme_bw()

  • 5
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
KEGG富集分析是一种常用的生物信息学分析方法,用于研究基因或蛋白质的功能和通路富集情况。在R语言中,可以使用Bioconductor中的包来进行KEGG富集分析,其中包括KEGGREST和clusterProfiler。 首先,你需要安装并加载这些包。可以使用下面的代码来完成这一步骤: install.packages("BiocManager") BiocManager::install("KEGGREST") BiocManager::install("clusterProfiler") library(KEGGREST) library(clusterProfiler) 接下来,你需要准备好你的基因列表,并使用KEGGREST包中的函数获取基因对应的KEGG通路信息。下面是一个示例代码,你可以根据你的需要进行修改: gene_list <- c("gene1", "gene2", "gene3") # 替换为你的基因列表 kegg_pathways <- keggGet("pathway", "hsa", "list") # 获取KEGG通路列表 gene_pathway <- enrichKEGG(gene = gene_list, organism = "hsa", pvalueCutoff = 0.05) # 进行KEGG富集分析 最后,你可以使用clusterProfiler包中的函数来可视化KEGG富集分析结果,比如绘制富集通路的柱状、网络等。以下是一个绘制柱状的示例代码: barplot(gene_pathway, showCategory = 10) # 显示前10个富集通路 通过以上步骤,你就可以在R语言中进行KEGG富集分析了。请注意,根据你的具体需求,你可能还需要进行一些参数的调整和结果的解释。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* [R语言clusProfiler进行GO与KEGG富集分析](https://blog.csdn.net/Joey_Liu666/article/details/124988292)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] - *3* [生信分析论文套路R语言代码](https://download.csdn.net/download/thtfhtfhvyyy/87244940)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_2"}}] [.reference_item style="max-width: 50%"] [ .reference_list ]

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值