工具 | 用R快速制作heatmap图,一学就会!

浏览: 2149

摘要: 这个heatmap图是不是很炫?知道是怎么做出来的吗?想学跟我来吧!

Clipboard Image.png

1、首先下载Rhttps://www.r-project.org/

Clipboard Image.png

根据自己的系统选择适合的版本安装

2、安装pheapmap程序包

由于我的是mac OS X版本,下面以Mac OS X版本示范安装pheapmap程序包。

Clipboard Image.png

Clipboard Image.png

Clipboard Image.png

Clipboard Image.png

打开R,菜单里面选择Packages & Data 下面的Package Installer,如果没有任何下面的Package选择,选择Get List,选择离我们最近的站点。在左侧Package Search里面搜索pheatmap,点击install selected。

3、准备数据文件

在Excel里面输入所需要画heatmap图数据,第一列为名称,后面的为均一化之后的数值。二代测序的表达量差异的结果,芯片表达差异的结果,甚至是qPCR的结果都可以。

Clipboard Image.png

Clipboard Image.png

然后将整理好的Excel文件另存为制表符分割的文本.txt,命名为all.txt。

Clipboard Image.png


打开all.txt确定一下,是否正确。

Clipboard Image.png

4 、制作heatmap图

启动R,先设置工作目录。选择菜单Misc里面的Change Working Directory,选择存放要做图的数据文件all.txt所在的文件夹。

复制输入以下代码,注意all.txt与文件名一致,miRNA_name与Excel里面第一列的名称一致,输出pdf的名称可以随意

data<-read.table("all.txt",head=T) #读入数据文件

attach=data

row.names(data)<-data$miRNA_name #提取行名

data<-data[,-1] #去除数据中的行名

data<-data.matrix(data) #转换为矩阵

library(pheatmap) #加载pheatmap函数

pdf("test.pdf", height=10, width=10) #将绘图输出到PDF

pheatmap(data,cluster_cols=FALSE,clustering_distance_row="correlation",clustering_method="complete",color=colorRampPalette(c("green","black","red")(100),revC=FALSE,scale="row",margins=c(5,10),fontsize_row=8,cellheight=10,cellwidth=30)  #默认参数

dev.off()

Clipboard Image.png

回车,在设置的工作目录里找test.pdf文件,即是生成的heatmap图。

其中输出的pdf("test.pdf", height=10, width=10),height, width可以根据图片大小自己调整

Clipboard Image.png

Clipboard Image.png

推荐 0
本文由 大数据人 创作,采用 知识共享署名-相同方式共享 3.0 中国大陆许可协议 进行许可。
转载、引用前需联系作者,并署名作者且注明文章出处。
本站文章版权归原作者及原出处所有 。内容为作者个人观点, 并不代表本站赞同其观点和对其真实性负责。本站是一个个人学习交流的平台,并不用于任何商业目的,如果有任何问题,请及时联系我们,我们将根据著作权人的要求,立即更正或者删除有关内容。本站拥有对此声明的最终解释权。

0 个评论

要回复文章请先登录注册