0
|
1 <tool id="phyloseq_ordinate" name="Phyloseq Ordination Plot" version="1.22.3.2" hidden="false">
|
|
2 <description>ordination plotting</description>
|
|
3 <requirements>
|
|
4 <requirement type="package" version="1.22.3">bioconductor-phyloseq</requirement>
|
|
5 <requirement type="package" version="1.20.0">r-getopt</requirement>
|
|
6 <requirement type="package" version="1.0.11">r-doparallel</requirement>
|
|
7 <requirement type="package" version="9.18">ghostscript</requirement>
|
|
8 </requirements>
|
|
9 <version_command><![CDATA[
|
|
10 echo $(R --version | grep version | grep -v GNU)", phyloseq version" $(R --vanilla --slave -e "library(phyloseq); cat(sessionInfo()\$otherPkgs\$phyloseq\$Version)" 2> /dev/null | grep -v -i "WARNING: ")
|
|
11 ]]></version_command>
|
|
12 <command detect_errors="exit_code"><![CDATA[
|
|
13 Rscript '${__tool_directory__}/phyloseq_ordinate_plot.R'
|
|
14 #if str($file_source.file_source_selector) == "set_biom":
|
|
15 --biom='$file_source.input'
|
|
16 --subset='$file_source.subsetColumn'
|
|
17 #else:
|
|
18 --otu_table='$OTU_TABLE'
|
|
19 --tax_table='$TAX_TABLE'
|
|
20 --meta_table='$META_TABLE'
|
|
21 --category='$file_source.category'
|
|
22 #end if
|
1
|
23 #if str($plot_selection.plot_type) == "2":
|
|
24 --plottype='$plot_selection.plot_type'
|
|
25 --kingdom='null'
|
|
26 #else:
|
|
27 --plottype='$plot_selection.plot_type'
|
|
28 --kingdom='$plot_selection.kingdom_field'
|
|
29 #end if
|
0
|
30 --method='$ord_method'
|
|
31 --distance='$distance'
|
1
|
32 --log='$log'
|
0
|
33 --outdir="$htmlfile.files_path"
|
|
34 --htmlfile='$htmlfile'
|
|
35 ]]>
|
|
36 </command>
|
|
37
|
|
38 <inputs>
|
|
39 <conditional name="file_source">
|
|
40 <param name="file_source_selector" type="select" label="Choose an input file type">
|
|
41 <option value="set_biom" selected="True">BIOM File</option>
|
|
42 <option value="set_table">TABULAR File</option>
|
|
43 </param>
|
|
44 <when value="set_biom">
|
|
45 <param format="biom1" name="input" type="data" label="Input File"/>
|
1
|
46 <param name="subsetColumn" type="select" label="Classify plot by">
|
0
|
47 <options>
|
|
48 <filter type="data_meta" ref="input" key="table_column_metadata_headers" />
|
|
49 </options>
|
|
50 </param>
|
|
51 </when>
|
|
52 <when value="set_table">
|
|
53 <param format="tabular" name="OTU_TABLE" type="data" label="OTU table"/>
|
|
54 <param format="tabular" name="TAX_TABLE" type="data" label="Taxonomy table"/>
|
|
55 <param format="tabular" name="META_TABLE" type="data" label="Metadata table"/>
|
1
|
56 <param name="category" type="data_column" data_ref="META_TABLE" use_header_names="True" label="Classify plot by"/>
|
0
|
57 </when>
|
|
58 </conditional>
|
|
59 <param name="ord_method" type="select" display="radio" label="Select Ordination method">
|
|
60 <option value="NMDS" selected="true">NMDS(Non-metric Multidimensional Scalling)</option>
|
|
61 <option value="DCA">DCA (Detrended Correspondence Analysis)</option>
|
|
62 <option value="CCA">CCA (Constrained Correspondence Analysis)</option>
|
|
63 <option value="RDA">RDA (Redundancy Analysis)</option>
|
|
64 <option value="DPCoA">DPCoA (Double Principle Coordinate Analysis)</option>
|
|
65 <option value="MDS">MDS (Multidimensional Scaling)</option>
|
|
66 <option value="PCoA">PCoA (Principal Coordinate Analysis)</option>
|
|
67 </param>
|
|
68 <param name="distance" type="select" display="radio" label="Select Distance method">
|
|
69 <option value="bray" selected="true">bray (Bray-Curtis distance)</option>
|
|
70 <option value="gower">gower (Gower's distance)</option>
|
|
71 <option value="jsd">jds (Jensen-Shannon Divergence)</option>
|
|
72 <option value="unifrac">unifrac (unweighted UniFrac distance)</option>
|
|
73 <option value="wunifrac">wunifrac (weighted-UniFrac distance)</option>
|
|
74 </param>
|
1
|
75 <!-- <param name="kingdom_field" type="select" display="radio" label="Select a taxonomic rank for the analysis">
|
0
|
76 <option value="Kingdom">Kingdom</option>
|
|
77 <option value="Phylum" selected="true">Phylum</option>
|
|
78 <option value="Class">Class</option>
|
|
79 <option value="Order">Order</option>
|
|
80 <option value="Family">Family</option>
|
|
81 <option value="Genus">Genus</option>
|
|
82 <option value="Species">Species</option>
|
|
83 </param>
|
1
|
84 -->
|
|
85 <conditional name="plot_selection">
|
|
86 <param name="plot_type" type="select" label="Type of Plot" help="Ordination plot type (Plot OTUs, Plot SAMPLE, Biplot graphic, Split plot graphic)">
|
|
87 <option value="1">OTU</option>
|
|
88 <option value="2" selected="true">SAMPLE</option>
|
|
89 <option value="3">BIPLOT</option>
|
|
90 <option value="4">SPLIT PLOT</option>
|
|
91 </param>
|
|
92 <when value="1">
|
|
93 <param name="kingdom_field" type="select" display="radio" label="Select a taxonomic rank for the analysis">
|
|
94 <option value="Kingdom">Kingdom</option>
|
|
95 <option value="Phylum" selected="true">Phylum</option>
|
|
96 <option value="Class">Class</option>
|
|
97 <option value="Order">Order</option>
|
|
98 <option value="Family">Family</option>
|
|
99 <option value="Genus">Genus</option>
|
|
100 <option value="Species">Species</option>
|
|
101 </param>
|
|
102 </when>
|
|
103 <when value="2">
|
|
104 </when>
|
|
105 <when value="3">
|
|
106 <param name="kingdom_field" type="select" display="radio" label="Select a taxonomic rank for the analysis">
|
|
107 <option value="Kingdom">Kingdom</option>
|
|
108 <option value="Phylum" selected="true">Phylum</option>
|
|
109 <option value="Class">Class</option>
|
|
110 <option value="Order">Order</option>
|
|
111 <option value="Family">Family</option>
|
|
112 <option value="Genus">Genus</option>
|
|
113 <option value="Species">Species</option>
|
|
114 </param>
|
|
115 </when>
|
|
116 <when value="4">
|
|
117 <param name="kingdom_field" type="select" display="radio" label="Select a taxonomic rank for the analysis">
|
|
118 <option value="Kingdom">Kingdom</option>
|
|
119 <option value="Phylum" selected="true">Phylum</option>
|
|
120 <option value="Class">Class</option>
|
|
121 <option value="Order">Order</option>
|
|
122 <option value="Family">Family</option>
|
|
123 <option value="Genus">Genus</option>
|
|
124 <option value="Species">Species</option>
|
|
125 </param>
|
|
126 </when>
|
|
127 </conditional>
|
0
|
128 </inputs>
|
|
129 <outputs>
|
1
|
130 <data format="txt" name="log" label="${tool.name} - stress value - distance matrix - sample column sum"/>
|
0
|
131 <data format="html" name="htmlfile" label="${tool.name} - ${ord_method.value_label}_${distance.value_label}.html"/>
|
|
132 </outputs>
|
|
133
|
|
134 <tests>
|
|
135 <test> <!-- Test #1: Test BIOM format input -->
|
|
136 <!-- Equivalent command (replace + with double dash): Rscript phyloseq_nmds.R +biom=test-data/GP.biom +subset=6 +method=NMDS
|
|
137 +distance=bray +kingdom=Phylum +cutoff=5 +keep=5 +outdir=outputdir +htmlfile=biom_out.html -->
|
|
138 <conditional name="file_source">
|
|
139 <param name="file_source_selector" value="set_biom"/>
|
|
140 <param name="input" value="GP.biom" ftype="biom1"/>
|
|
141 <param name="subsetColumn" value="Primer" />
|
|
142 </conditional>
|
|
143 <output name="htmlfile" ftype="html" file="biom_out.html" />
|
|
144 </test>
|
|
145 <test> <!-- Test #2: Test TABULAR format inputs -->
|
|
146 <!-- Equivalent command (replace + with double dash): Rscript phyloseq_nmds.R +otu_table=test-data/GP_OTU_TABLE.txt
|
|
147 +tax_table=test-data/GP_TAX_TABLE.txt +meta_table=test-data/GP_SAMPLE_TABLE.txt
|
|
148 +method="bray" +kingdom=Phylum +cutoff=5 + keep=5 +category=6 +outdir=outputdir +htmlfile=test.html -->
|
|
149 <conditional name="file_source">
|
|
150 <param name="file_source_selector" value="set_table" />
|
|
151 <param name="OTU_TABLE" value="GP_OTU_TABLE.txt" />
|
|
152 <param name="TAX_TABLE" value="GP_TAX_TABLE.txt" />
|
|
153 <param name="META_TABLE" value="GP_SAMPLE_TABLE.txt" />
|
|
154 <param name="category" value="6" />
|
|
155 </conditional>
|
|
156 <output name="htmlfile" ftype="html" file="test.html" />
|
|
157 </test>
|
|
158 </tests>
|
|
159
|
|
160 <help>
|
|
161 **What it does**
|
1
|
162 Creates various ordination plot using an R_package_called_phyloseq_,specifically according to the methods outlined in the Phyloseq_Ordination_Plot_Tutorial_.
|
0
|
163
|
1
|
164 .. _R_package_called_phyloseq: https://joey711.github.io/phyloseq/index.html
|
|
165 .. _Phyloseq_Ordination_Plot_Tutorial: https://joey711.github.io/phyloseq/plot_ordination-examples.html
|
0
|
166
|
|
167 -----
|
|
168
|
|
169 **Input**
|
|
170
|
1
|
171 - **Choose an input file type** - either 1 BIOM File or 3 TABULAR files (OTU, taxonomy and Metadata tables)
|
|
172 - **Classify plot by** - points on the plot will be colour coded according to their value of the selected feature
|
|
173 - **Select Ordination method** - several methods are supported
|
|
174 - **Select Distance method** - several methods are supported
|
0
|
175 - **Select a taxonomic rank for the analysis** - select a taxonomy for ordination plot
|
|
176
|
|
177 -----
|
|
178
|
|
179 =========
|
|
180 Resources
|
|
181 =========
|
|
182
|
|
183 **Wrapper Authors**
|
|
184
|
|
185 QFAB Bioinformatics (support@qfab.org)
|
|
186 Melbourne Bioinformatics
|
|
187 </help>
|
|
188 <citations>
|
|
189 <citation type="doi">10.18129/B9.bioc.phyloseq</citation>
|
|
190 </citations>
|
|
191 </tool>
|