pheatmap热图

# Create test matrix
test = matrix(rnorm(200), 20, 10)
test[1:10, seq(1, 10, 2)] = test[1:10, seq(1, 10, 2)] + 3
test[11:20, seq(2, 10, 2)] = test[11:20, seq(2, 10, 2)] + 2
test[15:20, seq(2, 10, 2)] = test[15:20, seq(2, 10, 2)] + 4
colnames(test) = paste("Test", 1:10, sep = "")
rownames(test) = paste("Gene", 1:20, sep = "")

Draw heatmaps

pheatmap(test) #默认行名与列名

pheatmap(test, kmeans_k = 2)

pheatmap(test, scale = "row", clustering_distance_rows = "correlation")

pheatmap(test, color = colorRampPalette(c("navy", "white", "firebrick3"))(50))

pheatmap(test, cluster_row = FALSE) #行聚类消失

pheatmap(test, legend = FALSE) #图例消失

Show text within cells

pheatmap(test, display_numbers = TRUE)

pheatmap(test, display_numbers = TRUE, number_format = "%.1e") #去掉斜杠

pheatmap(test, display_numbers = matrix(ifelse(test > 5, "*", ""), nrow(test)))

pheatmap(test, cluster_row = FALSE, legend_breaks = -1:4, legend_labels = c("0", "1e-4", "1e-3", "1e-2", "1e-1", "1"))

Fix cell sizes and save to file with correct size

画板出现问题：dev.off()

pheatmap(test, cellwidth = 15, cellheight = 12, main = "Example heatmap")

pheatmap(test, cellwidth = 15, cellheight = 12, fontsize = 8, filename = “test.pdf”)

为行和列生成注释

Display row and color annotations

annotation_col = data.frame(CellType = factor(rep(c("CT1", "CT2"), 5))

pheatmap(test, annotation_col = annotation_col) #annotation_col

annotation_col = data.frame(CellType = factor(rep(c("CT1", "CT2"), 5)), Time = 1:5
)

连续型的数据渐变色
离散型的数据：用不同的颜色

rownames(annotation_col) = paste("Test", 1:10, sep = "")annotation_row = data.frame(GeneClass = factor(rep(c("Path1", "Path2", "Path3"), c(10, 4, 6)))
)
rownames(annotation_row) = paste("Gene", 1:20, sep = "")

pheatmap(test, annotation_col = annotation_col, annotation_legend = FALSE)

pheatmap(test, annotation_col = annotation_col, annotation_row = annotation_row)

Change angle of text in the columns

pheatmap(test, annotation_col = annotation_col, annotation_row = annotation_row, angle_col = “45”)
pheatmap(test, annotation_col = annotation_col, angle_col = “0”)

Specify colors

ann_colors = list(
Time = c(“white”, “firebrick”),
CellType = c(CT1 = “#1B9E77”, CT2 = “#D95F02”),
GeneClass = c(Path1 = “#7570B3”, Path2 = “#E7298A”, Path3 = “#66A61E”)
)

pheatmap(test, annotation_col = annotation_col, annotation_colors = ann_colors, main = “Title”)

pheatmap(test, annotation_col = annotation_col, annotation_row = annotation_row,
annotation_colors = ann_colors)

pheatmap(test, annotation_col = annotation_col, annotation_colors = ann_colors[2])

Gaps in heatmaps

pheatmap(test, annotation_col = annotation_col, cluster_rows = FALSE, gaps_row = c(10, 14))

pheatmap(test, annotation_col = annotation_col, cluster_rows = FALSE, gaps_row = c(10, 14),
cutree_col = 2)

Show custom strings as row/col names

labels_row = c(“”, “”, “”, “”, “”, “”, “”, “”, “”, “”, “”, “”, “”, “”, “”,
“”, “”, “Il10”, “Il15”, “Il1b”)

pheatmap(test, annotation_col = annotation_col, labels_row = labels_row)

Specifying clustering from distance matrix

drows = dist(test, method = “minkowski”)
dcols = dist(t(test), method = “minkowski”)
pheatmap(test, clustering_distance_rows = drows, clustering_distance_cols = dcols)

Modify ordering of the clusters using clustering callback option

callback = function(hc, mat){
sv = svd(t(mat))$v[,1]
dend = reorder(as.dendrogram(hc), wts = sv)
as.hclust(dend)
}

pheatmap(test, clustering_callback = callback)

Not run:

Same using dendsort package

library(dendsort)

callback = function(hc, …){dendsort(hc)}
pheatmap(test, clustering_callback = callback)

End(Not run)

制作不易，麻烦点赞

pheatmap热图相关推荐

Pheatmap热图的绘制及如何调整图片
Pheatmap热图的绘制及如何调整图片 Pheatmap包是R语言绘制热图比较强大的软件包,当然现在也有很多资料介绍这个包的使用,但是今天我写的重点不是如何使用这个包绘制热图,而是如何绘制出更好看的 ...
r语言热图对列不进行聚类_R语言：手把手教你画pheatmap热图
R语言:手把手教你画pheatmap热图微生态导读: pheatmap默认会对输入矩阵数据的行和列同时进行聚类,但是也可以通过布尔型参数cluster_rows和cluster_cols设置是否对 ...
获取pheatmap热图聚类后和标准化后的结果
pheatmap是简单常用的热图绘制包,可以快速.简单.可定制的绘制漂亮热图.具体见R语言学习-热图简化和免费高颜值可定制在线绘图工具 ImageGP. 现在要解决的一个问题是图出来了,想看下转换后用 ...
热图pheatmap 热图3：热图行列分组信息注释
热图pheatmap library(pheatmap) getwd() setwd("D:\Win10 System\Documents\WeChat Files\wxid_f27yna0 ...
R：热图解释｜ pheatmap包参数及详细聚类图绘制流程（一篇解决热图绘制问题）
热图解释及pheatmap绘制热图一.热图绘制原理 1.1 热图介绍 1.2 热图绘制准备--均一化 1.3 热图绘制方式 1.4 热图数据查看示例二.pheatmap包简介 2.1 pheatm ...
245热图展示微生物组的物种和功能丰度或有无、距离矩阵
245热图展示微生物组的物种和功能丰度或有无本节作者:吴一磊中科院微生物所版本1.0.7,更新日期:2020年8月13日本项目永久地址:https://github.com/YongxinLi ...
R pdf大小_数据可视化 - 热图·R
热图(如下图所示),是一种展示样本与差异变量关系的可视化方式,可以通过R/python进行绘制.本文主要分别具体介绍如何使用R实现热图绘制. R绘制热图通常调用R的软件包"pheatmap ...
r语言聚类分析_技术贴 | R语言pheatmap聚类分析和热图
点击蓝字↑↑↑"微生态",轻松关注不迷路本文由阿童木根据实践经验而整理,希望对大家有帮助. 原创微文,欢迎转发转载. 导读 pheatmap默认会对输入矩阵数据的行和列同时进行聚 ...
服务器画热图显示无效的变量,使用pheatmap包绘制热图
加载所需R包 library(pheatmap) 设置工作路径 setwd("/Users/Davey/Desktop/VennDiagram/") # 清除当前环境中的变量 rm ...