Rで解析：同等性検定の実行「equivalence」パッケージ

t検定などでグループ間の平均に「意味のある差がある」ことを示すことがあると思います。でも、逆に「意味のある差がない」ことを統計的に証明するのは非常に難しいことです。

専門家で多くの議論がなされているようですが、個人的にはグループ間に統計的有意差がないから「グループは同じである」と結論付けて良いのか疑問です。まあ、いろいろな条件はありますが。

そんな問題を解決するかもしれない、同等性検定を実行できる「equivalence」パッケージからtoseコマンドを紹介します。多くのコマンドが収録されていますのでパッケージヘルプで確認してください。

パッケージのバージョンは0.7.2。実行コマンドはwindows 11のR version 4.1.3で確認しています。

パッケージのインストール

下記コマンドを実行してください。

#パッケージのインストール
install.packages("equivalence")

実行コマンドの紹介

詳細はコメント、パッケージヘルプを確認してください。

#&#12497;&#12483;&#12465;&#12540;&#12472;&#12398;&#35501;&#12415;&#36796;&#12415;
library("equivalence")

#&#12487;&#12540;&#12479;&#20363;&#12398;&#20316;&#25104;
set.seed(1234)
Data1 <- sample(1:5, 10, replace = TRUE)
Data2 <- sample(1:5, 10, replace = TRUE)

#&#24179;&#22343;&#20516;&#12391;&#21516;&#31561;&#24615;&#12434;&#26908;&#23450;:tost&#12467;&#12510;&#12531;&#12489;
Test <- tost(Data1, Data2, epsilon = 0.2)

#&#27161;&#28310;&#35492;&#24046;&#12398;&#24046;&#20998;
Test$se.diff
[1] 0.5744563

#tost&#12467;&#12510;&#12531;&#12489;&#12391;&#35373;&#23450;&#12375;&#12383;&#28961;&#35222;&#12377;&#12427;&#24046;
Test$epsilon
[1] 0.2

少しでも、あなたの解析が楽になりますように！！