R : Copyright 2003, The R Development Core Team
Version 1.6.2  (2003-01-10)

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))
> 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()
null device 
          1 
> 
> 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()
null device 
          1 
> 
> 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()
null device 
          1 
> 
> 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()
null device 
          1 
> 
> 
> 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()
null device 
          1 
> 
> 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()
null device 
          1 
> 
> 
> proc.time()
[1] 1.96 0.11 2.08 0.00 0.00
>