当您使用功能barplot2()
从库gplots
,将使用这种方法给出的例子。
首先,在帮助文件barplot2()
函数中给出了barplot。 ci.l
和ci.u
是假置信区间值。 Barplot应该保存为对象。
hh <- t(VADeaths)[1:2, 5:1]
mybarcol <- "gray20"
ci.l <- hh * 0.85
ci.u <- hh * 1.15
mp <- barplot2(hh, beside = TRUE,
col = c("grey12", "grey82"),
legend = colnames(VADeaths)[1:2], ylim = c(0, 100),
cex.names = 1.5, plot.ci = TRUE, ci.l = ci.l, ci.u = ci.u)
如果您看对象mp
,它包含所有酒吧的x坐标。
mp
[,1] [,2] [,3] [,4] [,5]
[1,] 1.5 4.5 7.5 10.5 13.5
[2,] 2.5 5.5 8.5 11.5 14.5
现在我使用上置信区间值来计算段的y值的坐标。分段将从比置信区间结束高1的位置开始。 y.cord
包含四行 - 第一行和第二行对应第一个栏,其他两行对应第二个栏。最高的y值是从每个柱对的置信区间的最大值计算出来的。 x.cord
值只是在mp
对象中重复相同的值,每个值为2次。
y.cord<-rbind(c(ci.u[1,]+1),c(apply(ci.u,2,max)+5),
c(apply(ci.u,2,max)+5),c(ci.u[2,]+1))
x.cord<-apply(mp,2,function(x) rep(x,each=2))
barplot由使用sapply()
后使5个线段使用计算出的坐标(因为此时有5组)。
sapply(1:5,function(x) lines(x.cord[,x],y.cord[,x]))
要绘制的线段上述文本计算x和y坐标,其中x是两个条形的x值的中间点和y值被从置信区间为每个杆对加上一些恒定的极大值来计算。然后使用功能text()
添加信息。
x.text<-colMeans(mp)
y.text<-apply(ci.u,2,max)+7
text(c("*","**","***","NS","***"),x=x.text,y=y.text)
multcomp中有一个plot.cld函数,您可以将字母放在您的酒吧上方指示重要性。 Perhabs这也适合你... – EDi 2013-03-20 22:28:07
还有'agricolae'包中的'bar.group',它为你打上字母。 – mnel 2013-03-20 22:35:05
如果您使用base R的'barplot',则可以存储像barstore < - barplot(1:3)'这样的条的中心点。为了验证这是否正常,请尝试'abline(v = barstore)'并注意垂直线都切断了条的中心。使用'段'可以使用这些存储点来绘制比较/交互线。 – thelatemail 2013-03-20 23:43:11