par(family="NimbusSan") postscript(file="fig1a.eps",width=8,height=6, horiz=F,onefile=F,pointsize=16,paper="special") Y <- c(22,5 ,5 ,9 ,1 ,2 ,6 ,13,9 ,21,35,13,6 ,5 ,7 ,2 ,4 ,3 ,3 ,5 ,24) barplot(height=Y,ylab="Number of subjects",ylim=c(0,38)) text(c(.7,12.7,24.7),c(23.5,36.5,25.5),c("22","35","24")) mtext(c(0,10,20,30,40,50,60,70,80,90,100),side=1,line=.5,at=2.4*0:10+.7) mtext("Percent of trials",side=1,line=2,at=12.7) text(1,35,"(a)") dev.off() postscript(file="fig1b.eps",width=8,height=6, horiz=F,onefile=F,pointsize=16,paper="special") Y <- c(3 ,1 ,2 ,5 ,0 ,1 ,4 ,5 ,4 ,9 ,14,11,3 ,4 ,2 ,2 ,3 ,3 ,2 ,4 ,18) barplot(height=Y,ylab="Number of subjects",ylim=c(0,38)) text(c(12.7,24.7),c(15.5,19.5),c("14","18")) mtext(c(0,10,20,30,40,50,60,70,80,90,100),side=1,line=.5,at=2.4*0:10+.7) mtext("Percent of trials",side=1,line=2,at=12.7) text(1,35,"(b)") dev.off() postscript(file="fig1c.eps",width=8,height=6, horiz=F,onefile=F,pointsize=16,paper="special") Y <- c(19,4,3,4,1,1,2,8,5,12,21,2,3,1,5,0,1,0,1,1,6) barplot(height=Y,ylab="Number of subjects",ylim=c(0,38)) text(c(.7,12.7),c(20.5,22.5),c("19","21")) mtext(c(0,10,20,30,40,50,60,70,80,90,100),side=1,line=.5,at=2.4*0:10+.7) mtext("Percent of trials",side=1,line=2,at=12.7) text(1,35,"(c)") dev.off() postscript(file="fig3a.eps",width=8,height=6, horiz=F,onefile=F,pointsize=16,paper="special") Y <- c(4,2,5,1,1,3,0,7,9,6,32,10,4,1,8,2,7,6,4,8,16) barplot(height=Y,ylab="Number of subjects",ylim=c(0,38)) text(c(.7,12.7,24.7),c(0,33.5,17.5),c("","32","16")) mtext(c(0,10,20,30,40,50,60,70,80,90,100),side=1,line=.5,at=2.4*0:10+.7) mtext("Percent of trials",side=1,line=2,at=12.7) text(1,35,"(a)") dev.off() postscript(file="fig3b.eps",width=8,height=6, horiz=F,onefile=F,pointsize=16,paper="special") Y <- c(1,0,2,0,0,2,0,4,3,3,6,6,2,0,8,2,5,5,2,4,13) barplot(height=Y,ylab="Number of subjects",ylim=c(0,38)) text(c(.7,12.7,24.7),c(0,0,14.5),c("","","13")) mtext(c(0,10,20,30,40,50,60,70,80,90,100),side=1,line=.5,at=2.4*0:10+.7) mtext("Percent of trials",side=1,line=2,at=12.7) text(1,35,"(b)") dev.off() postscript(file="fig3c.eps",width=8,height=6, horiz=F,onefile=F,pointsize=16,paper="special") Y <- c(3,2,3,1,1,1,0,3,6,3,26,4,2,1,0,0,2,1,2,4,3) barplot(height=Y,ylab="Number of subjects",ylim=c(0,38)) text(c(.7,12.7,24.7),c(23.5,27.5,25.5),c("","26","")) mtext(c(0,10,20,30,40,50,60,70,80,90,100),side=1,line=.5,at=2.4*0:10+.7) mtext("Percent of trials",side=1,line=2,at=12.7) text(1,35,"(c)") dev.off() ## Below are the data and other information for figures 2 and 4. The changes ## on the proof look okay, except for table 2 appearing before table 1. ## I would just like to throw in one more last-minute change, if possible, to ## the captions for figures 2 and 4. I would like the caption for figure 2 ## to read, "Percentage of trials in which participants chose the optimal ## over the competitive option as a function of a) reinforcement, b) gender, ## and c) reinforcement by gender in Experiment 2 (table 2 lists the ## distribution of men and women across the conditions displayed in 2c). ## Error bars are ±1 SEM. * p < .05. *** p < .0001."; and figure 4's caption ## to read, "Percentage of trials in which participants chose the optimal ## over the competitive option as a function of a) reinforcement, b) gender, ## and c) reinforcement by gender in Experiment 2 (there are 34 participants ## in each of the four conditions displayed in 4c). Error bars are ±1 SEM. * ## p < .05. *** p < .0001." library(gplots) h=c(38.6,58.8) ciu=c(2.8,2.9)+h cil=h-c(2.8,2.9) postscript(file="fig2a.eps",width=8,height=8, horiz=F,onefile=F,pointsize=19,paper="special") barplot2(height=h,ylim=c(0,100),ylab="Percent optimal",ci.u=ciu,ci.l=cil, names.arg=c("Non-monetary","Monetary"),plot.ci=T,space=.5,xlim=c(0,3.5)) text(c(1.6,3.2),c(42,62),c("38.6%","58.8% ***")) text(0.25,95,"(a)") dev.off() h=c(45.7,57.6) ciu=c(2.3,4.7)+h cil=h-c(2.3,4.7) postscript(width=8,height=8,pointsize=19, file="fig2b.eps",horiz=F,onefile=F,paper="special") # x11(width=8,height=8,pointsize=19) barplot2(height=h,ylim=c(0,100),ylab="Percent optimal",ci.u=ciu,ci.l=cil, names.arg=c("Women (n=150)","Men (n=50)"),plot.ci=T,space=.5,xlim=c(0,3.5)) text(c(1.6,3.2),c(50.1,62),c("47.5%","57.6% *")) text(0.25,95,"(b)") dev.off() h=matrix(c(38.0,39.1,53.6,73.5),2) ciu=h+c(3.0,6.9,3.3,4.7) cil=h-c(3.0,6.9,3.3,4.7) postscript(width=8,height=8,pointsize=19, file="fig2c.eps",horiz=F,onefile=F,paper="special") #x11(width=8,height=8,pointsize=19) barplot2(height=h,ylim=c(0,100),ylab="Percent optimal",ci.u=ciu,ci.l=cil,legend=T, names.arg=c("Non-monetary","Monetary"),plot.ci=T,beside=T,col=c("gray16","gray84")) legend(2,90,legend=c("Women (n=150)","Men (n=50)"),fill=c("gray16","gray84")) text(c(1.5,2.5,4.5,5.5),c(h)+10,c("38.0%","39.1%","53.6%","73.5%")) text(1.25,95,"(c)") dev.off() h=c(49.9,67.8) ciu=h+c(3.1,3.2) cil=h-c(3.1,3.2) postscript(file="fig4a.eps",width=8,height=8, horiz=F,onefile=F,pointsize=19,paper="special") barplot2(height=h,ylim=c(0,100),ylab="Percent optimal",ci.u=ciu,ci.l=cil, names.arg=c("Non-monetary","Monetary"),plot.ci=T,space=.5,xlim=c(0,3.5)) text(c(1.6,3.2),h+3,c("49.9%","67.8% ***")) text(0.25,95,"(a)") dev.off() h=c(53.9,63.8) ciu=h+c(3.3,3.2) cil=h-c(3.3,3.2) postscript(width=8,height=8,pointsize=19, file="fig4b.eps",horiz=F,onefile=F,paper="special") # x11(width=8,height=8,pointsize=19) barplot2(height=h,ylim=c(0,100),ylab="Percent optimal",ci.u=ciu,ci.l=cil, names.arg=c("Women (n=68)","Men (n=68)"),plot.ci=T,space=.5,xlim=c(0,3.5)) text(c(1.6,3.2),h+3,c("53.9%","63.8% *")) text(0.25,95,"(b)") dev.off() h=matrix(c(45.9,53.8,61.9,73.7),2) ciu=h+c(4.1,4.5,4.8,3.9) cil=h-c(4.1,4.5,4.8,3.9) postscript(width=8,height=8,pointsize=19, file="fig4c.eps",horiz=F,onefile=F,paper="special") #x11(width=8,height=8,pointsize=19) barplot2(height=h,ylim=c(0,100),ylab="Percent optimal",ci.u=ciu,ci.l=cil,legend=T, names.arg=c("Non-monetary","Monetary"),plot.ci=T,beside=T,col=c("gray16","gray84")) legend(2,95,legend=c("Women (n=68)","Men (n=68)"),fill=c("gray16","gray84")) text(c(1.5,2.5,4.5,5.5),c(h)+10,c("45.9%","53.8%","61.9%","73.7%")) text(1.25,95,"(c)") dev.off()