argsall <- commandArgs() input <- argsall[6] output <- argsall[7] tabl <- read.table(file=input, header=TRUE, sep="\t") fac <- levels(tabl[,2]) #fac #length(fac) ini <- rep(0,length(fac)) #for (q in 1:length(fac)) { # ini[q] <- 0 #} #ini for (i in 1:length(tabl[,2])) { for( j in 1:length(fac) ){ if(tabl[,2][i] == fac[j]){ ini[j] <- ini[j] + tabl[,3][i] } } } figureformat <- "fi.pdf" pdf(figureformat) par(mar=c(5,12,4,2)) barplot(ini, names.arg=fac,col="red", width=0.2, axis.lty=0, las=2, xlab="Frequence of ARGs", main="Bar Plot of ARGs Class", cex.axis=1, cex=0.7,offset=0, font=2, adj=0, horiz=T,border=NA, density=NA) dev.off() file.copy(figureformat, output, overwrite=TRUE) quit()