pnbinom のバックアップソース(No.1) - アールメカブ

アールメカブ


pnbinom のバックアップソース(No.1)

[[Rの備忘録]]

3 回表が出るまでに 9 回以上裏が出る確率

 > 1 - pnbinom(8, size=3, prob = .5)
 [1] 0.03271484

 > sum(dnbinom(9:100, size=3, prob = .5))
 [1] 0.03271484

繁桝『ベイズ統計入門』p.42 では以下の計算結果となっている。

 > sum(pnbinom(8:10, size=3, prob = .5, lower = F))
 [1] 0.06323242

 > pnbinom(9, size=3, prob = .5)
 [1] 0.9807129
 > pnbinom(10, size=3, prob = .5)
 [1] 0.9887695

 > dnbinom(10, size=3, prob = .5)
 [1] 0.00805664


 > pnbinom(9, size=3, prob = .5, low =T ) + 
         dnbinom(10, si = 3, pro = .5)
 [1] 0.9887695