#"Figure 2 Data"    
# "C" "UP" "UF" "UPF" 
#"N" 98.3 57.6 61.6 40
#"B" 97.7 63.2 69 47.5
#"B/O" 97.2 55.7 80.4 54.8
m2 <- t(matrix(c(
98.3,57.6,61.6,40,
97.7,63.2,69,47.5,
97.2,55.7,80.4,54.8),4,3))
colnames(m2) <- c("C","UP","UF","UPF")
par(mar=c(3,5,2,1,))
barplot(m2,beside=T,cex.names=1.5,cex.axis=1.5,ylim=c(0,100),cex.lab=1.5,
        ylab="Percent prudent decisions",col=c("gray50","gray70","gray90"))
text(c(1,2,3,5,6,7,9,10,11,13,14,15)+.5,10,labels=format(c(m2)),cex=1.5,srt=90)
legend(4.5,100,legend=c("No feedback","Behavioral","Behavioral/Outcome"),
              bty="n",cex=1.5,fill=c("gray50","gray70","gray90"))

#"Figure,3,Data",
m3 <- matrix(c("One",86.7,77.7,76.6,
"Two",97.5,94.3,96.6,
"Three",98.6,97.8,98.1,
"Four",98.3,97.7,97.2),4,4)
par(mar=c(5,5,1,1,))
plot(as.numeric(m3[2,]),axes=F,frame.plot=T,ylim=c(0,100),type="b",cex.lab=1.5,lty=1,
     xlab="Fifty-trial block",ylab="Percent prudent decisions")
lines(as.numeric(m3[3,]),type="b",lty=2)
lines(as.numeric(m3[4,]),type="b",lty=3)
axis(1,at=1:4,labels=m3[1,],cex.axis=1.5)
axis(2,at=10*0:10,labels=c(0,"",20,"",40,"",60,"",80,"",100),,cex.axis=1.5)
legend(2,50,lty=1:3,legend=c("No feedback","Behavioral","Behavioral/Outcome"),
              bty="n",cex=1.4)

#"Figure,4,Data",
m4 <- matrix(c("One",40.6,49.7,45.7,
"Two",44.6,53.8,53.4,
"Three",50,56,52.4,
"Four",57.6,63.2,55.7),4,4)
par(mar=c(5,5,1,1,))
plot(as.numeric(m4[2,]),axes=F,frame.plot=T,ylim=c(0,100),type="b",cex.lab=1.5,lty=1,
     xlab="Fifty-trial block",ylab="Percent prudent decisions")
lines(as.numeric(m4[3,]),type="b",lty=2)
lines(as.numeric(m4[4,]),type="b",lty=3)
axis(1,at=1:4,labels=m4[1,],cex.axis=1.5)
axis(2,at=10*0:10,labels=c(0,"",20,"",40,"",60,"",80,"",100),,cex.axis=1.5)
legend(1.5,90,lty=1:3,legend=c("No feedback","Behavioral","Behavioral/Outcome"),
              bty="n",cex=1.4)
dev.copy2eps(file="fig4.eps")

#"Figure,5,Data",
m5 <- matrix(c("One",46.1,39.1,53.8,
"Two",52.5,57,70.9,
"Three",60.2,63.1,75,
"Four",61.6,69,80.4),4,4)
par(mar=c(5,5,1,1,))
plot(as.numeric(m5[2,]),axes=F,frame.plot=T,ylim=c(0,100),type="b",cex.lab=1.5,lty=1,
     xlab="Fifty-trial block",ylab="Percent prudent decisions")
lines(as.numeric(m5[3,]),type="b",lty=2)
lines(as.numeric(m5[4,]),type="b",lty=3)
axis(1,at=1:4,labels=m5[1,],cex.axis=1.5)
axis(2,at=10*0:10,labels=c(0,"",20,"",40,"",60,"",80,"",100),,cex.axis=1.5)
legend(2,50,lty=1:3,legend=c("No feedback","Behavioral","Behavioral/Outcome"),
              bty="n",cex=1.4)
dev.copy2eps(file="fig5.eps")

#"Figure,6,Data",
m6 <- matrix(c("One",42.8,43.4,45.4,
"Two",41.7,46.8,49.3,
"Three",37.5,49.5,51.1,
"Four",40,47.5,54.8),4,4)
par(mar=c(5,5,1,1,))
plot(as.numeric(m6[2,]),axes=F,frame.plot=T,ylim=c(0,100),type="b",cex.lab=1.5,lty=1,
     xlab="Fifty-trial block",ylab="Percent prudent decisions")
lines(as.numeric(m6[3,]),type="b",lty=2)
lines(as.numeric(m6[4,]),type="b",lty=3)
axis(1,at=1:4,labels=m6[1,],cex.axis=1.5)
axis(2,at=10*0:10,labels=c(0,"",20,"",40,"",60,"",80,"",100),,cex.axis=1.5)
legend(1.5,90,lty=1:3,legend=c("No feedback","Behavioral","Behavioral/Outcome"),
              bty="n",cex=1.4)
dev.copy2eps(file="fig6.eps")