効果の大きさのペアごとの比較を取得しようとしています。でこれを行うことができますがcoh_d、繰り返し比較が行われます。たとえば、次のコードでは、setosa と versicolor は versicolor と setosa と同じです (反転した負/正符号を除いて)。
library(esvis)
iris<- iris
coh_d(Sepal.Length ~ Species, data=iris)
1.) 両方の組み合わせは報告する必要がありますか? そうでない場合、重複した比較を削除する方法はありますか?
2.) 各比較の出力に識別列を追加する方法はありますか。上記のコードでは、Sepal.Length というラベルが付けられた 6 つの行が追加されるため、複数の列の効果サイズを見つけると、比較がどの変数に属しているかがわかります。最終的な出力を次のようにしたいと思います。
a <- colnames(iris)
a <- a[1:4]
effect_fun<- function(y,x){
form2<- as.formula(paste0(y, "~", x))
res2<- effect<- rstatix::kruskal_effsize(data = iris, form2, conf.level = 0.95)
return(res2)
}
effect.sizes<- lapply(a, FUN = effect_fun, x="Species")
effect.sizes<- do.call(rbind, effect.sizes)
View(effect.sizes)
を使用するkruskal_effsizeと、必要な出力の形式が得られますが、ペアごとの比較が表示されないため、に切り替えましたcoh_d。
どんな助けでも大歓迎です。