##### R-code for making Figure 3.7

xbar = c(0.14,0.06,0.20,0.02,0.16,0.10,0.08,0.18,0.06,0.18,0.10,0.14,
         0.04,0.20,0.08,0.12,0.18,0.06,0.22,0.10)
i <- seq(1,20)

U=0.259
C=0.121
L=0

postscript("fig37.ps",width=4.5,height=4,horizontal=F)

plot(i,xbar,type="o",lty=1,pch=16,xlab="i",ylab=expression(p[i]),
     mgp=c(2,1,0),xlim=c(0,20), ylim=c(-0.01,0.27),cex=0.8)
lines(i,rep(U,length(i)),lty=2,cex=0.8)
lines(i,rep(C,length(i)),lty=2,cex=0.8)
lines(i,rep(L,length(i)),lty=2,cex=0.8)
axis(4,at=c(L,C,U),labels=c("L","C","U"),cex=0.8,las=1)

graphics.off()