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ちゃんを見にみなとみらい パシフィコ横浜まで。 今はバルキリー（実物大）もみなとみらい駅に。 | トップページ | 沼津港深海水族館（シーラカンスミュージアム）へ行ってきた。 »
「学問・資格」カテゴリの記事
- アメリカの空港で5G C-bandをを使うなという話で揉めてますが、どう考えても電波高度計(4.2-4.4GHz）が無茶苦茶で、衛星通信(3.7-4.2GHz)を邪魔しているのにC-Band(3.7–3.98GHz)が入ってくるとクレームとは、、、ジャイアニズムか！(2022.01.22)
- 不定方程式a*x+b*y=gcd(a,b) （ディオファントス方程式、ベズーの等式）の拡張ユークリッド互除法を使った計算をカシオの高精度計算サイトkeisan.casio.jpにUP!2022年共通テストの数学I・数学Aの問題4に出てても結果に自信ないので機械に頼る。。。(2022.01.19)
- 2022年度共通テストの数学I・数学Aの問題、例の防衛省のイージス・アシュアの縮尺違い地図で角度間違ってる話をそのまま持ってきてる！秋田県男鹿市の本山か。。。(2022.01.17)
- 高周波（RF・マイクロ波・ミリ波・5G）関連ニュース2022年1月18日 IEEE Microwave Magazineは女性高周波研究者特集、Microwave Journalはレーダで月面が見える。ドコモの誘電体線路アンテナ、3.45GHzのオークション結果、hBNを使ったスイッチ、Scikit-rfの使い方など。(2022.01.18)
- 犬や猫などのペットにマイクロチップを埋め込むことが義務化されるが、どんな無線通信技術を使っているか調べた。ISO11784/11785に準拠した134.2kHzを使うRFIDで、フェライトロッドアンテナ、キャパシタ、チップ内蔵。15桁の個体識別番号が読み取れる。(2022.01.10)
« 冷凍マンモス YUKAちゃんを見にみなとみらい パシフィコ横浜まで。 今はバルキリー（実物大）もみなとみらい駅に。 | トップページ | 沼津港深海水族館（シーラカンスミュージアム）へ行ってきた。 »
コメント