annotate ggplotscriptviolin.R @ 3:996c274d942d draft default tip

Uploaded
author moheydarian
date Wed, 22 Feb 2017 22:37:21 -0500
parents c7db13bd4c85
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
1
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
1 # Setup R error handling to go to stderr
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
2 options(show.error.messages=F, error=function(){cat(geterrmessage(),file=stderr());q("no",1,F)})
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
3
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
4 # We need to not crash galaxy with an UTF8 error on German LC settings.
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
5 loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8")
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
6
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
7
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
8 # Import library
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
9 library("getopt")
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
10 library("reshape2")
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
11 library("ggplot2")
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
12 options(stringAsfactors = FALSE, useFancyQuotes = FALSE)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
13 # Take in trailing command line arguments
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
14 args <- commandArgs(trailingOnly = TRUE)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
15
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
16
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
17 # get options, using the spec as defined by the enclosed list.
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
18 # we read the options from the default: commandArgs(TRUE).
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
19 option_specification = matrix(c(
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
20 'input', 'i', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
21 'title', 't',2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
22 'ggfill', 's', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
23 'ggcolor', 'z', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
24 'xlab', 'x', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
25 'ylab', 'y', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
26 'drawquartiles', 'a', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
27 'transform', 'w', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
28 'xaxismin', 'e', 2, 'integer',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
29 'xaxismax', 'f', 2, 'integer',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
30 'yaxismin', 'g', 2, 'integer',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
31 'yaxismax', 'h', 2, 'integer',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
32 'scaling', 'b', 2, 'character',
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
33 'output', 'o', 2, 'character'
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
34 ), byrow=TRUE, ncol=4);
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
35
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
36 # Parse options
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
37 options = getopt(option_specification);
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
38
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
39
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
40
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
41 # Print options to see what is going on
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
42 cat("\n input: ",options$input)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
43 cat("\n title: ",options$title)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
44 #cat("\n size: ",options$size)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
45 cat("\n xlab: ",options$xlab)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
46 cat("\n ylab: ",options$ylab)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
47 cat("\n output: ",options$output)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
48
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
49
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
50 if(options$scaling == "Automatic"){
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
51 gg_scaley = NULL
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
52 } else {
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
53 gg_scaley = ylim(options$yaxismin,options$yaxismax)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
54 cat("\n yaxismin: ",options$yaxismin)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
55 cat("\n yaxismax: ",options$yaxismax)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
56 }
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
57
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
58
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
59 integrated <- read.csv(options$input,sep='\t',header=TRUE)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
60 input <- melt(integrated)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
61
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
62
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
63
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
64 if(options$transform == "log2"){
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
65 input["value"] <- log2(input["value"])
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
66 }else if(options$transform == "log2plus1"){
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
67 input["value"] <- log2(input["value"]+1)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
68 }else if(options$transform == "log10"){
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
69 input["value"] <- log10(input["value"])
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
70 }else if(options$transform == "log10plus1"){
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
71 input["value"] <- log10(input["value"]+1)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
72 }else{
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
73 }
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
74
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
75 if(options$drawquartiles == "none"){
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
76 gg_quartile = NULL
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
77 } else {
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
78 gg_quartile = c(0.25, 0.5, 0.75)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
79 }
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
80
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
81
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
82
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
83 ggplot(input,aes(variable,value)) +geom_violin(scale = "area", colour = options$ggcolor, fill = options$ggfill, draw_quantiles = gg_quartile) +
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
84 gg_scaley+theme_bw()+xlab(options$xlab)+ylab(options$ylab)+ggtitle(options$title)
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
85
c7db13bd4c85 Uploaded
moheydarian
parents:
diff changeset
86 ggsave(file=options$output)