「統計学関連なんでもあり」の過去ログ---024

★ JavaScript版その他による正確確率検定の結果の不一致 ★

　374　JavaScript版その他による正確確率検定の結果の不一致　　わたなべ　　2003/05/12 (月) 15:04
　　375　Re: JavaScript版その他による正確確率検定の結果の不一致　　青木繁伸　　2003/05/12 (月) 15:21
　　　377　Re^2: JavaScript版その他による正確確率検定の結果の不一致　　わたなべ　　2003/05/12 (月) 15:38

374.　JavaScript版その他による正確確率検定の結果の不一致　　わたなべ　　2003/05/12 (月) 15:04

お伺いいたします.
こちらのサイトの統計「補助道具」を大変便利に使わせていただいています.
（既出でなければよいのですが）RxC分割表に関するFisherの正確確率検定の結果が，JavaScript，Exact testのページ，そしてRで異なりました.
（Peason, Fisherの各方法の違いというわけではありません.）
ちなみに，後二者は一致してました（JavaScript版が異なる）.
データは

1 101 2 42 0 102 （上の方法で順にp=0.021/p=0.082/p=0.082）
という，幾分問題が生じそうなデータですが，用いられているアルゴリズムが異なるのでしょうか?
端的に言って，どれが正確なのでしょうか?
ご回答，よろしくお願いします.

　　　　　[このページのトップへ]

375.　Re: JavaScript版その他による正確確率検定の結果の不一致　　青木繁伸　　2003/05/12 (月) 15:21

私がやってみたところでは，0.082 になりました。

>criterion=11.4772
>Fisher の方法による正確な P 値 = 0.08211360878
>査察した分割表の個数は 10 個

となりました。
使用したのは，Macintosh版のNetscape 7.02 です。
生起確率を加えるときの判定に実数値の比較をするのですが，JavaScript の実装上，同じであるべき数値が同じになっていない（マシンエプシロンレベルの違い）のでしょう。
これは，解決の難しい問題です。R を使うのが一番無難でしょうね。

　　　　　[このページのトップへ]

377.　Re^2: JavaScript版その他による正確確率検定の結果の不一致　　わたなべ　　2003/05/12 (月) 15:38

早速ありがとうございました.

> 私がやってみたところでは，0.082 になりました。
...
> 使用したのは，Macintosh版のNetscape 7.02 です。
> 生起確率を加えるときの判定に実数値の比較をするのですが，JavaScript の実装上，同じであるべき数値が同じになっていない（マシンエプシロンレベルの違い）のでしょう。
> これは，解決の難しい問題です。R を使うのが一番無難でしょうね。

なるほど.
MacのSafariで動かしていました.
確認してみましたところ，Explore5.2(Mac版)でも同じ値（0.082）となりました.
Safariがまだベータ版だからか，それ以前の問題なのか，詳しいことは分かりませんが，
注意したいと思います.
ありがとうございました!

　　　　　[このページのトップへ]