pythonic生物人

审稿人:就是这样,必须给minor revisions!

一图胜千言,好的图表让审稿人如沐春风!

本期主题-Heatmap美化。

Image

数据

使用“1.4.9 gene数据集”,

Image

Heatmap美化

美化前,
美化前,

这是一张聚类热力图(Cluster heatmap),它利用聚类算法对图中的行和列进行重新排序,将相似的数据点或变量聚集在一起,在行列边界区域添加了聚类的谱系图。通过聚类热图,既可以观察数据矩阵的分布情况,也可以分析聚类结果。

但是,仅仅通过颜色深浅来判断数据关联,太不直观,尝试拯救这张图:

## 17.2_01
# 未经过允许禁止转载!
# @Author  :  公众号:pythonic生物人
# @Desc    :  17.2 个性化Cluster heatmap(数值、聚类树、文本)

library(pheatmap)

# 绘制热图
pheatmap(
  gene_data, # 数据矩阵
# 🔸标准化设置
  scale = "row", # 按行标准化(可选,突出基因间的相对差异)
# 🔸聚类设置
  cluster_rows = TRUE, # 对基因(行)进行聚类
  cluster_cols = TRUE, # 对样本(列)进行聚类
  clustering_distance_rows = "euclidean", # 行聚类距离:欧几里得距离
  clustering_distance_cols = "euclidean", # 列聚类距离:欧几里得距离
# 🔸轴文本设置
  main = "Gene Expression Heatmap", # 添加标题
  fontsize = 10, # 标题大小设置
  fontsize_row = 10, # 调整行名(基因名)字体大小
  fontsize_col = 10, # 调整列名(样本名)字体大小
# 🔸颜色设置
  color = colorRampPalette(brewer.pal(11, "Spectral"))(50),
  border_color = "grey60", # 设置单元格边框颜色
# 🔸图中文本设置
  display_numbers = TRUE,
  number_format = "%.2f"
)
Image

主要在以下方面进行美化,使其符合顶刊美学:

  • 🔸按行标准化,突出Gene间的相对差异;
  • 🔸聚类树设置,聚类距离、行/列聚类、聚类算法、聚类树属性;
  • 🔸轴文本设置,标题、行/列名属性;
  • 🔸颜色设置,设置单元格边框颜色;
  • 🔸图中文本设置。

挑选一款喜欢的配色,

Image

获取绘图代码+测试数据+数据处理方法+依赖包安装方法,

👇(备注:299)

图片