#"Cooperation (y)"
#"Mean""S.E."
#"CI (x)","High Mean","Low Mean","High Mean","Low Mean"
f2a <- matrix(c(1,0.25,0.05,0.11,0.02,
2,0.19,0.12,0.09,0.05,
3,0.25,0.11,0.08,0.04,
4,0.31,0.16,0.08,0.07,
5,0.38,0.22,0.11,0.07,
6,0.44,0.28,0.14,0.1,
7,0.36,0.34,0.11,0.09,
8,0.37,0.25,0.11,0.09,
9,0.39,0.34,0.12,0.1),9,5,byrow=T)

#,"Prediction (y)"
#"Mean""S.E.",
#"CI (x)","High Mean","Low Mean","High Mean","Low Mean"
f2b <- matrix(c(1,0.59,0.16,0.1,0.04,
2,0.57,0.24,0.09,0.03,
3,0.6,0.27,0.1,0.05,
4,0.61,0.29,0.08,0.05,
5,0.71,0.43,0.09,0.07,
6,0.79,0.49,0.06,0.06,
7,0.75,0.53,0.05,0.04,
8,0.8,0.47,0.04,0.11,
9,0.8,0.66,0.04,0.08),9.5,byrow=T)

#,"Cooperation (y)"
#"Mean""S.E.",
#"CI (x)","High Range","Low Range","High Range","Low Range",
f3a <- matrix(c(1,0.2,NA,0.04,NA,
2,NA,NA,NA,NA,
3,NA,NA,NA,NA,
4,NA,0.38,NA,0.08,
5,0.37,0.46,0.07,0.09,
6,0.42,NA,0.1,NA,
7,NA,NA,NA,NA,
8,NA,NA,NA,NA,
9,NA,0.66,NA,0.11),9,5,byrow=T)

#,"Prediction (y)"
#"Mean""S.E.",
#"CI (x)","High Range","Low Range","High Range","Low Range"
f3b <- matrix(c(1,0.46,NA,0.09,NA,
2,NA,NA,NA,NA,
3,NA,NA,NA,NA,
4,NA,0.56,NA,0.07,
5,0.5,0.65,0.08,0.07,
6,0.5,NA,0.09,NA,
7,NA,NA,NA,NA,
8,NA,NA,NA,NA,
9,NA,0.82,NA,0.05),9,5,byrow=T)

#,"Cooperation (y)"
#"Mean""S.E.",
#"CI (x)","High Rank","Low Rank","High Rank","Low Rank"
f4a <- matrix(c(1,0.07,0.23,0.06,0.13,
2,0.17,NA,0.05,NA,
3,0.18,NA,0.03,NA,
4,0.27,NA,0.07,NA,
5,0.24,0.27,0.04,0.09,
6,NA,0.29,NA,0.11,
7,NA,0.33,NA,0.13,
8,NA,0.43,NA,0.12,
9,0.33,0.37,0.07,0.13),9,5,byrow=T)

#,"Prediction (y)"
#"Mean""S.E.",
#"CI (x)","High Rank","Low Rank","High Rank","Low Rank"
f4b <- matrix(c(1,0.32,0.48,0.06,0.06,
2,0.36,NA,0.05,NA,
3,0.45,NA,0.03,NA,
4,0.51,NA,0.05,NA,
5,0.54,0.53,0.04,0.03,
6,NA,0.56,NA,0.05,
7,NA,0.62,NA,0.06,
8,NA,0.65,NA,0.05,
9,0.69,0.71,0.04,0.07),9,5,byrow=T)

rownames(f2a) <- rownames(f2b) <- rownames(f3a) <- rownames(f3b) <-
  rownames(f4a) <- rownames(f4b) <- 1:9

library(gplots)

postscript(file="fig2a.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=18,paper="special")
plotCI(x=1:9-.05,y=f2a[,2],uiw=f2a[,4],xaxt="n",ylim=c(0,1),type="b",
       xlab="Cooperation Index",ylab="Proportion",main="Cooperation")
mtext(1:9,side=1,at=1:9,line=1)
plotCI(x=1:9+.05,y=f2a[,3],uiw=f2a[,5]+c(.018,0,0,0,0,0,0,0,0),ylim=c(0,1),lty=2,add=T)
lines(x=1:9+.05,y=f2a[,3],ylim=c(0,1),lty=2,type="b")
dev.off()

postscript(file="fig2b.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=18,paper="special")
plotCI(x=1:9-.05,y=f2b[,2],uiw=f2b[,4],xaxt="n",ylim=c(0,1),type="b",
       xlab="Cooperation Index",ylab="Proportion",main="Prediction")
mtext(1:9,side=1,at=1:9,line=1)
plotCI(x=1:9+.05,y=f2b[,3],uiw=f2b[,5]+c(.005,.008,0,0,0,0,0,0,0),ylim=c(0,1),lty=2,add=T)
lines(x=1:9+.05,y=f2b[,3],ylim=c(0,1),lty=2,type="b")
dev.off()

postscript(file="fig3a.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=18,paper="special")
barplot2(height=t(f3a[,2:3]),beside=T,xlab="Cooperation Index",ylab="Proportion",
         ylim=c(0,1),legend=c("Low Rank","High Rank"),main="Cooperation",
         plot.ci=T,ci.u=t(f3a[,2:3]+f3a[,4:5]),ci.l=t(f3a[,2:3]-f3a[,4:5]))
dev.off()

postscript(file="fig3b.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=18,paper="special")
barplot2(height=t(f3b[,2:3]),beside=T,axes=F,xlab="Cooperation Index",ylim=c(0,1),
         plot.ci=T,ci.u=t(f3b[,2:3]+f3b[,4:5]),ci.l=t(f3b[,2:3]-f3b[,4:5]),
         main="Prediction")
dev.off()

postscript(file="fig4a.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=18,paper="special")
barplot2(height=t(f4a[,2:3]),beside=T,xlab="Cooperation Index",ylab="Proportion",
         ylim=c(0,1),legend=c("Low Rank","High Rank"),main="Cooperation",
         plot.ci=T,ci.u=t(f4a[,2:3]+f4a[,4:5]),ci.l=t(f4a[,2:3]-f4a[,4:5]))
# smartlegend(x="left",y="top",c("Low Rank","High Rank"),fill=c("red","yellow"))
dev.off()

postscript(file="fig4b.eps",family="NimbusSan",width=8,height=8,
horiz=F,onefile=F,pointsize=18,paper="special")
barplot2(height=t(f4b[,2:3]),beside=T,axes=F,xlab="Cooperation Index",ylim=c(0,1),
         plot.ci=T,ci.u=t(f4b[,2:3]+f4b[,4:5]),ci.l=t(f4b[,2:3]-f4b[,4:5]),
         main="Prediction")
dev.off()