在GGPlot2直方图中为X值以上的任何值创建一个bin

时间:2012-07-23 17:25:41

标签: r ggplot2 histogram

使用ggplot2,我想创建一个直方图,其中X以上的任何内容都被分组到最终的bin中。例如,如果我的大多数发行版都在100到200之间,并且我希望以10分为单位,那么我希望将200以上的任何内容分类到" 200 +"中。

# create some fake data    
id <- sample(1:100000, 10000, rep=T)
visits <- sample(1:1200,10000, rep=T)

#merge to create a dataframe
df <- data.frame(cbind(id,visits))

#plot the data
hist <- ggplot(df, aes(x=visits)) + geom_histogram(binwidth=50)

如何限制X轴,同时仍然表示我想要限制的数据?

2 个答案:

答案 0 :(得分:5)

也许您正在寻找breaks的{​​{1}}参数:

geom_histogram

这看起来像这样(警告说这里假数据看起来很糟糕,轴也需要调整以匹配断点):

manual breaks on histogram

修改

也许其他人可以在这里权衡:

# create some fake data    
id <- sample(1:100000, 10000, rep=T)
visits <- sample(1:1200,10000, rep=T)

#merge to create a dataframe
df <- data.frame(cbind(id,visits))

#plot the data
require(ggplot2)
ggplot(df, aes(x=visits)) +
  geom_histogram(breaks=c(seq(0, 200, by=10), max(visits)), position = "identity") +
  coord_cartesian(xlim=c(0,210))

情节错误:

# create breaks and labels
brks <- c(seq(0, 200, by=10), max(visits))
lbls <- c(as.character(seq(0, 190, by=10)), "200+", "")
# true
length(brks)==length(lbls)

# hmmm
ggplot(df, aes(x=visits)) +
  geom_histogram(breaks=brks, position = "identity") +
  coord_cartesian(xlim=c(0,220)) +
  scale_x_continuous(labels=lbls)

看起来像this,但是在8个月前修复了。

答案 1 :(得分:3)

如果你想稍微捏一下bin标签的问题,那么只需将数据子集化并在新的牺牲数据框中创建分箱值:

id <- sample(1:100000, 10000, rep=T)
visits <- sample(1:1200,10000, rep=T)

#merge to create a dataframe
df <- data.frame(cbind(id,visits))
#create sacrificical data frame
dfsac <- df
dfsac$visits[dfsac$visits > 200 ] <- 200

然后使用breaks中的scale_x_continuous命令轻松定义bin标签:

ggplot(data=dfsac, aes(dfsac$visits)) + 
  geom_histogram(breaks=c(seq(0, 200, by=10)), 
                 col="black", 
                 fill="red") +
  labs(x="Visits", y="Count")+
  scale_x_continuous(limits=c(0, 200), breaks=c(seq(0, 200, by=10)), labels=c(seq(0,190, by=10), "200+"))

enter image description here