library("gplots") m.pg <- c(-2.6400, 3.6000, 6.0000, 3.6800, 5.4400) se.pg <- c(1.71938, 1.86548, 1.74738, 1.94484, 1.83492) m.pl <- c(-4.9600, -3.7600, -2.3200, -.1600, 6.5600) se.pl <- c(1.47024, 1.72170, 1.79139, 1.36587, 1.56852) postscript(file="fig1.eps",family="NimbusSan",width=8,height=8,horiz=F,pointsize=18,paper="special") plotCI(y=c(m.pg,m.pl),x=c(c(1:5)-.01,c(1:5)+.01),uiw=c(se.pg,se.pl),ylim=c(-6,8), ylab="Net IGT score",xlab="Block",lty=rep(c(1,2),c(5,5))) lines(y=m.pg,x=c(1:5)-.01,lty=1) lines(y=m.pl,x=c(1:5)+.01,lty=2) legend(3.6,-3.7,legend=c("Prior gain","Prior loss"),lty=c(1,2)) dev.off()