r - ggplot2 因子直方图显示概率质量而不是计数

Question

我正在尝试使用优秀的 ggplot2，使用 bar geom 来绘制概率质量而不是计数。但是，使用 aes(y=..密度..) 分布的总和不等于 1(但很接近)。我认为问题可能是由于因子的默认 binwidth 造成的。这是问题的示例，

age <- c(rep(0,4), rep(1,4))
mppf <- c(1,1,1,0,1,1,0,0)
data.test <- as.data.frame(cbind(age,mppf))
data.test$age <- as.factor(data.test$age)
data.test$mppf <- as.factor(data.test$mppf)
p.test.density <- ggplot(data.test, aes(mppf, group=age, fill=age)) +
geom_bar(aes(y=..density..), position='dodge') +
scale_y_continuous(limits=c(0,1))
dev.new()
print(p.test.density)

我可以通过保持 x 变量连续并设置 binwidth=1 来解决这个问题，但这看起来不太优雅。

data.test$mppf.numeric <- as.numeric(data.test$mppf)
p.test.density.numeric <- ggplot(data.test, aes(mppf.numeric, group=age, fill=age)) + 
geom_histogram(aes(y=..density..), position='dodge', binwidth=1)+ 
scale_y_continuous(limits=c(0,1))
dev.new()
print(p.test.density.numeric)

Answer 1

我认为您几乎已经弄清楚了，并且一旦您意识到您需要一个条形图而不是直方图，您就会明白了。

包含分类数据的条形图的默认宽度为 0.9(请参阅 ?stat_bin。geom_bar 的帮助页面不会提供默认的条形图宽度，但会向您发送到stat_bin以供进一步阅读。)。鉴于此，您的图显示了条形宽度为 0.9 的正确密度。只需将宽度更改为 1，您就会看到您期望看到的密度值。

ggplot(data.test, aes(x = mppf, group = age, fill = age)) +
  geom_bar(aes(y=..density..), position = "dodge", width = 1) +
  scale_y_continuous(limits=c(0,1))