裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

js-STAR

2017年11月12日 | 統計学

js-STAR version 8.9.6j(β版)
ちゃんとプログラム実行結果の検証をしているのか?

みっともないぞ

 

[直接確率計算2×2]で,どうやったら,「両側検定  :  p=0.0000   ns (.10<p)」なんて結果になるんだ?

「カイ二乗検定の結果」で,4つのセルが全く同じ場合(例えば,10,10,10,10)でもカイ二乗値が0ではない( 0.1 になる)とは,なんたる惨状!

[直接確率計算2×2]

    観測値1    観測値2
-----------------------------------------
群1    10    9
群2    7    8
-----------------------------------------

両側検定  :  p=0.0000   ns (.10<p)

以下略

=============================

「カイ二乗検定の結果」
(上段実測値,下段期待値)
----------------
10        10        
 10.000     10.000    
----------------
10        10        
 10.000     10.000    

x2(1)=   0.100   ,  ns

 

=============================

他にも p 値の表示エラーはいくつかのプログラムにある。

[直接確率計算2×2]

40  2
 5 30

の検定結果が 「両側検定  :  p=6.0000   ** (p<.01)」と表示される

どうもこれは,

> fisher.test(matrix(c(40, 5, 2, 30), 2))

    Fisher's Exact Test for Count Data

data:  matrix(c(40, 5, 2, 30), 2)
p-value =
6.108e-14

の,指数表示されたときの p 値の最上位桁の 6 を拾っているようだ。

 

40  2
 5 31

の検定結果が「両側検定  :  p=3e-14.0000   ** (p<.01)」と表示されてしまうのが笑えるのだが。

> fisher.test(matrix(c(40, 5, 2, 31), 2))

    Fisher's Exact Test for Count Data

data:  matrix(c(40, 5, 2, 31), 2)
p-value =
3e-14

この例から見ると,「指数表示されたときの上位から小数点の手前までと ".0000"」を表示しているようだ。

コメント
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする

PVアクセスランキング にほんブログ村

PVアクセスランキング にほんブログ村