2021年3月21日日曜日

ベイズ推定(東工大2021年入試・総合型選抜の問題)

 取りうる事象が2種類である試行を独立に n 回行ったときに、一方の事象が何回起こるかの確率分布を二項分布という。二項分布に従う典型的な確率変数としては、コイン投げによって表(もしくは裏)が出る回数が挙げられる。
 このように、我々の生活の中には様々な確率分布が存在する。ただ、残念ながら確率がわかっているという状況は現実では極めて稀なケースであり、現実には確率がわかっておらず、データから分布を推定しなければならないケースがある。
 今、形が歪んでおり、表が出る確率がわからないコインを考える。そのわからない確率をθとおく(ただし 0<θ<1)。コインを投げる回数を n、表が出る回数を k とする。このとき「このコインを使用する」という前提のもとで表が k 回出る条件付き確率は
で表される。ここで、0!=1 , nC0=1 とする。さて、いまθの値がわからないので、限られた回数だけコインを実際に投げて得たデータから確率θを推定したい。以下の問いに答えよ。ただし、各問ごとにコインの歪み具合は異なるとする。

問1 確率θの値を推定する方法の一つに、尤度(ゆうど)関数と呼ばれる概念を使った最尤推定法という方法がある。まず、尤度関数とは、上記の条件付き確率 P をθの関数とみなし、n や k には実際に投げて得たデータの値を代入して評価したものである。いま、ある歪んだコインを3回投げて3回とも表が出た。このときの尤度関数を L(θ) とおくとき、L(θ) を求めよ。

問2 最尤推定法とは、尤度関数が最大値をとるようなθをその推定値とするような推定法である。ある歪んだコインを5回投げたところ、表が出た回数について以下の表のような結果を得た。

     表/裏(コイン1)│ 表 裏
    ──────────┼─────
        回数    │ 4 1

このとき、尤度関数を求めたうえで、上で説明した最尤推定法に従ってθの推定値θ^ を求めよ。

問3 一般に、n 回コインを投げたときに k 回表が出たとする。このとき、最尤推定法を用いて得たθの推定値θ^ は k/n となることを証明せよ。

 次の問4以降は、歪み方の異なる2つのコイン(コイン1とコイン2)を考える。それぞれ表の出る確率がθ1 , θ2 であるが、前問と同様にそれらの値はわからないとする。Aさんがコイン1については20回、コイン2については50回投げたところ、表が出た回数について以下の表のような結果を得た。

     表/裏(コイン1)│ 表 裏
    ──────────┼─────
        回数    │ 11 9

     表/裏(コイン2)│ 表 裏
    ──────────┼─────
        回数    │ 24 26

 次に、別のBさんがAさんから1 , 2いずれかのコインを1つ受け取ったが、見た目が同じなためにどちらのコインを受け取ったかはわからず、したがって当初Bさんは「1/2 の確率でコイン1、1/2 の確率でコイン2」だと考えたとする。また、Bさんは、Aさんがそれぞれのコインを投げて得た上記の結果は知っているとする。
 その後、Bさんが自らそのコインを10回投げたところ、出た回数について以下の表のような結果を得た。

     表/裏(コイン1 or 2) │ 表 裏
    ─────────────┼─────
        回数       │ 5 5

さて、この結果からBさんは受け取ったコインが1か2かをどのように判断するべきだろうか。以下ではこれを考察しよう。

問4 一般的に、X が起こる確率を P(X)、Y が起こる確率を P(Y)、そして X と Y が同時に起こる確率を P(X,Y) とするとき、「X が実現した、という条件の下で Y が実現する確率」は
      P(Y|X)=P(X,Y)÷P(X)
で与えられ、これを Y の X に関する条件付き確率という。P(Y|X) は以下のようにも表現できることを示せ。
      P(Y|X)=P(X|Y)P(Y)÷P(X)
これをベイズの定理という。

問5 Bさんが得た結果「表が5、裏が5」を条件としたときの、Bさんが受けとったコインが1である条件付き確率と2である条件付き確率を比較したとき、もし前者が大きければBさんが受け取ったのはコイン1、後者が大きければコイン2と判断するとしよう。この場合、Bさんのコインは1,2のどちらと判断するのが妥当か、根拠も含めて答えよ。



問1  L(θ)=θ3

問2  L(θ)=5θ4(1ーθ) → L'(θ)=20θ3ー25θ4=0 → θ^=4/5

問3  θk(1ーθ)nーk

問4  P(X,Y)=P(Y|X)P(X)=P(X|Y)P(Y)

問5  

0 件のコメント:

コメントを投稿