检查heatmap.2 代码设计是不可能的。 heatmap.2 调用 image 来绘制热图,并将额外的关键字 args (...) 传递给它,因此调用 heatmap.2(..., cex.lab = n, ...) 会成功,而 image 会考虑 par,所以也可以这样可能的。但是这里xlab和ylab被替换为空字符串,并且刻度标签和轴标签都在后面绘制:
# from gplots heatmap.2 source:
image(1:nc, 1:nr, x, xlim = 0.5+ c(0, nc), ylim = 0.5+ c(0, nr),
axes = FALSE, xlab = '', ylab = '', col=col, breaks=breaks,
...)
然后绘制轴标签它使用mtext,它对par参数不敏感,但接受cex参数,这里没有提供:
# from gplots heatmap.2 source:
if(!is.null(xlab)) mtext(xlab, side = 1, line = margins[1] - 1.25)
if(!is.null(ylab)) mtext(ylab, side = 4, line = margins[2] - 1.25)
这意味着这些标签的大小始终相同,而热图有时会很大,所以这可能是个问题。也可以缩小图形设备,以及所有其他具有自定义选项的文本元素,然后轴标签的相对大小会更大。但是我们怎样才能以更灵活的方式解决这个问题呢?我修改了这样的代码:
# gplots heatmap.2 code modified:
## add row and column headings (xlab, ylab)
cex.lab <- ifelse(invalid(match.call()$cex.lab), 1.0, match.call()$cex.lab)
if(!is.null(xlab)) mtext(xlab, side = 1,
line = margins[1] - 1.25 + cex.lab / 5,
cex = cex.lab)
if(!is.null(ylab)) mtext(ylab, side = 4,
line = margins[2] - 1.25 + cex.lab / 5,
cex = cex.lab)
唯一的问题是定位参数现在是cex.lab / 5,这样可以避免在合理范围内与ticklabels重叠,但是可能需要在函数中添加另一个参数,如果不适合用户可以修改。
您可以将heatmap.2.R 从gplots source 复制到您的工作目录,修改并获取它:
source('heatmap.2.R')
但是,您还需要将这些行添加到开头:
require(gtools)
plot.dendrogram = stats:::plot.dendrogram