热图6：系列终章之“奇形怪状”的热图

走到这里，也代表热图这个系列到最终篇了。最后我们说说一些奇形怪状的热图。其实小编也发现，生物学现在整的不太好的一个方面就是发论文大家太注重图了，反倒忽略了数据的价值。我是认为用合理的图将数据价值表现出来即可，稍加修饰，无需太过。

但是很多小伙伴总是看到别人文章新奇的图，比较感兴趣，所以这里也示范一下一些新奇的热图，仅限于小编自己见过的！其他的小伙伴们可以留言分享。

一、环状热图

顾名思义，环状热图就是圆形的热图，将长形的热图卷一下就可以了。它唯一的好处可能就是热图太长的时候节省空间！

用circlize函数画图（ggplot2也可以实现，这里我们不再展示，有兴趣可以去网上查查）

首先加载数据并处理：

setwd("D:/物信息学/环形热图")
A <- read.csv("行列注释.csv", header = T, row.names = 1)
A <- as.matrix(A) #将表达矩阵转化为matrix
for (i in 1:nrow(A)) A[i, ] <- scale(log(unlist(A[i, ] + 1), 2)) #对数据进行标准化处理
annotation_row = data.frame(Pathway = c(rep("Wnt",20), rep("Inflammatory",32),rep("HIF",34)))#对行进行注释，用于后续的热图分裂
row.names(annotation_row) <- rownames(A)
annotation_row <- as.matrix(annotation_row)#在circlize函数中，需要为matrix

加载需要的包：

library(circlize)
library(dendextend)

设置颜色，画图：

col_color = colorRamp2(c(-2, 0, 2), c("#003399", "white", "#cccc00"))#设置legend颜色，范围
circos.heatmap(A, col = col_color, dend.side = "inside",#聚类放在环形内测rownames.side = "outside",#基因名放在环形外侧split = annotation_row,#用行注释分裂热图bg.border = "black")#背景边缘颜色
circos.clear()#画完图结束。一定要运行这个，不然后续画图会叠加

但是这样没有行名，也没有分组信息，也没有legend，需要额外添加。


circos.par(gap.after = c(2, 2,30))#让分裂的一个口大一点，可以添加行信息
circos.heatmap(A, col = col_color, dend.side = "inside",rownames.side = "outside",split = annotation_row,bg.border = "black",show.sector.labels = T)
#以下代码用于显示分组，但是很复杂
circos.track(track.index = get.current.track.index(), panel.fun = function(x, y) {if(CELL_META$sector.numeric.index == 3) { # the last sectorcn = colnames(A)n = length(cn)circos.text(rep(CELL_META$cell.xlim[2], n) + convert_x(1, "mm"), 1:n - 0.5, cn, cex = 0.5, adj = c(0, 0.5), facing = "inside")}
}, bg.border = NA)
#添加legend
library(ComplexHeatmap)
lgd = Legend(title = "Scale", col_fun = col_color)
grid.draw(lgd)

分组行名还是需要手动调整。

二、三角热图

总是在某些文章中（可能高分文章作图比较有性格吧）出现不一样的热图，除了上述讲的环状，还有三角热图，可以理解为热图只显示斜对角一半。

利用ComplexHeatmap正常加载数据处理，画热图加一个函数和改变参数就可以实现。


setwd("D:/生物信息学/环形热图")
A <- read.csv("行列注释.csv", header = T, row.names = 1)
A <- as.matrix(A) #将表达矩阵转化为matrix
for (i in 1:nrow(A)) A[i, ] <- scale(log(unlist(A[i, ] + 1), 2)) #对数据进行标准化处理
Heatmap(A, rect_gp = gpar(type = "none"), column_dend_side = "bottom",cell_fun = function(j, i, x, y, w, h, fill) {if(as.numeric(x) <= 1 - as.numeric(y)) {grid.rect(x, y, w, h, gp = gpar(fill = fill, col = fill))}})

这个热图基因表达量的时候真的没有什么意义。不过这种热图我们应该熟悉，一般用在相关系数的表达上，下三角是热图，上三角是数字，不过不是用这个函数实现的。

我们不妨做一下相关性分析及其可视化，或许在这里不会做无用功。

数据还是表达量数据，不过需要原始数据，不在经过做热图标准化那一步。

加载需要的数据包，并计算相关性和P值,画图：


library(ggcorrplot)
library(ggthemes)
#计算这几个样本的相关性
sample_cor <- round(cor(A),3)
#计算P值，看显著性
sample_P <- round(cor_pmat(A),3）ggcorrplot(sample_cor,#计算的相关系数矩阵method = 'circle',#热图显示为圆圈hc.order = T,#聚类hc.method = "ward.D",#聚类方法outline.color = "white",ggtheme = theme_bw(),type = "upper",#显示上三角colors = c("#6D9EC1","white","#E46726"),#热图颜色lab = T,#热图上显示相关系数lab_size = 2.5)

为了让这个图更加有意义，我们只显示P≤0.05的相关，也就是显著的。

ggcorrplot(sample_cor,#计算的相关系数矩阵method = 'circle',#热图显示为圆圈hc.order = T,#聚类hc.method = "ward.D",#聚类方法outline.color = "white",ggtheme = theme_bw(),type = "upper",#显示上三角colors = c("#6D9EC1","white","#E46726"),#热图颜色lab = T,#热图上显示相关系数lab_size = 2.5,p.mat = sample_P,#相关insig = "pch")#不显著的显示为X或者空白

走到这里，我想热图系列已经很圆满很丰富了，我们更多的还是要关注数据的实际意义，以及在适当的时候用合适的热图尽可能的表现数据意义！

系列结束总是给人一种轻松，祝大家学习愉快！

请大家多多支持关注和分享，更多精彩内容请关注我的公众号《KS科研分享与服务》

热图6：系列终章之“奇形怪状”的热图相关推荐

python画笔宽度_python教程系列终章：画笔控制
原标题:python教程系列终章:画笔控制 3. 画笔控制通过前面两节的学习,大家可能会感觉到实际上的编程相当于我们手握一直笔,通过各种指令控制笔的路径. 没错,我们平时写字或画画的各种操作都能够用 ...
《银河护卫队》系列终章IMAX中国首周末劲收3840万票房，全球斩获约2500万美元...
美通社消息,IMAX China宣布由漫威/迪士尼影业出品的<银河护卫队3>在中国上映首周末劲收票房3840万,IMAX票房占比则高达20%.作为该超级英雄系列的终章,<银河护卫队3 ...
python控制画笔尺寸_python教程系列终章：画笔控制
Beetle.penDown() Beetle.move(10) Beetle.penUp() Beetle.move(10) Beetle.penDown() Beetle.move(10) Bee ...
系列终章 - 放眼未来！进击吧， Blazor ！
Blazor 是一个 Web UI 框架,可通过 WebAssembly 在任意浏览器中运行 .Net . Blazor 旨在简化快速的单页面 .Net 浏览器应用的构建过程,它虽然使用了诸如 CSS ...
AD（二十二）本系列终章—PCB从制图到加工的各类注意事项
软件:Altium Designer 16 往期博客: AD(一) AD的工程创建 AD(二)电阻.电容模型的创建 AD(三)IC类元件模型创建 AD(四)排针类元件模型的创建 AD(五)光耦元件模 ...
一篇文章一张思维导图看懂Android学习最佳路线
一篇文章一张思维导图看懂Android学习最佳路线先上一张android开发知识点学习路线图思维导图 Android学习路线从4个阶段来对Android的学习过程做一个全面的分析:Android初级 ...
终章 - 软件工程实践总结作业
终章 - 软件工程实践总结作业一.请回望暑假时的第一次作业,你对于软件工程课程的想象当时正值酷暑,快开学了,在家里也是无所事事,就花费了一些心思写下了那份第一次作业.现在看来,当时的想象还是蛮符合 ...
怎么用matlab画无差别曲线,MATLAB系列第二章初等模型.ppt
<MATLAB系列第二章初等模型.ppt>由会员分享,可在线阅读,更多相关<MATLAB系列第二章初等模型.ppt(68页珍藏版)>请在人人文库网上搜索. 1.第二章初等模型 ...
SLAM导航机器人零基础实战系列：（六）SLAM建图与自主避障导航——2.google-cartographer机器人SLAM建图...
SLAM导航机器人零基础实战系列:(六)SLAM建图与自主避障导航--2.google-cartographer机器人SLAM建图摘要通过前面的基础学习,本章进入最为激动的机器人自主导航的学习.在 ...

热图6：系列终章之“奇形怪状”的热图

热图6：系列终章之“奇形怪状”的热图相关推荐

最新文章

热门文章