如何利用想要基因(感兴趣的GO条目)富集的BP分类画一张气泡图

软件:Rstudio

文件:GO富集结果,感兴趣的GO条目或者基因(本文以基因为例)

1,准备一个基因,将其在注释文件中所含有所有GOs号找出来。

 2,准备好GO富集结果

3,打开Rstudio

library(readr)
library(tidyverse)
library(org.Ji.eg.db)
library(clusterProfiler)
#建立索引,AVG是找到的所有的GO号
index_IL_BP <- which(AVG$ID %in% avgvscon$GO)
#选出15个GO条目,作为showCategory窗口的文件
description_IL_BP <- AVG_CONgo@result$Description[index_IL_BP][1:15]
#dotplot画图,label_format,可以使得较长的term一行展示50个字。
dotplot(AVG_CONgo,showCategory = description_IL_BP,title ="AVGvsCON_IL_BP",label_format=50)
barplot(AVG_CONgo,showCategory = description_IL_BP,title ="AVGvscon_IL_BP",label_format=50)

4,结果:

 

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值