カタバミさんのプログラミングノート

日曜プログラマーがプログラミング関係のメモを記録するブログです。

Rで二項分布とパーセント点の垂線を描画する

 X\sim \mathrm B(1000, 0.5), 0 \leq X \leq 1000 として、二項分布の曲線と95%点の垂線を描画します。

x <- dbinom(0:1000, 1000, 0.5)
plot(x, type="l", xlim=c(400, 600), yaxs="i", xlab="X", ylab="Prob.")
x1 <- qbinom(1-0.05, 1000, 0.5) #1-0.05 = 0.95 = 95%
segments(x1, 0, x1, x[x1])

f:id:by-unknown:20191015204238p:plain