2

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 引数を取り込めないようです。基本的な軸ラベル情報なしで、これらの全体的に見栄えの良いコントラスト グラフがあるのはイライラします... この問題の解決策/回避策を知っていますか?

4

2 に答える 2