Try   HackMD

ベルヌーイ試行の繰り返し

tags: probability-theory

ベルヌーイ試行を繰り返す

ベルヌーイ試行を繰り返すことを考える。1回目のベルヌーイ試行を

X1、2回目のベルヌーイ試行を
X2
と表すと、繰り返しに応じて添え字が増えていく。ベルヌーイ試行を、試行の順序通りに並べると
X1,X2,,Xi,

となる。
i
番目のベルヌーイ試行を
Xi
と記す。数列の一般項のようなものだが、添え字がある以外はすべて、同じ試行が並んでいる。

これらの試行は、互いに独立に同一の試行が繰り返されているものとする。コイン投げを繰り返す時には、投げ終えたコインを拾い、同じポジションに戻ってからまた投げる。受験を繰り返すのであれば、同じ学力のまま、テストだけが新しくなる想定でまた受験する。天気であれば、同じ気圧配置、同じ雲の配置が朝に繰り返される想定で、また予想する。

ベルヌーイ試行の繰り返しを止める

試行の繰り返しはいつか、止めなければならない。確率論の序盤で興味のある止め方は、次の3通りである。

  1. あらかじめ定めた
    n
    回の試行を終えたら、繰り返しを止める。
  2. 1
    回成功したら、試行の繰り返しを止める。
  3. r
    回成功したら、試行の繰り返しを止める。

総試行回数
n
を固定する場合(二項試行)

一つ目の止め方では、試行回数が予め決まっているので、成功回数

X1+X2++Xn
を数える。
Y=X1+X2++Xn

のモーメントは、それぞれのモーメントの和になる。
E[Y]=E[X1]+E[X2]++E[Xn]=np

分散も共分散が
0
なので、分散の和になる。
V[Y]=V[X1+X2++Xn]=V[X1]+V[X2]++V[Xn]+2i=1n1j=i+1nCov[Xi,Xj]=V[X1]+V[X2]++V[Xn]=np(1p)

モーメント母関数はそれぞれのモーメント母関数の積になる。
MY(t)=MX1(t)MX2(t)MXn(t)=(1p+pet)n

総試行回数

nを固定してベルヌーイ試行を互いに独立に繰り返したうちの、成功回数
Y
の確率関数は次のように導かれる。成功回数が
y
のすべての場合の確率の和を求める。

x1+x2++xn=yi=1npxi(1p)1xi=x1+x2++xn=ypixi(1p)nixi=x1+x2++xn=ypy(1p)ny=nCypy(1p)ny=n!y!(ny)!py(1p)ny

この確率関数を持つ確率分布を、二項分布という。二項分布に関する計算は、別の機会に譲る。

成功回数を
1
に固定する場合(幾何試行)

二つ目の止め方では、成功回数が決まっているので、試行回数

Nを数える。
N1=arginfn{X1+X2++Xn=1}

確率変数

N1は、
1
回成功するまでの総試行回数とする。この確率変数については、確率関数を先に導く。総試行回数が
k
回のときに1回目の成功があるという事象
{N1=k}

は、
k1
回連続して失敗し、
k
回目に初めて成功したことを意味する。各試行で表すと
{X1==Xk1=0,Xk=1}

となる。試行が互いに独立ならば、この事象の確率は
p(k)=Pr[N1=k]=(1p)k1p

である。この確率関数を持つ確率分布を、幾何分布という。確率が
p,p(1p),p(1p)2,

のように初項が
p
、公比が
1p
の幾何数列となっていることに由来する。幾何数列は等比数列ともいう。

幾何分布に関する計算も、別の機会に譲る。

成功回数を
r
に固定する場合(負の二項試行)

三つ目の止め方でも、成功回数が決まっているので、試行回数

Nを数える。
Nr=arginfn{X1+X2++Xn=r}

確率変数

Nrは、
r
回成功するまでの総試行回数とする。この確率変数については、確率関数を先に導く。総試行回数が
k
回のときに
1
回目の成功があるという事象
{Nr=k}

は、
k1
回目までに
r1
回成功していて、
k
回目に初めて
k
回目の成功があったことを意味する。

i回目の試行までの成功回数を
Yi
と置くと、これはそこまでの試行の総和
Yi=X1+X2++Xi

に等しい。
{Nr=k}
という事象を
Yi
を用いて表すと、
{Y1<r,Y2<r,,Yk1<r,Yk=r}

となる。
r=1
を代入してみると、
{Y1<1,Y2<1,,Yk1<1,Yk=1}

のように、
k
回目の試行で初めて
Y
1
となっていて、ここで初めて成功したことを表現できている。

ところで、

r回成功するまでの総試行は、
1
回成功するまでの試行の繰り返しを、さらに
r
回繰り返すことと等しい。
ここで、ベルヌーイ試行を
1
回成功するまで繰り返して成功したら止める試行を定義する。

N1=arginfn{X1+X2++Xn=1}

この試行を幾何試行と呼ぶ。幾何試行の

r回の繰り返しを
N1,1,N1,2,,N1,r

と置くと、これらは互いに独立に幾何分布に従う確率変数となる。そして、
Nr=N1,1+N1,2++N1,r

となる。これより、
Nr
が従う確率分布は同一の幾何分布に従う確率変数の和の分布に等しい。

確率変数

Nrの確率分布は、この事実を用いて導かれる。その確率分布を負の二項分布という。

負の二項分布に関する確率計算の詳細も、別の機会に譲る。

数える回数

上の説明では、二つ目の止め方と三つ目の止め方で、総試行回数を数えた。幾何分布の確率関数

p(k)=Pr[N1=k]=(1p)k1p
が、
k1
回連続して失敗してから、
k
回目に成功する確率であることから、失敗した回数
M
を数えても良い。

幾何試行について、両者は

M1=N11
の関係にある。また負の二項試行についても、両者は
Mr=Nrr

の関係にある。これらより、幾何試行では両者の平均は成功回数分の
1
だけ差がある
E[M1]=E[N1]1

が、分散は等しい。
V[M1]=V[N1]

負の二項試行でも、平均は成功回数分の
r
だけ差がある
E[Mr]=E[Nr]r

が、分散は同様に等しい。
V[Mr]=V[Nr]