母比率 0.5,試行回数 10 の二項分布で,x が 9 以上になる確率を求めるのに,
> pbinom(9, 10, 0.5, lower.tail=FALSE)
[1] 0.0009765625
としてはだめだ。
> pbinom(8, 10, 0.5, lower.tail=FALSE)
[1] 0.01074219
でなくてはならない。lower.tail=TRUE と,以下を考え合わせれば分かるが,うっかりすると最初に挙げたようにしてしまいそう。
> cbind(0:10, x, cumsum(x))
x
[1,] 0 0.0009765625 0.0009765625
[2,] 1 0.0097656250 0.0107421875
[3,] 2 0.0439453125 0.0546875000
[4,] 3 0.1171875000 0.1718750000
[5,] 4 0.2050781250 0.3769531250
[6,] 5 0.2460937500 0.6230468750
[7,] 6 0.2050781250 0.8281250000
[8,] 7 0.1171875000 0.9453125000
[9,] 8 0.0439453125 0.9892578125
[10,] 9 0.0097656250 0.9990234375
[11,] 10 0.0009765625 1.0000000000
> pbinom(9, 10, 0.5, lower.tail=FALSE)
[1] 0.0009765625
としてはだめだ。
> pbinom(8, 10, 0.5, lower.tail=FALSE)
[1] 0.01074219
でなくてはならない。lower.tail=TRUE と,以下を考え合わせれば分かるが,うっかりすると最初に挙げたようにしてしまいそう。
> cbind(0:10, x, cumsum(x))
x
[1,] 0 0.0009765625 0.0009765625
[2,] 1 0.0097656250 0.0107421875
[3,] 2 0.0439453125 0.0546875000
[4,] 3 0.1171875000 0.1718750000
[5,] 4 0.2050781250 0.3769531250
[6,] 5 0.2460937500 0.6230468750
[7,] 6 0.2050781250 0.8281250000
[8,] 7 0.1171875000 0.9453125000
[9,] 8 0.0439453125 0.9892578125
[10,] 9 0.0097656250 0.9990234375
[11,] 10 0.0009765625 1.0000000000
※コメント投稿者のブログIDはブログ作成者のみに通知されます