Я пытался найти ответы, как добавить несколько осей Y к гистограмме с различными способами иллюстрации моих данных, но не смог найти то, что хочу сделать. Пока мне удалось создать красивый штриховой график 1 . 
Теперь я хочу добавить больше оси Y, чтобы добавить другие переменные в виде точек над моими полосами, как на этом графике 2 .
Мой код пока выглядит так:
dataH <- read.delim("~/Speciale/Data - Narsarsuaq/barplot_H.txt")
bilanH <- aggregate(cbind(X189,X197,X205,X212)~treatment , data=dataH , mean)
rownames(bilanH) <- bilanH[,1]
bilanH <- as.matrix(bilanH[,-1])
lim <- 30
error.bar <- function(x, y, upper, lower=upper, length=0.1,...){
arrows(x,y+upper, x, y-lower, angle=90, code=3, length=length, ...)
}
stdevH <- aggregate(cbind(X189,X197,X205,X212)~treatment , data=dataH , sd)
rownames(stdevH) <- stdevH[,1]
stdevH <- as.matrix(stdevH[,-1]) * 1.96 / 10
ze_barplot <- barplot(bilanH , beside = T ,legend.text=T,density=c(0,15,15,15) , angle=c(0,45,0,135) , col="black" ,main="Total emission - High elevation", ylim=c(0,lim) , ylab="Emission rate (ng cm-2 leaf area h-1)")
barplotH <- error.bar(ze_barplot,bilanH, stdevH)
Есть ли простой способ добавить точки к моим графикам?