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()