R 中各种成功概率的二项分布图

Question

是否可以在 R 中为具有不同成功概率的二项分布随机变量绘制类似直方图的条/线？

试验次数 (n) 和样本空间保持不变。只有成功的概率 (p) 不同。为了让条形图彼此相邻，R 代码的外观如何？

这是我的意思的一个小例子：

yval <- 0:10 # sample space
plot(dbinom(yval, 10, 0.5) ~ yval, type = "h", col = "black", ylim = c(0, 0.35))
lines(dbinom(yval, 10, 1/6) ~ yval, type = "h", col = "red")
lines(dbinom(yval, 10, 0.6) ~ yval, type = "h", col = "green")
legend("topright", legend = c("p = 0.5", "p = 1/6", "p = 0.6"),
col = c("black", "red", "green"), lty = 1, cex = 0.7)

使用此代码，线条可以相互绘制。我需要在这里改变什么？

提前谢谢你。

Answer 1

你的意思是如下吗？

您可以在其他直方图的 x 值上添加一个小偏移量：

yval <- 0:10 # sample space
plot(dbinom(yval, 10, 0.5) ~ yval, type = "h", col = "black", ylim = c(0, 0.35))
lines(dbinom(yval, 10, 1/6) ~ I(yval + 0.1), type = "h", col = "red") # + 0.1
lines(dbinom(yval, 10, 0.6) ~ I(yval + 0.2), type = "h", col = "green") # + 0.2
legend("topright", legend = c("p = 0.5", "p = 1/6", "p = 0.6"),
col = c("black", "red", "green"), lty = 1, cex = 0.7)

Answer 2

为了结束我的评论，下面是我们可以用barplot 做的事情：

prob <- c(0.5, 1/6, 0.6)
yval <- 0:10
Y <- t(outer(yval, prob, dbinom, size = max(yval)))
barplot(Y, names.arg = yval, beside = TRUE, col = 1:3, border = 1:3,
        legend.text = paste0("p = ", format(prob, digits = 2)))

把它当作替代品。

---

备注1

注意上面outer的使用。为什么我们要多做一个t()？我们不能这样做吗？

Y <- outer(prob, yval, dbinom, size = max(yval))

没有。这会给你很多NaN 警告。检查args(dbinom)。功能dbinom预计yval在prob。

---

备注2

使用barplot可以轻松为您制作一些“副作用”，如下（使您更容易比较它们的形状）。

prob <- c(0.5, 1/6, 0.6)
yval <- 0:10
Y <- outer(yval, prob, dbinom, size = max(yval))  ## no `t()` now
barplot(Y, names.arg = paste0("p = ", format(prob, digits = 2)), beside = TRUE)

使用col 587654338 @未设置，使用灰色颜色，黑暗与杆的高度成比例。