2

我一直在使用包 multcompview 来直观地显示基于 anova 和 HSD tukey 测试的组之间的显着差异:

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') ) )

在此处输入图像描述

这是我的一个实际图表的一个示例: 在此处输入图像描述 该方法(我在 SO 中发布相关问题后发现)效果非常好,但我还没有找到能够添加 y 轴标签的方法(我需要标记Y 变量“排名漏洞”)。multcomp 函数似乎不接受 ylab 参数。在没有基本轴标签信息的情况下拥有这些整体好看的对比图令人沮丧……你知道这个问题的解决方案/解决方法吗?

4

2 回答 2

2

我也试过这个不错的,但不是自我解释的包。我发现当你想要标记轴或主标题时,你必须在 multcompBoxplot() 之后启动 title()。例如:

title(ylab = 'Response', main = 'Title')
于 2015-01-15T11:58:16.907 回答
1

在 SO 的任何其他人都没有回复并四处寻找各种复杂性不同的解决方案后,我发现了一个效果很好且非常简单的解决方案(使用 mtext):

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') ) )
mtext(side = 2, "Response", line = 2.3, cex=1.5)

在此处输入图像描述

于 2013-06-13T00:19:29.913 回答