Loading [MathJax]/jax/output/CommonHTML/jax.js

2018年12月2日日曜日

平均一定の二項分布たち

War-Games Advent Calendar 2018 https://adventar.org/calendars/3461 の記事です。すみません今年もゲームなしです。

平均一定の二項分布たち、試行回数を増やすと少数の法則が成立して、ポアソン分布になります。滅多に起きない事象についての法則で、製品の不良率の評価などに使われるそうです。

というのは、二項分布
Prob(X=k)=nCkpk(1p)nk
の平均がnpだから、平均N=npを一定に保って試行回数nを大きくすると確率pp=N/nに決まって、この確率自体は小さくなります。じゃあk回当たる確率
Probn(X=k)=n!(nk)!k!(N/n)k(1N/n)nk
nが増えるとどうなっていくのか、ということ。

これはつまりN=2n=3,4,6,12の場合をはねはねさんがMAの日記で提案したものです。
 2018年11月20日12:41
1対1の差しでの勝負、打撃の多いほうが勝ちというルールで先手だとして。みなさん、どれを選択しますか?
1.サイコロ振らずに2打撃
2.サイコロ3個振って、1-4で1打撃
3.サイコロ4個振って、1-3で1打撃
4.サイコロ6個振って、1-2で1打撃
5.サイコロ12個振って、1で1打撃

分散がnp(1p)=2(12/n)なので、 nが増えると増加しますから、どんどんブラッディになります、という指摘はコメントに見られます。

結論は、
limnProbn(X=k)=e22kk!
となってこれはポアソン分布に他ならないのですが、これだけだとあれなので

2018年11月20日17:24 さかみちさんのコメント 
計算が違っていなければ、選択肢5に向かって期待値である2以上の確率は下がり、逆に3以上の確率は上がるので。
N=2で確かめます。

負ける確率は、打撃0か1です。nが増えるにつれて1以下の確率が上がればよく、

まず
Probn(X=0)+Probn(X=1)=pn+np(1p)n1=(12/n)n+2(12/n)n1
です。

自然対数の底を定義するときに現れる数列(1+1/n)は単調増加なので、同じようにして(12/n)nも単調増加。

というわけで nが増えると2以上の確率は下がります。収束先は3e2=3/e2です。

2になる確率を計算すると,
Probn(X=2)=n(n1)2(2/n)2(12/n)n2=2(11/n)(12/n)n2=211/n(12/n)2(12/n)n
これはnが増えると減少し、2/e2に収束します。結局、打撃が3以上になる確率はnが増えると増えます。増えますけど、打撃が0か1になる確率もそれなりにあるので、選択は難しい。大きな打撃を与える確率は非常に小さくなりますから、少数の法則そのものです。

0 件のコメント:

コメントを投稿