tx <-read.table("../regr/lif_x_net_val.dat",header=T,colClasses="numeric")
ty <-read.table("../regr/lif_y_net_val.dat",header=T,colClasses="numeric")

xn.val <- tx[,1]
yn.val <- ty[,1]

tx <-read.table("../regr/lif_x_gross_val.dat",header=T,colClasses="numeric")
ty <-read.table("../regr/lif_y_gross_val.dat",header=T,colClasses="numeric")

xg.val <- tx[,1]
yg.val <- ty[,1]

rm(tx,ty)

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_intr.dat",header=T,colClasses="numeric")
target.val <- read.table("yhat_val_intr.dat",header=T,colClasses="numeric")
target.tst <- read.table("yhat_tst_intr.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

source("psopts.r");

idx <- seq(1,n.lrn,length=500)

y0.lrn <- c(0,sum(y.lrn-0.68))
x0.lrn <- c(0,n.lrn)

idx.yhat <- sort.list(yhat.lrn,decreasing=TRUE)
y1.lrn <- y.lrn[idx.yhat]-0.68 
y1.lrn <- cumsum(y1.lrn)
y1.lrn <- y1.lrn[idx]
x1.lrn <- 1:n.lrn
x1.lrn <- x1.lrn[idx]

postscript(file="cty_lif_intr_lrn.eps")

plot(x=c(x0.lrn,x1.lrn),y=c(y0.lrn,y1.lrn),
  ylab="dollars",xlab="size of mailing",type="n")
lines(x=x0.lrn,y=y0.lrn,col="green")
lines(x=x1.lrn,y=y1.lrn,col="red")

dev.off()

idx <- seq(1,n.val,length=500)

y0.val <- c(0,sum(y.val-0.68))
x0.val <- c(0,n.val)

idx.yhat <- sort.list(yhat.val,decreasing=TRUE)
y1.val <- y.val[idx.yhat]-0.68 
y1.val <- cumsum(y1.val)
y1.val <- y1.val[idx]
x1.val <- 1:n.val
x1.val <- x1.val[idx]

postscript(file="cty_lif_intr_val.eps")

plot(x=c(x0.val,x1.val),y=c(y0.val,y1.val),
  ylab="dollars",xlab="size of mailing",type="n")
lines(x=x0.val,y=y0.val,col="green")
lines(x=x1.val,y=y1.val,col="blue")

dev.off()

y0.lrn <- (y0.lrn/y0.lrn[2])*100
x0.lrn <- (x0.lrn/x0.lrn[2])*100

y1.lrn <- (y1.lrn/y1.lrn[500])*100
x1.lrn <- (x1.lrn/x1.lrn[500])*100

y1.val <- (y1.val/y1.val[500])*100
x1.val <- (x1.val/x1.val[500])*100

postscript(file="cty_lif_intr_both.eps")

plot(x=c(x0.lrn,x1.lrn,x1.val,xn.val),y=c(y0.lrn,y1.lrn,y1.val,yn.val),
  ylab="percent",xlab="percent",type="n")
lines(x=x0.lrn,y=y0.lrn,col="green")
lines(x=x1.lrn,y=y1.lrn,col="red")
lines(x=x1.val,y=y1.val,col="blue")
#lines(x=xn.val,y=yn.val,col="orange")

dev.off()

postscript(file="cty_lif_intr_rn.eps")

plot(x=c(x0.lrn,x1.lrn,x1.val,xn.val),y=c(y0.lrn,y1.lrn,y1.val,yn.val),
  ylab="percent",xlab="percent",type="n")
lines(x=x0.lrn,y=y0.lrn,col="green")
#lines(x=x1.lrn,y=y1.lrn,col="red")
lines(x=x1.val,y=y1.val,col="blue")
lines(x=xn.val,y=yn.val,col="orange")

dev.off()


idx <- seq(1,n.lrn,length=500)

y0.lrn <- c(0,sum(y.lrn))
x0.lrn <- c(0,n.lrn)

idx.yhat <- sort.list(yhat.lrn,decreasing=TRUE)
y1.lrn <- y.lrn[idx.yhat]
y1.lrn <- cumsum(y1.lrn)
y1.lrn <- y1.lrn[idx]
x1.lrn <- 1:n.lrn
x1.lrn <- x1.lrn[idx]

idx <- seq(1,n.val,length=500)

y0.val <- c(0,sum(y.val))
x0.val <- c(0,n.val)

idx.yhat <- sort.list(yhat.val,decreasing=TRUE)
y1.val <- y.val[idx.yhat]
y1.val <- cumsum(y1.val)
y1.val <- y1.val[idx]
x1.val <- 1:n.val
x1.val <- x1.val[idx]

y0.lrn <- (y0.lrn/y0.lrn[2])*100
x0.lrn <- (x0.lrn/x0.lrn[2])*100

y1.lrn <- (y1.lrn/y1.lrn[500])*100
x1.lrn <- (x1.lrn/x1.lrn[500])*100

y1.val <- (y1.val/y1.val[500])*100
x1.val <- (x1.val/x1.val[500])*100

postscript(file="cty_lif_intr_conv.eps")

plot(x=c(x0.lrn,x1.lrn,x1.val,xg.val),y=c(y0.lrn,y1.lrn,y1.val,yg.val),
  ylab="percent",xlab="percent",type="n")
lines(x=x0.lrn,y=y0.lrn,col="green")
lines(x=x1.lrn,y=y1.lrn,col="red")
lines(x=x1.val,y=y1.val,col="blue")
#lines(x=xg.val,y=yg.val,col="orange")

dev.off()

postscript(file="cty_lif_intr_rg.eps")

plot(x=c(x0.lrn,x1.lrn,x1.val,xg.val),y=c(y0.lrn,y1.lrn,y1.val,yg.val),
  ylab="percent",xlab="percent",type="n")
lines(x=x0.lrn,y=y0.lrn,col="green")
#lines(x=x1.lrn,y=y1.lrn,col="red")
lines(x=x1.val,y=y1.val,col="blue")
lines(x=xg.val,y=yg.val,col="orange")

dev.off()