今天复现一幅PNAS文章中的图，是个热图，重点不是这个热图，而是得到做热图的数据(数据代码已上传qq群)。原文及图片如下：
（reference：Convergent molecular, cellular, and cortical neuroimaging signatures of major depressive disorder）
构建数据
其实得到数据才是这一篇文章的重点内容，从原始的GO、KEGG分析结果csv文件中，我们挑选需要的通路，还有基因，构建做这个离散热图的文件。这里涉及到一个批量读入文件的操作，其实我们只有三个文件，这样做实属操作过度，不过这里只是一个例子，当你面对几十个文件要读入的时候，想起来这篇文章，获取有用！
setwd("D:/KS项目/公众号文章/复现PNAS")
#批量读入csv文件，这只是一种方法
library(dplyr)
data <-list.files(path = "D:/KS项目/公众号文章/复现PNAS",
pattern = "*.csv",
full.names = T)
data <- data%>%lapply(read.csv)
Enrich <- data[[1]]
gene_sel <- data[[2]]
pathway_sel <- data[[3]]
#挑选pathway，原图作者提供了20条，我可能输入有误吧，只有16，无所谓能做图即可
Enrich_sel <- Enrich[Enrich$Name %in% pathway_sel$pathway, ]
接下来就是鉴定我们选定的基因是否存在于通路基因中，通过循环构建数据。这样的数据是TRUE和FALSE的一个矩阵。由于后期无法聚类，所以将其转化为0，1矩阵。
#构建数据，就是看我们选择的基因是否存在于选定通路
enrich_pathway <- Enrich_sel$Name
Enrich_df <- gene_sel
for (i in 1:length(enrich_pathway)){
a <- Enrich_sel[i,]$Hit.in.Query.List
a <- as.character(unlist(strsplit(a, split = ",")))
b <- gene_sel$gene_sel %in% a
b <- as.data.frame(b)
colnames(b) <- enrich_pathway[i]
Enrich_df <- cbind(Enrich_df, b)
}
#构建矩阵
rownames(Enrich_df) <-Enrich_df[,1]
Enrich_df <- Enrich_df[,-1]
head(Enrich_df)
#浅尝一下，作个热图，做热图
Enrich_df[Enrich_df == "TRUE"] = 1
Enrich_df[Enrich_df == "FALSE"] = 0
Complexheatmap作图
library(ComplexHeatmap)
Heatmap(t(as.matrix(Enrich_df)),
cluster_rows = T,
cluster_columns = F,
show_column_names = T,
show_row_names = T,
row_names_side = 'left',
column_title = NULL,
heatmap_legend_param = list(
title=' ',
labels=c('Present','Absent'),
labels_gp = gpar(fontsize = 10),
border='black'),
col = c('#F2F2F0','#5A8FCA'),
rect_gp = gpar(col = "grey", lwd = 1),
row_names_gp = gpar(fontsize = 10),
column_names_gp = gpar(fontsize = 10))
我们看到有个缺点就是聚类和行名无法用不同颜色表示，我感觉Heatmap应该是无法直接实现了。后期还发现一个从没用过的做热图函数，还挺有意思的，它就可以实现聚类标注不同颜色，功能挺强大，缺点是参数名称太不友好，不像heatmap让人一看名字就知道是干嘛的，如下。
heatmaply作图
install.packages("heatmaply")
library(heatmaply)
heatmaply(t(as.matrix(Enrich_df)),
k_row=3,
colors=c('#F2F2F0','#5A8FCA'),
show_dendrogram=c(TRUE,FALSE),
column_text_angle=90,
row_dend_left=F,
grid_color ='grey',
hide_colorbar =T,
branches_lwd =0.5)
具体参数感兴趣可自行研究。其实还没有达到完全复现，ggplot是可以很轻松实现的，但这不是我的目的，不弄了。我想很多人和我一样，并不是专业计算机出身的人，所以不要追求代码写的漂亮，只要达到目的就可以。而且，写的过程中要自己思考，代码不是冷冰冰的，那都是逻辑啊！有错误了也是需要有逻辑的去思考去寻找，而不是什么也不看！
觉得分享对你有用的，点个赞、分享一下呗！更多精彩请至我的公众号《KS科研分享与服务》！
