我一直在使用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变量标记为"ranked“)。multcomp函数似乎不接受ylab参数。没有基本的轴标签信息…的这些整体好看的对比度图是令人沮丧的您是否知道此问题的解决方案/解决方法?
发布于 2015-01-15 19:58:17
我也尝试过这个很好的包,但不是自我解释的包。我发现,当你想要标记轴或主标题时,你必须在multcompBoxplot()之后启动title()。例如:
title(ylab = 'Response', main = 'Title')发布于 2013-06-13 08:19:30
在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)

https://stackoverflow.com/questions/16427175
复制相似问题