我们通过TCGA或者GEO数据库下载基因的表达谱数据,我们需要看一下某些基因在肿瘤组织和正常组织中的表达情况,我们可以通过差异分析后绘制火山图。获得的差异基因,结合下载的临床数据,我们可以进行生存分析【后台回复:TCGA-DEG,查看文章】。我们和还可以绘制热图

下面是我们通过之前从TCGA数据库下载LUAD的表达矩阵(标准化后的FPKM数据),数据格式如下:

我们做的基因有:

METTL3METTL14WTAPKIAA1429RBM15ZC3H13YTHDC1YTHDC2YTHDF1YTHDF2HNRNPCFTOALKBH5

绘制一些基因的热图代码:

## 微信公众号:MedBioInfoCloudoptions(stringsAsFactors = F)## 读入标准化的表达矩阵文件,TCGA直接下载好的HTSeq - FPKM,就可以用exp ".\\ExpMatrix.NT-TP.txt",header = T,                  sep = "\t",check.names = F)## 读入我们要分析的基因,也可以是来自你进行差异表达分析后具有显著差异表达的基因gene ".\\gene.txt",header = F,sep = geneNM 1],rownames(expData exp[geneNM,] conNum 59treatNum #####热图的绘制#install.packages("pheatmap")HeatMapData 1)library(pheatmap)Type "Normal",conNum),rep(names(Type) Type pdf("heatmap.pdf",height=5,width=20) # 图片宽高,根据样本数大小和基因数目修改pheatmap(HeatMapData,          annotation = Type,          color = colorRampPalette(c("green", "white", "red"))(50),#设置颜色         cluster_cols =F,         show_colnames = F,         scale="row",         fontsize = 10,         fontsize_row=10,         fontsize_col=3)dev.off()

差异表达基因热图怎么看_基因表达谱热图绘制相关推荐

  1. 三位数的茎叶图怎么看_如何用茎叶图表示三位数,什么是茎叶图

    茎叶图(Stem-and-Leaf Diagrams)又称"枝叶图",由统计学家约翰托奇设计,它的思路是将数组中的数按位数进行比较,将数的大小基本不变或变化不大的位作为一个主干(茎 ...

  2. 差异表达基因热图怎么看_多变的热图1(新手专用)

    热图(heatmap)用不同的颜色和颜色的深浅来直观的展示数据之间的差异.在测序类的文章里,几乎必有一幅热图用来展示差异表达基因.很多工具都可以完成热图的制作,今天这篇文章主要介绍利用R语言的 phe ...

  3. 差异表达基因热图怎么看_画图专题 | 热图函数pheatmap()

    画图专题 | 热图函数pheatamp() 百迈客医学 2018-10-31 先看一眼这个函数的参数,这么多,而且最后还有省略号.那么我们应该怎么合理使用这些参数让你的热图看起来更加高大上呢? phe ...

  4. sip 时序图_时序图怎么看_教你如何看懂时序图 - 什么是时序图_时序图怎么看_教你如何看懂时序图...

    时序图怎么看_教你如何看懂时序图 操作时序永远使用是任何一片IC芯片的最主要的内容.一个芯片的所有使用细节都会在它的官方器件手册上包含.所以使用一个器件事情,要充分做好的第一件事就是要把它的器件手册上 ...

  5. 差异表达基因热图怎么看_差异基因热图绘制:heatmap.2

    在RNA-seq数据分析中,差异表达基因分析是一项基本的技能,其中热图又是一种特别常见的用来展示差异表达基因分析结果的方式,今天分享一个非常好用的绘制热图的R函数:heatmap.2.该函数来自gpl ...

  6. 差异表达基因热图怎么看_为什么我代码里面选择top1000的sd基因绘制热图呢

    实际上写完了这个全网最好的差异分析代码:免费的数据分析付费的成品代码 我就可以收工了,但是永远不能低估粉丝的疑惑数量,任何一个细节都会被拿出来剖析. 比如代码里面我挑选了top1000的sd基因绘制热 ...

  7. 差异表达基因热图怎么看_热图的解读及边界聚类的意义

    热图(heatmap)是分子生物学文章里(尤其是microarray, RNA-seq相关论文)经常出现的图片.但热图一般有哪些用途,具体涉及哪些不同的参数画法,需要清楚了解. 1. 关于热图的用途 ...

  8. 三位数的茎叶图怎么看_茎叶图怎么看,茎叶图的特点(简要概括)

    茎叶图(Stem-and-Leaf Diagrams)又称"枝叶图",由统计学家约翰托奇设计,它的思路是将数组中的数按位数进行比较,将e68a84e8a2ade799bee5baa ...

  9. 图神经网络代码_第一篇:图神经网络(GNN)计算框架绪论

    写在开头: 这个专栏是为了总结我本科毕业设计中所设计的题目<基于GPU的图神经网络算法库的设计钰实现>.这半年来一直在这个方向上啃代码,读论文,真的学到了很多东西.尤其是阅读了大佬团队写的 ...

  10. 机械动作时序图怎么画_快速学习时序图:时序图简介、画法及实例

    原标题:快速学习时序图:时序图简介.画法及实例 时序图作为常用的UML交互图,可以直观的传达系统内外之间的交互过程,经常用在详细设计文档中.下面本文综合参考了多篇时序图的教程,根据作者的思路将时序图做 ...

最新文章

  1. 万物皆可文本时代来临?如何搞定NLP最强模型GPT
  2. 数据库:MySQL相关知识整理,值得收藏!
  3. 红魔1android版本能升到多少,红魔5G:实力全开,比快更快
  4. java_log_02
  5. python游戏程序中游戏对象是什么_Python游戏编程入门
  6. 罗技G500游戏鼠标
  7. 反向代理服务器域名解析
  8. 直观解释Gamma分布和Poisson分布
  9. Python非线性拟合自定义函数参数(对标MATLAB-nlinfit函数)
  10. go mysql 异步_GO操作MYSQL
  11. 下载json文件,解决浏览器对JSON文件链接直接打开问题
  12. DMCTF writeup
  13. 7-20 | 打印九九口诀表
  14. excel批量文件改名批量加后缀
  15. 构建TCP套接字(socket)的概念及具体步骤
  16. 安防工程商前期收集客户信息15种方法与详细流程
  17. CVE-2022-0847 Linux内核提权漏洞分析
  18. 美国北亚利桑那大学计算机专业排名,北亚利桑那大学排名 综合排名和专业排名介绍...
  19. 笔记:RT-Thread Studio 快速上手
  20. Java SE学习之【网络编程JavaSocket】

热门文章

  1. python中reduce是什么意思_python-reduce函数
  2. 筋膜枪方案-无刷马达方波运用1
  3. 倩女手游怎么查看服务器信息,倩女幽魂手游新增互通服务器 互通服务器信息介绍[图]...
  4. C++/MFC修行之路(6)停靠窗口
  5. kindle3nbsp;注册地址,美国地址…
  6. 第八周、第九周学习总结
  7. 计算机磁盘图标变了怎么改过来,本地磁盘图标变了怎么办 本地磁盘图标变解决方法...
  8. SUV世家广汽三菱,牵手微信企业号打造学习型企业
  9. ueditor php上传word,ueditor百度编辑器上传PDF并显示
  10. 2019,焉知非福,决胜新城