対立仮説の下での正しい数の推測の分布は、非中心超幾何分布に従います。これは、オッズ比の観点からパラメーター化されています。つまり、女性が実際にミルクが最初に(またはその逆に)追加されたときとは対照的に、実際には最初にお茶が最初に追加されました。オッズ比が1の場合、中央の超幾何分布が得られます。
これがうまくいくか見てみましょう。(非中心)超幾何分布の密度を計算するMCMCpack
機能を持つパッケージを使用して、説明のためにRを使用しますdnoncenhypergeom()
。これは、引数がありx
推測の正確な数について(注意してください:茶は本当に最初に追加されたときに、この2つの条件のいずれかの下で推測の正確な数は、例えば、ある)、引数をn1
、n2
とm1
4つの余白の3のために、とpsi
のために真のオッズ比。x
真のオッズ比が1である場合、0から4(すべてのマージンが4に等しい)の密度を計算してみましょう。
install.packages("MCMCpack")
library(MCMCpack)
sapply(0:4, function(x) dnoncenhypergeom(x, n1=4, n2=4, m1=4, psi=1))
これにより、
[1] 0.01428571 0.22857143 0.51428571 0.22857143 0.01428571
したがって、帰無仮説の下で、女性が8回の正しい推測を行う確率は1.43%です(つまり、お茶が最初に追加された場所では4カップすべてが正しく推測され、したがって、ミルクが最初に追加された場所では4カップすべてが正しく推測されます)。これは実際、フィッシャーが帰無仮説を棄却するのに十分であると考えた証拠の量です。
(.90/(1−.90))/(.10/(1−.10))=81odds(guessA|trueA)/odds(guessA|trueB))。女性が8カップすべてを正しく推測する可能性はどのくらいですか(つまり、彼女はお茶が最初に追加された場所で4カップすべてを正しく推測し、したがって牛乳が最初に追加された場所で4カップも正しく推測します)。
dnoncenhypergeom(4, n1=4, n2=4, m1=4, psi=81)
これにより、
[1] 0.8312221
したがって、電力は約83%です。