anova および HSD tukey 検定に基づいて、グループ間の有意差を視覚的に表示できるように、パッケージ multcompview を使用しています。
Group=c("G1","G1","G1","G1","G2","G2","G2","G2","G3","G3","G3","G3")
set.seed(0)
Vals=c(runif(4),runif(4)+0.7,runif(4)-0.7)
data=data.frame(Group)
data=cbind(data, Vals)
library(multcompView)
xzx <-multcompBoxplot(Vals~Group,data=data,sortFn=median, decreasing=FALSE,
horizontal=FALSE,
plotList=list(
boxplot=list(fig=c(0, 1, 0, 1), las=3,
cex.axis=1.5),
multcompLetters=list(
fig=c(0.87, 0.97, 0.115, 0.923), #0.1108, 0.9432 Top of
#page 18 manual for very convoluted explanation (c(y bottom, y top,x L, x R))
type='Letters') ) )
私の実際のグラフの 1 つの例を次に示します。
アプローチ (関連する質問を SO に投稿した後に見つけたもの) は非常にうまく機能しますが、y 軸ラベルを追加できる方法が見つかりませんでした ( Y変数「ランク付けされた脆弱性」)。multcomp 関数は ylab 引数を取り込めないようです。基本的な軸ラベル情報なしで、これらの全体的に見栄えの良いコントラスト グラフがあるのはイライラします... この問題の解決策/回避策を知っていますか?