複数の比較箱ひげ図の有意性バー (アスタリスク *、 *、* )を描画するためのループまたは関数を作成した人はいますか?
複数の変数 (19) と次のコードを持つデータ フレームを使用しています。
ボックスプロット:
boxplot$STATUS は、4 つのレベル (A、B、C、D) を持つ変数です。
par(mfrow=c(5,4))
for (i in c(1:19)){
boxplot(log(cito[,i]+1)~cito$STATUS)
}
TukeyHDT 事後テスト
TukeyHSD(aovp(log(cito$CRP+1)~cito$STATUS,center=T))
(cito$STATUS)そこで、 TukeyHSD で得られたグループ間の有意差の結果を boxplot のループにプロットしたいと思います。
助けていただければ幸いです。しばらくの間、私はオプションを探し続けます。