r - 一个方面的每个图的不同 Binwidth

Question

尝试了解为 geom_histogram 中的每个因子水平分配唯一的 binwidth。到目前为止虽然失败了。

这是可重现的数据

a <- rnorm(10,7,0.1)
b <- rnorm(10,13,5)
df <- data.frame(data = c(a,b),group=c(rep("a",10),rep("b",10)))
kk <- df%>%
  group_by(group)%>%
  mutate(bin=density(data)$bw)

binns <- round(unique(kk$bin),digits = 2)  # to get each binwidth for each group

ggplot()+
  geom_histogram(data=kk,aes(x=data, fill=group),binwidth=binss)+
  facet_wrap(~group,scales=c("free_y"))

Error in seq.default(round_any(range[1], size, floor),    round_any(range[2],  : 
  'from' must be of length 1
Error in seq.default(round_any(range[1], size, floor), round_any(range[2],  : 
  'from' must be of length 1
Error in exists(name, envir = env, mode = mode) : 
  argument "env" is missing, with no default

然后我试了一下

ggplot()+
  geom_histogram(data=kk,aes(x=data, fill=group),binwidth=c(binns[1],binns[2]))+
  facet_wrap(~group,scales=c("free_y"))

同样的错误发生了。我不明白为什么会出现同样的错误。

Answer 1

您可以迭代 bin 以创建图层

library(dplyr)
a <- rnorm(10,7,0.1)
b <- rnorm(10,13,5)
df <- data.frame(data = c(a,b),group=c(rep("a",10),rep("b",10)))
kk <- df %>%
  group_by(group) %>%
  mutate(bin=round(density(data)$bw, 2))
binns <- unique(kk$bin)

附上ggplot2

library(ggplot2)

在列表中为每个 bin 值创建一个直方图图层，其中仅包含此 bin 的数据。如果你有 30 个级别，你将有一个包含 30 个直方图图层的列表

lp_hist <- plyr::llply(binns, function(b) {
  geom_histogram(data = kk %>% 
                   filter(bin == b), 
                 aes(x = data, fill=group), 
                 binwidth = b)
  })

合并这些图层，将它们全部添加到 ggplot() 对象

p_hist <- Reduce("+", lp_hist, init = ggplot())

随心所欲的分面和缩放

p_hist + facet_grid(. ~ group, scales = "free_y")

您将获得所需的图形，即按面表示不同的 binwidth。

当心，因为 30 个级别将给出 30 个方面......很多。

注意:在 R 3.2 上使用 dplyr 0.4.3、ggplot2 2.1.0 和 plyr 1.8.3 .3