带有常见箱/图例中断的hexbin图的矩阵

时间:2013-01-06 21:00:36

标签: r plot

假设我有多次双变量实验的复制,我希望同时显示六角形分级图,具有常见的细胞计数。是否存在执行此操作的代码?有没有一种简单的方法可以修改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)

此代码不显示具有常见细胞计数的hexbin图的网格,但我希望如此。

1 个答案:

答案 0 :(得分:2)

使用网格图形绘制hexbin对象,因此par(mfrow = c(3,3))不会执行任何操作。每个图表都绘制在一个单独的页面上。要获得绘图选项的详细信息:

?gplot.hexbin

在这种情况下,我们要将maxcnt设置为最大的单元格数:

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

这会将相同的图例应用于每个图表。