par(family="NimbusSan")

f1 <- as.matrix(read.csv("fig1a.csv",header=F))
postscript(file="fig1a.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
plot(f1[1,3:9],f1[2,3:9],ty="b",xlab="Confidence",ylab="Bet acceptance",
     ylim=c(0,1))
points(f1[1,3:9],f1[3,3:9],ty="b",lty=2)
legend(.8,.3,lty=c(1,2),legend=c("Answers","Random"))
text(.55,.9,"(a)")
dev.off()

f1 <- as.matrix(read.csv("fig1b.csv",header=F))
postscript(file="fig1b.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
plot(f1[1,3:9],f1[2,3:9],ty="b",xlab="Confidence",ylab="Bet acceptance",
     ylim=c(0,1))
points(f1[1,3:9],f1[3,3:9],ty="b",lty=2)
text(.55,.9,"(b)")
dev.off()

f1 <- as.matrix(read.csv("fig1c.csv",header=F))
postscript(file="fig1c.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
plot(f1[1,3:9],f1[2,3:9],ty="b",xlab="Confidence",ylab="Bet acceptance",
     ylim=c(0,1))
points(f1[1,3:9],f1[3,3:9],ty="b",lty=2)
text(.55,.9,"(c)")
dev.off()

f1 <- as.matrix(read.csv("fig1d.csv",header=F))
postscript(file="fig1d.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
plot(f1[1,3:9],f1[2,3:9],ty="b",xlab="Confidence",ylab="Bet acceptance",
     ylim=c(0,1))
points(f1[1,3:9],f1[3,3:9],ty="b",lty=2)
text(.55,.9,"(d)")
dev.off()

b1 <- as.matrix(read.csv("fig2a.csv"))
colnames(b1) <- c("City\ndistance","State\narea","State\npopulation","State\ pop.\ndensity","City\npopulation")
postscript(file="fig2a.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
barplot(b1,beside=T,ylim=c(.5,1),xpd=F,legend=T)
abline(h=.5)
text(1.9,.95,"(a)")
dev.off()

b1 <- as.matrix(read.csv("fig2b.csv"))
colnames(b1) <- c("City\ndistance","State\narea","State\npopulation","State\ pop.\ndensity","City\npopulation")
postscript(file="fig2b.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
barplot(b1,beside=T,ylim=c(.5,1),xpd=F,legend=T)
abline(h=.5)
text(1.9,.95,"(b)")
dev.off()

f1 <- as.matrix(read.csv("fig3a.csv",header=F))
postscript(file="fig3a.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
plot(f1[1,3:9],f1[2,3:9],ty="b",xlab="Confidence",ylab="Bet acceptance",
     ylim=c(.25,1),lty=3,cex=.75) # st pop
points(f1[1,3:9],f1[3,3:9],ty="b",lty=2,cex=1) # st area
points(f1[1,3:9],f1[4,3:9],ty="b",lty=4,cex=.5) # st dens
points(f1[1,3:9],f1[5,3:9],ty="b",lty=5,cex=.25) # cit pop
points(f1[1,3:9],f1[6,3:9],ty="b",lty=1,cex=1.25) # cit dens
legend(.8,.6,lty=1:5,pt.cex=c(1.25,1,.75,.5,.25),merge=T,pch=1,
       legend=c("City dist.","State area","State pop.","State dens.","City pop."))
text(.55,.9,"(a)")
dev.off()

f1 <- as.matrix(read.csv("fig3b.csv",header=F))
postscript(file="fig3b.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
#x11(width=8,height=6,pointsize=16)
plot(f1[1,3:9],f1[2,3:9],ty="b",xlab="Confidence",ylab="Bet acceptance",
     ylim=c(.25,1),lty=3,cex=.75) # st pop
points(f1[1,3:9],f1[3,3:9],ty="b",lty=2,cex=1) # st area
points(f1[1,3:9],f1[4,3:9],ty="b",lty=4,cex=.5) # st dens
points(f1[1,3:9],f1[5,3:9],ty="b",lty=5,cex=.25) # cit pop
points(f1[1,3:9],f1[6,3:9],ty="b",lty=1,cex=1.25) # cit dist
legend(.8,.6,lty=1:5,pt.cex=c(1.25,1,.75,.5,.25),merge=T,pch=1,
       legend=c("City dist.","State area","State pop.","State dens.","City pop."))
text(.55,.9,"(b)")
dev.off()

f1 <- as.matrix(read.csv("fig4.csv",header=F))
postscript(file="fig4.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
plot(f1[1,3:9],f1[2,3:9],ty="b",xlab="Confidence",ylab="Bet acceptance",
     ylim=c(0,1))
points(f1[1,3:9],f1[3,3:9],ty="b",lty=2)
legend(.8,.3,lty=c(1,2),legend=c("Answers","Random"))
dev.off()

b1 <- as.matrix(read.csv("fig5.csv"))
colnames(b1) <- NULL
postscript(file="fig5.eps",family="NimbusSan",width=8,height=6,
           horiz=F,onefile=F,pointsize=16,paper="special")
barplot(b1,beside=T,ylim=c(2,5.2),xpd=F,col=c("grey24","grey76"))
mtext(c("Question 1\nCompetence","\nQuestion 2\nRelative to others",
                  "Question 3\nControl","\nQuestion 4\nProbability alterability"),
      side=1,at=c(2,5,8,11),line=c(1.5,1.5,1.5,1.5))
legend(3.6,5,legend=c("Answers","Random"),fill=c("grey24","grey76"))
mtext("*",side=3,at=c(2,11),line=c(-2,-.6),cex=2)
abline(h=2)
dev.off()