|
78
|
1 <tool id="QDNAseq-plot" name="QDNAseq-plot" version="0.0.1" force_history_refresh="True">
|
|
|
2
|
|
|
3 <requirements>
|
|
|
4
|
|
|
5 <!-- R 3.1.0 dependency will be used instead when available, now default R is used, see command -->
|
|
|
6 <!-- <requirement type="package" version="3.1.0">R</requirement> -->
|
|
|
7 <!-- <requirement type="package" version="1.2.2">qdnaseq</requirement> -->
|
|
|
8
|
|
|
9 </requirements>
|
|
|
10
|
|
|
11 <description>Plot QDNAseq copy-number/segments/calls profile</description>
|
|
|
12
|
|
|
13 <!-- change to /full/path/to/Rscript if required (eg /ccagc/lib/R/R-3.1.0/bin/Rscript) -->
|
|
|
14 <command interpreter="Rscript">
|
|
|
15 QDNAseq-plot.R
|
|
|
16 $cfg
|
|
|
17 </command>
|
|
|
18
|
|
|
19 <stdio>
|
|
|
20 <!-- Anything higher than 0 means the R script didnt finish (correctly) -->
|
|
|
21 <!-- Because different R packages deal with err/warn differently unable to waterproof this -->
|
|
|
22 <exit_code range="1:" level="fatal" description="R script finished too early, check log" />
|
|
|
23 </stdio>
|
|
|
24
|
|
|
25 <inputs>
|
|
|
26
|
|
|
27 <!-- ==================== -->
|
|
|
28 <!-- General inputs -->
|
|
|
29 <!-- ==================== -->
|
|
|
30
|
|
|
31 <!-- Job name: must contain non-whitespace chars -->
|
|
|
32 <param name="jobName" type="text" optional="false" label="Analysis/ouput name" help="Supply a name for the outputs to remind you what they contain" value="TEST">
|
|
|
33 <!-- <validator type="empty_field" /> -->
|
|
|
34 <validator type="regex" message="No whitespace characters allowed">^[^\s\\]*$</validator>
|
|
|
35 </param>
|
|
|
36
|
|
|
37 <!-- ==================== -->
|
|
|
38 <!-- Input RDS -->
|
|
|
39 <!-- ==================== -->
|
|
|
40 <param name="rdsFile" type="data" optional="False" format="rds" label="Input RDS file (QDNAseq R object)" help="RDS file should contain a QDNAseq/CGHcall object. If segment and/or call values are present these are added to the plot." />
|
|
|
41
|
|
|
42 <!-- ==================== -->
|
|
|
43 <!-- One pdf with all samples or a png of one chosen sample -->
|
|
|
44 <!-- ==================== -->
|
|
|
45 <!--
|
|
|
46 <conditional name="all_or_one">
|
|
|
47 <param name="show" type="select" label="Plot all samples in RDS object or choose one" help="">
|
|
|
48 <option value="all">Plot all samples</option>
|
|
|
49 <option value="one">Choose which one</option>
|
|
|
50 </param>
|
|
|
51 <when value="one">
|
|
|
52 <param name="sample_nr" type="integer" required="True" value="1" label="sample-nr (integer)" help="" />
|
|
|
53 </when>
|
|
|
54 <when value="all">
|
|
|
55 <param name="plotAll" type="hidden" value="" />
|
|
|
56 </when>
|
|
|
57
|
|
|
58 </conditional>
|
|
|
59 -->
|
|
|
60 <param name="sample_index" type="integer" required="True" value="1" label="Sample-index (integer)" help="The RDS input object can contain data from multiple individuals, this number tells the script which one to plot. Plotting multiple at the same time is not supported." />
|
|
|
61
|
|
|
62 <!-- ==================== -->
|
|
|
63 <!-- Optional advanced options -->
|
|
|
64 <!-- ==================== -->
|
|
|
65 <conditional name="advanced">
|
|
|
66 <param name="show" type="select" label="Extra options" help="Select yes to show extra options">
|
|
|
67 <option value="no">no</option>
|
|
|
68 <option value="yes">yes</option>
|
|
|
69 </param>
|
|
|
70 <when value="yes">
|
|
|
71 <!-- ==================== -->
|
|
|
72 <!-- Optional graphical/plotting options -->
|
|
|
73 <!-- ==================== -->
|
|
|
74 <param name="plot_width" size="3" type="integer" value="1440" label="Width of the png image produced" />
|
|
|
75 <param name="plot_height" size="3" type="integer" value="720" label="Height of the png image produced" />
|
|
|
76 </when>
|
|
|
77
|
|
|
78 <!-- need to set defaults because params are passed to R regardless of conditional opened/closed -->
|
|
|
79 <when value="no">
|
|
|
80 <param name="plot_width" type="hidden" value="1440" />
|
|
|
81 <param name="plot_height" type="hidden" value="720" />
|
|
|
82 </when>
|
|
|
83 </conditional>
|
|
|
84
|
|
|
85 </inputs>
|
|
|
86 <!-- ==================== -->
|
|
|
87 <!-- Config file to pass params to R script -->
|
|
|
88 <!-- ==================== -->
|
|
|
89 <configfiles>
|
|
|
90 <configfile name="cfg">
|
|
|
91 ## Desc: this file is sourced in QDNAseq-plot.R wrapper script
|
|
|
92 ## as means to pass all galaxy params to R
|
|
|
93
|
|
|
94 ## -----
|
|
|
95 ## required params
|
|
|
96 ## -----
|
|
|
97 "${jobName}" -> outputName
|
|
|
98 "${outputPng}" -> outputPngPath
|
|
|
99 "${rdsFile}" -> rdsFilePath
|
|
|
100 as.integer( "${sample_index}" ) -> SAMPLE_INDEX
|
|
|
101
|
|
|
102 ## -----
|
|
|
103 ## extra options
|
|
|
104 ## -----
|
|
|
105 as.integer( "${advanced.plot_width}" ) -> PLOT_WIDTH
|
|
|
106 as.integer( "${advanced.plot_height}" ) -> PLOT_HEIGHT
|
|
|
107
|
|
|
108 </configfile>
|
|
|
109 </configfiles>
|
|
|
110
|
|
|
111 <!-- ==================== -->
|
|
|
112 <!-- One image as output -->
|
|
|
113 <!-- ==================== -->
|
|
|
114 <outputs>
|
|
|
115 <data format="png" name="outputPng" label="QDNAseq: ${jobName} PNG plot" />
|
|
|
116 <!--
|
|
|
117 <data format="png" name="outputPng" label="QDNAseq: ${jobName} PNG plot">
|
|
|
118 <filter>(all_or_one['show'] == 'one')</filter>
|
|
|
119 </data>
|
|
|
120 -->
|
|
|
121 </outputs>
|
|
|
122
|
|
|
123 <help>
|
|
|
124 .. class:: infomark
|
|
|
125
|
|
|
126 **Introduction**
|
|
|
127
|
|
|
128 This tool is a wrapper for the plot function of the R Bioconductor package QDNAseq_
|
|
|
129
|
|
|
130 .. _QDNAseq: http://www.bioconductor.org/packages/release/bioc/html/QDNAseq.html
|
|
|
131
|
|
|
132 -----
|
|
|
133
|
|
|
134 .. class:: warningmark
|
|
|
135
|
|
|
136 As there is no R 3.1.0 galaxy-package yet (a requirement for QDNAseq), the **dependencies** need to be installed by hand and available to the user under which galaxy runs: R (3.1.0) and bioconductor package QDNAseq (>= 1.2.2). In case the default R is not 3.1.0, also the wrapper xml must be updated to include the correct path during installation of this tool.
|
|
|
137
|
|
|
138 -----
|
|
|
139
|
|
|
140 **Citation**
|
|
|
141
|
|
|
142 For the underlying QDNAseq R package please cite:
|
|
|
143 Scheinin I, Sie D, Bengtsson H, van de Wiel MA, Olshen AB, van Thuijl HF, van Essen HF, Eijk PP, Rustenburg F, Meijer GA, Reijneveld JC, Wesseling P, Pinkel D, Albertson DG and Ylstra B (2014). “DNA copy number analysis of fresh and formalin-fixed specimens by shallow whole-genome sequencing with identification and exclusion of problematic regions in the genome assembly.” Genome Research. doi:10.1101/gr.175141.114.
|
|
|
144
|
|
|
145 See also the bioconductor package_ documentation.
|
|
|
146
|
|
|
147 .. _package: http://www.bioconductor.org/packages/release/bioc/html/QDNAseq.html
|
|
|
148
|
|
|
149 .. image:: LGG150_copynumber.png
|
|
|
150 .. image:: LGG150_copynumberSegmented.png
|
|
|
151
|
|
|
152 </help>
|
|
|
153
|
|
|
154 </tool>
|
