r - 具有常见 bin/legend 中断的 hexbin 图矩阵

Question

假设我有许多双变量实验的重复，我希望在六边形分箱图中同时显示这些实验，并具有常见的细胞计数。是否有现有的代码可以做到这一点？有没有一种简单的方法来修改 hexbin 包来为我做到这一点？

例如：

library(hexbin)
x <- replicate(9, rnorm(10000), simplify=FALSE)
y <- replicate(9, rnorm(10000), simplify=FALSE)
h <- mapply(hexbin, x, y)
par(mfrow=c(3,3))
lapply(h, plot)

此代码不显示具有常见单元格计数的十六进制图网格，但我喜欢它。

Answer 1

hexbin 对象是使用网格图形绘制的，因此您的 par(mfrow=c(3,3)) 不会做任何事情。每个图表都绘制在单独的页面上。要获取绘图选项的详细信息：

?gplot.hexbin

在这种情况下，我们希望将 maxcnt 设置为最大的单元格计数：

lapply(h, plot, maxcnt=max(unlist(lapply(h, function(x) max(x@count)))))

这将对每个图形应用相同的图例。