`OneWayANOVATest`, `LeveneTest`, `FlignerKilleenTest` が使っている anova() 関数にバグがある。
全体の平均値を求めるときに,各群の平均値の平均を求めてしまっている。
Z̄ = mean(Z̄ᵢ)
これでは,各群のサンプルサイズが等しくない場合には正しい答えにならないのは自明である。
Z̄ = sum(Z̄ᵢ .* Nᵢ) / sum(Nᵢ).
としなければならない。
したがって `anova()` を使う関数 `OneWayANOVATest()`, `LeveneTest()`, `FlignerKilleenTest()` は,すべて誤った答えを出してしまう。
`OneWayANOVATest()` については,少なくとも 2021/08/02 にバグの指摘があったようだが,作者には伝わっっていないようである。
今回の指摘は,伝わるだろうか?一応 GitHub ではコメントしたのだけど。
# Julia 伝道師を気取っていたが,このようなことが度重なると,ちょっと困るなあ。
# Julia の熱心な使用者が少ないのが根本原因だろうなあ。
# 使用者が多ければ,おかしいと気づく人もそれなりに出てくるだろうから。
# 見逃されているというのは,見ている人が少ないからだけだろう。
※コメント投稿者のブログIDはブログ作成者のみに通知されます