r - 統計 : 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

一般に、視覚化に関する質問は、次のような場合に発生します。

成功の確率が0.8であるとすると、20回の操作のうち「少なくとも」16回成功する確率はどれくらいですか。

これは、次の2項式を使用して実行できます。

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

k = 16..2020のうち少なくとも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)

これにより、次のようになります。

お役に立てれば。