我的目标是使用 ggplot2 生成二项分布的两个重叠 PMF,根据我指定的颜色进行颜色编码,底部有一个图例。
到目前为止,我认为我已经正确设置了数据框。
successes <- c(seq(0,10,1),seq(0,10,1))
freq <- c(dbinom(seq(0,10,1),10,0.2),dbinom(seq(0,10,1),10,0.8))
class <- c(rep(' A ',11),rep(' B ',11))
df1 <- data.frame(cbind(successes,freq,class))
然而,这给出了错误的结果。
library(ggplot2)
g <- ggplot(df1, aes(successes),y=freq)
g + geom_bar(aes(fill = class))
我觉得我在效仿一个例子,却得到了完全不同的结果。这(几乎)做了我想要的:如果它给出相对频率,那将是准确的。
g <- ggplot(mpg, aes(class))
g + geom_bar(aes(fill = drv))
几个问题:
1)我的代码块哪里出错了?
2) 有没有更好的方法在一张图中向 PMF 显示?我没有决心使用直方图或条形图。
3) 我该如何设置才能让我能够选择颜色?
4) 如何对 x 轴上的值进行排序?它们不是类别。它们是数字 0-10,并且具有我想保留的自然顺序。
谢谢!
更新
以下两个块起作用。
successes <- c(seq(0,10,1),seq(0,10,1))
freq <- c(dbinom(seq(0,10,1),10,0.2),dbinom(seq(0,10,1),10,0.8))
class <- c(rep(' A ',11),rep(' B ',11))
df1 <- data.frame(successes,freq,class)
ggplot(df1, aes(successes ,y=freq, fill = class)) +
geom_bar(stat = "identity") +
scale_x_continuous(breaks = seq(0,10,1)) +
scale_fill_manual(values = c("blue", "green")) + theme_bw()
和
successes <- c(seq(0,10,1),seq(0,10,1))
freq <- c(dbinom(seq(0,10,1),10,0.2),dbinom(seq(0,10,1),10,0.8))
class <- c(rep(' A ',11),rep(' B ',11))
df1 <- data.frame(successes,freq,class)
ggplot(df1, aes(x=successes,y=freq),y=freq) +
geom_col(aes(fill = class)) +
scale_x_continuous(breaks = seq(0,10,1)) +
scale_fill_manual(values = c("blue", "green")) + theme_bw()