尝试了解为binwidth
中的每个要素级别分配唯一geom_histogram
。
到目前为止失败了。
这是可重复的数据
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)
#获取每个组的每个binwidth
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"))
发生了同样的错误。我无法理解为什么会出现同样的错误?
答案 0 :(得分:3)
你可以迭代垃圾箱来创建图层
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个方面......这很多。
Nota:在dplyr 0.4.3
上使用ggplot2 2.1.0
,plyr 1.8.3
和R 3.2.3