ベイズ統計の仮説検定〜頻度論の考え方に基づく検定〜【第3回】2017/10/102020/04/14ベイズ統計の仮説検定〜頻度論の考え方に基づく検定〜【第3回】ベイズ統計ライター:古澤嘉啓TwitterFacebookベイズ統計の仮説検定は全6回にわたって説明しています。他のページもぜひご覧ください。前回の前説でもお話しした通り、「やむを得ない理由がない限り、帰無仮説を棄却すべきではない」という従来の考え方に基づいたベイズ流の仮説検定の説明を行っていきます。前回の【第2回】ベイズ統計の仮説検定〜基本的な検定〜に引き続き、このグラフを例に説明していきます。従来の仮説検定の考え方は、ネイマン-ピアソン流の考え方で、1第1種の誤り確率を一定以下(主に5%など)に抑えながら2第2種の誤り確率を小さくするというものです。第1種の誤り確率は、「帰無仮説が正しいのにも関わらず、帰無仮説を棄却してしまう確率」であり、1を満たすことが従来の仮説検定の前提条件になります。この前提条件から、「あまりにもおかしなデータが得られない限り、帰無仮説を棄却するのは避けよう」という考えが見えてきます。これが初めに書いた「やむを得ない理由がない限り、帰無仮説を棄却すべきではない」という考え方につながっています。(参照→仮説検定とは?初心者にもわかりやすく解説!)しかし、上記のグラフで行われたベイズ流の仮説検定では、この前提条件に従っていません。第1種の誤り確率ではありませんが、「帰無仮説を満たす確率」が40%もあるのです。これは「やむを得ない理由がない限り、帰無仮説を棄却すべきではない」という考え方に反しています。そこで、「帰無仮説を満たす確率」を一定以下に抑えることを考えます。これが、頻度論の考え方に基づくベイズ流仮説検定になります。満たす確率」を一定以下に抑える事後オッズ比を用いて検定する従来の仮説検定では有意水準を5%とし検定を行うのがメジャーであるので、今回のベイズ流仮説検定でも、「帰無仮説を満たす確率」を5%以下に抑えることを考えていきます。ここで、事後オッズ比という概念を導入します。事後オッズ比は、事後分布において、「対立仮説を満たす確率よりも、帰無仮説を満たす確率がどのくらい大きいか」つまり、$$\frac{P(H_0|X)}{P(H_1|X)}$$で表されるつまり先の例でいうと、事後オッズ比は\(\frac{P(H_0|X)}{P(H_1|X)}=\frac{40\%}{60\%}=\frac{0.4}{0.6}=\frac{2}{3}\)と計算できます。この事後オッズ比が、\(\frac{1}{19}\)以下の場合帰無仮説を棄却し、\(\frac{1}{19}\)より上の場合帰無仮説を受容します。\(\frac{1}{19}\)が基準なのでしょうか。それはグラフを描いてみればすぐわかります。これが\(\frac{1}{19}\)の様子です。\(100\%×\frac{1}{1+19}=5\%\)で、赤い部分の確率が5%になっているとわかります。事後オッズ比が\(\frac{1}{19}\)以下なら、帰無仮説を満たす確率が5%以下であると言えるのです。実際に例題を用いて検定する〜頻度論に考え方に基づくベイズ流仮説検定〜では、実際に例題を用いて計算して見ましょう。成人男性の平均身長が170cmより大きいかどうかを調べたい。いま、成人男性の身長は\(N(\mu,10^2)\)に従っているとする。ひとりのデータ取ってきたところ、177cmであった。成人男性の平均身長が170cm以下であると言えるか。\(\mu\)は事前に\(N(172,20^2)\)に従っているとし、しきい値を\(\frac{1}{19}\)とする。仮説は次のようになります。170\)さて、まずは事後分布を計算しましょう。事後分布の平均、分散は以下のinfoから次のように計算できます。info母平均\(\mu\)、母分散\(\sigma^2\)に従う正規母集団から大きさ\(n\)の標本を抽出し、標本平均\(\overline{x}\)を得たとする。170|X=177)=1-0.251=0.749\)が得られます。よって、事後オッズ比は\(\frac{P(H_0|X)}{P(H_1|X)}=\frac{0.251}{0.749}=0.335\)となります。0.335は\(\frac{1}{19}\)よりも大きいので、帰無仮説を受容し、成人男性の平均身長は170cmより大きいとは言えない、と判断します。まとめいかがでしたでしょうか。このように、頻度論に考え方に基づくベイズ流仮説検定では、事後オッズ比を求めて検定していきます。次回は、検定の証拠の強さを表すベイズファクターについて説明していきます。G検定を取得するメリットを解説!最弱オセロを初めて攻略した天才オセロ高校生。TrendCOMMENTコメントを残すコメントをキャンセルメールアドレスが公開されることはありません。コメント*お名前*メールアドレス*プライバシーポリシーはこちら上に表示された文字を入力してください。