「統計学とRのはなし」にて
スライド38で,普通に起こること95%を中央に設定したのが間違いの元。左側のみが5%になる-λを設定しないといけない。
プログラムで言えば,qnorm(0.025) ではなくて qnorm(0.05) とすべし。
また,r <- sample(1:100, 1) は r <- sample(0:100, 1) とすべし。
根本的に,こういう類の問題はシミュレーションで解くものではない(プログラムもへたくそでまずいけど)。
z=qnorm(0.05) として (0.1-p)/sqrt{p(1-p)/100}=z を p について解けばよい。p=0.1603555 となる。
また,正規近似するのではなく,二項分布で正確に求めようとするならば,uniroot(function(p) pbinom(9, 100, prob=p)-0.05, c(0, 1)) で求める。この場合の答えは p=0.1518244
プログラムは,改造しようという気にもならない。
コメントを付けようと思ったけど,Facebook にログインしろだの SlideShare のアカウントは持っているかだの,小うるさいこというので,そこまでして教えてやる必要もないなあと思った次第。特殊なポリシーを持っていると,損するよね。
最新の画像[もっと見る]
- 算額(その1146) 6時間前
- 算額(その1144) 10時間前
- 算額(その1143) 1日前
- 算額(その1142) 2日前
- 算額(その1141) 2日前
- 算額(その1140) 2日前
- 算額(その1139) 3日前
- 算額(その1138) 3日前
- 算額(その1137) 4日前
- 算額(その1136) 4日前
※コメント投稿者のブログIDはブログ作成者のみに通知されます