Mercurial > repos > iuc > isoformswitchanalyzer

annotate IsoformSwitchAnalyzeR.R @ 2:db30020fee5b draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
author iuc
date Wed, 24 May 2023 06:06:36 +0000
parents 5ae218cee629
children bc68526d0e69
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
2
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
517 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
518 SwitchList <- analyzeNovelIsoformORF(
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
519 SwitchList,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
520 analysisAllIsoformsWithoutORF = TRUE,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
521 minORFlength = args$minORFlength,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
522 orfMethod = args$orfMethod,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
523 PTCDistance = args$PTCDistance,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
524 startCodons = "ATG",
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
525 stopCodons = c("TAA", "TAG", "TGA"),
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
526 showProgress = TRUE,
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
527 )
db30020fee5b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit d00628bfc2ac36914885e0941a2b4453ca0bff9a
iuc
parents: 1
diff changeset
528 }
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
529
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
530 ### Extract Sequences
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
531 SwitchList <- extractSequence(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
532 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
533 onlySwitchingGenes = args$onlySwitchingGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
534 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
535 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
536 extractNTseq = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
537 extractAAseq = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
538 removeShortAAseq = args$removeShortAAseq,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
539 removeLongAAseq = args$removeLongAAseq,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
540 removeORFwithStop = args$removeORFwithStop,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
541 addToSwitchAnalyzeRlist = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
542 writeToFile = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
543 pathToOutput = getwd(),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
544 outputPrefix = "isoformSwitchAnalyzeR_isoform",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
545 forceReExtraction = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
546 quiet = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
549 ### Summary
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
550 switchSummary <- extractSwitchSummary(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
551 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
552 filterForConsequences = args$filterForConsequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
553 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
554 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
555 onlySigIsoforms = args$onlySigIsoforms,
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
558 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
559 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
560 switchSummary,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
561 file = "switchSummary.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
562 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
563 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
564 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
565 row.names = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
570 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
571
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
572 # 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
573 #############################
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
574
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
575 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
576
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
577 ### Add annotation
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
578 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
579 SwitchList <- analyzeCPAT(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
580 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
581 pathToCPATresultFile = args$pathToCPATresultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
582 codingCutoff = args$codingCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
583 removeNoncodinORFs = args$removeNoncodingORFs
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
587 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
588 SwitchList <- analyzeCPC2(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
589 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
590 pathToCPC2resultFile = args$pathToCPC2resultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
591 removeNoncodinORFs = args$removeNoncodingORFs
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
595 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
596 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
597 full.names = TRUE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
598
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
599 SwitchList <- analyzePFAM(SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
600 pathToPFAMresultFile = pfamFiles,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
601 showProgress = FALSE)
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
604 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
605 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
606 full.names = TRUE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
607
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
608 SwitchList <- analyzeNetSurfP2(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
609 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
610 pathToNetSurfP2resultFile = netsurfFiles,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
611 smoothingWindowSize = args$smoothingWindowSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
612 probabilityCutoff = args$probabilityCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
613 minIdrSize = args$minIdrSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
614 showProgress = TRUE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
618 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
619 SwitchList <- analyzeIUPred2A(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
620 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
621 pathToIUPred2AresultFile = args$pathToIUPred2AresultFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
622 smoothingWindowSize = args$smoothingWindowSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
623 probabilityCutoff = args$probabilityCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
624 minIdrSize = args$minIdrSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
625 annotateBindingSites = args$annotateBindingSites,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
626 minIdrBindingSize = args$minIdrBindingSize,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
627 minIdrBindingOverlapFrac = args$minIdrBindingOverlapFrac,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
628 showProgress = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
629 quiet = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
633 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
634 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
635 full.names = TRUE)
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
636
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
637 SwitchList <- analyzeSignalP(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
638 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
639 pathToSignalPresultFile = signalpFiles,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
640 minSignalPeptideProbability = args$minSignalPeptideProbability
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
644 SwitchList <- analyzeAlternativeSplicing(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
645 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
646 onlySwitchingGenes = args$onlySwitchingGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
647 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
648 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
649 showProgress = TRUE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
652 SwitchList <- analyzeIntronRetention(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
653 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
654 onlySwitchingGenes = args$onlySwitchingGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
655 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
656 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
657 showProgress = TRUE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
660 consequences <- c(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
661 "intron_retention",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
662 "NMD_status",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
663 "isoform_seq_similarity",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
664 "ORF_genomic",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
665 "tss",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
666 "tts"
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
669 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
670 !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
671 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
672 consequences <- updatedConsequences
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
675 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
676 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
677 consequences <- updatedConsequences
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
680 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
681 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
682 consequences <- updatedConsequences
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
685 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
686 !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
687 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
688 consequences <- updatedConsequences
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
691 SwitchList <- analyzeSwitchConsequences(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
692 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
693 consequencesToAnalyze = consequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
694 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
695 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
696 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
697 ntCutoff = args$ntCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
698 ntJCsimCutoff = args$ntJCsimCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
699 AaCutoff = args$AaCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
700 AaFracCutoff = args$AaFracCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
701 AaJCsimCutoff = args$AaJCsimCutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
702 removeNonConseqSwitches = args$removeNonConseqSwitches,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
703 showProgress = TRUE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
707 ### Visual analysis
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
708 # Top genes
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
709
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
710 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
711
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
712 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
713
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
714 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
715 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
716 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
717 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
718 width = 9
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
721 switchPlot(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
722 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
723 gene = args$gene,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
724 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
725 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
726 IFcutoff = args$IFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
727 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
728 rescaleTranscripts = args$rescaleTranscripts,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
729 reverseMinus = args$reverseMinus,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
730 addErrorbars = args$addErrorbars,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
731 logYaxis = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
732 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
733 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
734 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
735
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
736 } else {
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
737 mostSwitchingGene <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
738 extractTopSwitches(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
739 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
740 n = Inf,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
741 filterForConsequences = args$filterForConsequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
742 extractGenes = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
743 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
744 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
745 inEachComparison = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
746 sortByQvals = args$sortByQvals
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
749 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
750 mostSwitchingGene,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
751 file = "mostSwitchingGene.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
752 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
753 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
754 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
755 row.names = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
759 switchPlotTopSwitches(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
760 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
761 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
762 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
763 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
764 n = args$genesToPlot,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
765 sortByQvals = args$sortByQvals,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
766 pathToOutput = getwd(),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
767 fileType = "pdf"
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
770 outputFile <-
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
771 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
772
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
773 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
774 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
775 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
776 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
777 width = 12
0
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
780 consequenceSummary <- extractConsequenceSummary(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
781 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
782 consequencesToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
783 includeCombined = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
784 asFractionTotal = args$asFractionTotal,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
785 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
786 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
787 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
788 plotGenes = args$plotGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
789 simplifyLocation = args$simplifyLocation,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
790 removeEmptyConsequences = args$removeEmptyConsequences,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
791 returnResult = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
792 localTheme = theme_bw()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
793 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
794 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
795
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
796 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
797 consequenceSummary,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
798 file = "consequencesSummary.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
799 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
800 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
801 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
802 row.names = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
806 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
807 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
808 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
809 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
810 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
811 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
812 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
813 consequenceEnrichment <- extractConsequenceEnrichment(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
814 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
815 consequencesToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
816 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
817 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
818 countGenes = args$countGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
819 analysisOppositeConsequence = args$analysisOppositeConsequence,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
820 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
821 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
822 minEventsForPlotting = 10,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
823 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
824 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
825 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
826
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
827 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
828 consequenceEnrichment,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
829 file = "consequencesEnrichment.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
830 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
831 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
832 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
833 row.names = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
837 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
838 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
839 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
840 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
841 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
842 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
843 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
844 splicingEnrichment <- extractSplicingEnrichment(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
845 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
846 splicingToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
847 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
848 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
849 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
850 countGenes = args$countGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
851 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
852 minEventsForPlotting = 10,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
853 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
854 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
855 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
856
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
857 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
858 splicingEnrichment,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
859 file = "splicingEnrichment.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
860 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
861 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
862 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
863 row.names = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
867 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
868 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
869 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
870 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
871 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
872 width = 12
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
873 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
874 splicingSummary <- extractSplicingSummary(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
875 SwitchList,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
876 splicingToAnalyze = "all",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
877 asFractionTotal = args$asFractionTotal,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
878 alpha = args$alpha,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
879 dIFcutoff = args$dIFcutoff,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
880 onlySigIsoforms = args$onlySigIsoforms,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
881 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
882 plotGenes = args$plotGenes,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
883 localTheme = theme_bw(),
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
884 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
885 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
886 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
887
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
888 write.table(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
889 splicingSummary,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
890 file = "splicingSummary.tsv",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
891 quote = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
892 sep = "\t",
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
893 col.names = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
894 row.names = FALSE
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
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
898 ### Volcano like plot:
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
899 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
900
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
901 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
902 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
903 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
904 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
905 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
906 )
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
907
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
908 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
909 geom_point(
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
910 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
911 size = 1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
912 ) +
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
913 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
914 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
915 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
916 #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
917 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
918 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
919 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
920 print(p)
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
921 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
922
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
923 ### 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
924 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
925 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
926 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
927 height = 6,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
928 width = 9
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
929 )
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
930 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
931 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
932 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
933 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
934 size = 1) +
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
935 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
936 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
937 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
938 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
939 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
940 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
941 print(p)
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
942 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
943
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
944 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
945 pdf(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
946 file = outputFile,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
947 onefile = FALSE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
948 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
949 width = 14
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
950 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
951 splicingGenomeWide <- extractSplicingGenomeWide(
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
952 SwitchList,
1
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
953 featureToExtract = "isoformUsage",
5ae218cee629 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 7d278967f9c8d2c6ce0a8f83be2c444822746bbf
iuc
parents: 0
diff changeset
954 splicingToAnalyze = "all",
0
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
955 plot = TRUE,
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
956 returnResult = TRUE
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
957 )
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
958 dev.off()
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
959 }
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
960 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
961
bb611fa3bc3b planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/isoformswitchanalyzer commit 2c61e4c6151000201dd9a8323722a380bc235380
iuc
parents:
diff changeset
962 }