R の記事も多くなったけど，問題のあるものもある

> My Notes

> 統計学とかR(R言語)とかPython3とかプログラミングとかの覚え書きとか走り書きとか。 座右の銘にしたい: All work and no play makes Jack a dull boy.

> 2017-09-14

> R(R言語) ノンパラメトリック検定(適合度検定、カイ二乗検定、基準値との比較、一様性の検定)、chisq.test()

http://my-notes.hatenablog.com/entry/2017/09/14/001050

で，観測値が理論比に従うかの検定についての記事がある。

＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝

コードは冗長であるが，要するに，

あるカテゴリー変数について母比率が p = c(0.4, 0.2, 0.3, 0.1) であるとき，

標本での観察値が x = c(152, 102, 58, 188) のとき，標本比率が母比率と同じか？という検定をすること。

筆者は，

> x = c(152, 102, 58, 188)
> chisq.test(x)
    Chi-squared test for given probabilities

data:  x
X-squared = 77.728, df = 3, p-value < 2.2e-16

としてしまっている。これでは，母比率が一様である，つまり  p= c(0.25, 0.25, 0.25, 0.25) の検定をしていることになる。

> chisq.test(x, p= c(0.25, 0.25, 0.25, 0.25))

    Chi-squared test for given probabilities

data:  x
X-squared = 77.728, df = 3, p-value < 2.2e-16 # X-squared が同じ値

しかし，母比率は p = c(0.4, 0.2, 0.3, 0.1) なのだから，これは間違っている。

> chisq.test(x, p= c(0.4, 0.2, 0.3, 0.1))

    Chi-squared test for given probabilities

data:  x
X-squared = 448.87, df = 3, p-value < 2.2e-16

とすべき。

まあ，all or none の結論から言えばどちらも，帰無仮説を棄却することには成るんだが，間違いは間違い。

R を自習するということなら問題ないが，間違った自習結果を公表するのは問題がある。

些細なことではあるが，大きな問題は，「間違いを指摘しようにも，コメントを受け付ける態勢にない」ということ。

天上天下唯我独尊というのはいただけないなあ。（前には，コメントを受け付けていて，コメントしたこともあると思うのだが，それが不愉快だったのでコメント拒否にしたのかな）