library(gplots) #m1 <- matrix(10-(c(261,250,219,235,269,265)-241)/19.1,2,3) #e1 <- matrix(10-(c(253,241,204,207,252,251)-241)/19.1,2,3) #f1 <- 2*m1-e1 # lower f1 <- rbind(c(9.35, 11.89, 9.41), c(10.7, 11.75, 9.42)) e1 <- rbind(c(8.45, 10.3, 7.59),c(8.3, 8.85, 7.98)) m1 <- rbind(c(8.9, 11.1, 8.5),c(9.5, 10.3, 8.7)) postscript(file="fig2.eps",family="NimbusSan",width=12,height=8, horiz=F,onefile=F,pointsize=18,paper="special") barplot2(m1,plot.ci=T,ci.u=e1,ci.l=f1,beside=T,col=(c("gray40","gray80")), ylab="Seconds",xlab="",ylim=c(0,14),grid.lwd = 2, plot.grid=TRUE) mtext(c("Emotion prevalent","Evaluation prevalent","Neutral in prevalence"), side=1,line=.5,at=c(2,5,8)) legend(6.25,14.1,bty="n",fill=c("gray40","gray80"),c("Emotion question","Moral question")) dev.off() #grid.col="white", border="white")