画图专题 | 热图函数pheatamp()

百迈客医学 2018-10-31

先看一眼这个函数的参数，这么多，而且最后还有省略号。那么我们应该怎么合理使用这些参数让你的热图看起来更加高大上呢？

pheatmap(mat, color = colorRampPalette(rev(brewer.pal(n = 7, name =

"RdYlBu")))(100), kmeans_k = NA, breaks = NA, border_color = "grey60",

cellwidth = NA, cellheight = NA, scale = "none", cluster_rows = TRUE,

cluster_cols = TRUE, clustering_distance_rows = "euclidean",

clustering_distance_cols = "euclidean", clustering_method = "complete",

clustering_callback = identity2, cutree_rows = NA, cutree_cols = NA,

treeheight_row = ifelse((class(cluster_rows) == "hclust") || cluster_rows,

50, 0), treeheight_col = ifelse((class(cluster_cols) == "hclust") ||

cluster_cols, 50, 0), legend = TRUE, legend_breaks = NA,

legend_labels = NA, annotation_row = NA, annotation_col = NA,

annotation = NA, annotation_colors = NA, annotation_legend = TRUE,

annotation_names_row = TRUE, annotation_names_col = TRUE,

drop_levels = TRUE, show_rownames = T, show_colnames = T, main = NA,

fontsize = 10, fontsize_row = fontsize, fontsize_col = fontsize,

display_numbers = F, number_format = "%.2f", number_color = "grey30",

fontsize_number = 0.8 * fontsize, gaps_row = NULL, gaps_col = NULL,

labels_row = NULL, labels_col = NULL, filename = NA, width = NA,

height = NA, silent = FALSE, na_col = "#DDDDDD", ...)

1、数据准备

#用来画热图的数据

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

部分数据显示如下：

2、画图

如使用默认参数画出来的热图，是不是很不好看？基因名和样本名乱成一堆，也看不出来那些样本聚类到了一起...

3、参数调整

#颜色参数

我们现在将其色调改为红白蓝，看下图，是不是立马好看多了？不知道该选什么颜色？网址颜色大全https://www.color-hex.com/color-names.html送给你。

#数据变换参数

我们现在将数据进行行均一化，聚类方法选择complete，进行均一化之后的热图还没有前面的好看，查看数据表达情况，我们发现最小的表达值3.149和最大的表达值13.533之间并没有相差很大，因此这套数据其实不适合做均一化处理，画出来的热图更好看。

#色度条

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

我们对其进行字符修改，可以看到色度条上相应的数字变成了表示颜色高低表达的字符。

#注释条

这里我们以列注释条为例。在注释条中，这里可以看到正常样本的大类中混有几个癌症样本，癌症大类中也混入了几个正常样本？你有没有想过是为什么呢？当然，你也可以选择不对列进行聚类，这样样本的顺序可以随你自己决定怎样放在一起。

#其他修改参数

我们修改了字体大小，不显示样本名，以及设置图片标题，到此时，经过一系列参数更改，是不是已经比一开始使用默认参数好看多了呢？(可以增加一个对比图，GSE19804_tumor_vs_normal_heatmap6.png与GSE19804_tumor_vs_normal_heatmap1.png的 对比)。

当然还有一些其他的用到不多的参数，留给读者自己去实验一下吧...

#小格子参数设置--热图是由一个个的小四方格子组成的，每一个小格子代表一个基因在一个样本内的表达情况。

#输出文件参数设置，一般可以直接将画好的热图以png格式或者pdf格式进行写出。