# program arrl.prg (No.2: final plots) # # par(mfrow=c(1,1)) par(mar=c(7,10,7,10)) nage<-length(unique(xs)) avasmat<-matrix(0,nage,4) avasmat[,1]<-unique(xs) avasmat[,2]<-unique(exp(lowx)-cons) # :1 avasmat[,3]<-unique(exp(mea)-cons) # :2 avasmat[,4]<-unique(exp(upp)-cons) # ;3 # plot(x,exp(yval)-cons,ylim=c(ylow,yup),xlab=xname,ylab=yname,pch=1) #: 4 plot(x,exp(yval)-cons,ylim=c(ylow,yup),xlab=xname,ylab=yname,pch=1) #: 4 # # ++++ When untransformed data are used, the above parts 1-4 should be # ++++ replaced by the following parts 1-4. # # avasmat[,2]<-unique(lowx) # :1 # avasmat[,3]<-unique(mea) # :2 # avasmat[,4]<-unique(upp) # :3 # plot(x,yval,ylim=c(ylow,yup),xlab=xname,ylab=yname,pch=1) # :4 # lines(avasmat[,1],avasmat[,2],col=5) lines(avasmat[,1],avasmat[,3],col=18) lines(avasmat[,1],avasmat[,4],col=5) # # end of program No.2