我一直在使用包 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 参数。在没有基本轴标签信息的情况下拥有这些整体好看的对比图令人沮丧……你知道这个问题的解决方案/解决方法吗?