annotate IsoformSwitchAnalyzeR.R @ 3:bc68526d0e69 draft

planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
author iuc
date Wed, 24 May 2023 13:33:37 +0000
parents db30020fee5b
children 8575def9219d
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
1 # Load the IsoformSwitchAnalyzeR library
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
2 library(IsoformSwitchAnalyzeR,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
3 quietly = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
4 warn.conflicts = FALSE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
5 library(argparse, quietly = TRUE, warn.conflicts = FALSE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
6 library(dplyr, quietly = TRUE, warn.conflicts = FALSE)
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
7 library(ggplot2, quietly = TRUE, warn.conflicts = FALSE)
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
8
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
10 # setup R error handling to go to stderr
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
11 options(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
12 show.error.messages = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
13 error = function() {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
14 cat(geterrmessage(), file = stderr())
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
15 q("no", 1, FALSE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
16 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
17 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
18
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
19 # we need that to not crash galaxy with an UTF8 error on German LC settings.
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
20 loc <- Sys.setlocale("LC_MESSAGES", "en_US.UTF-8")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
21
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
22 ################################################################################
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
23 ### Input Processing
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
24 ################################################################################
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
25
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
26
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
27 # Collect arguments from command line
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
28 parser <- ArgumentParser(description = "IsoformSwitcheR R script")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
29
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
30 parser$add_argument("--modeSelector")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
31 parser$add_argument("--parentDir", required = FALSE, help = "Parent directory")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
32 parser$add_argument("--readLength",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
33 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
34 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
35 help = "Read length (required for stringtie)")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
36 parser$add_argument("--annotation", required = FALSE, help = "Annotation")
2
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
37 parser$add_argument("--stringtieAnnotation", required = FALSE, help = "Stringtie annotation")
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
38 parser$add_argument("--transcriptome", required = FALSE, help = "Transcriptome")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
39 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
40 "--fixStringTieAnnotationProblem",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
41 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
42 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
43 help = "Fix StringTie annotation problem"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
44 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
45 parser$add_argument("--countFiles", required = FALSE, help = "Count files")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
46 parser$add_argument("--toolSource", required = FALSE, help = "Tool source")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
47 parser$add_argument("--rObject", required = FALSE, help = "R object")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
48 parser$add_argument("--IFcutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
49 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
50 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
51 help = "IFcutoff")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
52 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
53 "--geneExpressionCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
54 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
55 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
56 help = "Gene expression cutoff"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
57 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
58 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
59 "--isoformExpressionCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
60 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
61 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
62 help = "Isoform expression cutoff"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
63 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
64 parser$add_argument("--alpha",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
65 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
66 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
67 help = "")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
68 parser$add_argument("--dIFcutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
69 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
70 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
71 help = "dIF cutoff")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
72 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
73 "--onlySigIsoforms",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
74 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
75 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
76 help = "Only significative isoforms"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
77 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
78 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
79 "--filterForConsequences",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
80 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
81 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
82 help = "Filter for consequences"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
83 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
84 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
85 "--removeSingleIsformGenes",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
86 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
87 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
88 help = "Remove single isoform genes"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
89 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
90 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
91 "--keepIsoformInAllConditions",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
92 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
93 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
94 help = "Keep isoform in all conditions"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
95 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
96 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
97 "--correctForConfoundingFactors",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
98 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
99 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
100 help = "Correct for confunding factors"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
101 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
102 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
103 "--overwriteIFvalues",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
104 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
105 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
106 help = "Overwrite IF values"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
107 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
108 parser$add_argument(
2
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
109 "--removeNonConvensionalChr",
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
110 required = FALSE,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
111 action = "store_true",
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
112 help = "Remove non-conventional chromosomes"
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
113 )
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
114 parser$add_argument(
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
115 "--reduceToSwitchingGenes",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
116 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
117 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
118 help = "Reduce to switching genes"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
119 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
120 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
121 "--reduceFurtherToGenesWithConsequencePotential",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
122 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
123 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
124 help = "Reduce further to genes with consequence potential"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
125 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
126 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
127 "--keepIsoformInAllConditions2",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
128 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
129 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
130 help = "Keep isoform in ll conditions"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
131 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
132 parser$add_argument("--minORFlength",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
133 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
134 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
135 help = "")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
136 parser$add_argument("--orfMethod", required = FALSE, help = "ORF methods")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
137 parser$add_argument("--PTCDistance",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
138 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
139 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
140 help = "")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
141 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
142 "--removeShortAAseq",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
143 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
144 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
145 help = "Remove short aminoacid sequences"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
146 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
147 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
148 "--removeLongAAseq",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
149 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
150 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
151 help = "Remove long aminoacid sequences"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
152 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
153 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
154 "--removeORFwithStop",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
155 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
156 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
157 help = "Remove ORF with stop codon"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
158 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
159 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
160 "--onlySwitchingGenes",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
161 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
162 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
163 help = "Only switching genes"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
164 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
165 parser$add_argument("--analysisMode", required = FALSE, help = "Analyze all isoforms with differential usage or single genes")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
166 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
167 "--genesToPlot",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
168 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
169 default = 10,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
170 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
171 help = "Number of genes to plot"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
172 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
173 parser$add_argument("--gene", required = FALSE, help = "Gene ID to analyze")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
174 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
175 "--sortByQvals",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
176 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
177 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
178 help = "Sort genes by Q-val values"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
179 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
180 parser$add_argument("--countGenes",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
181 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
182 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
183 help = "Count genes")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
184 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
185 "--asFractionTotal",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
186 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
187 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
188 help = "Plot gene expresson as fraction of total"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
189 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
190 parser$add_argument("--plotGenes",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
191 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
192 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
193 help = "Plot genes instead of isoforms")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
194 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
195 "--simplifyLocation",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
196 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
197 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
198 help = "Simplify localtion"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
199 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
200 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
201 "--removeEmptyConsequences",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
202 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
203 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
204 help = "Remove empty consequences"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
205 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
206 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
207 "--analysisOppositeConsequence",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
208 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
209 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
210 help = "Analysi opposite consequences"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
211 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
212 parser$add_argument("--pathToCPATresultFile",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
213 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
214 help = "Path to CPAT result file")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
215 parser$add_argument("--pathToCPC2resultFile",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
216 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
217 help = "Path to CPC2 result file")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
218 parser$add_argument("--pathToPFAMresultFile",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
219 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
220 help = "Path to PFAM result file")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
221 parser$add_argument("--pathToNetSurfP2resultFile",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
222 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
223 help = "Path to NetSurfP2 result file")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
224 parser$add_argument("--pathToSignalPresultFile",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
225 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
226 help = "Path to signalP result file")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
227 parser$add_argument("--pathToIUPred2AresultFile",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
228 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
229 help = "Path to IUPred2A result file")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
230 parser$add_argument("--codingCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
231 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
232 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
233 help = "Codding cutoff")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
234 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
235 "--removeNoncodingORFs",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
236 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
237 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
238 help = "Remove non-coding ORFs"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
239 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
240 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
241 "--minSignalPeptideProbability",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
242 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
243 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
244 help = "Minimul signal peptide probability"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
245 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
246 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
247 "--smoothingWindowSize",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
248 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
249 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
250 help = "Smoothing windows size"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
251 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
252 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
253 "--probabilityCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
254 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
255 type = "double",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
256 help = "Probability cutoff"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
257 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
258 parser$add_argument("--minIdrSize",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
259 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
260 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
261 help = "Min Idr size")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
262 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
263 "--annotateBindingSites",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
264 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
265 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
266 help = "Annotate binding sites"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
267 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
268 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
269 "--minIdrBindingSize",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
270 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
271 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
272 help = "Minimun Idr binding size"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
273 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
274 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
275 "--minIdrBindingOverlapFrac",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
276 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
277 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
278 help = ""
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
279 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
280 parser$add_argument("--ntCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
281 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
282 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
283 help = "Nucleotide cutoff")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
284 parser$add_argument("--ntFracCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
285 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
286 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
287 help = "Nucleotide fraction cutoff")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
288 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
289 "--ntJCsimCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
290 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
291 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
292 help = "Nucleotide Jaccard simmilarity cutoff"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
293 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
294 parser$add_argument("--AaCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
295 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
296 type = "integer",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
297 help = "Aminoacid cutoff")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
298 parser$add_argument("--AaFracCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
299 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
300 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
301 help = "Aminoacid fraction cutoff")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
302 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
303 "--AaJCsimCutoff",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
304 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
305 type = "numeric",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
306 help = "Aminoacid Jaccard similarity cutoff"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
307 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
308 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
309 "--removeNonConseqSwitches",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
310 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
311 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
312 help = "Remove switches without consequences"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
313 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
314 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
315 "--rescaleTranscripts",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
316 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
317 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
318 help = "Rescale transcripts"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
319 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
320 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
321 "--reverseMinus",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
322 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
323 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
324 help = "Reverse minus"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
325 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
326 parser$add_argument(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
327 "--addErrorbars",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
328 action = "store_true",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
329 required = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
330 help = "Add error bars"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
331 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
332
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
333
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
334 args <- parser$parse_args()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
335
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
336 # Data import
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
337 ###################
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
338
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
339 if (args$modeSelector == "data_import") {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
340
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
341 quantificationData <- importIsoformExpression(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
342 parentDir = args$parentDir,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
343 addIsofomIdAsColumn = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
344 readLength = args$readLength
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
345 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
346
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
347 ### Make design matrix
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
348 myDesign <- data.frame(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
349 sampleID = colnames(quantificationData$abundance)[-1],
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
350 condition = gsub(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
351 "[[:digit:]]+",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
352 "",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
353 colnames(quantificationData$abundance)[-1]
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
354 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
355 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
356
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
357 if (args$toolSource == "stringtie") {
2
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
358 if (!is.null(args$stringtieAnnotation)) {
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
359 SwitchList <- importRdata(
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
360 isoformCountMatrix = quantificationData$counts,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
361 isoformRepExpression = quantificationData$abundance,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
362 designMatrix = myDesign,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
363 removeNonConvensionalChr = args$removeNonConvensionalChr,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
364 isoformExonAnnoation = args$stringtieAnnotation,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
365 isoformNtFasta = args$transcriptome,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
366 addAnnotatedORFs = FALSE,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
367 showProgress = TRUE,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
368 fixStringTieAnnotationProblem = args$fixStringTieAnnotationProblem
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
369 )
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
370
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
371 SwitchList <- addORFfromGTF(
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
372 SwitchList,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
373 removeNonConvensionalChr = args$removeNonConvensionalChr,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
374 pathToGTF = args$annotation
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
375 )
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
376
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
377 } else {
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
378 SwitchList <- importRdata(
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
379 isoformCountMatrix = quantificationData$counts,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
380 isoformRepExpression = quantificationData$abundance,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
381 designMatrix = myDesign,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
382 removeNonConvensionalChr = args$removeNonConvensionalChr,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
383 isoformNtFasta = args$transcriptome,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
384 isoformExonAnnoation = args$annotation,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
385 showProgress = TRUE,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
386 fixStringTieAnnotationProblem = args$fixStringTieAnnotationProblem
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
387 )
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
388 }
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
389
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
390 } else {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
391 SwitchList <- importRdata(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
392 isoformCountMatrix = quantificationData$counts,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
393 isoformRepExpression = quantificationData$abundance,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
394 designMatrix = myDesign,
2
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
395 removeNonConvensionalChr = args$removeNonConvensionalChr,
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
396 isoformExonAnnoation = args$annotation,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
397 isoformNtFasta = args$transcriptome,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
398 showProgress = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
399 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
400 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
401
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
402 geneCountMatrix <- extractGeneExpression(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
403 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
404 extractCounts = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
405 addGeneNames = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
406 addIdsAsColumns = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
407 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
408
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
409 if (args$countFiles == "collection") {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
410
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
411 expressionDF <- data.frame(geneCountMatrix)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
412
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
413 myDesign$condition[length(myDesign$condition)]
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
414
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
415 dataframe_factor1 <- expressionDF %>% select(matches(myDesign$condition[1]))
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
416 dataframe_factor2 <- expressionDF %>% select(matches(myDesign$condition[length(myDesign$condition)]))
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
417
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
418
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
419 lf1 <- as.list(as.data.frame(dataframe_factor1))
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
420 sampleNames1 <- colnames(as.data.frame(dataframe_factor1))
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
421
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
422 lf2 <- as.list(as.data.frame(dataframe_factor2))
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
423 sampleNames2 <- colnames(as.data.frame(dataframe_factor2))
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
424
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
425 geneNames <- row.names(as.data.frame(expressionDF))
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
426
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
427
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
428 for (index in seq_along(lf1)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
429 tabular_expression <- data.frame(geneNames, lf1[index])
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
430 colnames(tabular_expression) <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
431 c("Geneid", sampleNames1[index])
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
432 filename <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
433 paste(sampleNames1[index], "dataset.tabular", sep = "_")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
434 output_path <- paste("./count_files/factor1/", filename, sep = "")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
435 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
436 tabular_expression,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
437 output_path,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
438 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
439 row.names = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
440 quote = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
441 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
442 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
443 for (index in seq_along(lf2)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
444 tabular_expression <- data.frame(geneNames, lf2[index])
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
445 colnames(tabular_expression) <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
446 c("Geneid", sampleNames2[index])
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
447 filename <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
448 paste(sampleNames2[index], "dataset.tabular", sep = "_")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
449 output_path <- paste("./count_files/factor2/", filename, sep = "")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
450 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
451 tabular_expression,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
452 output_path,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
453 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
454 row.names = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
455 quote = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
456 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
457 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
458 } else if (args$countFiles == "matrix") {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
459 expressionDF <- data.frame(geneCountMatrix)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
460 geneNames <- row.names(expressionDF)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
461
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
462 expressionDF <- cbind(geneNames, expressionDF)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
463 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
464 as.data.frame(expressionDF),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
465 "./count_files/matrix.tabular",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
466 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
467 row.names = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
468 quote = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
469 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
470 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
471 as.data.frame(myDesign),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
472 "./count_files/samples.tabular",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
473 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
474 row.names = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
475 quote = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
476 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
477 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
478
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
479 save(SwitchList, file = "SwitchList.Rda")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
480
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
481 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
482
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
483 if (args$modeSelector == "first_step") {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
484
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
485 # First part of the analysis
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
486 #############################
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
487
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
488 load(file = args$rObject)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
489
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
490 ### Filter
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
491 SwitchList <- preFilter(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
492 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
493 IFcutoff = args$IFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
494 geneExpressionCutoff = args$geneExpressionCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
495 isoformExpressionCutoff = args$isoformExpressionCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
496 removeSingleIsoformGenes = args$removeSingleIsformGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
497 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
498 keepIsoformInAllConditions = args$keepIsoformInAllConditions,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
499 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
500 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
501 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
502
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
503 ### Test for isoform switches
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
504 SwitchList <- isoformSwitchTestDEXSeq(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
505 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
506 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
507 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
508 correctForConfoundingFactors = args$correctForConfoundingFactors,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
509 overwriteIFvalues = args$overwriteIFvalues,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
510 reduceToSwitchingGenes = args$reduceToSwitchingGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
511 reduceFurtherToGenesWithConsequencePotential = args$reduceFurtherToGenesWithConsequencePotential,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
512 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
513 keepIsoformInAllConditions = args$keepIsoformInAllConditions2,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
514 showProgress = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
515 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
516
3
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
517 # Analyze missing annotated isoforms by default
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
518 SwitchList <- analyzeNovelIsoformORF(
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
519 SwitchList,
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
520 analysisAllIsoformsWithoutORF = TRUE,
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
521 minORFlength = args$minORFlength,
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
522 orfMethod = args$orfMethod,
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
523 PTCDistance = args$PTCDistance,
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
524 startCodons = "ATG",
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
525 stopCodons = c("TAA", "TAG", "TGA"),
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
526 showProgress = TRUE,
bc68526d0e69 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 44cb89c86b7a48993b6987a3a40129782fd7a093
iuc
parents: 2
diff changeset
527 )
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
528
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
529 ### Extract Sequences
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
530 SwitchList <- extractSequence(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
531 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
532 onlySwitchingGenes = args$onlySwitchingGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
533 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
534 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
535 extractNTseq = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
536 extractAAseq = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
537 removeShortAAseq = args$removeShortAAseq,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
538 removeLongAAseq = args$removeLongAAseq,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
539 removeORFwithStop = args$removeORFwithStop,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
540 addToSwitchAnalyzeRlist = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
541 writeToFile = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
542 pathToOutput = getwd(),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
543 outputPrefix = "isoformSwitchAnalyzeR_isoform",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
544 forceReExtraction = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
545 quiet = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
546 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
547
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
548 ### Summary
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
549 switchSummary <- extractSwitchSummary(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
550 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
551 filterForConsequences = args$filterForConsequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
552 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
553 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
554 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
555 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
556
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
557 save(SwitchList, file = "SwitchList.Rda")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
558 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
559 switchSummary,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
560 file = "switchSummary.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
561 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
562 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
563 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
564 row.names = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
565 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
566
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
567 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
568
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
569 if (args$modeSelector == "second_step") {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
570
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
571 # Second part of the analysis
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
572 #############################
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
573
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
574 load(file = args$rObject)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
575
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
576 ### Add annotation
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
577 if (!is.null(args$pathToCPATresultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
578 SwitchList <- analyzeCPAT(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
579 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
580 pathToCPATresultFile = args$pathToCPATresultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
581 codingCutoff = args$codingCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
582 removeNoncodinORFs = args$removeNoncodingORFs
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
583 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
584 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
585
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
586 if (!is.null(args$pathToCPC2resultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
587 SwitchList <- analyzeCPC2(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
588 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
589 pathToCPC2resultFile = args$pathToCPC2resultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
590 removeNoncodinORFs = args$removeNoncodingORFs
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
591 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
592 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
593
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
594 if (!is.null(args$pathToPFAMresultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
595 pfamFiles <- list.files(path = args$pathToPFAMresultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
596 full.names = TRUE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
597
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
598 SwitchList <- analyzePFAM(SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
599 pathToPFAMresultFile = pfamFiles,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
600 showProgress = FALSE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
601 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
602
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
603 if (!is.null(args$pathToNetSurfP2resultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
604 netsurfFiles <- list.files(path = args$pathToNetSurfP2resultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
605 full.names = TRUE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
606
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
607 SwitchList <- analyzeNetSurfP2(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
608 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
609 pathToNetSurfP2resultFile = netsurfFiles,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
610 smoothingWindowSize = args$smoothingWindowSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
611 probabilityCutoff = args$probabilityCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
612 minIdrSize = args$minIdrSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
613 showProgress = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
614 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
615 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
616
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
617 if (!is.null(args$pathToIUPred2AresultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
618 SwitchList <- analyzeIUPred2A(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
619 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
620 pathToIUPred2AresultFile = args$pathToIUPred2AresultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
621 smoothingWindowSize = args$smoothingWindowSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
622 probabilityCutoff = args$probabilityCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
623 minIdrSize = args$minIdrSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
624 annotateBindingSites = args$annotateBindingSites,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
625 minIdrBindingSize = args$minIdrBindingSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
626 minIdrBindingOverlapFrac = args$minIdrBindingOverlapFrac,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
627 showProgress = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
628 quiet = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
629 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
630 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
631
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
632 if (!is.null(args$pathToSignalPresultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
633 signalpFiles <- list.files(path = args$pathToSignalPresultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
634 full.names = TRUE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
635
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
636 SwitchList <- analyzeSignalP(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
637 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
638 pathToSignalPresultFile = signalpFiles,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
639 minSignalPeptideProbability = args$minSignalPeptideProbability
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
640 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
641 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
642
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
643 SwitchList <- analyzeAlternativeSplicing(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
644 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
645 onlySwitchingGenes = args$onlySwitchingGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
646 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
647 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
648 showProgress = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
649 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
650
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
651 SwitchList <- analyzeIntronRetention(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
652 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
653 onlySwitchingGenes = args$onlySwitchingGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
654 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
655 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
656 showProgress = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
657 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
658
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
659 consequences <- c(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
660 "intron_retention",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
661 "NMD_status",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
662 "isoform_seq_similarity",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
663 "ORF_genomic",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
664 "tss",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
665 "tts"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
666 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
667
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
668 if (!is.null(args$pathToCPATresultFile) ||
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
669 !is.null(args$pathToCPC2resultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
670 updatedConsequences <- c(consequences, "coding_potential")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
671 consequences <- updatedConsequences
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
672 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
673
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
674 if (!is.null(args$pathToPFAMresultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
675 updatedConsequences <- c(consequences, "domains_identified")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
676 consequences <- updatedConsequences
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
677 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
678
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
679 if (!is.null(args$pathToSignalPresultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
680 updatedConsequences <- c(consequences, "signal_peptide_identified")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
681 consequences <- updatedConsequences
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
682 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
683
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
684 if (!is.null(args$pathToNetSurfP2resultFile) ||
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
685 !is.null(args$pathToIUPred2AresultFile)) {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
686 updatedConsequences <- c(consequences, "IDR_identified", "IDR_type")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
687 consequences <- updatedConsequences
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
688 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
689
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
690 SwitchList <- analyzeSwitchConsequences(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
691 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
692 consequencesToAnalyze = consequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
693 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
694 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
695 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
696 ntCutoff = args$ntCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
697 ntJCsimCutoff = args$ntJCsimCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
698 AaCutoff = args$AaCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
699 AaFracCutoff = args$AaFracCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
700 AaJCsimCutoff = args$AaJCsimCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
701 removeNonConseqSwitches = args$removeNonConseqSwitches,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
702 showProgress = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
703 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
704
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
705
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
706 ### Visual analysis
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
707 # Top genes
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
708
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
709 if (args$analysisMode == "single") {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
710
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
711 outputFile <- file.path(getwd(), "single_gene.pdf")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
712
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
713 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
714 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
715 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
716 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
717 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
718 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
719
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
720 switchPlot(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
721 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
722 gene = args$gene,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
723 condition1 = myDesign$condition[1],
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
724 condition2 = myDesign$condition[length(myDesign$condition)],
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
725 IFcutoff = args$IFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
726 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
727 rescaleTranscripts = args$rescaleTranscripts,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
728 reverseMinus = args$reverseMinus,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
729 addErrorbars = args$addErrorbars,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
730 logYaxis = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
731 localTheme = theme_bw(base_size = 8)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
732 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
733 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
734
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
735 } else {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
736 mostSwitchingGene <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
737 extractTopSwitches(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
738 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
739 n = Inf,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
740 filterForConsequences = args$filterForConsequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
741 extractGenes = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
742 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
743 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
744 inEachComparison = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
745 sortByQvals = args$sortByQvals
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
746 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
747
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
748 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
749 mostSwitchingGene,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
750 file = "mostSwitchingGene.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
751 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
752 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
753 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
754 row.names = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
755 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
756
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
757
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
758 switchPlotTopSwitches(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
759 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
760 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
761 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
762 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
763 n = args$genesToPlot,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
764 sortByQvals = args$sortByQvals,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
765 pathToOutput = getwd(),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
766 fileType = "pdf"
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
767 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
768
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
769 outputFile <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
770 file.path(getwd(), "extractConsequencesSummary.pdf")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
771
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
772 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
773 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
774 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
775 height = 6,
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
776 width = 12
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
777 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
778
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
779 consequenceSummary <- extractConsequenceSummary(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
780 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
781 consequencesToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
782 includeCombined = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
783 asFractionTotal = args$asFractionTotal,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
784 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
785 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
786 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
787 plotGenes = args$plotGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
788 simplifyLocation = args$simplifyLocation,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
789 removeEmptyConsequences = args$removeEmptyConsequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
790 returnResult = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
791 localTheme = theme_bw()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
792 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
793 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
794
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
795 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
796 consequenceSummary,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
797 file = "consequencesSummary.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
798 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
799 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
800 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
801 row.names = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
802 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
803
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
804
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
805 outputFile <- file.path(getwd(), "consequencesEnrichment.pdf")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
806 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
807 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
808 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
809 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
810 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
811 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
812 consequenceEnrichment <- extractConsequenceEnrichment(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
813 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
814 consequencesToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
815 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
816 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
817 countGenes = args$countGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
818 analysisOppositeConsequence = args$analysisOppositeConsequence,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
819 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
820 localTheme = theme_bw(base_size = 12),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
821 minEventsForPlotting = 10,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
822 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
823 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
824 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
825
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
826 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
827 consequenceEnrichment,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
828 file = "consequencesEnrichment.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
829 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
830 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
831 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
832 row.names = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
833 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
834
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
835
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
836 outputFile <- file.path(getwd(), "splicingEnrichment.pdf")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
837 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
838 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
839 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
840 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
841 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
842 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
843 splicingEnrichment <- extractSplicingEnrichment(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
844 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
845 splicingToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
846 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
847 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
848 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
849 countGenes = args$countGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
850 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
851 minEventsForPlotting = 10,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
852 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
853 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
854 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
855
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
856 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
857 splicingEnrichment,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
858 file = "splicingEnrichment.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
859 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
860 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
861 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
862 row.names = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
863 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
864
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
865
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
866 outputFile <- file.path(getwd(), "splicingSummary.pdf")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
867 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
868 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
869 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
870 height = 6,
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
871 width = 12
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
872 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
873 splicingSummary <- extractSplicingSummary(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
874 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
875 splicingToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
876 asFractionTotal = args$asFractionTotal,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
877 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
878 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
879 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
880 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
881 plotGenes = args$plotGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
882 localTheme = theme_bw(),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
883 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
884 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
885 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
886
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
887 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
888 splicingSummary,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
889 file = "splicingSummary.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
890 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
891 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
892 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
893 row.names = FALSE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
894 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
895
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
896
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
897 ### Volcano like plot:
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
898 outputFile <- file.path(getwd(), "volcanoPlot.pdf")
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
899
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
900 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
901 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
902 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
903 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
904 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
905 )
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
906
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
907 p <- ggplot(data = SwitchList$isoformFeatures, aes(x = dIF, y = -log10(isoform_switch_q_value))) +
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
908 geom_point(
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
909 aes(color = abs(dIF) > 0.1 & isoform_switch_q_value < 0.05), # default cutoff
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
910 size = 1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
911 ) +
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
912 geom_hline(yintercept = -log10(0.05), linetype = "dashed") + # default cutoff
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
913 geom_vline(xintercept = c(-0.1, 0.1), linetype = "dashed") + # default cutoff
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
914 facet_wrap(~ condition_2) +
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
915 #facet_grid(condition_1 ~ condition_2) + # alternative to facet_wrap if you have overlapping conditions
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
916 scale_color_manual("Signficant\nIsoform Switch", values = c("black", "red")) +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
917 labs(x = "dIF", y = "-Log10 ( Isoform Switch Q Value )") +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
918 theme_bw()
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
919 print(p)
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
920 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
921
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
922 ### Switch vs Gene changes:
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
923 outputFile <- file.path(getwd(), "switchGene.pdf")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
924 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
925 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
926 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
927 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
928 )
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
929 p <- ggplot(data = SwitchList$isoformFeatures,
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
930 aes(x = gene_log2_fold_change, y = dIF)) +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
931 geom_point(aes(color = abs(dIF) > 0.1 &
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
932 isoform_switch_q_value < 0.05),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
933 size = 1) +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
934 facet_wrap(~ condition_2) +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
935 geom_hline(yintercept = 0, linetype = "dashed") +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
936 geom_vline(xintercept = 0, linetype = "dashed") +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
937 scale_color_manual("Signficant\nIsoform Switch", values = c("black", "red")) +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
938 labs(x = "Gene log2 fold change", y = "dIF") +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
939 theme_bw()
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
940 print(p)
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
941 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
942
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
943 outputFile <- file.path(getwd(), "splicingGenomewide.pdf")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
944 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
945 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
946 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
947 height = 6,
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
948 width = 14
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
949 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
950 splicingGenomeWide <- extractSplicingGenomeWide(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
951 SwitchList,
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
952 featureToExtract = "isoformUsage",
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
953 splicingToAnalyze = "all",
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
954 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
955 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
956 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
957 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
958 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
959 save(SwitchList, file = "SwitchList.Rda")
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
960
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
961 }