###getfigures.r
#
#  Generate Figures for the bkr paper.
#  Use results in directory ./output/*

##Draw boxplots for Circle 20 lambdas.
c20gam <- read.csv("output/Circle20.lambda.csv")
names(c20gam) <- c(1:20);
postscript("c20lambdaboxgam.ps");
boxplot(c20gam, ylab=expression(lambda), outline=T);
dev.off();

c20mix <- read.csv("output/Circle20.lambda.mix.csv")
names(c20mix) <- c(1:20);
postscript("c20lambdaboxmix.ps");
boxplot(c20mix, ylab=expression(lambda), outline=T);
dev.off();

##Draw boxplots for Boston Housing 20 lambdas.
bh <- read.csv("output/BostonHousing.lambda.csv");
bhname <- c("CRIM", "ZN", "INDUS", "CHAS", "NOX", "RM", "AGE",
            "DIS", "RAD", "TAX", "PTRATIO", "B", "LSTAT");
names(bh) <- rep("", length(bhname));
postscript("bhlambdabox.ps");
boxplot(bh, ylab=expression(lambda), outline=T);
mtext(bhname, side=1, line=1, at=1:length(bhname), las=3);
dev.off();

##Obtain Tables for the bkr paper
#
getmean <- function(fname){
  fname <- paste("output/", fname, ".comp.txt", sep="");
  tab <- as.matrix(read.table(fname, header=T, sep=","));
  mintab <- apply(tab, 1, min);
  for(i in 1:dim(tab)[1]){
    tab[i, ] <- tab[i, ]/mintab[i];
  }
  return(apply(tab, 2, mean));
}

getrate <- function(fname){
  fname <- paste("output/", fname, ".comp.txt", sep="");
  tab <- as.matrix(read.table(fname, header=T, sep=","));
  return(apply(tab, 2, mean));
}


tabreg <- rbind(getmean("Friedman1"),
                getmean("Friedman2"),
                getmean("Friedman3"),
                getmean("BostonHousing"),
                getmean("BodyFat"),
                getmean("Basketball"));

tabclass <- rbind(getrate("Circle2"),
                  getrate("Circle5"),
                  getrate("Circle20"),
                  getrate("Ionosphere"),
                  getrate("Sonar"),
                  getrate("WDBC"));
                  

library(xtable);
xtable(tabreg, digits=3);
xtable(tabclass, digits=4);