r - 统计 : binomial distribution and simple graph in R

Question

我正在尝试为 R 中的二项式分布制作一个简单的图表。

所以问题是“有 20 名患者，对 4 名患者成功手术的概率是多少，(假设概率 = 0.8)。”

我做的是

x <- rbinom(4, size= 20, prob = 0.8)
y <- pbinom(x, size = 20, prob = 0.8)
plot(y, type="l")

但我不确定这是否是正确的绘图方式..

Answer 1

一般来说，当您的问题类似于以下内容时，您的可视化问题就会出现:

What's the probability of having "at least" 16 successes out of 20 operations given that the probability of success is 0.8?

这可以通过使用二项式公式来完成:

p(x=k) = choose(n, k) * .8^k * .2^(n-k) # equivalent to dbinom(k, n, prob=0.8)

对于 k = 16..20，我们需要相同的值并将所有这些值相加，以获得 20 次中至少有 16 次成功的概率。这是使用 dbinom< 完成的作为:

sum(dbinom(16:20, 20, prob=0.8)) # 0.6296483

请注意，在此成功率 (0.8) 下，至少 4 次成功的概率仅为 1。也就是说，我们肯定至少会成功 4 次。这就是为什么我选择相对较高的成功率。

绘制此图(使用 ggplot2):

df <- data.frame(x=1:20, prob=dbinom(1:20, 20, prob=0.8))
require(ggplot2)
ggplot(data=dd, aes(x=x,y=prob)) + geom_line() +
      geom_ribbon(data=subset(dd,x>=16 & x<=20),aes(ymax=prob),ymin=0,
                  fill="red", colour = NA, alpha = 0.5)

这给你类似的东西:

希望这对您有所帮助。