##### R-code for making Figure 3.5. This code needs to
##### read data from "example31.summary" generated by 
##### the R-code "Example31.r"

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

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

postscript("fig35.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]
sbar <- dat$V3[1:24]
i <- seq(1,24)

U=84.636
C=79.533
L=74.430

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.43,79.533,84.636),labels=c("L","C","U"),cex=0.8,las=1)
title(xlab="(a)",cex=0.9)

U=7.468
C=3.575
L=0

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

graphics.off()