裏 RjpWiki

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

Julia の `OneWayANOVATest`, `LeveneTest`, `FlignerKilleenTest` は,全滅!!!

2022年05月14日 | Julia

`OneWayANOVATest`, `LeveneTest`, `FlignerKilleenTest`   が使っている anova() 関数にバグがある。

全体の平均値を求めるときに,各群の平均値の平均を求めてしまっている。

    Z̄ = mean(Z̄ᵢ)

これでは,各群のサンプルサイズが等しくない場合には正しい答えにならないのは自明である。

    Z̄ = sum(Z̄ᵢ .* Nᵢ) / sum(Nᵢ).

としなければならない。

したがって `anova()` を使う関数 `OneWayANOVATest()`, `LeveneTest()`, `FlignerKilleenTest()` は,すべて誤った答えを出してしまう。

`OneWayANOVATest()` については,少なくとも 2021/08/02 にバグの指摘があったようだが,作者には伝わっっていないようである。

今回の指摘は,伝わるだろうか?一応 GitHub ではコメントしたのだけど。

# Julia 伝道師を気取っていたが,このようなことが度重なると,ちょっと困るなあ。
# Julia の熱心な使用者が少ないのが根本原因だろうなあ。
# 使用者が多ければ,おかしいと気づく人もそれなりに出てくるだろうから。
# 見逃されているというのは,見ている人が少ないからだけだろう。

コメント    この記事についてブログを書く
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする
« 一応,技術文書だと思うのに | トップ | BrownForsytheTest よ,お前... »
最新の画像もっと見る

コメントを投稿

Julia」カテゴリの最新記事