##### R-code for making Figure 3.1. It reads data from
##### "example31.summary" generated by the code 
##### "example31.r"

dat <- read.table("example31.summary")

xbar <- dat$V1[1:20]
rbar <- dat$V2[1:20]
i <- seq(1,20)

postscript("fig31.ps",width=7.2,height=3.6,horizontal=F)

par(mfrow=c(1,2), mar=c(4,4,1,2))

xbar <- dat$V1[1:24]
rbar <- dat$V2[1:24]
i <- seq(1,24)

U=84.577
C=79.533
L=74.489

plot(i[1:20],xbar[1:20],type="o",lty=1,pch=16,xlab="i",
     ylab=expression(bar(X)[i]),mgp=c(2,1,0),xlim=c(0,25), 
     ylim=c(70,86), cex=0.8)
lines(i[21:24],xbar[21:24],type="o",lty=3,pch=1,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(74.489,79.533,84.577),labels=c("L","C","U"),cex=0.8,las=1)
title(xlab="(a)",cex=0.9)

U=18.490
C=8.745
L=0

plot(i[1:20],rbar[1:20],type="o",lty=1,pch=16,xlab="i",
     ylab=expression(bar(R)[i]),mgp=c(2,1,0),xlim=c(0,25), 
     ylim=c(-1,25), cex=0.8)
lines(i[21:24],rbar[21:24],type="o",lty=3,pch=1,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(0,8.745,18.49),labels=c("L","C","U"),cex=0.8,las=1)
title(xlab="(b)",cex=0.9)

graphics.off()