Rによるフィッシャーの直接確率検定

Rにてフィッシャーの直接確率検定 (Fisher exact test) を行う．直接確率検定は正確確率検定と訳されることも多い．推定統計学の父，Ronald Aylmer Fisher によって開発された．本検定法は分割表 (クロス集計表) における2つ以上のカテゴリーが独立であるかどうかを調べる手法である．すなわち，観測された2つ以上の分類が互いに無関係であるか，または何らかの関係性を有しているかを検定する方法である．同様の検定を行うものにカイ二乗検定 (適合度検定・独立性の検定) があるが，分割表のセルの期待値に10未満のものがある場合はカイ二乗検定は避け，フィッシャーの直接確率検定を用いるほうが良いとされる．フィッシャーの直接確率検定のコンセプトは非常にシンプルである．本検定法では，観測された分割表に対し，観測され得る値の全出現パターンの生起確率をそれぞれ超幾何分布にて計算し，それらの中で観測された分割表の生起確率よりも小さい生起確率を有する分割表の生起確率の和をp値と定義する．その名の通り，データの出現パターンを直接計算して確率を割り出すナイーブな手法であるが，それ故にデータ数が大きい場合は非常に大きな計算量が必要になり，計算機なしでは計算できなくなる．R ではデフォルトでインストールされている関数 fisher.test() を用いることで実行できる．

以下のように2x2の分割表が得られたとき，行および列で示される2要因が互いに独立であるかどうかは以下のように検定する．帰無仮説 (H₀) は2要因間に独立性が成り立つこと，すなわち，2要因の比率に差がない，または，2要因が互いに無関係であることである．有意水準5%にて検定する．

p値の計算

最初に，Rを起動させ，以下のコマンドにて，上の分割表を変数 mx に格納する．

$ mx=matrix(c(2, 10, 15, 3), nrow=2, byrow=T)

フィッシャーの直接確率検定はコマンド fisher.test() にて実行する．この関数は R にデフォルトでインストールされているパッケージに含まれている．検定は，最も簡単には，以下のようにコマンドを打つ．

$ fisher.test(mx)

これを実行した結果は以下のようになる．

 
        Fisher's Exact Test for Count Data

data:  mx
p-value = 0.0005367
alternative hypothesis: true odds ratio is not equal to 1
95 percent confidence interval:
 0.003325764 0.363182271
sample estimates:
odds ratio 
0.04693661 

以上の結果では，p値が0.0005367となっており，有意水準を5%で帰無仮説が棄却され，行および列の2要因の間には何らかの関連性があることが示唆された．alternative hypothesis とある行には，対立仮説が明記されている．上の場合の対立仮説は，オッズ比が1とは有意に異なること，である．オッズ比とは上の2x2の分割表における ad/cd を計算した値であり，要因間の関連の強さをあらわす指標である．独立性の検定とは超幾何分布によるオッズ比の検定ともいえる．この対立仮説のオッズ比の値は以下の or コマンドで別の値に指定することが可能である．結果の次の行には，オッズ比の95%信頼区間が表示される．計算させる信頼区間の幅等は以下のオプションで変更できる．最後の行の値は，超幾何分布に基づいたオッズ比の最尤推定量である．

また，上は2x2の分割表における検定であるが，以下のような2x2より大きい分割表の検定も実行できる．

以下のように実行する．

$ mx=matrix(c(2, 10, 8, 15, 3, 2), nrow=2, byrow=T)

$ fisher.test(mx)

これを実行した結果は以下のようになる．

 
        Fisher's Exact Test for Count Data

data:  mx
p-value = 0.0001025
alternative hypothesis: two.sided

このように，2x2より大きい分割表を対象にした検定では，信頼区間やオッズ比は計算されない．結果の p-value と書かれている値がp値である．その下には，検定を両側検定で行ったことが明記されている．両側検定を行うか片側検定を行うかは以下のオプションで選択できる．

オプション一覧

オプションには以下のようなものがある．

正規分布

標準化

パラメトリックとノンパラメトリック