« 冷凍マンモス YUKAちゃんを見にみなとみらい パシフィコ横浜まで。 今はバルキリー(実物大)もみなとみらい駅に。 | トップページ | 沼津港深海水族館(シーラカンスミュージアム)へ行ってきた。 »

2013年7月22日 (月)

N人中k人の誕生日が一致する確率は?(R言語のpbirthday, qbirthday)【再掲】

Twitterで誕生日のパラドックスの話が流れていたので、昔書いたものを再掲してみる。

ある集まりの中で、2人の誕生日が一致する確率が案外高いということがある。ちょっと直感に反するところなので誕生日のパラドックスと呼ばれたりする。

では、これを拡張して、、、

N人中、k人の誕生日が一致する確率は?というのは結構難しい(というかちゃんと厳密な式が書いてある論文がリンク切れになっていた、、、)。

しかし、非常によく合う近似式が知られている。

”Methods for Studying Coincidences” 、PERSl DlACONlS and FREDERICK MOSTELLER,

Journal of the American Statistical Association
December 1989, Vol. 84, No. 408, Applications & Case Studies

の中の、、、

Birthday_problem2013


という式。

R言語のpbirthday, qbirthdayという関数に使われていたりする。

これ:http://stat.ethz.ch/R-manual/R-patched/library/stats/html/birthday.html

では図示すると、、、

50%を超えるのは

Birthday_problem20132

2人一致:23人、 3人一致:88人、4人一致:187人、5人一致:313人、6人一致:459人とかになる。

本当にこの近似式いいの?という疑問があるので、メルセンヌツイスタでN=1000000でモンテカルロシミュレーションをしたのがこれ。

Birthday_problem20133

なかなかの精度。
--

でこの式をお手軽に計算したいときは、カシオの高精度計算サイトkeisan.casio.jpにアップしている式をどうぞ。

誕生日の一致する確率(N人中k人以上)

誕生日が一致する確率(N人中2~5人、グラフ表示)

« 冷凍マンモス YUKAちゃんを見にみなとみらい パシフィコ横浜まで。 今はバルキリー(実物大)もみなとみらい駅に。 | トップページ | 沼津港深海水族館(シーラカンスミュージアム)へ行ってきた。 »

学問・資格」カテゴリの記事

コメント

コメントを書く

(ウェブ上には掲載しません)

トラックバック

この記事のトラックバックURL:
http://app.cocolog-nifty.com/t/trackback/512682/57844620

この記事へのトラックバック一覧です: N人中k人の誕生日が一致する確率は?(R言語のpbirthday, qbirthday)【再掲】:

« 冷凍マンモス YUKAちゃんを見にみなとみらい パシフィコ横浜まで。 今はバルキリー(実物大)もみなとみらい駅に。 | トップページ | 沼津港深海水族館(シーラカンスミュージアム)へ行ってきた。 »

最近の記事

最近のコメント

2017年11月
      1 2 3 4
5 6 7 8 9 10 11
12 13 14 15 16 17 18
19 20 21 22 23 24 25
26 27 28 29 30    
フォト
無料ブログはココログ