裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

二項分布の上側確率

2010年07月14日 | ブログラミング
母比率 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
コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« R 2,7 の新機能など | トップ | sd(list())は? »
最新の画像もっと見る

コメントを投稿

ブログラミング」カテゴリの最新記事