# fig1
Ch=14
postscript(file="fig1.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=16,paper="special")
plot(c(0,110),c(0,100),type="n",axes=F,xlab="",ylab="")
rect(0,80,20,100,col="gray80")
rect(0+Ch,80-Ch,20+Ch,100-Ch,col="gray80")
rect(0+2*Ch,80-2*Ch,20+2*Ch,100-2*Ch,col="gray80")
rect(0+3*Ch,80-3*Ch,20+3*Ch,100-3*Ch,col="gray80")
rect(0+4*Ch,80-4*Ch,20+4*Ch,100-4*Ch,col="gray80")
text(20,98,pos=4,labels="fixation cue appears")
text(20+Ch,98-Ch,pos=4,labels="saccad targets appear")
text(20+2*Ch,98-2*Ch,pos=4,labels="go cue presented")
text(20+3*Ch,98-3*Ch,pos=4,labels="saccade executed")
text(20+4*Ch,98-4*Ch,pos=4,labels="reward delivered")
points(c(10, 10+Ch,10+Ch-6,10+Ch+6, 10+2*Ch-6,10+2*Ch+6, 10+3*Ch-6),
       c(90, 90-Ch,90-Ch-6,90-Ch+6, 90-2*Ch-6,90-2*Ch+6, 90-3*Ch-6),
       pch=20)
arrows(0,70,3.5*Ch,70-3.5*Ch,lwd=2)
text(25,80-3*Ch,labels="time")
arrows(10+3*Ch,90-3*Ch,10+3*Ch-5,90-3*Ch-5,length=0.1)
par(new=TRUE)
xspline(87+c(0,2,0,-2),33+c(4,0,-2,0),open=F,shape=c(0,1,1,1))
dev.off()

## Figure 2:

## Probability of choosing risky over alternating, population: 0.858 +/- 0.007
## Probability of choosing alternating over risky, population: 0.142 +/- 0.007

## Probability of choosing risky over alternating, monkey N: 0.797 +/- 0.014
## Probability of choosing alternating over risky, monkey N: 0.203 +/- 0.014

## Probability of choosing risky over alternating, monkey B: 0.957 +/- 0.005
## Probability of choosing alternating over risky, monkey B: 0.043 +/- 0.005

## Probability of choosing risky over alternating, monkey E: 0.875 +/- 0.010
## Probability of choosing alternating over risky, monkey E: 0.125 +/- 0.010

library(gplots)
postscript(file="fig2.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=16,paper="special")
m1=matrix(c(.858,.142,.797,.203,.957,.043,.875,.125),2,)
c1=matrix(c(.007,.007,.014,.014,.005,.005,.010,.010),2,)
barplot2(m1,ci.u=m1+c1,ci.l=m1-c1,plot.ci=T,beside=T,space=c(0,3),
         ylab="choice frequency",
         col=c("white","gray90"),xlim=c(2,22),ylim=c(0,1.05),
         names.arg=c("population","monkey N","monkey B","monkey E"))
text(c(2,3.6,7,8.6,12,13.6,17,18.6)+1.5,m1+c1,labels=c("risky","alternating"),
     pos=3)
dev.off()

## Figure 3:

## Change in large reward, -35 : 0.571 +/- 0.010
## Change in large reward, +0 : 0.624 +/- 0.006
## Change in large reward, +35 : 0.654 +/- 0.009

## Change in small reward, -35 : 0.637 +/- 0.009
## Change in small reward, +0 : 0.612 +/- 0.006
## Change in small reward, +35 : 0.619 +/- 0.012

postscript(file="fig3.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=16,paper="special")
plotCI(c(.571,.624,.654),uiw=c(.01,.006,.009)+.0005,type="b",
       ylab="probability of choosing risky option",xlab="",
       ylim=c(.5,.7),xaxt="n")
plotCI(c(.637,.612,.619),uiw=c(.009,.006,.012)+.0005,add=T)
points(c(.637,.612,.619),type="b",lty=2)
mtext(c("-35","0","+35"),side=1,at=1:3)
mtext(expression(paste("relative expected value of risky option (",
    symbol("\155"),"L juice)",sep="")),side=1,line=1.5)
text(c(1.1,1.1),c(.57,.64),c("change in large reward","change in small reward"),
     pos=4)
dev.off()