File: plot_ExN50_statistic.Rscript

package info (click to toggle)
trinityrnaseq 2.2.0%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: stretch
  • size: 212,452 kB
  • ctags: 5,067
  • sloc: perl: 45,552; cpp: 19,678; java: 11,865; sh: 1,485; makefile: 613; ansic: 427; python: 313; xml: 83
file content (14 lines) | stat: -rwxr-xr-x 399 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
#!/usr/bin/env Rscript

args<-commandArgs(TRUE)

dat_filename = args[1]
pdf_filename = paste(dat_filename, ".plot.pdf", sep='')
pdf(pdf_filename)
data = read.table(dat_filename, com='', header=T, row.names=1)
plot(data$E.N50, xlab="Pct. Ex.", ylab="N50 contig length", col='blue', t='b')

write(cat("ExN50 data plotted as:", pdf_filename), stderr())

quit(save = "no", status = 0, runLast = FALSE)