假设我有多次双变量实验的复制,我希望同时显示六角形分级图,具有常见的细胞计数。是否存在执行此操作的代码?有没有一种简单的方法可以修改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图的网格,但我希望如此。
答案 0 :(得分:2)
使用网格图形绘制hexbin对象,因此par(mfrow = c(3,3))不会执行任何操作。每个图表都绘制在一个单独的页面上。要获得绘图选项的详细信息:
?gplot.hexbin
在这种情况下,我们要将maxcnt设置为最大的单元格数:
lapply(h, plot, maxcnt=max(unlist(lapply(h, function(x) max(x@count)))))
这会将相同的图例应用于每个图表。