ポスト

#統計 Rのfisher.testの信頼区間は本質的に上のP値関数から計算されているが、P値は別の定義で計算されている。(ゆえに信頼区間とP値が整合的でない。)

メニューを開く

黒木玄 Gen Kuroki@genkuroki

みんなのコメント

メニューを開く

#統計 fnh = FNH(a+c, b+d, a+b, OR)とおくと、Rのfisher.testのP値は pvalue(a,b,c,d|OR) = sum(pdf(fnh, x) for x in support(fnh) if pdf(fnh, x)⪅pdf(fnh, a)) のOR=1の場合を使って計算されている。離散分布の数値計算で≤を使うと失敗するので、⪅と書いている。純粋数学的には≤でよい。

黒木玄 Gen Kuroki@genkuroki

人気ポスト

もっと見る
Yahoo!リアルタイム検索アプリ