多分组热图不用愁,Pheatmap来帮忙

【Pheatmap 绘制多分组热图】

热图作为实验数据分布情况的直观展示方法,已经成为高分文章的不错选择,它不仅可以对数据质量进行具像化展示,还可以对数据和样品进行聚类。在R中有多个包均可绘制热图,今天小锐给大家介绍一个可以绘制多个分组热图的R包, libarary(pheatmap)。

福特锐界PLUS 豪横实验室

广告

首先是准备工作

install.packages('pheatmap')#安装R包

library(pheatmap)#加载R包

展开剩余90%

data=read.table('for_plot_filter.txt',header=TRUE,row.names=1,sep='t',check.names=F,quote='')#加载数据。

数据文件的示例图

画个基础图形

pheatmap(data,filename = 'test.png') #绘制图形,并保存文件。

聚类调整

在默认情况下,生成的热图是含有行列聚类,如果不想显示聚类,可加入参数cluster_row = FALSE, cluster_col = FALSE,样品会按照输入顺序进行绘制。

pheatmap(data, cluster_row = TRUE, cluster_col = FALSE,filename = 'test.png')#行聚类,列不聚类

绘制单个分组信息

绘制单个分组信息,输入分组信息annotation_col添加列的注释信息。

colgroup=read.table('group.txt',sep='t',header=F,row.names=1,check.names=F,quote='')

colnames(colgroup)=c('Group')#分组文件设置列名

pheatmap(data,annotation_col = colgroup, filename = 'test.png')

左侧图为分组文件内容

绘制多个分组信息

绘制多个分组信息,annotation_col和annotation_row参数同时添加行和列的注释信息

colgroup=read.table('group.txt',sep='t',header=F,row.names=1,check.names=F,quote='')

rowgroup=read.table('taxon.txt',sep='t',header=F,row.names=1,check.names=F,quote='')

colnames(colgroup)=c('Group','Time','Other')#分组文件设置列名

colnames(rowgroup)=c('Class')

pheatmap(data,annotation_col = colgroup,

annotation_row=rowgroup, filename = 'test.png')

自定颜色列表

利用ann_colors自定注释信息的颜色列表

col = colorRampPalette(c('lightblue', 'yellow', 'orange', 'red'),bias=3)(300)#设置数据颜色

ann_colors = list( Time = c('white', 'firebrick'), Group = c(B ='#00447E',L= '#F34800'),

Class = c(taxon1 = '#64A10E', taxon2 = '#795EA2', taxon3 = '#3370CC') )#设置注释信息颜色

pheatmap(data,annotation_col = colgroup, annotation_row=rowgroup,

col=col,annotation_colors = ann_colors, filename = 'test.png')

结果分块

利用cutree_col,将列聚类结果分块要求队列聚类,gaps_row参数将行结果分块,要求对行不聚类

pheatmap(data, cluster_row = FALSE, cluster_col = TRUE, annotation_col = colgroup,

annotation_row=rowgroup, col=col, gaps_row = c(5, 10),#在5和10行添加分隔

annotation_colors = ann_colors, cutree_col=2,#聚类结果分成两类

filename = 'test.png')

其他参数

sacle: 按照'row”、 'column'或者'none'对数据进行归一化

legend_breaks:设定图例显示范围

legend_labels:设置图例标签

legend = FALSE:去掉图例

annotation_legend = FALSE:去掉注释图例

border_color:设定每个热图格子的边框色

display_numbers = TRUE:设定在每个热图格子中显示相应的数值,

number_color:设置数值字体的颜色

cellwidth和cellheight:设定每个热图格子的宽度和高度

main:添加主标题

fontsize_row 和fontsize_row:设置行名列名文字大小

下课铃响

今天的教学就到这咯

下期的生信课堂我们再见~

> 一款超级简单的热图绘制软件—HemI

> 还不会绘制好看的Circos热图?快来了解一下

> 代码放送 | 热图绘制函数heatmap.2()的参数解析

供稿 | 宋瑞雪

编辑 | 鲁淑妮

暗中观察

默默关注

首赞 +1
平台声明
(0)

相关推荐

  • R语言学习系列之“多变的热图”

    咱公众号也不能只做一个系列,所以经过深思熟虑,打算将来慢慢增加一些内容,主要有以下几个系列TCGA数据分析系列GEO数据分析系列"老板给一个基因,我该怎么办"系列文献阅读系列R语言 ...

  • 为什么我代码里面选择top1000的sd基因绘制热图呢

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

  • 16s分析之差异展示(热图)

    前两天我向大家推了16s做差异分析的两个包(没有看的请点击下面链接): 1.16s分析之差异分析(DESeq2) 2.16s分析之差异OTU 挑选(edgeR) 差异做出来了如何展示,也是一个值得思考 ...

  • R语言代码相关疑问标准提问

    关于如何提问,如何高效沟通,其实我们讲解了非常多了,比如我一直推崇的邮件交流:如果你希望我回答你的问题 ,然后也会随机抽取粉丝提问进行解答:答读者问第一弹:R里面差异分析的limma包用法细节 .也高 ...

  • 使用NMF代替层次聚类

    前面我们在教程:使用R包deconstructSigs根据已知的signature进行比例推断,顺利的把508个病人,根据11个signature进行了比例推断,得到的比例矩阵以普通的热图,以及phe ...

  • 快来使用ggheatmap强化你的热图吧!

    创作原因 用法 参数 可视化 综合示例 结语 创作原因 目前最为常见的热图绘制R包,主要包括pheatmap和ComplexHeatmap(仅个人使用习惯).它们强大的功能,基本可以满足所有科研人员的 ...

  • 技术贴 | R语言:构建一个转录代谢互作调控网络:(二)热图的美化以及大样本分组信息的快速注释

    本文由可爱的乔巴根据实践经验而整理,希望对大家有帮助. 原创微文,欢迎转发转载. 导读 上期介绍了利用WGCNA包中的Cor函数和corPvalueStudent函数计算两组小样本的相关性并进行热图可 ...

  • ComplexHeatmap绘制热图详细步骤

    pheatmap是一个非常受欢迎的绘制热图的R包.ComplexHeatmap包即是受之启发而来.你可以发现Heatmap()函数中很多参数都与pheatmap()相同.在pheatmap的时代(请允 ...

  • 技术贴 | R语言:手把手教你画pheatmap热图

    导读: pheatmap默认会对输入矩阵数据的行和列同时进行聚类,但是也可以通过布尔型参数cluster_rows和cluster_cols设置是否对行或列进行聚类,具体看分析需求.利用display ...

  • 天热,不用愁吃啥,8道凉菜,每天换着做,省事好做还好吃

    天热,就容易犯愁吃什么.和大家分享8道小凉菜做法,简单省事又爽口开胃,很适合夏天哦~ 第一道:菠菜粉丝 菠菜,在焯水之后,质地变得更为鲜嫩也更易入味,与粉丝搭配在一起,一绿一白,颜色一看就很清爽诱人, ...

  • 美食:天冷天热,不用愁吃啥,8道凉菜,每天换着做,省事好做还好吃

    第一道:菠菜粉丝 菠菜,在焯水之后,质地变得更为鲜嫩也更易入味,与粉丝搭配在一起,一绿一白,颜色一看就很清爽诱人,再拌一些调味料,每一口都特别开胃下饭. 具体制作方法 第一步:准备食材.菠菜1小把.干 ...

  • 微信表情包居然还能这么玩,斗图不用愁!

    将「雷科技Lite」收藏为我的小程序,不再错过精彩内容 小雷公司有个运营妹纸是表情包大神,无论你发什么文字都能用表情包回答你~ 这就算了,还时不时的发表情包整蛊我. 除了这种屏幕脏了,怎么擦都擦不掉的 ...

  • 美食:天冷心热,不用愁吃啥,8道凉菜,每天换着做,省事好做还好吃

    天冷.天热,就容易犯愁吃什么.和大家分享8道小凉菜做法,简单省事又爽口开胃,很适合夏天哦~ 第一道:菠菜粉丝 菠菜,在焯水之后,质地变得更为鲜嫩也更易入味,与粉丝搭配在一起,一绿一白,颜色一看就很清爽 ...

  • 老年人白发不用愁 偏方来帮忙

    中医认为与精血亏虚.气滞血瘀有关.1.当归12克.白芍12克.川芎10克.熟地12克.白果5个.首乌9克.桑叶10克,水煎服,每日一剂.2.生地15克.白薇15克.赤芍12克.丹皮12克.地骨皮12克 ...

  • 小学英语:背熟这16张图,3年词汇不用愁,5天牢记7000词!收藏

    小学英语背熟这16张单词分类表,低年级的词汇不用愁!全部都在这里了! 王老师推出了二年级数学下册所有必考的专项练习 只要9.9 赶紧扫一扫下单购买吧!

  • 技术贴 | R语言pheatmap聚类分析和热图

    本文由阿童木根据实践经验而整理,希望对大家有帮助. 原创微文,欢迎转发转载. 导读 pheatmap默认会对输入矩阵数据的行和列同时进行聚类,但是也可以通过布尔型参数cluster_rows和clus ...