热图pheatmap()函数

微点

修改于 2019-11-01 17:59:24

3.3K0

修改于 2019-11-01 17:59:24

文章被收录于专栏：sci666sci666

先看一眼这个函数的参数，这么多！而且最后还有省略号！

那么我们应该怎么合理使用这些参数让你的热图看起来更加高大上呢？

数据准备：

此次例子，我们选择了一套GEO数据库的肺癌数据，数据编号为GSE19804，120个样本，其中包含60个癌症样本和60个癌旁正常样本，前面我们使用t检验，并对p值进行BH校正，筛选fdr小于0.01的基因中前40个在癌症相对于正常样本中显著差异表达的基因进行热图绘制。

数据格式如下：

默认参数画图：

#绘图 pheatmap(heatmap_data)

是不是很不好看？基因名和样本名乱成一堆，也看不出来那些样本聚类到了一起…

参数调整：

#颜色参数：

color	表示颜色，用来画热图的颜色，可以自己定义，默认值为colorRampPalette(rev(brewer.pal(n = 7, name =”RdYlBu”)))(100)，RdYlBu也就是Rd红色，Yi黄色，Bu蓝色的过度，则主调色为红黄蓝。

#数据变换参数：

scale	是指对数值进行均一化处理，在基因表达量的数据中，有些基因表达量极低，有些基因表达量极高，因此把每个基因在不同处理和重复中的数据转换为平均值为0，方差为1的数据，可以看出每个基因在某个处理和重复中表达量是高还是低，一般选择做row均一化。
clustering_method	表示聚类方法，值可以是hclust的任何一种，如”ward.D”,”single”, “complete”, “average”, “mcquitty”, “median”, “centroid”, “ward.D2″。
cluster_rows	表示行是否聚类，值可以是FALSE或TRUE
clustering_distance_rows	行距离度量的方法，如欧氏距离
cutree_rows	行聚类数
treeheight_row	行聚类树的高度，默认为50
gaps_row	对行进行分割，就不应对相应的行进行聚类
cluster_cols	表示列是否聚类，值可以是FALSE或TRUE
clustering_distance_cols	列距离度量的方法
cutree_cols	列聚类数
treeheight_col	列聚类树的高度，默认为50
gaps_col	对列进行分割，就不应对相应的列进行聚类

#色度条参数

就是热图右上角那个小小的长方条

legend	逻辑值，是否显示色度条，默认为T
legend_breaks	显示多少个颜色数值段
legend_labels	对色度条上对应位置的字符进行修改

#注释条

annotation_colors	对标签的颜色进行修改
annotation_legend	是否显示标签注释条
annotation_row	数据框格式，用来定义热图所在行的注释条
annotation_names_row	逻辑值，是否显示行标签名称
annotation_col	数据框格式，用来定义热图所在列的注释条
annotation_names_col	逻辑值，是否显示列标签名称

#其他修改参数

main	设置图的标题
fontsize	是设置所有除主图以外的标签的大小
number_color	字体的颜色
show_rownames	是否显示行名
fontsize_row	行名的字体大小
labels_row	X轴坐标名设置
show_colnames	是否显示列名
fontsize_col	列名的字体大小
labels_col	y轴坐标名设置

经过小编一系列参数更改，修改如下：

#加载包

library(pheatmap) library(RColorBrewer)

#设置颜色

color.1 <- colorRampPalette(rev(c(“#ff0000”, “#000000”, “#00ff00”)))(100)

#注释条

#设置列的注释条

anno_col=data.frame(sampleType=factor(rep(c(“tumor”,”normal”),each=60)))

#与数据样本名对应 row.names(anno_col)=colnames(heatmap_data)#设置列的注释条的颜色 ann_color=list(sampleType=c(tumor=’#cd0000′,normal=’#3a5fcd’))

#绘图

pheatmap(heatmap_data,scale=”none”,color=color,annotation_col=anno_col,annotation_colors=ann_color,fontsize=8,show_colnames=F,main=”GSE19804:tumor vs normal heatmap”,legend_breaks=c(4,12),legend_labels=c(‘low’,’high’))

是不是已经比一开始使用默认参数好看多了呢？如下：

当然还有一些其他的用到不多的参数

留给读者自己去实验一下吧…

#小格子参数设置

热图是由一个个的小四方格子组成的，每一个小格子代表一个基因在一个样本内的表达情况

fontsize_number	小格子中数字大小
display_numbers	逻辑值，是否在小格子中显示数字
number_format	小格子中数字显示形式，但仅有在display_numbers=T时才能使用
na_col	设置小格子为缺失值时的颜色
cellwidth	表示每个小格子的宽度
cellheight	表示每个小格子的高度