#!/usr/bin/env Rscript
library(ggplot2)
args = commandArgs(trailingOnly=TRUE)
groups = read.csv(args[1],header=T)
data = read.csv(args[2])
for (i in c(2:10)) {
  name = names(groups)[i]
  cols = c(2:15)
  group = data[groups[i] == 1,cols]
  freq = 100*colSums(group,na.rm=TRUE)/colSums(!is.na(group))
  algos = toupper(names(data)[cols])
  algos = gsub("HIGH",'HC',algos)
  algos = gsub(".CONFIDENCE",'',algos)
  algos = gsub("\\.",'-',algos)
  algos <- factor(algos,levels=rev(c(
  "MP2D-LAZAR-ALL",
  "MP2D-LAZAR-HC",
  "MP2D-RF",
  "MP2D-LR",
  "MP2D-LR2",
  "MP2D-NN",
  "MP2D-SVM",
  "CDK-LAZAR-ALL",
  "CDK-LAZAR-HC",
  "CDK-RF",
  "CDK-LR",
  "CDK-LR2",
  "CDK-NN",
  "CDK-SVM"
  )))
  plot = ggplot(data.frame(freq),aes(x=freq,y=algos)) + geom_bar(stat="identity") + xlab("% mutagenic") + ylab(element_blank()) + xlim(c(0,100))# + theme(axis.text.x = element_text(angle=90))
  ggsave(paste("figures/",name,".png",sep=""))
}