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
の中の、、、
R言語のpbirthday, qbirthdayという関数に使われていたりする。
これ:http://stat.ethz.ch/R-manual/R-patched/library/stats/html/birthday.html
では図示すると、、、
50%を超えるのは
2人一致:23人、 3人一致:88人、4人一致:187人、5人一致:313人、6人一致:459人とかになる。
本当にこの近似式いいの?という疑問があるので、メルセンヌツイスタでN=1000000でモンテカルロシミュレーションをしたのがこれ。
なかなかの精度。
--
でこの式をお手軽に計算したいときは、カシオの高精度計算サイトkeisan.casio.jpにアップしている式をどうぞ。
« 冷凍マンモス YUKAちゃんを見にみなとみらい パシフィコ横浜まで。 今はバルキリー(実物大)もみなとみらい駅に。 | トップページ | 沼津港深海水族館(シーラカンスミュージアム)へ行ってきた。 »
「学問・資格」カテゴリの記事
- 高周波(RF・マイクロ波・ミリ波・5G)関連ニュース2021年2月16日 IEEE Microwave Magazineの特集はオールデジタルのRFID、Microwave JournalはEバンド ミリ波通信に衛星や気球を使う話、アメリカの半導体企業がバイデンに投資を迫る、(2021.02.17)
- カオスを生じる電気回路、Chua’s circuitをLTspiceで回路シミュレーションしてみる。(2021.02.19)
- Labyrinth Chaos(迷宮カオス)を生むThomas-Rössler方程式のパラメータbを色々変えて、Python+Scipyでルンゲクッタ8次のDOP853(Dormand&Prince)を使って計算してGIFアニメ(2021.02.16)
- フィッツヒュー・南雲 (FitzHugh-Nagumo) 方程式をPython+Scipyでルンゲクッタ8次のDOP853(Dormand Prince)で計算。(2021.02.23)
- 「水晶振動子の等価回路計算」をカシオの高精度計算サイトkeisan.casio.jpの自作式としてUP! インピーダンスの大きさと位相がグラフ化できる。(2021.02.12)
« 冷凍マンモス YUKAちゃんを見にみなとみらい パシフィコ横浜まで。 今はバルキリー(実物大)もみなとみらい駅に。 | トップページ | 沼津港深海水族館(シーラカンスミュージアム)へ行ってきた。 »
コメント