## |          |No        |External-M|Internal-M|No        |External- |Internal- |
## |          |event-Mean|ean       |ean       |event-Std.|Std. Error|Std. Error|
## |          |          |          |          |Error     |          |          |
## |Implausibl|6.792     |6.321     |7.033     |0.230     |0.175     |0.210     |
## |e         |          |          |          |          |          |          |
## |Plausible |6.415     |5.213     |5.975     |0.210     |0.180     |0.232     |

library(gplots)

d1 <- matrix(c(6.792, 6.321, 7.033, 6.415, 5.213, 5.975),3,2)
d2 <- matrix(c(0.230, 0.175, 0.210, 0.210, 0.180, 0.232),3,2)

postscript(file="fig1.eps",width=8,height=8,
horiz=F,onefile=F,pointsize=16,paper="special")
barplot2(height=d1,plot.ci=T,ci.l=d1-d2,ci.u=d1+d2,xpd=F,
         beside=T,names.arg=c("Implausible","Plausible"),
         legend.text=c("No event","External","Internal"),
         col=c("gray90","gray65","gray40"),ylim=c(1,9),
         ylab="Need-for-treatment judgments (1-9)")
dev.off()