source("mlp.r") b <- as.matrix (read.table("bdata",head=FALSE)) X <- b[1:400,1:2] y <- b[1:400,3] X.tst <- b[401:500,1:2] y.tst <- b[401:500,3] test.err <- function (y.tst, p) { mean (log (ifelse(y.tst==1,p,1-p))) } cat("Average test log probability with training mean", test.err(y.tst,mean(y)), "\n") r<- mlp.cv (y, X, 0.002, 0.002, 1:100, 10000, 10, TRUE) p1 <- mlp.predict (X.tst, r) cat("Average test log probabiity:", test.err(y.tst,p1), "\n")