差异表达基因热图怎么看_获得差异表达基因后-基因功能注释
基因功能注释一言以概之,就是symbol--->ENTREZID--->GO/KEGG/GSEA。
第一步:加载包
library(DOSE)
library(GO.db)
library(org.Hs.eg.db)
library(topGO)
library(GSEABase)
library(clusterProfiler)
没安装的请自行安装
第二步:symbol--->ENTREZID
###上调基因注释
gene <- resSig_up$symbol
gene = bitr(gene, fromType="SYMBOL", toType="ENTREZID", OrgDb="org.Hs.eg.db")
org.Hs.eg.db是symbol转换成ENTREZID所依据的数据库,ENTREZID关联着很多文献数据,为后面的功能注释提供线索,所以需要ID转换。
第三步:GO和KEGG注释
###细胞组分
ego_CC <- enrichGO(gene = gene$ENTREZID,
OrgDb= org.Hs.eg.db,
ont = "CC",
pAdjustMethod = "BH",
minGSSize = 1,
pvalueCutoff = 0.01,
qvalueCutoff = 0.01,
readable = TRUE)
###生物过程
ego_BP <- enrichGO(gene = gene$ENTREZID,
OrgDb= org.Hs.eg.db,
ont = "BP",
pAdjustMethod = "BH",
minGSSize = 1,
pvalueCutoff = 0.01,
qvalueCutoff = 0.01,
readable = TRUE)
###分子功能:
ego_MF <- enrichGO(gene = gene$ENTREZID,
OrgDb= org.Hs.eg.db,
ont = "MF",
pAdjustMethod = "BH",
minGSSize = 1,
pvalueCutoff = 0.01,
qvalueCutoff = 0.01,
readable = TRUE)
pdf(file="ego_MF.pdf")
barplot(ego_MF, showCategory=20,title="EnrichmentGO_MF", drop=T)#条状图,按p从小到大排的
dev.off()
pdf(file="ego_BP.pdf")
dotplot(ego_BP,title="EnrichmentGO_BP_dot")#点图,按富集的数从大到小的
dev.off()
也可以一步到位:
go <- enrichGO(gene = gene$ENTREZID, OrgDb = org.Hs.eg.db, ont='ALL',
pAdjustMethod = 'BH',
pvalueCutoff = 0.01,
qvalueCutoff = 0.01,keyType = 'ENTREZID'
)
dim(go)
dim(go[go$ONTOLOGY=='BP',])
dim(go[go$ONTOLOGY=='CC',])
dim(go[go$ONTOLOGY=='MF',])
gbarplot(go,showCategory=20,drop=T)
dotplot(go,showCategory=50)
###KEGG分析
kk <- enrichKEGG(gene = gene$ENTREZID,
organism ="human",
pvalueCutoff = 0.01,
qvalueCutoff = 0.01,
minGSSize = 1,
#readable = TRUE ,
use_internal_data =FALSE)
pdf(file="barplot_kk.pdf")
barplot(kk)
dev.off()
pdf(file="dotplot_kk.pdf")
dotplot(kk)
dev.off()
差异表达基因热图怎么看_获得差异表达基因后-基因功能注释相关推荐
- 差异表达基因热图怎么看_陈根:从基因层面看衰老本质,被误会的端粒差异
文/陈根 衰老令人难以理解. 自然界几十亿年演化出的身体,为何会迎来衰老和死亡?整个20世纪,科学家都试图发现衰老的秘密. 从衰老的本质来说,衰老的发生一般是从微小的基因层面上开始的.这与染色体和染色 ...
- 差异表达基因热图怎么看_多变的热图1(新手专用)
热图(heatmap)用不同的颜色和颜色的深浅来直观的展示数据之间的差异.在测序类的文章里,几乎必有一幅热图用来展示差异表达基因.很多工具都可以完成热图的制作,今天这篇文章主要介绍利用R语言的 phe ...
- 差异表达基因热图怎么看_画图专题 | 热图函数pheatmap()
画图专题 | 热图函数pheatamp() 百迈客医学 2018-10-31 先看一眼这个函数的参数,这么多,而且最后还有省略号.那么我们应该怎么合理使用这些参数让你的热图看起来更加高大上呢? phe ...
- 差异表达基因热图怎么看_差异基因热图绘制:heatmap.2
在RNA-seq数据分析中,差异表达基因分析是一项基本的技能,其中热图又是一种特别常见的用来展示差异表达基因分析结果的方式,今天分享一个非常好用的绘制热图的R函数:heatmap.2.该函数来自gpl ...
- 差异表达基因热图怎么看_为什么我代码里面选择top1000的sd基因绘制热图呢
实际上写完了这个全网最好的差异分析代码:免费的数据分析付费的成品代码 我就可以收工了,但是永远不能低估粉丝的疑惑数量,任何一个细节都会被拿出来剖析. 比如代码里面我挑选了top1000的sd基因绘制热 ...
- 差异表达基因热图怎么看_热图的解读及边界聚类的意义
热图(heatmap)是分子生物学文章里(尤其是microarray, RNA-seq相关论文)经常出现的图片.但热图一般有哪些用途,具体涉及哪些不同的参数画法,需要清楚了解. 1. 关于热图的用途 ...
- 差异表达基因热图怎么看_【生信】憨憨版分析转录组基因表达差异(二)
[简单热图] # 调用下载的pheatmap包.library(pheatmap)#设置当前工作目录,setwd("目标路径")即保存待测文件的目录,使R能够调用.setwd(&q ...
- 差异表达基因热图怎么看_四个基因的能量代谢模型构建与验证,轻松发4分sci!...
今天分享一篇2019年6月发表于J Cell Physiol(IF:4.522)的文章,标题是Identification of a four‐gene metabolic signature pre ...
- 差异表达基因热图怎么看_学徒作业我想看为什么这几个基因的表达量相关性非常高...
不知道多少人看完了全网第一个单细胞课程(免费基础课程),免费学习地址在B站:https://www.bilibili.com/video/av38741055 而且完成结业考核20题:https:/ ...
最新文章
- 【Alpha】开发日志Day4-0715
- css实验内容,12个令人惊叹的CSS实验项目
- 《Angular从零到一》导读
- 93后阿里P7晒出工资单:原来是狠补了这个~真香
- 【题意+解析】1041 Be Unique (20 分)_18行代码AC
- Linux系统编程(七)消息队列
- 使用HTML5实现刮刮卡效果
- 开源php ajax 图片编辑器 Open source PHP-based Ajax image editor
- 大数据workshop:《云数据·大计算:海量日志数据分析与应用》之《社交数据分析:好友推荐》篇...
- context c语言作用,理解 Go context
- 熊猫DataFrame apply()示例
- Android-JNI开发系列《五》局部引用全局引用全局弱引用缓存策略
- Windows Xp 优化文件 的一段BAT代码
- 扫雷游戏网页版_如何成为扫雷高手,数学老师教你玩转这款复古小游戏
- 真正的落雷(打印图形练习题)C语言
- 21秋期末考试建设工程法规10221k1
- 训练集、验证集、测试集划分
- 血栓清道夫机器人_血栓的清道夫终于找到了!每天吃一点,血管通畅,血栓消失了...
- linux虚拟桌面设置不同的背景,虚拟桌面增强器为Windows 10中的每个虚拟桌面设置不同的壁纸 | MOS86...
- c++面向对象的学生管理系统