R : Copyright 2002, The R Development Core Team
Version 1.5.1  (2002-06-17)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type `license()' or `licence()' for distribution details.

R is a collaborative project with many contributors.
Type `contributors()' for more information.

Type `demo()' for some demos, `help()' for on-line help, or
`help.start()' for a HTML browser interface to help.
Type `q()' to quit R.

> invisible(options(echo = TRUE))
> plt <- read.table("bv_plt.dat",
+          header=F,colClasses="numeric",col.names=c("cp","mse.lrn","mse.val"))
> 
> plt <- plt[1:11,]
> plt$cp <- -log(plt$cp)
> xlim <- c(plt$cp,plt$cp)
> ylim <- c(plt$mse.lrn,plt$mse.val)
> 
> source("psopts.r");
> postscript(file="bv_plt.eps")
> plot(x=xlim,y=ylim,type="n",xlab="-log(cp)",ylab="mse")
> points(x=plt$cp,y=plt$mse.lrn,pch="o",cex=1.0)
> points(x=plt$cp,y=plt$mse.val,pch="o",cex=1.0)
> lines(x=plt$cp,y=plt$mse.lrn,type="l",lty="solid",col="red",lwd=2)
> lines(x=plt$cp,y=plt$mse.val,type="l",lty="solid",col="blue",lwd=2)
> dev.off()
null device 
          1 
> 
> proc.time()
[1] 0.98 0.08 1.04 0.00 0.00
>