6

我有一个包含数百行和数十列的矩阵，并希望绘制热图。如果我使用原生 R 函数：

heatmap(matrix(sample(1:10000),nrow=1000,ncol=10))

我得到一个行标题难以辨认的数字。我假设生成的图像与当前绘图设备的规格相匹配。我想控制行的高度，即使这不能显示在我当前的设备上。简单地写入一个高 PNG 只会在图像上方/下方添加空白：

png( '/looks_the_same_with_whitespace.png', width=500, height=1500)
heatmap(matrix(sample(1:10000),nrow=1000,ncol=10))
dev.off()

有没有一种干净的方法可以做到这一点，也许是通过欺骗 R 以为我有一个非常高的显示器？来自受良好支持的库的函数也是一个很好的答案。

4

1 回答 1

8

前段时间，我也遇到了同样的问题。heatmap.2包中的函数gplots解决了这个问题。但是，我不太明白查看这么多基因（或 rowlabs）有什么好处。但正如你所问的那样，你会做这样的事情：关键是改变热图的layout（参见？heatmap.2，？layout），它在绘图设备上的 2x2 网格上绘制热图（？layout 中的示例更好地解释了这一点）。之后，您可能希望通过相应更改来更改cexof 。您的情况可能需要对这些值进行一些调整才能获得所需的结果。rowlabscexRow

library(gplots)

row.scaled.expr <- matrix(sample(1:10000),nrow=1000,ncol=10)

png( 'heatmap_without_whitespace_smaller_row_lab.png', width=500, height=1500)
heatmap.2(row.scaled.expr, dendrogram ='row',
                 Colv=FALSE, col=greenred(800), 
                 key=FALSE, keysize=1.0, symkey=FALSE, density.info='none',
                 trace='none', colsep=1:10,
                 sepcolor='white', sepwidth=0.05,
                 scale="none",cexRow=0.2,cexCol=2,
                 labCol = colnames(row.scaled.expr),                 
                 hclustfun=function(c){hclust(c, method='mcquitty')},
                 lmat=rbind( c(0, 3), c(2,1), c(0,4) ), lhei=c(0.25, 4, 0.25 ),                 
                 )
dev.off()
于 2011-02-14T06:34:36.140 回答