# HG changeset patch # User mb2013 # Date 1384950016 18000 # Node ID fe749b5b462e112ed2077a04ac534446b4816857 # Parent 0cef5c0382c8bf804fcb18826ee9372b0f5e93d2 Uploaded diff -r 0cef5c0382c8 -r fe749b5b462e barplot.R --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/barplot.R Wed Nov 20 07:20:16 2013 -0500 @@ -0,0 +1,26 @@ +#Plottool makes a graph barplot of variance +#developer: M.Baak + +#commands extracting of commandline +args <- commandArgs(TRUE) + +#input files and options +input <- args[1] +title <- args[2] +#output file +output <- args[3] + +#reading of input files +read <- read.csv(file <- input,header = TRUE) + +stdev<- read[,1] #standard deviation +pov = stdev^2/sum(stdev^2) #variance + +png(output) #output in png format + + +#creating the barplot with the variance of pca +suppressMessages(barplot(pov, main = title, names.arg = c(1:NROW(pov)), ylim = c(0,1), col = heat.colors(3))) + + +graphics.off()