Mercurial > repos > jasper > cluster_picker

changeset 0:597bff809239 draft default tip

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author jasper
date Thu, 13 Apr 2017 16:23:47 -0400
parents
children
files cluster_picker-661b74ba4cc8/cluster-picker.xml cluster_picker-661b74ba4cc8/test-data/._.DS_Store cluster_picker-661b74ba4cc8/test-data/dummy.fasta cluster_picker-661b74ba4cc8/test-data/dummy.figTree cluster_picker-661b74ba4cc8/test-data/dummy.nwk
diffstat 5 files changed, 83 insertions(+), 0 deletions(-) [+]
line wrap: on
line diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/cluster_picker-661b74ba4cc8/cluster-picker.xml	Thu Apr 13 16:23:47 2017 -0400
@@ -0,0 +1,73 @@
+<tool id="cluster-pick" name="Cluster Picker" version="1.2.3">
+  <description>cluster picker</description>
+<requirements>
+    <requirement type="package" version="1.2.3">cluster-picker</requirement>
+</requirements>
+<command> <![CDATA[ln -s '$sequence_file' seq.dat &&
+    ln -s '$tree_file' tree.dat &&
+    cluster-picker seq.dat tree.dat $init_threshold $main_threshold $dis_threshold $size $genetic_distance
+  ]]> </command>
+<inputs>
+  <param format="fasta" name="sequence_file" type="data" label="align"
+    help="Please enter fasta format sequences file name to process"/>
+  <param format="nhx" name="tree_file" type="data" label="tree"
+    help="Please enter newick format tree, with branch lengths and node support"/>
+  <param name="init_threshold" type="float" value="0.9" optional="true"
+    help="Please enter initial support threshold for these initial subtrees" label="initial threshold">
+  </param>
+  <param name="main_threshold" type="float" value="0.9" optional="true"
+    help="Please enter main support threshold for clusters" label="main threshold">
+  </param>
+  <param name="dis_threshold" type="float" value="0.045" optional="true"
+    help="Please enter genetic distance threshold for clusters in %" label="distance threshold">
+  </param>
+  <param name="size" type="integer" value="10" optional="true"
+    help="To output all cluster names for clusters of size >= X enter X" label="size">
+  </param>
+  <param name="genetic_distance" type="select" label="genetic distance"
+    help="Please enter scoring type for genetic distance\n
+          gap      	= disregard sites with -, ~, or n\n
+          abs      	= count absolute character differences\n
+          valid    	= only count differences for sites with nucleotides: a, c, t, g in both sequences\n
+          ambiguity	= disregard sites with -, ~, or n and do not count ambiguities as differences (e.g. a vs r is not a difference)">
+      <option value="gap">gap</option>
+      <option value="abs">abs</option>
+      <option value="valid">valid</option>
+      <option value="ambiguity">ambiguity</option>
+  </param>
+</inputs>
+<outputs>
+  <data name='fig' format="figTree" from_work_dir="tree_clusterPicks.nwk.figTree"/>
+  <data name='tr' format="nhx" from_work_dir="tree_clusterPicks.nwk"/>
+  <data name='log' format="txt" from_work_dir="tree_clusterPicks_log.txt"/>
+  <data name='seq' format="fasta" from_work_dir="seq.dat_tree_clusterPicks.fasta"/>
+</outputs>
+
+<tests>
+    <test>
+        <param name="sequence_file" value="dummy.fasta" />
+        <param name="tree_file" value="dummy.nhx" />
+        <param name="init_threshold" value="0.9" />
+        <param name="main_threshold" value="0.9" />
+        <param name="dis_threshold" value="0.045" />
+        <param name="size" value="10" />
+        <param name="genetic_distance" value="gap" />
+        <param name="main_threshold" file="0.9" />
+        <output name="fig" file="dummy.figTree" />
+    </test>
+</tests>
+
+<help>
+Cluster identification strategies differ between studies and as a consequence cluster definitions vary.
+</help>
+<citations>
+   <citation type="bibtex">@ARTICLE{cluster-picker,
+   author = {Ragonnet-Cronin, Manon, et al;},
+   title = {Automated analysis of phylogenetic clusters},
+   journal = {BMC bioinformatics},
+   year = {2013},
+   volume = {14.1},
+   pages = {317}}
+ </citation>
+</citations>
+</tool>
Binary file cluster_picker-661b74ba4cc8/test-data/._.DS_Store has changed
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/cluster_picker-661b74ba4cc8/test-data/dummy.fasta	Thu Apr 13 16:23:47 2017 -0400
@@ -0,0 +1,8 @@
+-e >seq1
+ATTA
+>seq2
+ATTA
+>seq3
+ATTA
+>seq4
+ATTA
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/cluster_picker-661b74ba4cc8/test-data/dummy.figTree	Thu Apr 13 16:23:47 2017 -0400
@@ -0,0 +1,1 @@
+#NEXUS
begin taxa;
dimensions ntax=4;
	taxlabels
	seq1
	seq2
	seq3
	seq4
	;
end;
begin trees;
	tree tree_1 = [&R] (((seq1:0.001,seq2:0.002)[&support=0.0]:0.001,seq3:0.005)[&support=0.0]:0.004,seq4:0.003)[&support=0.0]:0.001;
end;
begin figtree;
	set appearance.backgroundColorAttribute="User Selection";
	set appearance.backgroundColour=#-1;
	set appearance.branchColorAttribute="User Selection";
	set appearance.branchLineWidth=1.0;
	set appearance.foregroundColour=#-16777216;
	set appearance.selectionColour=#-2144520576;
	set branchLabels.colorAttribute="User Selection";
	set branchLabels.displayAttribute="Branch times";
	set branchLabels.fontName="sansserif";
	set branchLabels.fontSize=8;
	set branchLabels.fontStyle=0;
	set branchLabels.isShown=false;
	set branchLabels.significantDigits=4;
	set layout.expansion=0;
	set layout.layoutType="RECTILINEAR";
	set layout.zoom=0;
	set nodeBars.barWidth=4.0;
	set nodeLabels.colorAttribute="User Selection";
	set nodeLabels.displayAttribute="support";
	set nodeLabels.fontName="Arial";
	set nodeLabels.fontSize=10;
	set nodeLabels.fontStyle=0;
	set nodeLabels.isShown=true;
	set nodeLabels.significantDigits=3;
	set polarLayout.alignTipLabels=false;
	set polarLayout.angularRange=0;
	set polarLayout.rootAngle=0;
	set polarLayout.rootLength=100;
	set polarLayout.showRoot=true;
	set radialLayout.spread=0.0;
	set rectilinearLayout.alignTipLabels=false;
	set rectilinearLayout.curvature=0;
	set rectilinearLayout.rootLength=100;
	set scale.offsetAge=0.0;
	set scale.rootAge=1.0;
	set scale.scaleFactor=1.0;
	set scale.scaleRoot=false;
	set scaleAxis.automaticScale=true;
	set scaleAxis.fontSize=8.0;
	set scaleAxis.isShown=false;
	set scaleAxis.lineWidth=1.0;
	set scaleAxis.majorTicks=1.0;
	set scaleAxis.origin=0.0;
	set scaleAxis.reverseAxis=false;
	set scaleAxis.showGrid=true;
	set scaleAxis.significantDigits=4;
	set scaleBar.automaticScale=true;
	set scaleBar.fontSize=10.0;
	set scaleBar.isShown=true;
	set scaleBar.lineWidth=1.0;
	set scaleBar.scaleRange=0.0;
	set scaleBar.significantDigits=4;
	set tipLabels.colorAttribute="User Selection";
	set tipLabels.displayAttribute="Names";
	set tipLabels.fontName="Arial";
	set tipLabels.fontSize=10;
	set tipLabels.fontStyle=0;
	set tipLabels.isShown=true;
	set tipLabels.significantDigits=4;
	set trees.order=true;
	set trees.orderType="increasing";
	set trees.rooting=false;
	set trees.rootingType="User Selection";
	set trees.transform=false;
	set trees.transformType="cladogram";
end;
\ No newline at end of file
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/cluster_picker-661b74ba4cc8/test-data/dummy.nwk	Thu Apr 13 16:23:47 2017 -0400
@@ -0,0 +1,1 @@
+(((seq1:0.001,seq2:0.002):0.001,seq3:0.005):0.004,seq4:0.003):0.001;