target.lrn <- read.table("../lrn/num/472.dat",header=T,colClasses="numeric")
target.val <- read.table("../val/num/472.dat",header=T,colClasses="numeric")
target.tst <- read.table("../tst/num/472.dat",header=T,colClasses="numeric")
 
y.lrn <- target.lrn[,1]
y.val <- target.val[,1]
y.tst <- target.tst[,1]
y     <- c(y.lrn,y.val,y.tst)

n.lrn <- length(y.lrn)
n.val <- length(y.val)
n.tst <- length(y.tst)
n     <- length(y)

rm(target.lrn,target.val,target.tst)

target.lrn <- read.table("yhat_lrn_05.dat",header=T,colClasses="numeric")
target.val <- read.table("yhat_val_05.dat",header=T,colClasses="numeric")
target.tst <- read.table("yhat_tst_05.dat",header=T,colClasses="numeric")
 
yhat.lrn <- target.lrn[,1]
yhat.val <- target.val[,1]
yhat.tst <- target.tst[,1]
yhat     <- c(yhat.lrn,yhat.val,yhat.tst)

idx.lrn <- 1:n.lrn
idx.val <- (n.lrn+1):(n.lrn+n.val)
idx.tst <- (n.lrn+n.val+1):n

y1 <- yhat.lrn
y1 <- y1-0.68
y1 <- sort(y1)
y1 <- y1[n.lrn:1]
y1 <- cumsum(y1) 

idx <- 1:n.lrn

print(paste("maximum profit in learning sample is ", max(y1)))
print(paste("maximum occurs at ", idx[y1==max(y1)]))

y1 <- y.lrn
y1 <- y1-0.68
y1 <- sort(y1)
y1 <- y1[n.lrn:1]
y1 <- cumsum(y1) 

idx <- 1:n.lrn

print(paste("maximum profit in learning sample is ", max(y1)))
print(paste("maximum occurs at ", idx[y1==max(y1)]))

