Mercurial > repos > mb2013 > nepenthes_3dpca
comparison variance_plot.R @ 19:f3d9bb18110a draft default tip
Uploaded
| author | mb2013 |
|---|---|
| date | Tue, 20 May 2014 03:30:07 -0400 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| 18:a207e838f232 | 19:f3d9bb18110a |
|---|---|
| 1 #Plottool makes a graph barplot of the variance | |
| 2 #MB | |
| 3 | |
| 4 #commands extracting of commandline | |
| 5 args <- commandArgs(TRUE) | |
| 6 | |
| 7 #input files and options | |
| 8 input <- args[1] | |
| 9 title <- args[2] | |
| 10 #output file | |
| 11 output <- args[3] | |
| 12 | |
| 13 #reading of input files | |
| 14 read <- read.csv(file <- input,header = TRUE) | |
| 15 | |
| 16 stdev<- read[,1] #standard deviation | |
| 17 pov = stdev^2/sum(stdev^2) #variance | |
| 18 | |
| 19 png(output) #output in png format | |
| 20 | |
| 21 | |
| 22 #creating the barplot with the variance of pca | |
| 23 suppressMessages(barplot(pov, main = title, names.arg = c(1:NROW(pov)), ylim = c(0,1), col = heat.colors(3))) | |
| 24 | |
| 25 | |
| 26 graphics.off() |