是否可以制作200k行的热图?

时间:2015-10-15 15:13:15

标签: r heatmap pheatmap

我想知道是否可以制作200k行的热图?我有一个矩阵,其基因组坐标为行,每列代表每个患者中该区域的存在与否。所以我有15个病人(专栏)。现在,当我尝试使用heatmap.2或pheatmap时,我得到内存分配问题,我将如何使用整个矩阵生成热图。我的矩阵的值是jsut 0和1,我想根据热图绘制一些假设。我怎样才能使用它。我在我的笔记本电脑上尝试了一次,但它没有帮助,所以我试着在我们的Linux集群中非常强大。但仍显示以下错误。我该如何解决这个问题。我还添加了sessionInfo()。感谢任何解决方法

data<-read.delim("path/H3_marks_map.txt",sep="\t", row.names=1)
pdf(file="path/map/H3_maps1.pdf")
pheatmap(data,scale="none")
Error: cannot allocate vector of size 170.5 Gb
dev.off()

sessionInfo()                           
R version 3.1.2 (2014-10-31)
Platform: x86_64-unknown-linux-gnu (64-bit)

locale:
 [1] LC_CTYPE=C             LC_NUMERIC=C           LC_TIME=C.UTF-8       
 [4] LC_COLLATE=C.UTF-8     LC_MONETARY=C.UTF-8    LC_MESSAGES=C.UTF-8   
 [7] LC_PAPER=C.UTF-8       LC_NAME=C              LC_ADDRESS=C          
[10] LC_TELEPHONE=C         LC_MEASUREMENT=C.UTF-8 LC_IDENTIFICATION=C   

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] RColorBrewer_1.1-2 pheatmap_1.0.7    

loaded via a namespace (and not attached):
[1] colorspace_1.2-6 grid_3.1.2       gtable_0.1.2     munsell_0.4.2   
[5] plyr_1.8.3       Rcpp_0.12.1      scales_0.3.0     tools_3.1.2 

0 个答案:

没有答案